pith. sign in

arxiv: 2412.10142 · v3 · submitted 2024-12-13 · 🧮 math.AP

Nonlinear discrete Schr\"odinger equations with a point defect

Dirk Hennig This is my paper

Pith reviewed 2026-05-23 07:31 UTC · model grok-4.3

classification 🧮 math.AP
keywords discrete nonlinear Schrödinger equationdelta potentialpoint defectexcitation thresholdlocalized ground statesfocusing nonlinearityscattering
0
0 comments X

The pith

A delta potential imposes explicit l2-norm thresholds for exponentially localized ground states in the discrete nonlinear Schrödinger equation under focusing nonlinearity.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper studies the d-dimensional discrete nonlinear Schrödinger equation with general power nonlinearity and a single delta potential modeling a point defect. It proves existence of spatially exponentially localized time-periodic ground states for focusing nonlinearity and derives how an attractive or repulsive defect shifts the lower excitation threshold given by a supercritical l2 norm, supplying explicit formulas for those thresholds. A reader would care because the threshold quantifies when the defect and the nonlinearity together trap energy rather than letting it disperse across the lattice. The work further supplies upper l2-norm bounds under which defocusing nonlinearity preserves the linear bound states created by the defect and shows that solutions with subcritical norm scatter to the linear problem in higher integrability spaces.

Core claim

For focusing nonlinearity we prove the existence of a spatially exponentially localized and time-periodic ground state and investigate the impact of an attractive respectively repulsive delta potential on the existence of an excitation threshold, i.e. supercritical l2 norm, for the creation of such ground states. Explicit expressions for the lower excitation thresholds are given. Reciprocally, we discuss the influence of defocusing nonlinearity on the durability of the linear bound states and provide upper thresholds of the l2-norm for their preservation. Regarding the asymptotic behavior of the solutions we establish that for a l2-norm below the excitation threshold the solutions scatter to

What carries the argument

The delta potential as a point defect that breaks translational invariance to produce a linear bound state of negative or positive energy, combined with the focusing power nonlinearity that produces self-trapping once the l2 norm exceeds an explicit threshold.

If this is right

  • Explicit lower l2-norm thresholds for ground-state existence can be computed for both attractive and repulsive defects.
  • Linear bound states persist under defocusing nonlinearity only when the l2 norm stays below explicit upper thresholds.
  • Solutions whose l2 norm lies below the excitation threshold scatter to linear solutions in l^{p>2}.
  • The ground states remain exponentially localized in space for all times.
  • The results apply in any dimension d and for any power of the nonlinearity.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The explicit thresholds could be used to predict the minimal excitation strength needed to localize energy around a single impurity in a discrete lattice.
  • The same competition between defect-induced binding and nonlinearity may appear in other discrete nonlinear wave models with localized potentials.
  • Numerical continuation from the linear bound state could be used to track how the threshold changes with defect strength.
  • The scattering result implies that weak excitations around the defect behave essentially linearly at long times.

Load-bearing premise

The standard functional-analytic properties of the discrete Laplacian with the delta potential together with the specific power nonlinearity suffice to close the existence arguments without extra restrictions on dimension or exponent.

What would settle it

A direct numerical simulation of the time-dependent equation in one dimension with norm slightly above the claimed threshold that fails to converge to a time-periodic exponentially localized profile would falsify the existence statement.

read the original abstract

We study the $d$-dimensional discrete nonlinear Schr\"odinger equation with general power nonlinearity and a delta potential. Our interest lies in the interplay between two localization mechanisms. On the one hand, the attractive (repulsive) delta potential acting as a point defect breaks the translational invariance of the lattice so that a linear staggering (non-staggering) bound state is formed with negative (positive) energy. On the other hand, focusing nonlinearity may lead to self-trapping of excitation energy. For focusing nonlinearity we prove the existence of a spatially exponentially localized and time-periodic ground state and investigate the impact of an attractive respectively repulsive delta potential on the existence of an excitation threshold, i.e. supercritical $l^2$ norm, for the creation of such a ground states. Explicit expressions for the lower excitation thresholds are given. Reciprocally, we discuss the influence of defocusing nonlinearity on the durability of the linear bound states and provide upper thresholds of the $l^2-$norm for their preservation. Regarding the asymptotic behavior of the solutions we establish that for a $l^2-$norm below the excitation threshold the solutions scatter to a solution of the linear problem in $l^{p>2}$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript studies the d-dimensional discrete nonlinear Schrödinger equation with a delta point defect and power nonlinearity. For focusing cases it proves existence of exponentially localized time-periodic ground states via variational methods, derives explicit lower excitation thresholds (supercritical ℓ² norms) that depend on the sign and strength of the defect, examines how defocusing nonlinearity affects persistence of linear bound states (with upper thresholds), and shows scattering to linear solutions below the excitation threshold.

