Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2602.14150

Approximating the $S$ matrix for solving the Marchenko equation: the case of channels with different thresholds

N. A. Khokhlov

quant-phnucl-thSubmitted Feb 15, 2026• Announced Feb 17, 2026
arXivHTMLPDF

Abstract

This work extends previous results on the inverse scattering problem within the framework of Marchenko theory (fixed-$l$ inversion). In particular, I approximate an $n$-channel $S$-matrix as a function of the first-channel momentum $q$ by a sum of a rational term and a truncated sinc series for each matrix element. Relativistic kinematics are taken into account through the correct momentum-energy relation, and the necessary minor generalization of Marchenko theory is given. For energies where only a subset of scattering channels is open, the analytic structure of the $S$-matrix is analyzed. I demonstrate that the submatrix corresponding to closed channels, particularly near their thresholds, can be reconstructed from the experimentally accessible submatrix of open channels.The convergence of the proposed method is verified by applying it to data generated from a direct solution of the scattering problem for a known potential, and comparing the reconstructed potential with the original one. Finally, the method is applied to the analysis of $S_{31}$ $πN$ scattering data.

Comments
14 pages, 6 figures

Related papers

Inverse scattering using the Heitler equation

Anthony J. Devaney, Arthur B. Weglein

The problem of determining a scattering potential from scattered field data is addressed within the two integral equation formulation of scattering theory first proposed by Goldberger (1964). The lowest-order approximation within this formulation casts the inverse problem entirely in terms of the Heitler equation which relates the on-shell T matrix directly to the scattering potential. This lowest-order approximation is applied to the case of non-relativistic quantum mechanical inverse scattering from one-dimensional local potentials and is found to yield an exact solution for the case of a one-dimensional delta function potential.

Inverse Problems•Dec 1, 1989•8 citations

Two-body interaction. I. Inverse scattering problem, nonlocal potentials, and a conception of quantum measurement theory

V. M. Muzafarov

The present paper generalizes preceding papers of the author [1]–[6] and opens a cycle of works concerning the general quantum-mechanical inverse scattering problem (for a given partial channel) and its solution in analytic form in a nontraditional setting: the scattering wave function in the momentum representation is recovered immediately from the scattering data. The solution we find for the problem in question is given by a family of phase-equivalent wave functions with explicitly controlled functional-analytic arbitrariness. This arbitrariness corresponds to the arbitrariness in the choice of phase-equivalent nonlocal potentials and demands including some additional concepts of quantum measurement theory to a complete description of quantum-mechanical systems. We discuss several related problems of quantum-mechanical probability distributions and their statistical interpretation. The mathematical apparatus of the approach thus developed is illustrated by the example of describing the nucleon-nucleon interaction in uncoupled partial channels. In conclusion we announce the topics of forthcoming works of the series.

arXiv:nucl-th/0702053•Russian Journal of Mathematical Physics•Mar 1, 2007•3 citations

Kirchhoff's Rule for Quantum Wires. II: The Inverse Problem with Possible Applications to Quantum Computers

V. Kostrykin, R. Schrader

In this article we continue our investigations of one particle quantum scattering theory for Schrödinger operators on a set of connected (idealized one-dimensional) wires forming a graph with an arbitrary number of open ends. The Hamiltonian is given as minus the Laplace operator with suitable linear boundary conditions at the vertices (the local Kirchhoff law). In “Kirchhoff's rule for quantum wires” [J. Phys. A: Math. Gen. 32, 595–630 (1999)] we provided an explicit algebraic expression for the resulting (on-shell) S-matrix in terms of the boundary conditions and the lengths of the internal lines and we also proved its unitarity. Here we address the inverse problem in the simplest context with one vertex only but with an arbitrary number of open ends. We provide an explicit formula for the boundary conditions in terms of the S-matrix at a fixed, prescribed energy. We show that any unitary n × n matrix may be realized as the S-matrix at a given energy by choosing appropriate (unique) boundary conditions. This might possibly be used for the design of elementary gates in quantum computing. As an illustration we calculate the boundary conditions associated to the unitary operators of some elementary gates for quantum computers and raise the issue whether in general the unitary operators associated to quantum gates should rather be viewed as scattering operators instead of time evolution operators for a given time associated to a quantum mechanical Hamiltonian. We also suggest an approach by which the S-matrix in our context may be obtained from “scattering experiments”, another aspect of the inverse problem. Finally we extend our previous discussion, how our approach is related to von Neumann's theory of selfadjoint extensions.

arXiv:quant-ph/9910053•Fortschritte der Physik•Aug 1, 2000•110 citations

Design of Quantum Gates Using Quantum Scattering Theory

Harish Parthasarathy, Monika Aggarwal

Given two Hamiltonians \(H_0,H_1=H_0+V\) corresponding respectively to a free particle Hamiltonian and the free particle Hamiltonian plus its interaction energy with a scattering centre, we can compute the S-matrix, i.e., the matrix that determines the quantum mechanical amplitude of the process involving a particle coming from a free particle state at an infinite distance from the scattering centre at time \(t=-\infty\), entering into an initial scattered state, interacting with the scattering centre, entering into a final scattered state, and then escaping away to infinity at time \(t\rightarrow\infty\) to a final free particle state. Such amplitudes were first studied by Lippmann and Schwinger (Steven Weinberg, The Quantum Theory of Fields, Vol. 1, W.O. Amrein, Hilbert Space Methods in Quantum Mechanics) and then analyzed in a mathematically precise way by [T. Kato], [W. Amrein], and [K.B. Sinha]. In this paper, we first derive the Lippmann-Schwinger equation for the input and output scattered states in terms of the input and output free particle states using intuitive arguments originally due to Lippmann and Schwinger, and then, using rigorous operator theoretic arguments revolving around the spectral theorem for unbounded self-adjoint operators in a Hilbert space, we obtain an elegant formula for the scattering matrix elements in terms of the interaction potential and the free particle Hamiltonian between two free particle states corresponding to the same measure. We then discuss a computationally more efficient method based on the Dyson series expansion (Steven Weinberg, The Quantum Theory of Fields, Vol. 1), which works even in the case when the interaction Hamiltonian is time varying, in contrast to the Lippmann-Schwinger method, for calculating the S-matrix. We assume that the scattering potential, which can be time varying, can be controlled by incorporating control parameters and explain how the resulting unitary S-matrix can be made to approximate a given unitary gate in infinite-dimensional Hilbert space by optimizing the Hilbert-Schmidt/Frobenius norm distance between the given gate and the S-matrix gate, in which the latter is approximated by a truncated Dyson series. Of course, the S-matrix given by the Dyson series is a Taylor series functional expansion in the interaction potential and hence is a highly nonlinear function of the control parameters, and hence we propose numerical methods for such optimization. We then observe that if the control parameters are taken as random variables with joint moments, then the TPCP map involving statistical averaging of the adjoint action of the S-matrix acting on an initial mixed state w.r.t. the probability distribution of the control parameters will result in a final mixed state that is linear in the statistical moments of the control parameters, and hence optimization w.r.t. these moments can be easily carried out. Our method of TPCP map design is to optimize the sum of Frobenius distance squares between a sequence of desired output mixed states and the S-matrix based formula for the TPCP map acting on the corresponding input mixed states w.r.t. the statistical moments of the control parameters. We then generalize this method to include the case of TPCP maps obtained from the Hudson-Parthasarathy quantum noisy Schrödinger equation. In summary, we simulate the S-matrix gate for independent realizations of the control parameters having a probability distribution defined by their joint moments obtained by the above optimization procedure and apply the adjoint action of these independent S-matrices on a given input state, and take their ensemble average to obtain a good approximation to the action of the desired TPCP map on any given input state. This suggests that TPCP maps of arbitrarily large size can be realized using scattering experiments in the laboratory.

Qeios•Apr 16, 2024•0 citations

Relativistic J-matrix method of scattering in 1 + 1 space-time I. Theory *

Abdulaziz D Alhaidari

Abstract We make a relativistic extension of the one-dimensional J-matrix method of scattering. The relativistic potential matrix is a combination of vector, scalar, and pseudo-scalar components. These are non-singular short-range potential functions (not necessarily analytic like a square well) such that they are well represented by their matrix elements in a finite subset of a square integrable basis set. This set is chosen to support a tridiagonal symmetric matrix representation for the free Dirac operator. We derive the transmission and reflection coefficients. This work will be followed by another where we apply the theory to obtain scattering information for different potential coupling modes including those with evidence of Klein tunneling and of supercritical resonances.

arXiv:2001.06298•Journal of Physics Communications•Sep 1, 2020•0 citations

Low-energy positron scattering from metastable helium

Anonymous, Hui-Li Han, Ting-Yun Shi +1 more