Significance. The explicit threshold formulas constitute a concrete, falsifiable contribution that distinguishes the interplay between defect-induced and self-trapping localization. The compactness argument that exploits the delta potential to obtain Palais-Smale condition uniformly in dimension and exponent is standard yet cleanly applied here; the resulting statements hold without artificial restrictions on d or σ.

major comments (2)
  1. [§4.2, Theorem 4.3] §4.2, Theorem 4.3: the lower threshold formula is stated as explicit, yet the proof invokes a mountain-pass value whose dependence on the defect parameter is only shown to be monotone; an explicit closed-form expression would require solving the associated linear eigenvalue problem explicitly, which is only done for d=1. Clarify whether the formula is closed-form or merely monotone in the defect strength.
  2. [§5] §5, the scattering statement below threshold: the proof uses Strichartz-type estimates on the linear evolution with delta potential. The decay rates quoted appear to rely on the spectral gap of the linear operator; verify that the constant in the dispersive estimate remains uniform when the defect strength approaches the critical value at which the linear bound state disappears.
minor comments (2)
  1. [Notation] Notation: the symbol E_p for the excitation threshold is introduced in §3 but reused with different subscripts in §4; a single consistent definition table would help.
  2. [Figure 1] Figure 1 caption: the plotted profiles are labeled “ground states” but the caption does not indicate whether they correspond to the mountain-pass or to the linear eigenfunction; add a brief clarification.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and the constructive comments. We address each major comment point-by-point below. Minor revisions will be made to improve clarity on the variational characterization of the thresholds and the dependence of dispersive constants.

read point-by-point responses

  1. Referee: [§4.2, Theorem 4.3] §4.2, Theorem 4.3: the lower threshold formula is stated as explicit, yet the proof invokes a mountain-pass value whose dependence on the defect parameter is only shown to be monotone; an explicit closed-form expression would require solving the associated linear eigenvalue problem explicitly, which is only done for d=1. Clarify whether the formula is closed-form or merely monotone in the defect strength.

    Authors: We agree that the lower excitation thresholds are expressed via the mountain-pass value of the constrained functional, which is shown to depend monotonically on the defect strength λ. This value is explicit in the sense that it is given by a direct variational formula (the infimum of the energy on the Nehari manifold subject to the ℓ²-norm constraint), without additional implicit parameters. For d=1 the linear eigenvalue problem is solvable in closed form, yielding an algebraic expression for the threshold; for general d the expression remains variational but is still explicit and monotone. We will revise the statement of Theorem 4.3 and the surrounding discussion to clarify this distinction between closed-form (d=1) and variational-explicit (general d) expressions. revision: partial

  2. Referee: [§5] §5, the scattering statement below threshold: the proof uses Strichartz-type estimates on the linear evolution with delta potential. The decay rates quoted appear to rely on the spectral gap of the linear operator; verify that the constant in the dispersive estimate remains uniform when the defect strength approaches the critical value at which the linear bound state disappears.

    Authors: The Strichartz estimates employed in the scattering argument are the standard ones for the linear discrete Schrödinger operator with a fixed delta defect of strength λ, where |λ| lies strictly below the critical value at which the linear bound state ceases to exist; this guarantees a positive spectral gap. The constants in these estimates depend on the size of the gap and are therefore uniform for each fixed λ. As λ approaches the critical value the gap closes and the constants deteriorate, but the theorem is formulated for a given fixed defect strength in the regime where the bound state persists. We will add a short remark after the statement of the scattering result noting the dependence of the constants on the spectral gap and that uniformity holds when λ is bounded away from criticality. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained via standard analysis

full rationale

The paper establishes existence of ground states and excitation thresholds for the discrete NLS with delta defect via variational methods applied to the stationary profile equation obtained from the time-periodic ansatz. Compactness follows from the defect breaking translation invariance in l^2(Z^d), valid for arbitrary d and power without further restrictions. Linear bound states are constructed via the resolvent of the discrete Laplacian plus delta, a standard construction. No load-bearing self-citations, no fitted parameters renamed as predictions, and no self-definitional reductions appear in the claimed results; all thresholds are derived from the functional-analytic properties of the equation itself rather than by construction from the paper's inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard discrete Sobolev embeddings, variational methods for existence, and asymptotic analysis for scattering; no free parameters are fitted to data and no new entities are postulated.

axioms (2)
  • standard math Standard properties of the discrete Laplacian and embedding theorems in l^p spaces on Z^d
    Invoked implicitly to establish existence of bound states and to control the nonlinearity.
  • domain assumption Existence of linear bound states for the delta potential
    Used as the starting point for the interplay with nonlinearity.

pith-pipeline@v0.9.0 · 5734 in / 1367 out tokens · 29444 ms · 2026-05-23T07:31:49.292701+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

53 extracted references · 53 canonical work pages

  1. [1]

    Morsch, M

    O. Morsch, M. Oberthaler, Dynamics of Bose–Einstein condensates in optical lattice . Rev. Mod. Phys. 78, 179–215 (2006); M.A. Porter, R. Carretero-Gonz´ alez, P.G. Kevrekidis, and B. Malomed, Nonlinear lattice dynamics of Bose–Einstein condensates , Chaos 15, 015115 (2005)

    work page 2006

  2. [2]

    Skorobogatiy and J

    M. Skorobogatiy and J. Yang, Fundamentals of photonic crystal guiding Cambridge University Press: Cambridge, UK, (2008)

    work page 2008

  3. [3]

    Christodoulides and R.I

    D.N. Christodoulides and R.I. Joseph, Discrete self-focusing in nonlinear arrays of coupled wave guides, Opt. Lett. 13, 794–796 (1988); H.S. Eisenberg, Y. Silberberg, R. Morando tti, A.R. Boyd, and J.S. Aitchison, Discrete spatial optical solitons in waveguide arrays , Phys. Rev. Lett. 81, 3383–3386 (1988); D.N. Christodoulides, F. Lederer, Y. Silberberg,...

    work page 1988

  4. [4]

    F. Ye, D. Mihalache, B. Hu, and N.C. Panoiu, Subwavelength plasmonic lattice solitons in arrays of meta llic nanowires, Phys. Rev. Lett. 104, 106802 (2010)

    work page 2010

  5. [5]

    P. G. Kevrekidis, K. O. Rasmussen and A. R. Bishop, The discrete nonlinear Schr¨ odinger equation: A survey of recent results, Int. Journal of Modern Physics B 15, 2833–2900 (2001)

    work page 2001

  6. [6]

    Eilbeck and M

    J.C. Eilbeck and M. Johansson, The discrete nonlinear Sc hr¨ odinger equation-20 years on in: L. V´ azquez, R.S. MacKay, M.P. Zorzano (Eds.), Localization and Energy Transfer in Nonlinear Systems. World Scientific, Singapore, pp. 44–67 (2003)

    work page 2003

  7. [7]

    Kevredikis, The Nonlinear Discrete Schr¨ odinger Equation: Mathematical Analysis, Numerial Computations, and Physical Perspectives (Springer-Verlag, Berlin, Heidelberg, 2009)

    P.G. Kevredikis, The Nonlinear Discrete Schr¨ odinger Equation: Mathematical Analysis, Numerial Computations, and Physical Perspectives (Springer-Verlag, Berlin, Heidelberg, 2009)

    work page 2009

  8. [8]

    Cazenave and P.L

    T. Cazenave and P.L. Lions, Orbital stability of standing waves for some nonlinear Schr ¨ odinger equation, Comm. Math. Phys. 85, 549-561 (1982)

    work page 1982

  9. [9]

    Cazenave, An introduction to nonlinear Schr¨ odinger equations, IM-UFRJ, Rio de Janeiro (1989)

    T. Cazenave, An introduction to nonlinear Schr¨ odinger equations, IM-UFRJ, Rio de Janeiro (1989)

    work page 1989

  10. [10]

    Cazenave, Semilinear Schr¨ odinger equations, vol

    T. Cazenave, Semilinear Schr¨ odinger equations, vol. 10 Courant Lecture Notes in Mathematics, AMS, Provide nce (2003). NONLINEAR DISCRETE SCHR ¨ODINGER EQUATIONS WITH A POINT DEFECT 23

    work page 2003

  11. [11]

    Stein, Harmonic analysis, real-variable methods, orthogonality and oscillatory integrals , Princeton University Press, 1993

    E. Stein, Harmonic analysis, real-variable methods, orthogonality and oscillatory integrals , Princeton University Press, 1993

    work page 1993

  12. [12]

    Weinstein, Excitation thrsholds for nonlinear localized modes on latt ices, Nonlinearity 12

    M.I. Weinstein, Excitation thrsholds for nonlinear localized modes on latt ices, Nonlinearity 12. 673–691 (1999)

    work page 1999

  13. [13]

    Weinstein, Lyapunov stability of ground states of nonlinear dispersiv e evolution equations , Comm

    M.I. Weinstein, Lyapunov stability of ground states of nonlinear dispersiv e evolution equations , Comm. Pure Appl. Math. 39, 51 (1986)

    work page 1986

  14. [14]

    Cuevas, N.I

    J. Cuevas, N.I. Karachalios, and F. Palmero, Lower and upper estimates on the excitation threshold for br eathers in DNLS lattices , XXX

    work page

  15. [15]

    Caudrelier, M

    V. Caudrelier, M. Mintchev, E. Ragoucy, Solving the quantum nonlinear Schr¨ odinger equation with δ-type impurity, J. Math. Phys. 46 1-24 (2005)

    work page 2005

  16. [16]

    Bellazi, M

    B. Bellazi, M. Mintchev, Quantum field theory on star graphs , J. Phys. A, Math. Theor. 39 1101–1117 (2006)

    work page 2006

  17. [17]

    Le Coz, R

    S. Le Coz, R. Fukuizumi, G. Fibich, B. Ksherim, and Y. Siv an, Instability of bound states of a nonlinear Schr¨ odinger equation with a Dirac potential, Physica D 237 1103–1128 (2008)

    work page 2008

  18. [18]

    Fukuizumi, and L

    R. Fukuizumi, and L. Jeanjean, Stability of standing waves for a nonlinear Schr¨ odinger eq uation with a repulsive Dirac delta potential , Discrete Contin. Dyn. Syst. 21 121–136 (2008)

    work page 2008

  19. [19]

    Fukuizumi, M

    R. Fukuizumi, M. Ohta, T. Ozawa, Nonlinear Schr¨ odinger equation with a point defect, Ann. Inst. Henri Poincar´ e, Anal. Non Lin´ eaire25 837–845 (2008)

    work page 2008

  20. [20]

    Deift and J

    P. Deift and J. Park, Long-time asymptotic stability for solutions of the NLS equ ation with a delta potential amd even initial data , Int. Math. Res. Notices 24, 5505–5624 (2011)

    work page 2011

  21. [21]

    Belloni, R.W

    M. Belloni, R.W. Richard, The infinite well and Dirac delta function potentials as peda gogical, mathematical and physical models in quantum mechanics , Phys. Rep. 540 25-122 (2014)

    work page 2014

  22. [22]

    Banica and N

    V. Banica and N. Visciglia, Scattering for NLS with a delta potential , J. Diff. Eqs. 260, 4410–4439 (2016)

    work page 2016

  23. [23]

    Ikeda and T

    M. Ikeda and T. Inui, Global dynamics below the standing waves for the focusing se milinear Schr¨ odinger equation with a repulsive Dirac delta potential , Anal. PDE 10 481–512 (2017)

    work page 2017

  24. [24]

    Ardila and T

    A.H. Ardila and T. Inui, Threshold scattering for the focusing NLS with a repulsive D irac delta potential , J. Diff. Eqs. 313, 54-84 (2022)

    work page 2022

  25. [25]

    Merle, Determination of blow-up solutions with mini mal mass for nonlinear Schr¨ odinger equations with criticl power, Duke Math

    F. Merle, Determination of blow-up solutions with mini mal mass for nonlinear Schr¨ odinger equations with criticl power, Duke Math. J 69 427-454 (1993)

    work page 1993

  26. [26]

    Merle and P

    F. Merle and P. Raphael, On universality of blow-up profile for L2 critical nonlinear Schr¨ odinger equation, Invent. Math. 156 565–672 (2004)

    work page 2004

  27. [27]

    Merle and P

    F. Merle and P. Raphael, The blow-up dynamic and upper bound on the blow-up rate for cr itical nonlinear Schr¨ odinger equation, Ann. of Math. 161 157–222 (2005)

    work page 2005

  28. [28]

    Perelman, On the blow up phenomenon for the critical nonlinear Schr¨ od inger equation in 1D , Ann

    G. Perelman, On the blow up phenomenon for the critical nonlinear Schr¨ od inger equation in 1D , Ann. Henri. Poincar´ e2 605–673 (2001)

    work page 2001

  29. [29]

    Planchon and P

    F. Planchon and P. Raphael, Existence and stability of the log-log blow-up dynamics for the L2-critical nonlinear Schr¨ odinger equation in a domain, Ann. Henri Poincar´ e,8 1177–1219 (2007)

    work page 2007

  30. [30]

    Raphael, Stability of the log-log bound for blow up solutions to the cr itical non linear Schr¨ odinger equation, Math

    P. Raphael, Stability of the log-log bound for blow up solutions to the cr itical non linear Schr¨ odinger equation, Math. Ann. 331 577–609 (2005)

    work page 2005

  31. [31]

    Goodman, P.J

    R.H. Goodman, P.J. Holmes, and M.I. Weinstein, Strong NLS soliton-defect interactions, Physica D 192, 215–248 (1987)

    work page 1987

  32. [32]

    Dodson, Global well-posedness and scattering for the mass critical nonlinear Schr¨ odinger equation with mass below the mass of the ground state , Adv

    B. Dodson, Global well-posedness and scattering for the mass critical nonlinear Schr¨ odinger equation with mass below the mass of the ground state , Adv. Math. 285 1589-1618 (2015)

    work page 2015

  33. [33]

    Guerekata and A

    G.M. Guerekata and A. Pankov, Global well-posedness for discrete nonlinear Schr¨ odinger equation , Appl. Anal. 89 1513-1521 (2010)

    work page 2010

  34. [34]

    Hennig, N.I

    D. Hennig, N.I. Karachalios, amd J. Cuevas-Maraver, The closeness of the Ablowitz-Ladik lattice to the discrete nonlinear Schr¨ odinger equation, J. Diff. Eqs. 316 346-363 (2022)

    work page 2022

  35. [35]

    Wang, The well-posedness and scattering theory of nonlinear Schr ¨ odiner equations on lattice , arXiv: 2408.01174v1

    J. Wang, The well-posedness and scattering theory of nonlinear Schr ¨ odiner equations on lattice , arXiv: 2408.01174v1

    work page arXiv

  36. [36]

    Zakharov, S.V

    V.E. Zakharov, S.V. Manakov,S.P. Novikov, and L.P. Pit aevskii, Solitons: The Inverse Scattering Method Nauka Publishers: Moscow, Russia, 1980; English Translation: Co nsultants Bureau: New York, NY, USA, 1984

    work page 1980

  37. [37]

    Ablowitz and H

    M.J. Ablowitz and H. Segur, Solitons and Inverse Scattering Method SIAM: Philadelphia, PA, USA, 1981

    work page 1981

  38. [38]

    Calogero and A

    F. Calogero and A. Degasperis, Spectral Transform and Solitons: Tools to Solve and Investi gate Nonlinear Evolution Equations North-Holland: New York, NY, USA, 1982

    work page 1982

  39. [39]

    Newell, Solitons in Mathematics and Physics SIAM: Philadelphia, PA, USA, 1985

    A.C. Newell, Solitons in Mathematics and Physics SIAM: Philadelphia, PA, USA, 1985. 24 DIRK HENNIG

    work page 1985

  40. [40]

    Ablowitz and B.M

    M.J. Ablowitz and B.M. Herbst, On homoclinic structure and numerically induced chaos for t he nonlinear Schr¨ odinger equationSIAM J. Appl. Math. 50, 339–351 (1990)

    work page 1990

  41. [41]

    D. Levi, M. Petrera, and C. Scimiterna, On the integrability of the discrete nonlinear Schr¨ odinge r equation , Europhys. Lett. 84, 10003 (2008)

    work page 2008

  42. [42]

    Komech, E.A

    A.I. Komech, E.A. Kopylova, and M. Kunze, Dispersive estimates for the 1D discrete Schr¨ odinger and K lein- Gordon equationa, Appl. Anal. 85, 1487-1508 (2005)

    work page 2005

  43. [43]

    Stefanov and P.G

    A. Stefanov and P.G. Kevrekidis, Asymptotic behaviour of small solutions for the discrete Sc hr¨ odinger and Klein-Gordon equations, Nonlinearity 18, 1841–1857 (2005)

    work page 2005

  44. [44]

    Berestycki and T

    H. Berestycki and T. Cazenave, Instabile des etats stationaires dans le equations de Schr¨ odinger et Klein-Gordon nonlineares, C. R. Acad. Sci. Paris 293, 489–492 (1981)

    work page 1981

  45. [45]

    Sulem and P.L

    C. Sulem and P.L. Sulem, The nonlinear Schr¨ odinger equation, self-focusing and wa ve collapse , Applied Mathematical Sciences 139, Springer-Verlag, New York (199 9)

    work page

  46. [46]

    Weinstein, Nonlinear Schr¨ odinger equations and sharp interpolation estimates, Comm

    M.I. Weinstein, Nonlinear Schr¨ odinger equations and sharp interpolation estimates, Comm. Math, Phys. 87, 567–576 (1983)

    work page 1983

  47. [47]

    R. S. MacKay and S. Aubry, Proof of existence of breathers for time-reversible or Hami ltonian networks of weakly coupled oscillators , Nonlinearity 7 (1994), 1623–1643

    work page 1994

  48. [48]

    Lions, The concentration compactness principle in the calculus of variation: The locally compac case , Ann

    P.-L. Lions, The concentration compactness principle in the calculus of variation: The locally compac case , Ann. Inst. Henri Poincar´ e, Analysis Nonlineaire 1, 223 (1984)

    work page 1984

  49. [49]

    Teschl, Jacobi operatos and completely integrable lattices , Americam Mathematical Society, Providence, RI, 2000

    G. Teschl, Jacobi operatos and completely integrable lattices , Americam Mathematical Society, Providence, RI, 2000

    work page 2000

  50. [50]

    Douglas, Banach algebra techniques in operator theory , Graduate Texts in Mathematics, Springer New York (1998)

    R.G. Douglas, Banach algebra techniques in operator theory , Graduate Texts in Mathematics, Springer New York (1998)

    work page 1998

  51. [51]

    Hennig, Breather solutions in conservative and dissipative nonlin ear Klein-Gordon lattices , J

    D. Hennig, Breather solutions in conservative and dissipative nonlin ear Klein-Gordon lattices , J. Fixed Point Theory Appl. 26 (2024)

    work page 2024

  52. [52]

    Hichem and L

    H. Hichem and L. Song, Strict mononotonicity of the global branch of solutions in t he L2 norm and uniqueness of the normalized ground state for various classes of PDEs: T wo general methods with some examples , (2023) arXiv:2302.09681v1

    work page arXiv 2023

  53. [53]

    Zhang and J

    Z. Zhang and J. Sun, Normalized solutions of NLS equations with mixed nonlinear ities, Adv Nonl. Anal. 13 (204), 20240004. Email address, Email:: dirkhennig@uth.gr Department of Mathematics, University of Thessaly, Lamia G R35100, Greece

    work page