Low-energy positron scattering from singlet and triplet metastable He($1s2s$) is investigated using the $R$-matrix propagation method in hyperspherical coordinates. Elastic and positronium-formation cross sections are reported, and near-threshold resonance structures are analyzed in terms of eigenphase sums and the time-delay matrix. For the triplet target, the calculated cross sections are in agreement with available convergent-close-coupling results and display the expected threshold behavior. Beyond the known $S$-wave features, Feshbach resonance series in higher partial waves extending up to highly excited atomic thresholds are systematically identified. The time-delay matrix further uncovers hidden resonances that produce obvious structures in the positronium-formation cross sections.

arXiv:2606.29253•Physical Review A•Sep 17, 2026•0 citations
Approximating the SSS matrix for solving the Marchenko equation: the case of channels with different thresholds
This work extends previous results on the inverse scattering problem within the framework of Marchenko theory (fixed-lll inversion). In particular, I approximate an nnn-channel SSS-matrix as a function of the first-channel momentum qqq by a sum of a rational term and a truncated sinc series for each matrix element. Relativistic kinematics are taken into account through the correct momentum-energy relation, and the necessary minor generalization of Marchenko theory is given.
For energies where only a subset of scattering channels is open, the analytic structure of the SSS-matrix is analyzed. I demonstrate that the submatrix corresponding to closed channels, particularly near their thresholds, can be reconstructed from the experimentally accessible submatrix of open channels.The convergence of the proposed method is verified by applying it to data generated from a direct solution of the scattering problem for a known potential, and comparing the reconstructed potential with the original one. Finally, the method is applied to the analysis of S31S_{31}S31​ πNπNπN scattering data.
The problem of determining a scattering potential from scattered field data is addressed within the two integral equation formulation of scattering theory first proposed by Goldberger (1964). The lowest-order approximation within this formulation casts the inverse problem entirely in terms of the Heitler equation which relates the on-shell T matrix directly to the scattering potential. This lowest-order approximation is applied to the case of non-relativistic quantum mechanical inverse scattering from one-dimensional local potentials and is found to yield an exact solution for the case of a one-dimensional delta function potential.
The present paper generalizes preceding papers of the author [1]–[6] and opens a cycle of works concerning the general quantum-mechanical inverse scattering problem (for a given partial channel) and its solution in analytic form in a nontraditional setting: the scattering wave function in the momentum representation is recovered immediately from the scattering data. The solution we find for the problem in question is given by a family of phase-equivalent wave functions with explicitly controlled functional-analytic arbitrariness. This arbitrariness corresponds to the arbitrariness in the choice of phase-equivalent nonlocal potentials and demands including some additional concepts of quantum measurement theory to a complete description of quantum-mechanical systems. We discuss several related problems of quantum-mechanical probability distributions and their statistical interpretation. The mathematical apparatus of the approach thus developed is illustrated by the example of describing the nucleon-nucleon interaction in uncoupled partial channels. In conclusion we announce the topics of forthcoming works of the series.
In this article we continue our investigations of one particle quantum scattering theory for Schrödinger operators on a set of connected (idealized one-dimensional) wires forming a graph with an arbitrary number of open ends. The Hamiltonian is given as minus the Laplace operator with suitable linear boundary conditions at the vertices (the local Kirchhoff law). In “Kirchhoff's rule for quantum wires” [J. Phys. A: Math. Gen. 32, 595–630 (1999)] we provided an explicit algebraic expression for the resulting (on-shell) S-matrix in terms of the boundary conditions and the lengths of the internal lines and we also proved its unitarity. Here we address the inverse problem in the simplest context with one vertex only but with an arbitrary number of open ends. We provide an explicit formula for the boundary conditions in terms of the S-matrix at a fixed, prescribed energy. We show that any unitary n × n matrix may be realized as the S-matrix at a given energy by choosing appropriate (unique) boundary conditions. This might possibly be used for the design of elementary gates in quantum computing. As an illustration we calculate the boundary conditions associated to the unitary operators of some elementary gates for quantum computers and raise the issue whether in general the unitary operators associated to quantum gates should rather be viewed as scattering operators instead of time evolution operators for a given time associated to a quantum mechanical Hamiltonian. We also suggest an approach by which the S-matrix in our context may be obtained from “scattering experiments”, another aspect of the inverse problem. Finally we extend our previous discussion, how our approach is related to von Neumann's theory of selfadjoint extensions.
Given two Hamiltonians H0,H1=H0+VH_0,H_1=H_0+VH0​,H1​=H0​+V corresponding respectively to a free particle Hamiltonian and the free particle Hamiltonian plus its interaction energy with a scattering centre, we can compute the S-matrix, i.e., the matrix that determines the quantum mechanical amplitude of the process involving a particle coming from a free particle state at an infinite distance from the scattering centre at time t=−∞t=-\inftyt=−∞, entering into an initial scattered state, interacting with the scattering centre, entering into a final scattered state, and then escaping away to infinity at time t→∞t\rightarrow\inftyt→∞ to a final free particle state. Such amplitudes were first studied by Lippmann and Schwinger (Steven Weinberg, The Quantum Theory of Fields, Vol. 1, W.O. Amrein, Hilbert Space Methods in Quantum Mechanics) and then analyzed in a mathematically precise way by [T. Kato], [W. Amrein], and [K.B. Sinha]. In this paper, we first derive the Lippmann-Schwinger equation for the input and output scattered states in terms of the input and output free particle states using intuitive arguments originally due to Lippmann and Schwinger, and then, using rigorous operator theoretic arguments revolving around the spectral theorem for unbounded self-adjoint operators in a Hilbert space, we obtain an elegant formula for the scattering matrix elements in terms of the interaction potential and the free particle Hamiltonian between two free particle states corresponding to the same measure. We then discuss a computationally more efficient method based on the Dyson series expansion (Steven Weinberg, The Quantum Theory of Fields, Vol. 1), which works even in the case when the interaction Hamiltonian is time varying, in contrast to the Lippmann-Schwinger method, for calculating the S-matrix. We assume that the scattering potential, which can be time varying, can be controlled by incorporating control parameters and explain how the resulting unitary S-matrix can be made to approximate a given unitary gate in infinite-dimensional Hilbert space by optimizing the Hilbert-Schmidt/Frobenius norm distance between the given gate and the S-matrix gate, in which the latter is approximated by a truncated Dyson series. Of course, the S-matrix given by the Dyson series is a Taylor series functional expansion in the interaction potential and hence is a highly nonlinear function of the control parameters, and hence we propose numerical methods for such optimization. We then observe that if the control parameters are taken as random variables with joint moments, then the TPCP map involving statistical averaging of the adjoint action of the S-matrix acting on an initial mixed state w.r.t. the probability distribution of the control parameters will result in a final mixed state that is linear in the statistical moments of the control parameters, and hence optimization w.r.t. these moments can be easily carried out. Our method of TPCP map design is to optimize the sum of Frobenius distance squares between a sequence of desired output mixed states and the S-matrix based formula for the TPCP map acting on the corresponding input mixed states w.r.t. the statistical moments of the control parameters. We then generalize this method to include the case of TPCP maps obtained from the Hudson-Parthasarathy quantum noisy Schrödinger equation. In summary, we simulate the S-matrix gate for independent realizations of the control parameters having a probability distribution defined by their joint moments obtained by the above optimization procedure and apply the adjoint action of these independent S-matrices on a given input state, and take their ensemble average to obtain a good approximation to the action of the desired TPCP map on any given input state. This suggests that TPCP maps of arbitrarily large size can be realized using scattering experiments in the laboratory.
Abstract We make a relativistic extension of the one-dimensional J-matrix method of scattering. The relativistic potential matrix is a combination of vector, scalar, and pseudo-scalar components. These are non-singular short-range potential functions (not necessarily analytic like a square well) such that they are well represented by their matrix elements in a finite subset of a square integrable basis set. This set is chosen to support a tridiagonal symmetric matrix representation for the free Dirac operator. We derive the transmission and reflection coefficients. This work will be followed by another where we apply the theory to obtain scattering information for different potential coupling modes including those with evidence of Klein tunneling and of supercritical resonances.
Low-energy positron scattering from singlet and triplet metastable He(1s2s1s2s1s2s) is investigated using the RRR-matrix propagation method in hyperspherical coordinates. Elastic and positronium-formation cross sections are reported, and near-threshold resonance structures are analyzed in terms of eigenphase sums and the time-delay matrix. For the triplet target, the calculated cross sections are in agreement with available convergent-close-coupling results and display the expected threshold behavior. Beyond the known SSS-wave features, Feshbach resonance series in higher partial waves extending up to highly excited atomic thresholds are systematically identified. The time-delay matrix further uncovers hidden resonances that produce obvious structures in the positronium-formation cross sections.