Sigma set scattering equations in nuclear reaction theory

International Nuclear Information System (INIS)

Kowalski, K.L.; Picklesimer, A.

1982-01-01

The practical applications of partially summed versions of the Rosenberg equations involving only special subsets (sigma sets) of the physical amplitudes are investigated with special attention to the Pauli principle. The requisite properties of the transformations from the pair labels to the set of partitions labeling the sigma set of asymptotic channels are established. New, well-defined, scattering integral equations for the antisymmetrized transition operators are found which possess much less coupling among the physically distinct channels than hitherto expected for equations with kernels of equal complexity. In several cases of physical interest in nuclear physics, a single connected-kernel equation is obtained for the relevant antisymmetrized elastic scattering amplitude

Transformation properties of the integrable evolution equations

International Nuclear Information System (INIS)

Konopelchenko, B.G.

1981-01-01

Group-theoretical properties of partial differential equations integrable by the inverse scattering transform method are discussed. It is shown that nonlinear transformations typical to integrable equations (symmetry groups, Baecklund-transformations) and these equations themselves are contained in a certain universal nonlinear transformation group. (orig.)

Compactness and robustness: Applications in the solution of integral equations for chemical kinetics and electromagnetic scattering

Science.gov (United States)

Zhou, Yajun

This thesis employs the topological concept of compactness to deduce robust solutions to two integral equations arising from chemistry and physics: the inverse Laplace problem in chemical kinetics and the vector wave scattering problem in dielectric optics. The inverse Laplace problem occurs in the quantitative understanding of biological processes that exhibit complex kinetic behavior: different subpopulations of transition events from the "reactant" state to the "product" state follow distinct reaction rate constants, which results in a weighted superposition of exponential decay modes. Reconstruction of the rate constant distribution from kinetic data is often critical for mechanistic understandings of chemical reactions related to biological macromolecules. We devise a "phase function approach" to recover the probability distribution of rate constants from decay data in the time domain. The robustness (numerical stability) of this reconstruction algorithm builds upon the continuity of the transformations connecting the relevant function spaces that are compact metric spaces. The robust "phase function approach" not only is useful for the analysis of heterogeneous subpopulations of exponential decays within a single transition step, but also is generalizable to the kinetic analysis of complex chemical reactions that involve multiple intermediate steps. A quantitative characterization of the light scattering is central to many meteoro-logical, optical, and medical applications. We give a rigorous treatment to electromagnetic scattering on arbitrarily shaped dielectric media via the Born equation: an integral equation with a strongly singular convolution kernel that corresponds to a non-compact Green operator. By constructing a quadratic polynomial of the Green operator that cancels out the kernel singularity and satisfies the compactness criterion, we reveal the universality of a real resonance mode in dielectric optics. Meanwhile, exploiting the properties of

Integral equation methods for electromagnetics

CERN Document Server

Volakis, John

2012-01-01

This text/reference is a detailed look at the development and use of integral equation methods for electromagnetic analysis, specifically for antennas and radar scattering. Developers and practitioners will appreciate the broad-based approach to understanding and utilizing integral equation methods and the unique coverage of historical developments that led to the current state-of-the-art. In contrast to existing books, Integral Equation Methods for Electromagnetics lays the groundwork in the initial chapters so students and basic users can solve simple problems and work their way up to the mo

High-order integral equation methods for problems of scattering by bumps and cavities on half-planes.

Science.gov (United States)

Pérez-Arancibia, Carlos; Bruno, Oscar P

2014-08-01

This paper presents high-order integral equation methods for the evaluation of electromagnetic wave scattering by dielectric bumps and dielectric cavities on perfectly conducting or dielectric half-planes. In detail, the algorithms introduced in this paper apply to eight classical scattering problems, namely, scattering by a dielectric bump on a perfectly conducting or a dielectric half-plane, and scattering by a filled, overfilled, or void dielectric cavity on a perfectly conducting or a dielectric half-plane. In all cases field representations based on single-layer potentials for appropriately chosen Green functions are used. The numerical far fields and near fields exhibit excellent convergence as discretizations are refined-even at and around points where singular fields and infinite currents exist.

Roy–Steiner equations for πN scattering

Directory of Open Access Journals (Sweden)

Ruiz de Elvira J.

2014-06-01

Full Text Available In this talk, we present a coupled system of integral equations for the πN → πN (s-channel and ππ → N̅N (t-channel lowest partial waves, derived from Roy–Steiner equations for pion–nucleon scattering. After giving a brief overview of this system of equations, we present the solution of the t-channel sub-problem by means of Muskhelishvili–Omnès techniques, and solve the s-channel sub-problem after finding a set of phase shifts and subthreshold parameters which satisfy the Roy–Steiner equations.

A boundary integral equation method using auxiliary interior surface approach for acoustic radiation and scattering in two dimensions.

Science.gov (United States)

Yang, S A

2002-10-01

This paper presents an effective solution method for predicting acoustic radiation and scattering fields in two dimensions. The difficulty of the fictitious characteristic frequency is overcome by incorporating an auxiliary interior surface that satisfies certain boundary condition into the body surface. This process gives rise to a set of uniquely solvable boundary integral equations. Distributing monopoles with unknown strengths over the body and interior surfaces yields the simple source formulation. The modified boundary integral equations are further transformed to ordinary ones that contain nonsingular kernels only. This implementation allows direct application of standard quadrature formulas over the entire integration domain; that is, the collocation points are exactly the positions at which the integration points are located. Selecting the interior surface is an easy task. Moreover, only a few corresponding interior nodal points are sufficient for the computation. Numerical calculations consist of the acoustic radiation and scattering by acoustically hard elliptic and rectangular cylinders. Comparisons with analytical solutions are made. Numerical results demonstrate the efficiency and accuracy of the current solution method.

A calderón-preconditioned single source combined field integral equation for analyzing scattering from homogeneous penetrable objects

KAUST Repository

Valdés, Felipe

2011-06-01

A new regularized single source equation for analyzing scattering from homogeneous penetrable objects is presented. The proposed equation is a linear combination of a Calderón-preconditioned single source electric field integral equation and a single source magnetic field integral equation. The equation is immune to low-frequency and dense-mesh breakdown, and free from spurious resonances. Unlike dual source formulations, this equation involves operator products that cannot be discretized using standard procedures for discretizing standalone electric, magnetic, and combined field operators. Instead, the single source equation proposed here is discretized using a recently developed technique that achieves a well-conditioned mapping from div- to curl-conforming function spaces, thereby fully respecting the space mapping properties of the operators involved, and guaranteeing accuracy and stability. Numerical results show that the proposed equation and discretization technique give rise to rapidly convergent solutions. They also validate the equation\\'s resonant free character. © 2006 IEEE.

Transient analysis of electromagnetic wave interactions on high-contrast scatterers using volume electric field integral equation

KAUST Repository

Sayed, Sadeed Bin

2014-07-01

A marching on-in-time (MOT)-based time domain volume electric field integral equation (TD-VEFIE) solver is proposed for accurate and stable analysis of electromagnetic wave interactions on high-contrast scatterers. The stability is achieved using band-limited but two-sided (non-causal) temporal interpolation functions and an extrapolation scheme to cast the time marching into a causal form. The extrapolation scheme is designed to be highly accurate for oscillating and exponentially decaying fields, hence it accurately captures the physical behavior of the resonant modes that are excited inside the dielectric scatterer. Numerical results demonstrate that the resulting MOT scheme maintains its stability as the number of resonant modes increases with the contrast of the scatterer.

On a new series of integrable nonlinear evolution equations

International Nuclear Information System (INIS)

Ichikawa, Y.H.; Wadati, Miki; Konno, Kimiaki; Shimizu, Tohru.

1980-10-01

Recent results of our research are surveyed in this report. The derivative nonlinear Schroedinger equation for the circular polarized Alfven wave admits the spiky soliton solutions for the plane wave boundary condition. The nonlinear equation for complex amplitude associated with the carrier wave is shown to be a generalized nonlinear Schroedinger equation, having the ordinary cubic nonlinear term and the derivative of cubic nonlinear term. A generalized scheme of the inverse scattering transformation has confirmed that superposition of the A-K-N-S scheme and the K-N scheme for the component equations valids for the generalized nonlinear Schroedinger equation. Then, two types of new integrable nonlinear evolution equation have been derived from our scheme of the inverse scattering transformation. One is the type of nonlinear Schroedinger equation, while the other is the type of Korteweg-de Vries equation. Brief discussions are presented for physical phenomena, which could be accounted by the second type of the new integrable nonlinear evolution equation. Lastly, the stationary solitary wave solutions have been constructed for the integrable nonlinear evolution equation of the second type. These solutions have peculiar structure that they are singular and discrete. It is a new challenge to construct singular potentials by the inverse scattering transformation. (author)

A volume integral equation solver for quantum-corrected transient analysis of scattering from plasmonic nanostructures

KAUST Repository

Sayed, Sadeed Bin; Uysal, Ismail Enes; Bagci, Hakan; Ulku, H. Arda

2018-01-01

Quantum tunneling is observed between two nanostructures that are separated by a sub-nanometer gap. Electrons “jumping” from one structure to another create an additional current path. An auxiliary tunnel is introduced between the two structures as a support for this so that a classical electromagnetic solver can account for the effects of quantum tunneling. The dispersive permittivity of the tunnel is represented by a Drude model, whose parameters are obtained from the electron tunneling probability. The transient scattering from the connected nanostructures (i.e., nanostructures plus auxiliary tunnel) is analyzed using a time domain volume integral equation solver. Numerical results demonstrating the effect of quantum tunneling on the scattered fields are provided.

A Highly Stable Marching-on-in-Time Volume Integral Equation Solver for Analyzing Transient Wave Interactions on High-Contrast Scatterers

KAUST Repository

Bagci, Hakan

2014-01-06

Time domain integral equation (TDIE) solvers represent an attractive alternative to finite difference (FDTD) and finite element (FEM) schemes for analyzing transient electromagnetic interactions on composite scatterers. Current induced on a scatterer, in response to a transient incident field, generates a scattered field. First, the scattered field is expressed as a spatio-temporal convolution of the current and the Green function of the background medium. Then, a TDIE is obtained by enforcing boundary conditions and/or fundamental field relations. TDIEs are often solved for the unknown current using marching on-in-time (MOT) schemes. MOT-TDIE solvers expand the current using local spatio-temporal basis functions. Inserting this expansion into the TDIE and testing the resulting equation in space and time yields a lower triangular system of equations (termed MOT system), which can be solved by marching in time for the coefficients of the current expansion. Stability of the MOT scheme often depends on how accurately the spatio-temporal convolution of the current and the Green function is discretized. In this work, band-limited prolate-based interpolation functions are used as temporal bases in expanding the current and discretizing the spatio-temporal convolution. Unfortunately, these functions are two sided, i.e., they require ”future” current samples for interpolation, resulting in a non-causal MOT system. To alleviate the effect of non-causality and restore the ability to march in time, an extrapolation scheme can be used to estimate the future values of the currents from their past values. Here, an accurate, stable and band-limited extrapolation scheme is developed for this purpose. This extrapolation scheme uses complex exponents, rather than commonly used harmonics, so that propagating and decaying mode fields inside the dielectric scatterers are accurately modeled. The resulting MOT scheme is applied to solving the time domain volume integral equation (VIE

A Stable Marching on-in-time Scheme for Solving the Time Domain Electric Field Volume Integral Equation on High-contrast Scatterers

KAUST Repository

Sayed, Sadeed Bin

2015-05-05

A time domain electric field volume integral equation (TD-EFVIE) solver is proposed for characterizing transient electromagnetic wave interactions on high-contrast dielectric scatterers. The TD-EFVIE is discretized using the Schaubert- Wilton-Glisson (SWG) and approximate prolate spherical wave (APSW) functions in space and time, respectively. The resulting system of equations can not be solved by a straightforward application of the marching on-in-time (MOT) scheme since the two-sided APSW interpolation functions require the knowledge of unknown “future” field samples during time marching. Causality of the MOT scheme is restored using an extrapolation technique that predicts the future samples from known “past” ones. Unlike the extrapolation techniques developed for MOT schemes that are used in solving time domain surface integral equations, this scheme trains the extrapolation coefficients using samples of exponentials with exponents on the complex frequency plane. This increases the stability of the MOT-TD-EFVIE solver significantly, since the temporal behavior of decaying and oscillating electromagnetic modes induced inside the scatterers is very accurately taken into account by this new extrapolation scheme. Numerical results demonstrate that the proposed MOT solver maintains its stability even when applied to analyzing wave interactions on high-contrast scatterers.

A Stable Marching on-in-time Scheme for Solving the Time Domain Electric Field Volume Integral Equation on High-contrast Scatterers

KAUST Repository

Sayed, Sadeed Bin; Ulku, Huseyin; Bagci, Hakan

2015-01-01

A time domain electric field volume integral equation (TD-EFVIE) solver is proposed for characterizing transient electromagnetic wave interactions on high-contrast dielectric scatterers. The TD-EFVIE is discretized using the Schaubert- Wilton-Glisson (SWG) and approximate prolate spherical wave (APSW) functions in space and time, respectively. The resulting system of equations can not be solved by a straightforward application of the marching on-in-time (MOT) scheme since the two-sided APSW interpolation functions require the knowledge of unknown “future” field samples during time marching. Causality of the MOT scheme is restored using an extrapolation technique that predicts the future samples from known “past” ones. Unlike the extrapolation techniques developed for MOT schemes that are used in solving time domain surface integral equations, this scheme trains the extrapolation coefficients using samples of exponentials with exponents on the complex frequency plane. This increases the stability of the MOT-TD-EFVIE solver significantly, since the temporal behavior of decaying and oscillating electromagnetic modes induced inside the scatterers is very accurately taken into account by this new extrapolation scheme. Numerical results demonstrate that the proposed MOT solver maintains its stability even when applied to analyzing wave interactions on high-contrast scatterers.

Inverse scattering solution of the Chew-Low equation

International Nuclear Information System (INIS)

Nakano, K.

1985-01-01

Techniques for solving the inverse scattering problem are applied to the Chew-Low equation to obtain the nucleon form factor directly from the experimental phase shifts. A new dispersion relation is derived for the P 11 wave because of its sign-changing phase shift. A self-consistent solution for each channel is obtained, but the universality of form factor is not confirmed. Also, an iterative procedure based on Omnes' method is developed in order to solve coupled-channel, singular integral equations. (orig.)

Relativistic two-and three-particle scattering equations using instant and light-front dynamics

International Nuclear Information System (INIS)

Adhikari, S.K.; Tomio, L.; Frederico, T.

1992-01-01

Starting from the Bethe-Salpeter equation for two particles in the ladder approximation and integrating over the time component of momentum we derive three dimensional scattering integral equations satisfying constraints of unitarity and relativity, both employing the light-front and instant-form variables. The equations we arrive at are those first derived by Weinberg and by Blankenbecler and Sugar, and are shown to be related by a transformation of variables. Hence we show how to perform and relate identical dynamical calculation using these two equations. We extends this procedure to the case of three particles interacting via two-particle separable potentials. Using light-front and instant form variables we suggest a couple of three dimensional three-particle scattering equations satisfying constraints of two and three-particle unitarity and relativity. The three-particle light-front equation is shown to be approximately related by a transformation of variables to one of the instant-form three-particle equations. (author)

Solution of neutron slowing down equation including multiple inelastic scattering

International Nuclear Information System (INIS)

El-Wakil, S.A.; Saad, A.E.

1977-01-01

The present work is devoted the presentation of an analytical method for the calculation of elastically and inelastically slowed down neutrons in an infinite non absorbing homogeneous medium. On the basis of the Central limit theory (CLT) and the integral transform technique the slowing down equation including inelastic scattering in terms of the Green function of elastic scattering is solved. The Green function is decomposed according to the number of collisions. A formula for the flux at any lethargy O (u) after any number of collisions is derived. An equation for the asymptotic flux is also obtained

Low-frequency scaling of the standard and mixed magnetic field and Müller integral equations

KAUST Repository

Bogaert, Ignace; Cools, Kristof; Andriulli, Francesco P.; Bagci, Hakan

2014-01-01

The standard and mixed discretizations for the magnetic field integral equation (MFIE) and the Müller integral equation (MUIE) are investigated in the context of low-frequency (LF) scattering problems involving simply connected scatterers

Inverse scattering transform for the KPI equation on the background of a one-line soliton*Inverse scattering transform for the KPI equation on the background of a one-line soliton

Science.gov (United States)

Fokas, A. S.; Pogrebkov, A. K.

2003-03-01

We study the initial value problem of the Kadomtsev-Petviashvili I (KPI) equation with initial data u(x1,x2,0) = u1(x1)+u2(x1,x2), where u1(x1) is the one-soliton solution of the Korteweg-de Vries equation evaluated at zero time and u2(x1,x2) decays sufficiently rapidly on the (x1,x2)-plane. This involves the analysis of the nonstationary Schrödinger equation (with time replaced by x2) with potential u(x1,x2,0). We introduce an appropriate sectionally analytic eigenfunction in the complex k-plane where k is the spectral parameter. This eigenfunction has the novelty that in addition to the usual jump across the real k-axis, it also has a jump across a segment of the imaginary k-axis. We show that this eigenfunction can be reconstructed through a linear integral equation uniquely defined in terms of appropriate scattering data. In turn, these scattering data are uniquely constructed in terms of u1(x1) and u2(x1,x2). This result implies that the solution of the KPI equation can be obtained through the above linear integral equation where the scattering data have a simple t-dependence.

A boundary integral equation for boundary element applications in multigroup neutron diffusion theory

International Nuclear Information System (INIS)

Ozgener, B.

1998-01-01

A boundary integral equation (BIE) is developed for the application of the boundary element method to the multigroup neutron diffusion equations. The developed BIE contains no explicit scattering term; the scattering effects are taken into account by redefining the unknowns. Boundary elements of the linear and constant variety are utilised for validation of the developed boundary integral formulation

Light propagation in finite-sized photonic crystals: multiple scattering using an electric field integral equation

DEFF Research Database (Denmark)

Kristensen, Philip Trøst; Lodahl, Peter; Mørk, Jesper

2010-01-01

We present an accurate, stable, and efficient solution to the Lippmann–Schwinger equation for electromagnetic scattering in two dimensions. The method is well suited for multiple scattering problems and may be applied to problems with scatterers of arbitrary shape or non-homogenous background mat...

Two loop integrals and QCD scattering

International Nuclear Information System (INIS)

Anastasiou, C.

2001-04-01

We present the techniques for the calculation of one- and two-loop integrals contributing to the virtual corrections to 2→2 scattering of massless particles. First, tensor integrals are related to scalar integrals with extra powers of propagators and higher dimension using the Schwinger representation. Integration By Parts and Lorentz Invariance recurrence relations reduce the number of independent scalar integrals to a set of master integrals for which their expansion in ε = 2 - D/2 is calculated using a combination of Feynman parameters, the Negative Dimension Integration Method, the Differential Equations Method, and Mellin-Barnes integral representations. The two-loop matrix-elements for light-quark scattering are calculated in Conventional Dimensional Regularisation by direct evaluation of the Feynman diagrams. The ultraviolet divergences are removed by renormalising with the MS-bar scheme. Finally, the infrared singular behavior is shown to be in agreement with the one anticipated by the application of Catani's formalism for the infrared divergences of generic QCD two-loop amplitudes. (author)

Invariant imbedding equations for linear scattering problems

International Nuclear Information System (INIS)

Apresyan, L.

1988-01-01

A general form of the invariant imbedding equations is investigated for the linear problem of scattering by a bounded scattering volume. The conditions for the derivability of such equations are described. It is noted that the possibility of the explicit representation of these equations for a sphere and for a layer involves the separation of variables in the unperturbed wave equation

Study of scattering cross section of a plasma column using Green's function volume integral equation method

Science.gov (United States)

Soltanmoradi, Elmira; Shokri, Babak

2017-05-01

In this article, the electromagnetic wave scattering from plasma columns with inhomogeneous electron density distribution is studied by the Green's function volume integral equation method. Due to the ready production of such plasmas in the laboratories and their practical application in various technological fields, this study tries to find the effects of plasma parameters such as the electron density, radius, and pressure on the scattering cross-section of a plasma column. Moreover, the incident wave frequency influence of the scattering pattern is demonstrated. Furthermore, the scattering cross-section of a plasma column with an inhomogeneous collision frequency profile is calculated and the effect of this inhomogeneity is discussed first in this article. These results are especially used to determine the appropriate conditions for radar cross-section reduction purposes. It is shown that the radar cross-section of a plasma column reduces more for a larger collision frequency, for a relatively lower plasma frequency, and also for a smaller radius. Furthermore, it is found that the effect of the electron density on the scattering cross-section is more obvious in comparison with the effect of other plasma parameters. Also, the plasma column with homogenous collision frequency can be used as a better shielding in contrast to its inhomogeneous counterpart.

Analysis of electromagnetic wave interactions on nonlinear scatterers using time domain volume integral equations

KAUST Repository

Ulku, Huseyin Arda

2014-07-06

Effects of material nonlinearities on electromagnetic field interactions become dominant as field amplitudes increase. A typical example is observed in plasmonics, where highly localized fields “activate” Kerr nonlinearities. Naturally, time domain solvers are the method of choice when it comes simulating these nonlinear effects. Oftentimes, finite difference time domain (FDTD) method is used for this purpose. This is simply due to the fact that explicitness of the FDTD renders the implementation easier and the material nonlinearity can be easily accounted for using an auxiliary differential equation (J.H. Green and A. Taflove, Opt. Express, 14(18), 8305-8310, 2006). On the other hand, explicit marching on-in-time (MOT)-based time domain integral equation (TDIE) solvers have never been used for the same purpose even though they offer several advantages over FDTD (E. Michielssen, et al., ECCOMAS CFD, The Netherlands, Sep. 5-8, 2006). This is because explicit MOT solvers have never been stabilized until not so long ago. Recently an explicit but stable MOT scheme has been proposed for solving the time domain surface magnetic field integral equation (H.A. Ulku, et al., IEEE Trans. Antennas Propag., 61(8), 4120-4131, 2013) and later it has been extended for the time domain volume electric field integral equation (TDVEFIE) (S. B. Sayed, et al., Pr. Electromagn. Res. S., 378, Stockholm, 2013). This explicit MOT scheme uses predictor-corrector updates together with successive over relaxation during time marching to stabilize the solution even when time step is as large as in the implicit counterpart. In this work, an explicit MOT-TDVEFIE solver is proposed for analyzing electromagnetic wave interactions on scatterers exhibiting Kerr nonlinearity. Nonlinearity is accounted for using the constitutive relation between the electric field intensity and flux density. Then, this relation and the TDVEFIE are discretized together by expanding the intensity and flux - sing half

Adaptive integral equation methods in transport theory

International Nuclear Information System (INIS)

Kelley, C.T.

1992-01-01

In this paper, an adaptive multilevel algorithm for integral equations is described that has been developed with the Chandrasekhar H equation and its generalizations in mind. The algorithm maintains good performance when the Frechet derivative of the nonlinear map is singular at the solution, as happens in radiative transfer with conservative scattering and in critical neutron transport. Numerical examples that demonstrate the algorithm's effectiveness are presented

Radiative transfer equation accounting for rotational Raman scattering and its solution by the discrete-ordinates method

International Nuclear Information System (INIS)

Rozanov, Vladimir V.; Vountas, Marco

2014-01-01

Rotational Raman scattering of solar light in Earth's atmosphere leads to the filling-in of Fraunhofer and telluric lines observed in the reflected spectrum. The phenomenological derivation of the inelastic radiative transfer equation including rotational Raman scattering is presented. The different forms of the approximate radiative transfer equation with first-order rotational Raman scattering terms are obtained employing the Cabannes, Rayleigh, and Cabannes–Rayleigh scattering models. The solution of these equations is considered in the framework of the discrete-ordinates method using rigorous and approximate approaches to derive particular integrals. An alternative forward-adjoint technique is suggested as well. A detailed description of the model including the exact spectral matching and a binning scheme that significantly speeds up the calculations is given. The considered solution techniques are implemented in the radiative transfer software package SCIATRAN and a specified benchmark setup is presented to enable readers to compare with own results transparently. -- Highlights: • We derived the radiative transfer equation accounting for rotational Raman scattering. • Different approximate radiative transfer approaches with first order scattering were used. • Rigorous and approximate approaches are shown to derive particular integrals. • An alternative forward-adjoint technique is suggested as well. • An additional spectral binning scheme which speeds up the calculations is presented

Relativistic wave equations and compton scattering

International Nuclear Information System (INIS)

Sutanto, S.H.; Robson, B.A.

1998-01-01

Full text: Recently an eight-component relativistic wave equation for spin-1/2 particles was proposed.This equation was obtained from a four-component spin-1/2 wave equation (the KG1/2 equation), which contains second-order derivatives in both space and time, by a procedure involving a linearisation of the time derivative analogous to that introduced by Feshbach and Villars for the Klein-Gordon equation. This new eight-component equation gives the same bound-state energy eigenvalue spectra for hydrogenic atoms as the Dirac equation but has been shown to predict different radiative transition probabilities for the fine structure of both the Balmer and Lyman a-lines. Since it has been shown that the new theory does not always give the same results as the Dirac theory, it is important to consider the validity of the new equation in the case of other physical problems. One of the early crucial tests of the Dirac theory was its application to the scattering of a photon by a free electron: the so-called Compton scattering problem. In this paper we apply the new theory to the calculation of Compton scattering to order e 2 . It will be shown that in spite of the considerable difference in the structure of the new theory and that of Dirac the cross section is given by the Klein-Nishina formula

On integrability conditions of the equations of nonsymmetrical chiral field on SO(4)

International Nuclear Information System (INIS)

Tskhakaya, D.D.

1990-01-01

Possibility of integrating the equations of nonsymmetrical chiral field on SO(4) by means of the inverse scattering method is investigated. Maximal number of the motion integrals is found for the corresponding system of ordinary differential equations

Graphical analyses of connected-kernel scattering equations

International Nuclear Information System (INIS)

Picklesimer, A.

1983-01-01

Simple graphical techniques are employed to obtain a new (simultaneous) derivation of a large class of connected-kernel scattering equations. This class includes the Rosenberg, Bencze-Redish-Sloan, and connected-kernel multiple scattering equations as well as a host of generalizations of these and other equations. The basic result is the application of graphical methods to the derivation of interaction-set equations. This yields a new, simplified form for some members of the class and elucidates the general structural features of the entire class

Path integral approach to electron scattering in classical electromagnetic potential

International Nuclear Information System (INIS)

Xu Chuang; Feng Feng; Li Ying-Jun

2016-01-01

As is known to all, the electron scattering in classical electromagnetic potential is one of the most widespread applications of quantum theory. Nevertheless, many discussions about electron scattering are based upon single-particle Schrodinger equation or Dirac equation in quantum mechanics rather than the method of quantum field theory. In this paper, by using the path integral approach of quantum field theory, we perturbatively evaluate the scattering amplitude up to the second order for the electron scattering by the classical electromagnetic potential. The results we derive are convenient to apply to all sorts of potential forms. Furthermore, by means of the obtained results, we give explicit calculations for the one-dimensional electric potential. (paper)

Graphical analyses of connected-kernel scattering equations

International Nuclear Information System (INIS)

Picklesimer, A.

1982-10-01

Simple graphical techniques are employed to obtain a new (simultaneous) derivation of a large class of connected-kernel scattering equations. This class includes the Rosenberg, Bencze-Redish-Sloan, and connected-kernel multiple scattering equations as well as a host of generalizations of these and other equations. The graphical method also leads to a new, simplified form for some members of the class and elucidates the general structural features of the entire class

Solution of the Helmholtz-Poincare Wave Equation using the coupled boundary integral equations and optimal surface eigenfunctions

International Nuclear Information System (INIS)

Werby, M.F.; Broadhead, M.K.; Strayer, M.R.; Bottcher, C.

1992-01-01

The Helmholtz-Poincarf Wave Equation (H-PWE) arises in many areas of classical wave scattering theory. In particular it can be found for the cases of acoustical scattering from submerged bounded objects and electromagnetic scattering from objects. The extended boundary integral equations (EBIE) method is derived from considering both the exterior and interior solutions of the H-PWECs. This coupled set of expressions has the advantage of not only offering a prescription for obtaining a solution for the exterior scattering problem, but it also obviates the problem of irregular values corresponding to fictitious interior eigenvalues. Once the coupled equations are derived, they can be obtained in matrix form by expanding all relevant terms in partial wave expansions, including a bi-orthogonal expansion of the Green's function. However some freedom in the choice of the surface expansion is available since the unknown surface quantities may be expanded in a variety of ways so long as closure is obtained. Out of many possible choices, we develop an optimal method to obtain such expansions which is based on the optimum eigenfunctions related to the surface of the object. In effect, we convert part of the problem (that associated with the Fredholms integral equation of the first kind) an eigenvalue problem of a related Hermitian operator. The methodology will be explained in detail and examples will be presented

High-order boundary integral equation solution of high frequency wave scattering from obstacles in an unbounded linearly stratified medium

Science.gov (United States)

Barnett, Alex H.; Nelson, Bradley J.; Mahoney, J. Matthew

2015-09-01

We apply boundary integral equations for the first time to the two-dimensional scattering of time-harmonic waves from a smooth obstacle embedded in a continuously-graded unbounded medium. In the case we solve, the square of the wavenumber (refractive index) varies linearly in one coordinate, i.e. (Δ + E +x2) u (x1 ,x2) = 0 where E is a constant; this models quantum particles of fixed energy in a uniform gravitational field, and has broader applications to stratified media in acoustics, optics and seismology. We evaluate the fundamental solution efficiently with exponential accuracy via numerical saddle-point integration, using the truncated trapezoid rule with typically 102 nodes, with an effort that is independent of the frequency parameter E. By combining with a high-order Nyström quadrature, we are able to solve the scattering from obstacles 50 wavelengths across to 11 digits of accuracy in under a minute on a desktop or laptop.

Exact solutions to the Boltzmann equation by mapping the scattering integral into a differential operator

International Nuclear Information System (INIS)

Zabadal, Jorge; Borges, Volnei; Van der Laan, Flavio T.; Santos, Marcio G.

2015-01-01

This work presents a new analytical method for solving the Boltzmann equation. In this formulation, a linear differential operator is applied over the Boltzmann model, in order to produce a partial differential equation in which the scattering term is absent. This auxiliary equation is solved via reduction of order. The exact solution obtained is employed to define a precursor for the buildup factor. (author)

Exact solutions to the Boltzmann equation by mapping the scattering integral into a differential operator

Energy Technology Data Exchange (ETDEWEB)

Zabadal, Jorge; Borges, Volnei; Van der Laan, Flavio T., E-mail: jorge.zabadal@ufrgs.br, E-mail: borges@ufrgs.br, E-mail: ftvdl@ufrgs.br [Universidade Federal do Rio Grande do Sul (UFRGS), Porto Alegre, RS (Brazil). Departamento de Engenharia Mecanica. Grupo de Pesquisas Radiologicas; Ribeiro, Vinicius G., E-mail: vinicius_ribeiro@uniritter.edu.br [Centro Universitario Ritter dos Reis (UNIRITTER), Porto Alegre, RS (Brazil); Santos, Marcio G., E-mail: phd.marcio@gmail.com [Universidade Federal do Rio Grande do Sul (UFRGS), Tramandai, RS (Brazil). Departamento Interdisciplinar do Campus Litoral Norte

2015-07-01

This work presents a new analytical method for solving the Boltzmann equation. In this formulation, a linear differential operator is applied over the Boltzmann model, in order to produce a partial differential equation in which the scattering term is absent. This auxiliary equation is solved via reduction of order. The exact solution obtained is employed to define a precursor for the buildup factor. (author)

Volume integral equation for electromagnetic scattering: Rigorous derivation and analysis for a set of multilayered particles with piecewise-smooth boundaries in a passive host medium

Science.gov (United States)

Yurkin, Maxim A.; Mishchenko, Michael I.

2018-04-01

We present a general derivation of the frequency-domain volume integral equation (VIE) for the electric field inside a nonmagnetic scattering object from the differential Maxwell equations, transmission boundary conditions, radiation condition at infinity, and locally-finite-energy condition. The derivation applies to an arbitrary spatially finite group of particles made of isotropic materials and embedded in a passive host medium, including those with edges, corners, and intersecting internal interfaces. This is a substantially more general type of scatterer than in all previous derivations. We explicitly treat the strong singularity of the integral kernel, but keep the entire discussion accessible to the applied scattering community. We also consider the known results on the existence and uniqueness of VIE solution and conjecture a general sufficient condition for that. Finally, we discuss an alternative way of deriving the VIE for an arbitrary object by means of a continuous transformation of the everywhere smooth refractive-index function into a discontinuous one. Overall, the paper examines and pushes forward the state-of-the-art understanding of various analytical aspects of the VIE.

Extreme Scale FMM-Accelerated Boundary Integral Equation Solver for Wave Scattering

KAUST Repository

AbdulJabbar, Mustafa Abdulmajeed; Al Farhan, Mohammed; Al-Harthi, Noha A.; Chen, Rui; Yokota, Rio; Bagci, Hakan; Keyes, David E.

2018-01-01

scattering, which uses FMM as a matrix-vector multiplication inside the GMRES iterative method. Our FMM Helmholtz kernels treat nontrivial singular and near-field integration points. We implement highly optimized kernels for both shared and distributed memory

Direct time integration of Maxwell's equations in linear dispersive media with absorption for scattering and propagation of femtosecond electromagnetic pulses

Science.gov (United States)

Joseph, Rose M.; Hagness, Susan C.; Taflove, Allen

1991-01-01

The initial results for femtosecond pulse propagation and scattering interactions for a Lorentz medium obtained by a direct time integration of Maxwell's equations are reported. The computational approach provides reflection coefficients accurate to better than 6 parts in 10,000 over the frequency range of dc to 3 x 10 to the 16th Hz for a single 0.2-fs Gaussian pulse incident upon a Lorentz-medium half-space. New results for Sommerfeld and Brillouin precursors are shown and compared with previous analyses. The present approach is robust and permits 2D and 3D electromagnetic pulse propagation directly from the full-vector Maxwell's equations.

The kinematic algebras from the scattering equations

International Nuclear Information System (INIS)

Monteiro, Ricardo; O’Connell, Donal

2014-01-01

We study kinematic algebras associated to the recently proposed scattering equations, which arise in the description of the scattering of massless particles. In particular, we describe the role that these algebras play in the BCJ duality between colour and kinematics in gauge theory, and its relation to gravity. We find that the scattering equations are a consistency condition for a self-dual-type vertex which is associated to each solution of those equations. We also identify an extension of the anti-self-dual vertex, such that the two vertices are not conjugate in general. Both vertices correspond to the structure constants of Lie algebras. We give a prescription for the use of the generators of these Lie algebras in trivalent graphs that leads to a natural set of BCJ numerators. In particular, we write BCJ numerators for each contribution to the amplitude associated to a solution of the scattering equations. This leads to a decomposition of the determinant of a certain kinematic matrix, which appears naturally in the amplitudes, in terms of trivalent graphs. We also present the kinematic analogues of colour traces, according to these algebras, and the associated decomposition of that determinant

Time-independent integral equation for Maxwell's system. Application of radar cross section computation

International Nuclear Information System (INIS)

Pujols, Agnes

1991-01-01

We prove that the scattering operator for the wave equation in the exterior of an non-homogeneous obstacle exists. Its distribution kernel is represented by a time-dependent boundary integral equation. A space-time integral variational formulation is developed for determining the current induced by the scattering of an electromagnetic wave by an homogeneous object. The discrete approximation of the variational problem using a finite element method in both space and time leads to stable convergent schemes, giving a numerical code for perfectly conducting cylinders. (author) [fr

Integral equations

CERN Document Server

Moiseiwitsch, B L

2005-01-01

Two distinct but related approaches hold the solutions to many mathematical problems--the forms of expression known as differential and integral equations. The method employed by the integral equation approach specifically includes the boundary conditions, which confers a valuable advantage. In addition, the integral equation approach leads naturally to the solution of the problem--under suitable conditions--in the form of an infinite series.Geared toward upper-level undergraduate students, this text focuses chiefly upon linear integral equations. It begins with a straightforward account, acco

Complete integrability of the difference evolution equations

International Nuclear Information System (INIS)

Gerdjikov, V.S.; Ivanov, M.I.; Kulish, P.P.

1980-01-01

The class of exactly solvable nonlinear difference evolution equations (DEE) related to the discrete analog of the one-dimensional Dirac problem L is studied. For this starting from L we construct a special linear non-local operator Λ and obtain the expansions of w and σ 3 deltaw over its eigenfunctions, w being the potential in L. This allows us to obtain compact expressions for the integrals of motion and to prove that these DEE are completely integrable Hamiltonian systems. Moreover, it is shown that there exists a hierarchy of Hamiltonian structures, generated by Λ, and the action-angle variables are explicity calculated. As particular cases the difference analog of the non-linear Schroedinger equation and the modified Korteweg-de-Vries equation are considered. The quantization of these Hamiltonian system through the use of the quantum inverse scattering method is briefly discussed [ru

Splines and their reciprocal-bases in volume-integral equations

International Nuclear Information System (INIS)

Sabbagh, H.A.

1993-01-01

The authors briefly outline the use of higher-order splines and their reciprocal-bases in discretizing the volume-integral equations of electromagnetics. The discretization is carried out by means of the method of moments, in which the expansion functions are the higher-order splines, and the testing functions are the corresponding reciprocal-basis functions. These functions satisfy an orthogonality condition with respect to the spline expansion functions. Thus, the method is not Galerkin, but the structure of the resulting equations is quite regular, nevertheless. The theory is applied to the volume-integral equations for the unknown current density, or unknown electric field, within a scattering body, and to the equations for eddy-current nondestructive evaluation. Numerical techniques for computing the matrix elements are also given

Inverse random source scattering for the Helmholtz equation in inhomogeneous media

Science.gov (United States)

Li, Ming; Chen, Chuchu; Li, Peijun

2018-01-01

This paper is concerned with an inverse random source scattering problem in an inhomogeneous background medium. The wave propagation is modeled by the stochastic Helmholtz equation with the source driven by additive white noise. The goal is to reconstruct the statistical properties of the random source such as the mean and variance from the boundary measurement of the radiated random wave field at multiple frequencies. Both the direct and inverse problems are considered. We show that the direct problem has a unique mild solution by a constructive proof. For the inverse problem, we derive Fredholm integral equations, which connect the boundary measurement of the radiated wave field with the unknown source function. A regularized block Kaczmarz method is developed to solve the ill-posed integral equations. Numerical experiments are included to demonstrate the effectiveness of the proposed method.

Comments on the integrability of the loop-space chiral equations

International Nuclear Information System (INIS)

Gu, C.; Wang, L.L.C.

1980-01-01

A demonstration is given how the ordinary space chiral equations provide the existence conditions for the infinite number of conserved currents and how these currents are related to the so-called inverse-scattering equations, whose integrability is provided by the original chiral equations. Loop-space chiral equations are introduced. The integrability conditions of the non-local currents in two possible different situations are discussed. In the first case, the generating functions are functionals of the loop alone. The integrability conditions are not satisfied and higher order conserved non-local currents do not exist. In the second case, the generating functions are functionals of the loop as well as a parameter the integrability conditions at a restricted point of the parameter are satisfied, however there is an infinite fold of arbitrariness. It indicates that additional guiding principles are needed in addition to the original loop-space chiral equation in order to uniquely determine the infinite conserved non-local currents as functionals of the loop and the parameter

Roy-Steiner equations for πN scattering

Science.gov (United States)

de Elvira, J. Ruiz; Ditsche, C.; Hoferichter, M.; Kubis, B.; Meißner, U.-G.

2015-10-01

In this talk, we briefly review our ongoing collaboration to precisely determine the low-energy πN scattering amplitude by means of Roy-Steiner equations. After giving a brief overview of this system of dispersive equations and their application to πN scattering, we proceed to solve for the lower partial waves of the s-channel (πN → πN) and the t-channel l( {π π to bar NN} right) sub-problems.

A Comprehensive Review of Boundary Integral Formulations of Acoustic Scattering Problems

Directory of Open Access Journals (Sweden)

S.I. Zaman

2000-12-01

Full Text Available This is a review presenting an overview of the developments in boundary integral formulations of the acoustic scattering problems. Generally, the problem is formulated in one of two ways viz. Green’s representation formula, and the Layer-theoretic formulation utilizing either a simple-layer or a double-layer potential. The review presents and expounds the major contributions in this area over the last four decades. The need for a robust and improved formulation of the exterior scattering problem (Neumann or Dirichlet arose due to the fact that the classical formulation failed to yield a unique solution at (acoustic wave-numbers which correspond to eigenvalues (eigenfrequencies of the corresponding interior scattering problem. Moreover, this correlation becomes more pronounced as the wave-numbers become larger i.e. as the (acoustic frequency increases. The robust integral formulations which are discussed here yield Fredholms integral equations of the second kind which are more amenable to computation than the first kind. However, the integral equation involves a hypersingular kernel which creates ill-conditioning in the final matrix representation. This is circumvented by a regularisation technique. An extensive useful list of references is also presented here for researchers in this area.

Two-body Dirac equations for nucleon-nucleon scattering

International Nuclear Information System (INIS)

Liu Bin; Crater, Horace

2003-01-01

We investigate the nucleon-nucleon interaction by using the meson exchange model and the two-body Dirac equations of constraint dynamics. This approach to the two-body problem has been successfully tested for QED and QCD relativistic bound states. An important question we wish to address is whether or not the two-body nucleon-nucleon scattering problem can be reasonably described in this approach as well. This test involves a number of related problems. First we must reduce our two-body Dirac equations exactly to a Schroedinger-like equation in such a way that allows us to use techniques to solve them already developed for Schroedinger-like systems in nonrelativistic quantum mechanics. Related to this, we present a new derivation of Calogero's variable phase shift differential equation for coupled Schroedinger-like equations. Then we determine if the use of nine meson exchanges in our equations gives a reasonable fit to the experimental scattering phase shifts for n-p scattering. The data involve seven angular momentum states including the singlet states 1 S 0 , 1 P 1 , 1 D 2 and the triplet states 3 P 0 , 3 P 1 , 3 S 1 , 3 D 1 . Two models that we have tested give us a fairly good fit. The parameters obtained by fitting the n-p experimental scattering phase shift give a fairly good prediction for most of the p-p experimental scattering phase shifts examined (for the singlet states 1 S 0 , 1 D 2 and triplet states 3 P 0 , 3 P 1 ). Thus the two-body Dirac equations of constraint dynamics present us with a fit that encourages the exploration of a more realistic model. We outline generalizations of the meson exchange model for invariant potentials that may possibly improve the fit

Scattering equations, supergravity integrands, and pure spinors

Energy Technology Data Exchange (ETDEWEB)

Adamo, Tim; Casali, Eduardo [Department of Applied Mathematics & Theoretical Physics, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA (United Kingdom)

2015-05-25

The tree-level S-matrix of type II supergravity can be computed in scattering equation form by correlators in a worldsheet theory analogous to a chiral, infinite tension limit of the pure spinor formalism. By defining a non-minimal version of this theory, we give a prescription for computing correlators on higher genus worldsheets which manifest space-time supersymmetry. These correlators are conjectured to provide the loop integrands of supergravity scattering amplitudes, supported on the scattering equations. We give non-trivial evidence in support of this conjecture at genus one and two with four external states. Throughout, we find a close correspondence with the pure spinor formalism of superstring theory, particularly regarding regulators and zero-mode counting.

Scattering equations, supergravity integrands, and pure spinors

International Nuclear Information System (INIS)

Adamo, Tim; Casali, Eduardo

2015-01-01

The tree-level S-matrix of type II supergravity can be computed in scattering equation form by correlators in a worldsheet theory analogous to a chiral, infinite tension limit of the pure spinor formalism. By defining a non-minimal version of this theory, we give a prescription for computing correlators on higher genus worldsheets which manifest space-time supersymmetry. These correlators are conjectured to provide the loop integrands of supergravity scattering amplitudes, supported on the scattering equations. We give non-trivial evidence in support of this conjecture at genus one and two with four external states. Throughout, we find a close correspondence with the pure spinor formalism of superstring theory, particularly regarding regulators and zero-mode counting.

Single-site Green function of the Dirac equation for full-potential electron scattering

Energy Technology Data Exchange (ETDEWEB)

Kordt, Pascal

2012-05-30

I present an elaborated analytical examination of the Green function of an electron scattered at a single-site potential, for both the Schroedinger and the Dirac equation, followed by an efficient numerical solution, in both cases for potentials of arbitrary shape without an atomic sphere approximation. A numerically stable way to calculate the corresponding regular and irregular wave functions and the Green function is via the angular Lippmann-Schwinger integral equations. These are solved based on an expansion in Chebyshev polynomials and their recursion relations, allowing to rewrite the Lippmann-Schwinger equations into a system of algebraic linear equations. Gonzales et al. developed this method for the Schroedinger equation, where it gives a much higher accuracy compared to previous perturbation methods, with only modest increase in computational effort. In order to apply it to the Dirac equation, I developed relativistic Lippmann-Schwinger equations, based on a decomposition of the potential matrix into spin spherical harmonics, exploiting certain properties of this matrix. The resulting method was embedded into a Korringa-Kohn-Rostoker code for density functional calculations. As an example, the method is applied by calculating phase shifts and the Mott scattering of a tungsten impurity. (orig.)

Single-site Green function of the Dirac equation for full-potential electron scattering

International Nuclear Information System (INIS)

Kordt, Pascal

2012-01-01

I present an elaborated analytical examination of the Green function of an electron scattered at a single-site potential, for both the Schroedinger and the Dirac equation, followed by an efficient numerical solution, in both cases for potentials of arbitrary shape without an atomic sphere approximation. A numerically stable way to calculate the corresponding regular and irregular wave functions and the Green function is via the angular Lippmann-Schwinger integral equations. These are solved based on an expansion in Chebyshev polynomials and their recursion relations, allowing to rewrite the Lippmann-Schwinger equations into a system of algebraic linear equations. Gonzales et al. developed this method for the Schroedinger equation, where it gives a much higher accuracy compared to previous perturbation methods, with only modest increase in computational effort. In order to apply it to the Dirac equation, I developed relativistic Lippmann-Schwinger equations, based on a decomposition of the potential matrix into spin spherical harmonics, exploiting certain properties of this matrix. The resulting method was embedded into a Korringa-Kohn-Rostoker code for density functional calculations. As an example, the method is applied by calculating phase shifts and the Mott scattering of a tungsten impurity. (orig.)

Inverse scattering transform for the time dependent Schroedinger equation with applications to the KPI equation

Energy Technology Data Exchange (ETDEWEB)

Xin, Zhou [Wisconsin Univ., Madison (USA). Dept. of Mathematics

1990-03-01

For the direct-inverse scattering transform of the time dependent Schroedinger equation, rigorous results are obtained based on an operator-triangular-factorization approach. By viewing the equation as a first order operator equation, similar results as for the first order n x n matrix system are obtained. The nonlocal Riemann-Hilbert problem for inverse scattering is shown to have solution. (orig.).

Inverse scattering transform for the time dependent Schroedinger equation with applications to the KPI equation

International Nuclear Information System (INIS)

Zhou Xin

1990-01-01

For the direct-inverse scattering transform of the time dependent Schroedinger equation, rigorous results are obtained based on an operator-triangular-factorization approach. By viewing the equation as a first order operator equation, similar results as for the first order n x n matrix system are obtained. The nonlocal Riemann-Hilbert problem for inverse scattering is shown to have solution. (orig.)

Multidimensional singular integrals and integral equations

CERN Document Server

Mikhlin, Solomon Grigorievich; Stark, M; Ulam, S

1965-01-01

Multidimensional Singular Integrals and Integral Equations presents the results of the theory of multidimensional singular integrals and of equations containing such integrals. Emphasis is on singular integrals taken over Euclidean space or in the closed manifold of Liapounov and equations containing such integrals. This volume is comprised of eight chapters and begins with an overview of some theorems on linear equations in Banach spaces, followed by a discussion on the simplest properties of multidimensional singular integrals. Subsequent chapters deal with compounding of singular integrals

Integrability of a system of two nonlinear Schroedinger equations

International Nuclear Information System (INIS)

Zhukhunashvili, V.Z.

1989-01-01

In recent years the inverse scattering method has achieved significant successes in the integration of nonlinear models that arise in different branches of physics. However, its region of applicability is still restricted, i.e., not all nonlinear models can be integrated. In view of the great mathematical difficulties that arise in integration, it is clearly worth testing a model for integrability before turning to integration. Such a possibility is provided by the Zakharov-Schulman method. The question of the integrability of a system of two nonlinear Schroedinger equations is resolved. It is shown that the previously known cases exhaust all integrable variants

Extreme Scale FMM-Accelerated Boundary Integral Equation Solver for Wave Scattering

KAUST Repository

AbdulJabbar, Mustafa Abdulmajeed

2018-03-27

Algorithmic and architecture-oriented optimizations are essential for achieving performance worthy of anticipated energy-austere exascale systems. In this paper, we present an extreme scale FMM-accelerated boundary integral equation solver for wave scattering, which uses FMM as a matrix-vector multiplication inside the GMRES iterative method. Our FMM Helmholtz kernels treat nontrivial singular and near-field integration points. We implement highly optimized kernels for both shared and distributed memory, targeting emerging Intel extreme performance HPC architectures. We extract the potential thread- and data-level parallelism of the key Helmholtz kernels of FMM. Our application code is well optimized to exploit the AVX-512 SIMD units of Intel Skylake and Knights Landing architectures. We provide different performance models for tuning the task-based tree traversal implementation of FMM, and develop optimal architecture-specific and algorithm aware partitioning, load balancing, and communication reducing mechanisms to scale up to 6,144 compute nodes of a Cray XC40 with 196,608 hardware cores. With shared memory optimizations, we achieve roughly 77% of peak single precision floating point performance of a 56-core Skylake processor, and on average 60% of peak single precision floating point performance of a 72-core KNL. These numbers represent nearly 5.4x and 10x speedup on Skylake and KNL, respectively, compared to the baseline scalar code. With distributed memory optimizations, on the other hand, we report near-optimal efficiency in the weak scalability study with respect to both the logarithmic communication complexity as well as the theoretical scaling complexity of FMM. In addition, we exhibit up to 85% efficiency in strong scaling. We compute in excess of 2 billion DoF on the full-scale of the Cray XC40 supercomputer.

Inverse Scattering Problem For The Schrödinger Equation With An Additional Quadratic Potential On The Entire Axis

Science.gov (United States)

Guseinov, I. M.; Khanmamedov, A. Kh.; Mamedova, A. F.

2018-04-01

We consider the Schrödinger equation with an additional quadratic potential on the entire axis and use the transformation operator method to study the direct and inverse problems of the scattering theory. We obtain the main integral equations of the inverse problem and prove that the basic equations are uniquely solvable.

Fast Near-Field Calculation for Volume Integral Equations for Layered Media

DEFF Research Database (Denmark)

Kim, Oleksiy S.; Meincke, Peter; Breinbjerg, Olav

2005-01-01

. Afterwards, the scattered electric field can be easily computed at a regular rectangular grid on any horizontal plane us-ing a 2-dimensional FFT. This approach provides significant speedup in the near-field calculation in comparison to a straightforward numerical evaluation of the ra-diation integral since......An efficient technique based on the Fast Fourier Transform (FFT) for calculating near-field scattering by dielectric objects in layered media is presented. A higher or-der method of moments technique is employed to solve the volume integral equation for the unknown induced volume current density...

Approximate solutions for the two-dimensional integral transport equation. The critically mixed methods of resolution

International Nuclear Information System (INIS)

Sanchez, Richard.

1980-11-01

This work is divided into two part the first part (note CEA-N-2165) deals with the solution of complex two-dimensional transport problems, the second one treats the critically mixed methods of resolution. These methods are applied for one-dimensional geometries with highly anisotropic scattering. In order to simplify the set of integral equation provided by the integral transport equation, the integro-differential equation is used to obtain relations that allow to lower the number of integral equation to solve; a general mathematical and numerical study is presented [fr

Fringe integral equation method for a truncated grounded dielectric slab

DEFF Research Database (Denmark)

Jørgensen, Erik; Maci, S.; Toccafondi, A.

2001-01-01

The problem of scattering by a semi-infinite grounded dielectric slab illuminated by an arbitrary incident TMz polarized electric field is studied by solving a new set of “fringe” integral equations (F-IEs), whose functional unknowns are physically associated to the wave diffraction processes...

Remarks on the inverse scattering transform associated with toda equations

Science.gov (United States)

Ablowitz, Mark J.; Villorroel, J.

The Inverse Scattering Transforms used to solve both the 2+1 Toda equation and a novel reduction, the Toda differential-delay equations are outlined. There are a number of interesting features associated with these systems and the related scattering theory.

Solving very large scattering problems using a parallel PWTD-enhanced surface integral equation solver

KAUST Repository

Liu, Yang

2013-07-01

The computational complexity and memory requirements of multilevel plane wave time domain (PWTD)-accelerated marching-on-in-time (MOT)-based surface integral equation (SIE) solvers scale as O(NtNs(log 2)Ns) and O(Ns 1.5); here N t and Ns denote numbers of temporal and spatial basis functions discretizing the current [Shanker et al., IEEE Trans. Antennas Propag., 51, 628-641, 2003]. In the past, serial versions of these solvers have been successfully applied to the analysis of scattering from perfect electrically conducting as well as homogeneous penetrable targets involving up to Ns ≈ 0.5 × 106 and Nt ≈ 10 3. To solve larger problems, parallel PWTD-enhanced MOT solvers are called for. Even though a simple parallelization strategy was demonstrated in the context of electromagnetic compatibility analysis [M. Lu et al., in Proc. IEEE Int. Symp. AP-S, 4, 4212-4215, 2004], by and large, progress in this area has been slow. The lack of progress can be attributed wholesale to difficulties associated with the construction of a scalable PWTD kernel. © 2013 IEEE.

Handbook of integral equations

CERN Document Server

Polyanin, Andrei D

2008-01-01

This handbook contains over 2,500 integral equations with solutions as well as analytical and numerical methods for solving linear and nonlinear equations. It explores Volterra, Fredholm, WienerHopf, Hammerstein, Uryson, and other equations that arise in mathematics, physics, engineering, the sciences, and economics. This second edition includes new chapters on mixed multidimensional equations and methods of integral equations for ODEs and PDEs, along with over 400 new equations with exact solutions. With many examples added for illustrative purposes, it presents new material on Volterra, Fredholm, singular, hypersingular, dual, and nonlinear integral equations, integral transforms, and special functions.

Numerical Integration of the Transport Equation For Infinite Homogeneous Media

Energy Technology Data Exchange (ETDEWEB)

Haakansson, Rune

1962-01-15

The transport equation for neutrons in infinite homogeneous media is solved by direct numerical integration. Accounts are taken to the anisotropy and the inelastic scattering. The integration has been performed by means of the trapezoidal rule and the length of the energy intervals are constant in lethargy scale. The machine used is a Ferranti Mercury computer. Results are given for water, heavy water, aluminium water mixture and iron-aluminium-water mixture.

Inverse scattering transform for the time dependent Schrödinger equation with applications to the KPI equation

Science.gov (United States)

Zhou, Xin

1990-03-01

For the direct-inverse scattering transform of the time dependent Schrödinger equation, rigorous results are obtained based on an opertor-triangular-factorization approach. By viewing the equation as a first order operator equation, similar results as for the first order n x n matrix system are obtained. The nonlocal Riemann-Hilbert problem for inverse scattering is shown to have solution.

An algebraic approach to the scattering equations

Energy Technology Data Exchange (ETDEWEB)

Huang, Rijun; Rao, Junjie [Zhejiang Institute of Modern Physics, Zhejiang University,Hangzhou, 310027 (China); Feng, Bo [Zhejiang Institute of Modern Physics, Zhejiang University,Hangzhou, 310027 (China); Center of Mathematical Science, Zhejiang University,Hangzhou, 310027 (China); He, Yang-Hui [School of Physics, NanKai University,Tianjin, 300071 (China); Department of Mathematics, City University,London, EC1V 0HB (United Kingdom); Merton College, University of Oxford,Oxford, OX14JD (United Kingdom)

2015-12-10

We employ the so-called companion matrix method from computational algebraic geometry, tailored for zero-dimensional ideals, to study the scattering equations. The method renders the CHY-integrand of scattering amplitudes computable using simple linear algebra and is amenable to an algorithmic approach. Certain identities in the amplitudes as well as rationality of the final integrand become immediate in this formalism.

An algebraic approach to the scattering equations

International Nuclear Information System (INIS)

Huang, Rijun; Rao, Junjie; Feng, Bo; He, Yang-Hui

2015-01-01

We employ the so-called companion matrix method from computational algebraic geometry, tailored for zero-dimensional ideals, to study the scattering equations. The method renders the CHY-integrand of scattering amplitudes computable using simple linear algebra and is amenable to an algorithmic approach. Certain identities in the amplitudes as well as rationality of the final integrand become immediate in this formalism.

Angular finite volume method for solving the multigroup transport equation with piecewise average scattering cross sections

International Nuclear Information System (INIS)

Calloo, A.; Vidal, J.F.; Le Tellier, R.; Rimpault, G.

2011-01-01

This paper deals with the solving of the multigroup integro-differential form of the transport equation for fine energy group structure. In that case, multigroup transfer cross sections display strongly peaked shape for light scatterers and the current Legendre polynomial expansion is not well-suited to represent them. Furthermore, even if considering an exact scattering cross sections representation, the scattering source in the discrete ordinates method (also known as the Sn method) being calculated by sampling the angular flux at given directions, may be wrongly computed due to lack of angular support for the angular flux. Hence, following the work of Gerts and Matthews, an angular finite volume solver has been developed for 2D Cartesian geometries. It integrates the multigroup transport equation over discrete volume elements obtained by meshing the unit sphere with a product grid over the polar and azimuthal coordinates and by considering the integrated flux per solid angle element. The convergence of this method has been compared to the S_n method for a highly anisotropic benchmark. Besides, piecewise-average scattering cross sections have been produced for non-bound Hydrogen atoms using a free gas model for thermal neutrons. LWR lattice calculations comparing Legendre representations of the Hydrogen scattering multigroup cross section at various orders and piecewise-average cross sections for this same atom are carried out (while keeping a Legendre representation for all other isotopes). (author)

On the asymptotic solution to a class of linear integral equations

International Nuclear Information System (INIS)

Gautesen, A.K.

1988-01-01

The authors consider Fredholm integral equations of the first kind whose kernels are a function of the difference between two points times a large parameter. Conditions on the kernel are stated in terms of a function corresponding to a Wiener-Hopf factorization of the Fourier transform of the kernel. They give the complete asymptotic expansions of the solution to the integral equations. As applications of the author's results, the author considers the steady-state, acoustical scattering of a plane wave by both a hard strip and a soft strip. The author's results are uniform with respect to the direction of incidence

MOT solution of the PMCHWT equation for analyzing transient scattering from conductive dielectrics

KAUST Repository

Uysal, Ismail Enes

2015-01-01

Transient electromagnetic interactions on conductive dielectric scatterers are analyzed by solving the Poggio-Miller-Chan-Harrington-Wu-Tsai (PMCHWT) surface integral equation with a marching on-in-time (MOT) scheme. The proposed scheme, unlike the previously developed ones, permits the analysis on scatterers with multiple volumes of different conductivity. This is achieved by maintaining an extra temporal convolution that only depends on permittivity and conductivity of these volumes. Its discretization and computation come at almost no additional cost and do not change the computational complexity of the resulting MOT solver. Accuracy and applicability of the MOT-PMCHWT solver are demonstrated by numerical examples.

A Wavelet-Enhanced PWTD-Accelerated Time-Domain Integral Equation Solver for Analysis of Transient Scattering from Electrically Large Conducting Objects

KAUST Repository

Liu, Yang

2018-02-26

A wavelet-enhanced plane-wave time-domain (PWTD) algorithm for efficiently and accurately solving time-domain surface integral equations (TD-SIEs) on electrically large conducting objects is presented. The proposed scheme reduces the memory requirement and computational cost of the PWTD algorithm by representing the PWTD ray data using local cosine wavelet bases (LCBs) and performing PWTD operations in the wavelet domain. The memory requirement and computational cost of the LCB-enhanced PWTD-accelerated TD-SIE solver, when applied to the analysis of transient scattering from smooth quasi-planar objects with near-normal incident pulses, scale nearly as O(Ns log Ns) and O(Ns 1.5 ), respectively. Here, Ns denotes the number of spatial unknowns. The efficiency and accuracy of the proposed scheme are demonstrated through its applications to the analysis of transient scattering from a 185 wave-length-long NASA almond and a 123-wavelength long Air-bus-A320 model.

PREFACE: Symmetries and integrability of difference equations Symmetries and integrability of difference equations

Science.gov (United States)

Levi, Decio; Olver, Peter; Thomova, Zora; Winternitz, Pavel

2009-11-01

The concept of integrability was introduced in classical mechanics in the 19th century for finite dimensional continuous Hamiltonian systems. It was extended to certain classes of nonlinear differential equations in the second half of the 20th century with the discovery of the inverse scattering transform and the birth of soliton theory. Also at the end of the 19th century Lie group theory was invented as a powerful tool for obtaining exact analytical solutions of large classes of differential equations. Together, Lie group theory and integrability theory in its most general sense provide the main tools for solving nonlinear differential equations. Like differential equations, difference equations play an important role in physics and other sciences. They occur very naturally in the description of phenomena that are genuinely discrete. Indeed, they may actually be more fundamental than differential equations if space-time is actually discrete at very short distances. On the other hand, even when treating continuous phenomena described by differential equations it is very often necessary to resort to numerical methods. This involves a discretization of the differential equation, i.e. a replacement of the differential equation by a difference one. Given the well developed and understood techniques of symmetry and integrability for differential equations a natural question to ask is whether it is possible to develop similar techniques for difference equations. The aim is, on one hand, to obtain powerful methods for solving `integrable' difference equations and to establish practical integrability criteria, telling us when the methods are applicable. On the other hand, Lie group methods can be adapted to solve difference equations analytically. Finally, integrability and symmetry methods can be combined with numerical methods to obtain improved numerical solutions of differential equations. The origin of the SIDE meetings goes back to the early 1990s and the first

One-way spatial integration of hyperbolic equations

Science.gov (United States)

Towne, Aaron; Colonius, Tim

2015-11-01

In this paper, we develop and demonstrate a method for constructing well-posed one-way approximations of linear hyperbolic systems. We use a semi-discrete approach that allows the method to be applied to a wider class of problems than existing methods based on analytical factorization of idealized dispersion relations. After establishing the existence of an exact one-way equation for systems whose coefficients do not vary along the axis of integration, efficient approximations of the one-way operator are constructed by generalizing techniques previously used to create nonreflecting boundary conditions. When physically justified, the method can be applied to systems with slowly varying coefficients in the direction of integration. To demonstrate the accuracy and computational efficiency of the approach, the method is applied to model problems in acoustics and fluid dynamics via the linearized Euler equations; in particular we consider the scattering of sound waves from a vortex and the evolution of hydrodynamic wavepackets in a spatially evolving jet. The latter problem shows the potential of the method to offer a systematic, convergent alternative to ad hoc regularizations such as the parabolized stability equations.

Parallel, explicit, and PWTD-enhanced time domain volume integral equation solver

KAUST Repository

Liu, Yang

2013-07-01

Time domain volume integral equations (TDVIEs) are useful for analyzing transient scattering from inhomogeneous dielectric objects in applications as varied as photonics, optoelectronics, and bioelectromagnetics. TDVIEs typically are solved by implicit marching-on-in-time (MOT) schemes [N. T. Gres et al., Radio Sci., 36, 379-386, 2001], requiring the solution of a system of equations at each and every time step. To reduce the computational cost associated with such schemes, [A. Al-Jarro et al., IEEE Trans. Antennas Propagat., 60, 5203-5215, 2012] introduced an explicit MOT-TDVIE method that uses a predictor-corrector technique to stably update field values throughout the scatterer. By leveraging memory-efficient nodal spatial discretization and scalable parallelization schemes [A. Al-Jarro et al., in 28th Int. Rev. Progress Appl. Computat. Electromagn., 2012], this solver has been successfully applied to the analysis of scattering phenomena involving 0.5 million spatial unknowns. © 2013 IEEE.

Conformational effect on small angle neutron scattering behavior of interacting polyelectrolyte solutions: a perspective of integral equation theory.

Science.gov (United States)

Shew, Chwen-Yang; Do, Changwoo; Hong, Kunlun; Liu, Yun; Porcar, Lionel; Smith, Gregory S; Chen, Wei-Ren

2012-07-14

We present small angle neutron scattering (SANS) measurements of deuterium oxide (D(2)O) solutions of linear and star sodium poly(styrene sulfonate) (NaPSS) as a function of polyelectrolyte concentration. Emphasis is on understanding the dependence of their SANS coherent scattering cross section I(Q) on the molecular architecture of single polyelectrolyte. The key finding is that for a given concentration, star polyelectrolytes exhibit more pronounced characteristic peaks in I(Q), and the position of the first peak occurs at a smaller Q compared to their linear counterparts. Based on a model of integral equation theory, we first compare the SANS experimental I(Q) of salt-free polyelectrolyte solutions with that predicted theoretically. Having seen their satisfactory qualitative agreement, the dependence of counterion association behavior on polyelectrolyte geometry and concentration is further explored. Our predictions reveal that the ionic environment of polyelectrolyte exhibits a strong dependence on polyelectrolyte geometry at lower polyelectrolyte concentration. However, when both linear and star polyelectrolytes exceed their overlap concentrations, the spatial distribution of counterion is found to be essentially insensitive to polyelectrolyte geometry due to the steric effect.

Solution of volume-surface integral equations using higher-order hierarchical Legendre basis functions

DEFF Research Database (Denmark)

Kim, Oleksiy S.; Meincke, Peter; Breinbjerg, Olav

2007-01-01

The problem of electromagnetic scattering by composite metallic and dielectric objects is solved using the coupled volume-surface integral equation (VSIE). The method of moments (MoM) based on higher-order hierarchical Legendre basis functions and higher-order curvilinear geometrical elements...... with the analytical Mie series solution. Scattering by more complex metal-dielectric objects are also considered to compare the presented technique with other numerical methods....

Angular finite volume method for solving the multigroup transport equation with piecewise average scattering cross sections

Energy Technology Data Exchange (ETDEWEB)

Calloo, A.; Vidal, J.F.; Le Tellier, R.; Rimpault, G., E-mail: ansar.calloo@cea.fr, E-mail: jean-francois.vidal@cea.fr, E-mail: romain.le-tellier@cea.fr, E-mail: gerald.rimpault@cea.fr [CEA, DEN, DER/SPRC/LEPh, Saint-Paul-lez-Durance (France)

2011-07-01

This paper deals with the solving of the multigroup integro-differential form of the transport equation for fine energy group structure. In that case, multigroup transfer cross sections display strongly peaked shape for light scatterers and the current Legendre polynomial expansion is not well-suited to represent them. Furthermore, even if considering an exact scattering cross sections representation, the scattering source in the discrete ordinates method (also known as the Sn method) being calculated by sampling the angular flux at given directions, may be wrongly computed due to lack of angular support for the angular flux. Hence, following the work of Gerts and Matthews, an angular finite volume solver has been developed for 2D Cartesian geometries. It integrates the multigroup transport equation over discrete volume elements obtained by meshing the unit sphere with a product grid over the polar and azimuthal coordinates and by considering the integrated flux per solid angle element. The convergence of this method has been compared to the S{sub n} method for a highly anisotropic benchmark. Besides, piecewise-average scattering cross sections have been produced for non-bound Hydrogen atoms using a free gas model for thermal neutrons. LWR lattice calculations comparing Legendre representations of the Hydrogen scattering multigroup cross section at various orders and piecewise-average cross sections for this same atom are carried out (while keeping a Legendre representation for all other isotopes). (author)

Retarded potentials and time domain boundary integral equations a road map

CERN Document Server

Sayas, Francisco-Javier

2016-01-01

This book offers a thorough and self-contained exposition of the mathematics of time-domain boundary integral equations associated to the wave equation, including applications to scattering of acoustic and elastic waves. The book offers two different approaches for the analysis of these integral equations, including a systematic treatment of their numerical discretization using Galerkin (Boundary Element) methods in the space variables and Convolution Quadrature in the time variable. The first approach follows classical work started in the late eighties, based on Laplace transforms estimates. This approach has been refined and made more accessible by tailoring the necessary mathematical tools, avoiding an excess of generality. A second approach contains a novel point of view that the author and some of his collaborators have been developing in recent years, using the semigroup theory of evolution equations to obtain improved results. The extension to electromagnetic waves is explained in one of the appendices...

Steady-state equations of even flux and scattering

International Nuclear Information System (INIS)

Verwaerde, D.

1985-11-01

Some mathematical properties of steady-state equation of even flux are shown in variational formalism. This theoretical frame allows to study the existence of a solution and its asymptotical behavior in opaque media (i.e. the relation with scattering equation). At last it allows to qualify the convergence velocity of resolution iterative processes used practically [fr

Integral equations with difference kernels on finite intervals

CERN Document Server

Sakhnovich, Lev A

2015-01-01

This book focuses on solving integral equations with difference kernels on finite intervals. The corresponding problem on the semiaxis was previously solved by N. Wiener–E. Hopf and by M.G. Krein. The problem on finite intervals, though significantly more difficult, may be solved using our method of operator identities. This method is also actively employed in inverse spectral problems, operator factorization and nonlinear integral equations. Applications of the obtained results to optimal synthesis, light scattering, diffraction, and hydrodynamics problems are discussed in this book, which also describes how the theory of operators with difference kernels is applied to stable processes and used to solve the famous M. Kac problems on stable processes. In this second edition these results are extensively generalized and include the case of all Levy processes. We present the convolution expression for the well-known Ito formula of the generator operator, a convolution expression that has proven to be fruitful...

Application of the Galerkin's method to the solution of the one-dimensional integral transport equation: generalized collision probabilities taken in account the flux gradient and the linearly anisotropic scattering

International Nuclear Information System (INIS)

Sanchez, Richard.

1975-04-01

For the one-dimensional geometries, the transport equation with linearly anisotropic scattering can be reduced to a single integral equation; this is a singular-kernel FREDHOLM equation of the second kind. When applying a conventional projective method that of GALERKIN, to the solution of this equation the well-known collision probability algorithm is obtained. Piecewise polynomial expansions are used to represent the flux. In the ANILINE code, the flux is supposed to be linear in plane geometry and parabolic in both cylindrical and spherical geometries. An integral relationship was found between the one-dimensional isotropic and anisotropic kernels; this allows to reduce the new matrix elements (issuing from the anisotropic kernel) to classic collision probabilities of the isotropic scattering equation. For cylindrical and spherical geometries used an approximate representation of the current was used to avoid an additional numerical integration. Reflective boundary conditions were considered; in plane geometry the reflection is supposed specular, for the other geometries the isotropic reflection hypothesis has been adopted. Further, the ANILINE code enables to deal with an incoming isotropic current. Numerous checks were performed in monokinetic theory. Critical radii and albedos were calculated for homogeneous slabs, cylinders and spheres. For heterogeneous media, the thermal utilization factor obtained by this method was compared with the theoretical result based upon a formula by BENOIST. Finally, ANILINE was incorporated into the multigroup APOLLO code, which enabled to analyse the MINERVA experimental reactor in transport theory with 99 groups. The ANILINE method is particularly suited to the treatment of strongly anisotropic media with considerable flux gradients. It is also well adapted to the calculation of reflectors, and in general, to the exact analysis of anisotropic effects in large-sized media [fr

Scattering Behavior of Waveguide Channels of a New Coupled Integrable Dispersionless System

International Nuclear Information System (INIS)

Souleymanou, Abbagari; Kuetche, Victor K.; Bouetou, Thomas B.; Kofane, Timoleon C.

2011-01-01

Based upon the powerful Hirota method for unearthing soliton solutions to nonlinear partial differential evolution equations, we investigate the scattering properties of a new coupled integrable dispersionless system while surveying the interactions between its self-confined travelling wave solutions. As a result, we ascertain three types of scattering features depending strongly upon a characteristic parameter. Using such findings to depict soliton solutions with nonzero angular momenta, we derive an extended form of the dispersionless system, which is valuable for further physical applications. (general)

On the initial condition problem of the time domain PMCHWT surface integral equation

KAUST Repository

Uysal, Ismail Enes

2017-05-13

Non-physical, linearly increasing and constant current components are induced in marching on-in-time solution of time domain surface integral equations when initial conditions on time derivatives of (unknown) equivalent currents are not enforced properly. This problem can be remedied by solving the time integral of the surface integral for auxiliary currents that are defined to be the time derivatives of the equivalent currents. Then the equivalent currents are obtained by numerically differentiating the auxiliary ones. In this work, this approach is applied to the marching on-in-time solution of the time domain Poggio-Miller-Chan-Harrington-Wu-Tsai surface integral equation enforced on dispersive/plasmonic scatterers. Accuracy of the proposed method is demonstrated by a numerical example.

Combined Helmholtz Integral Equation - Fourier series formulation of acoustical radiation and scattering problems

CSIR Research Space (South Africa)

Fedotov, I

2006-07-01

Full Text Available The Combined Helmholtz Integral Equation – Fourier series Formulation (CHIEFF) is based on representation of a velocity potential in terms of Fourier series and finding the Fourier coefficients of this expansion. The solution could be substantially...

Recent advances in marching-on-in-time schemes for solving time domain volume integral equations

KAUST Repository

Sayed, Sadeed Bin; Ulku, Huseyin Arda; Bagci, Hakan

2015-01-01

Transient electromagnetic field interactions on inhomogeneous penetrable scatterers can be analyzed by solving time domain volume integral equations (TDVIEs). TDVIEs are constructed by setting the summation of the incident and scattered field intensities to the total field intensity on the volumetric support of the scatterer. The unknown can be the field intensity or flux/current density. Representing the total field intensity in terms of the unknown using the relevant constitutive relation and the scattered field intensity in terms of the spatiotemporal convolution of the unknown with the Green function yield the final form of the TDVIE. The unknown is expanded in terms of local spatial and temporal basis functions. Inserting this expansion into the TDVIE and testing the resulting equation at discrete times yield a system of equations that is solved by the marching on-in-time (MOT) scheme. At each time step, a smaller system of equations, termed MOT system is solved for the coefficients of the expansion. The right-hand side of this system consists of the tested incident field and discretized spatio-temporal convolution of the unknown samples computed at the previous time steps with the Green function.

Recent advances in marching-on-in-time schemes for solving time domain volume integral equations

KAUST Repository

Sayed, Sadeed Bin

2015-05-16

Transient electromagnetic field interactions on inhomogeneous penetrable scatterers can be analyzed by solving time domain volume integral equations (TDVIEs). TDVIEs are constructed by setting the summation of the incident and scattered field intensities to the total field intensity on the volumetric support of the scatterer. The unknown can be the field intensity or flux/current density. Representing the total field intensity in terms of the unknown using the relevant constitutive relation and the scattered field intensity in terms of the spatiotemporal convolution of the unknown with the Green function yield the final form of the TDVIE. The unknown is expanded in terms of local spatial and temporal basis functions. Inserting this expansion into the TDVIE and testing the resulting equation at discrete times yield a system of equations that is solved by the marching on-in-time (MOT) scheme. At each time step, a smaller system of equations, termed MOT system is solved for the coefficients of the expansion. The right-hand side of this system consists of the tested incident field and discretized spatio-temporal convolution of the unknown samples computed at the previous time steps with the Green function.

Low-frequency scaling of the standard and mixed magnetic field and Müller integral equations

KAUST Repository

Bogaert, Ignace

2014-02-01

The standard and mixed discretizations for the magnetic field integral equation (MFIE) and the Müller integral equation (MUIE) are investigated in the context of low-frequency (LF) scattering problems involving simply connected scatterers. It is proved that, at low frequencies, the frequency scaling of the nonsolenoidal part of the solution current can be incorrect for the standard discretization. In addition, it is proved that the frequency scaling obtained with the mixed discretization is correct. The reason for this problem in the standard discretization scheme is the absence of exact solenoidal currents in the rotated RWG finite element space. The adoption of the mixed discretization scheme eliminates this problem and leads to a well-conditioned system of linear equations that remains accurate at low frequencies. Numerical results confirm these theoretical predictions and also show that, when the frequency is lowered, a finer and finer mesh is required to keep the accuracy constant with the standard discretization. © 1963-2012 IEEE.

Scattering Of Nonplanar Acoustic Waves

Science.gov (United States)

Gillman, Judith M.; Farassat, F.; Myers, M. K.

1995-01-01

Report presents theoretical study of scattering of nonplanar acoustic waves by rigid bodies. Study performed as part of effort to develop means of predicting scattering, from aircraft fuselages, of noise made by rotating blades. Basic approach was to model acoustic scattering by use of boundary integral equation to solve equation by the Galerkin method.

Numerical method for solving integral equations of neutron transport. II

International Nuclear Information System (INIS)

Loyalka, S.K.; Tsai, R.W.

1975-01-01

In a recent paper it was pointed out that the weakly singular integral equations of neutron transport can be quite conveniently solved by a method based on subtraction of singularity. This previous paper was devoted entirely to the consideration of simple one-dimensional isotropic-scattering and one-group problems. The present paper constitutes interesting extensions of the previous work in that in addition to a typical two-group anisotropic-scattering albedo problem in the slab geometry, the method is also applied to an isotropic-scattering problem in the x-y geometry. These results are compared with discrete S/sub N/ (ANISN or TWOTRAN-II) results, and for the problems considered here, the proposed method is found to be quite effective. Thus, the method appears to hold considerable potential for future applications. (auth)

Inverse scattering scheme for the Dirac equation at fixed energy

International Nuclear Information System (INIS)

Leeb, H.; Lehninger, H.; Schilder, C.

2001-01-01

Full text: Based on the concept of generalized transformation operators a new hierarchy of Dirac equations with spherical symmetric scalar and fourth component vector potentials is presented. Within this hierarchy closed form expressions for the solutions, the potentials and the S-matrix can be given in terms of solutions of the original Dirac equation. Using these transformations an inverse scattering scheme has been constructed for the Dirac equation which is the analog to the rational scheme in the non-relativistic case. The given method provides for the first time an inversion scheme with closed form expressions for the S-matrix for non-relativistic scattering problems with central and spin-orbit potentials. (author)

The ICVSIE: A General Purpose Integral Equation Method for Bio-Electromagnetic Analysis.

Science.gov (United States)

Gomez, Luis J; Yucel, Abdulkadir C; Michielssen, Eric

2018-03-01

An internally combined volume surface integral equation (ICVSIE) for analyzing electromagnetic (EM) interactions with biological tissue and wide ranging diagnostic, therapeutic, and research applications, is proposed. The ICVSIE is a system of integral equations in terms of volume and surface equivalent currents in biological tissue subject to fields produced by externally or internally positioned devices. The system is created by using equivalence principles and solved numerically; the resulting current values are used to evaluate scattered and total electric fields, specific absorption rates, and related quantities. The validity, applicability, and efficiency of the ICVSIE are demonstrated by EM analysis of transcranial magnetic stimulation, magnetic resonance imaging, and neuromuscular electrical stimulation. Unlike previous integral equations, the ICVSIE is stable regardless of the electric permittivities of the tissue or frequency of operation, providing an application-agnostic computational framework for EM-biomedical analysis. Use of the general purpose and robust ICVSIE permits streamlining the development, deployment, and safety analysis of EM-biomedical technologies.

Regge-pole description of potential scattering by means of the phase-integral method

International Nuclear Information System (INIS)

Amaha, A.

1992-01-01

The application of Regge-pole theory to different atomic and molecular scattering has shown to have promising interpretational power in the differential cross sections. Differential cross sections can be analysed in terms of interference between the 'background' amplitude and a few Regge-pole positions of the scattering matrix (S matrix) representing surface waves around the interaction region. By the analytic continuation of the radial Schroedinger differential equation into the complex plane of angular momentum one can determine the analytic properties of the S matrix which contains the physical information in the scattering processes. For interaction potentials fulfilling certain properties, the study of the S matrix leads to the study of the F matrix introduced by Froeman and Froeman for the treatment of connection problems for phase-integral solutions of the differential equation. In this thesis the quantum mechanical scattering problem is analysed in the framework of Regge-pole theory with the use of the complex-angular-momentum formalism. To determine the S matrix, the relevant F matrix elements which give the stokes constants are derived and their properties are studied. The poles of the S matrix for particular complex values of the angular momentum quantum number are the Regge-poles. Using the Regge-pole positions and residues together with the background integral, the differential cross sections are calculated and compared with corresponding partial-wave representations

Relativistic three-particle dynamical equations: I. Theoretical development

International Nuclear Information System (INIS)

Adhikari, S.K.; Tomio, L.; Frederico, T.

1993-11-01

Starting from the two-particle Bethe-Salpeter equation in the ladder approximation and integrating over the time component of momentum, three dimensional scattering integral equations satisfying constrains of relativistic unitarity and covariance are rederived. These equations were first derived by Weinberg and by Blankenbecler and Sugar. These two-particle equations are shown to be related by a transformation of variables. Hence it is shown to perform and relate dynamical calculation using these two equations. Similarly, starting from the Bethe-Salpeter-Faddeev equation for the three-particle system and integrating over the time component of momentum, several three dimensional three-particle scattering equations satisfying constraints of relativistic unitary and covariance are derived. Two of these three-particle equations are related by a transformation of variables as in the two-particle case. The three-particle equations obtained are very practical and suitable for performing relativistic scattering calculations. (author)

Integral equations and their applications

CERN Document Server

Rahman, M

2007-01-01

For many years, the subject of functional equations has held a prominent place in the attention of mathematicians. In more recent years this attention has been directed to a particular kind of functional equation, an integral equation, wherein the unknown function occurs under the integral sign. The study of this kind of equation is sometimes referred to as the inversion of a definite integral. While scientists and engineers can already choose from a number of books on integral equations, this new book encompasses recent developments including some preliminary backgrounds of formulations of integral equations governing the physical situation of the problems. It also contains elegant analytical and numerical methods, and an important topic of the variational principles. Primarily intended for senior undergraduate students and first year postgraduate students of engineering and science courses, students of mathematical and physical sciences will also find many sections of direct relevance. The book contains eig...

General algebraic theory of identical particle scattering

International Nuclear Information System (INIS)

Bencze, G.; Redish, E.F.

1978-01-01

We consider the nonrelativistic N-body scattering problem for a system of particles in which some subsets of the particles are identical. We demonstrate how the particle identity can be included in a general class of linear integral equations for scattering operators or components of scattering operators. The Yakubovskii, Yakubovskii--Narodestkii, Rosenberg, and Bencze--Redish--Sloan equations are included in this class. Algebraic methods are used which rely on the properties of the symmetry group of the system. Operators depending only on physically distinguishable labels are introduced and linear integral equations for them are derived. This procedure maximally reduces the number of coupled equations while retaining the connectivity properties of the original equations

Inclusion of exact exchange in the noniterative partial-differential-equation method of electron-molecule scattering - Application to e-N2

Science.gov (United States)

Weatherford, C. A.; Onda, K.; Temkin, A.

1985-01-01

The noniterative partial-differential-equation (PDE) approach to electron-molecule scattering of Onda and Temkin (1983) is modified to account for the effects of exchange explicitly. The exchange equation is reduced to a set of inhomogeneous equations containing no integral terms and solved noniteratively in a difference form; a method for propagating the solution to large values of r is described; the changes in the polarization potential of the original PDE method required by the inclusion of exact static exchange are indicated; and the results of computations for e-N2 scattering in the fixed-nuclei approximation are presented in tables and graphs and compared with previous calculations and experimental data. Better agreement is obtained using the modified PDE method.

Transport equation solving methods

International Nuclear Information System (INIS)

Granjean, P.M.

1984-06-01

This work is mainly devoted to Csub(N) and Fsub(N) methods. CN method: starting from a lemma stated by Placzek, an equivalence is established between two problems: the first one is defined in a finite medium bounded by a surface S, the second one is defined in the whole space. In the first problem the angular flux on the surface S is shown to be the solution of an integral equation. This equation is solved by Galerkin's method. The Csub(N) method is applied here to one-velocity problems: in plane geometry, slab albedo and transmission with Rayleigh scattering, calculation of the extrapolation length; in cylindrical geometry, albedo and extrapolation length calculation with linear scattering. Fsub(N) method: the basic integral transport equation of the Csub(N) method is integrated on Case's elementary distributions; another integral transport equation is obtained: this equation is solved by a collocation method. The plane problems solved by the Csub(N) method are also solved by the Fsub(N) method. The Fsub(N) method is extended to any polynomial scattering law. Some simple spherical problems are also studied. Chandrasekhar's method, collision probability method, Case's method are presented for comparison with Csub(N) and Fsub(N) methods. This comparison shows the respective advantages of the two methods: a) fast convergence and possible extension to various geometries for Csub(N) method; b) easy calculations and easy extension to polynomial scattering for Fsub(N) method [fr

Application of the Radiative Transfer Equation (RTE) to Scattering by ...

African Journals Online (AJOL)

Application of the Radiative Transfer Equation (RTE) to Scattering by a Dust Aerosol Layer. ... Incident radiation in its journey through the atmosphere before reaching the earth surface encounters particles of different sizes and composition such as dust aerosols resulting in interactions that lead to absorption and scattering.

Coupled radiative transfer equation and diffusion approximation model for photon migration in turbid medium with low-scattering and non-scattering regions

International Nuclear Information System (INIS)

Tarvainen, Tanja; Vauhkonen, Marko; Kolehmainen, Ville; Arridge, Simon R; Kaipio, Jari P

2005-01-01

In this paper, a coupled radiative transfer equation and diffusion approximation model is extended for light propagation in turbid medium with low-scattering and non-scattering regions. The light propagation is modelled with the radiative transfer equation in sub-domains in which the assumptions of the diffusion approximation are not valid. The diffusion approximation is used elsewhere in the domain. The two equations are coupled through their boundary conditions and they are solved simultaneously using the finite element method. The streamline diffusion modification is used to avoid the ray-effect problem in the finite element solution of the radiative transfer equation. The proposed method is tested with simulations. The results of the coupled model are compared with the finite element solutions of the radiative transfer equation and the diffusion approximation and with results of Monte Carlo simulation. The results show that the coupled model can be used to describe photon migration in turbid medium with low-scattering and non-scattering regions more accurately than the conventional diffusion model

Bound states and scattering in four-body systems

International Nuclear Information System (INIS)

Narodetsky, I.M.

1979-01-01

It is the purpose of this review to provide the clear and elementary introduction in the integral equation method and to demonstrate explicitely its usefulness for the physical applications. The existing results concerning the application of the integral equation technique for the four-nucleon bound states and scattering are reviewed.The treatment is based on the quasiparticle approach that permits the simple interpretation of the equations in terms of quasiparticle scattering. The mathematical basis for the quasiparticle approach is the Hilbert-Schmidt theorem of the Fredholm integral equation theory. This paper contains the detailed discussion of the Hilbert-Schmidt expansion as applied to the 2-particle amplitudes and to the 3 + 1 and 2 + 2 amplitudes which are the kernels of the four-body equations. The review contains essentially the discussion of the four-body quasiparticle equations and results obtained for bound states and scattering

Roy-Steiner equations for pion-nucleon scattering

Science.gov (United States)

Ditsche, C.; Hoferichter, M.; Kubis, B.; Meißner, U.-G.

2012-06-01

Starting from hyperbolic dispersion relations, we derive a closed system of Roy-Steiner equations for pion-nucleon scattering that respects analyticity, unitarity, and crossing symmetry. We work out analytically all kernel functions and unitarity relations required for the lowest partial waves. In order to suppress the dependence on the high energy regime we also consider once- and twice-subtracted versions of the equations, where we identify the subtraction constants with subthreshold parameters. Assuming Mandelstam analyticity we determine the maximal range of validity of these equations. As a first step towards the solution of the full system we cast the equations for the π π to overline N N partial waves into the form of a Muskhelishvili-Omnès problem with finite matching point, which we solve numerically in the single-channel approximation. We investigate in detail the role of individual contributions to our solutions and discuss some consequences for the spectral functions of the nucleon electromagnetic form factors.

On integrability of the Killing equation

Science.gov (United States)

Houri, Tsuyoshi; Tomoda, Kentaro; Yasui, Yukinori

2018-04-01

Killing tensor fields have been thought of as describing the hidden symmetry of space(-time) since they are in one-to-one correspondence with polynomial first integrals of geodesic equations. Since many problems in classical mechanics can be formulated as geodesic problems in curved space and spacetime, solving the defining equation for Killing tensor fields (the Killing equation) is a powerful way to integrate equations of motion. Thus it has been desirable to formulate the integrability conditions of the Killing equation, which serve to determine the number of linearly independent solutions and also to restrict the possible forms of solutions tightly. In this paper, we show the prolongation for the Killing equation in a manner that uses Young symmetrizers. Using the prolonged equations, we provide the integrability conditions explicitly.

Modeling of Graphene Planar Grating in the THz Range by the Method of Singular Integral Equations

Science.gov (United States)

Kaliberda, Mstislav E.; Lytvynenko, Leonid M.; Pogarsky, Sergey A.

2018-04-01

Diffraction of the H-polarized electromagnetic wave by the planar graphene grating in the THz range is considered. The scattering and absorption characteristics are studied. The scattered field is represented in the spectral domain via unknown spectral function. The mathematical model is based on the graphene surface impedance and the method of singular integral equations. The numerical solution is obtained by the Nystrom-type method of discrete singularities.

Linear integral equations and soliton systems

International Nuclear Information System (INIS)

Quispel, G.R.W.

1983-01-01

A study is presented of classical integrable dynamical systems in one temporal and one spatial dimension. The direct linearizations are given of several nonlinear partial differential equations, for example the Korteweg-de Vries equation, the modified Korteweg-de Vries equation, the sine-Gordon equation, the nonlinear Schroedinger equation, and the equation of motion for the isotropic Heisenberg spin chain; the author also discusses several relations between these equations. The Baecklund transformations of these partial differential equations are treated on the basis of a singular transformation of the measure (or equivalently of the plane-wave factor) occurring in the corresponding linear integral equations, and the Baecklund transformations are used to derive the direct linearization of a chain of so-called modified partial differential equations. Finally it is shown that the singular linear integral equations lead in a natural way to the direct linearizations of various nonlinear difference-difference equations. (Auth.)

Three-dimensional integral equation approach to light scattering, extinction cross sections, local density of states, and quasi-normal modes

DEFF Research Database (Denmark)

de Lasson, Jakob Rosenkrantz; Mørk, Jesper; Kristensen, Philip Trøst

2013-01-01

We present a numerical formalism for solving the Lippmann–Schwinger equation for the electric field in three dimensions. The formalism may be applied to scatterers of different shapes and embedded in different background media, and we develop it in detail for the specific case of spherical scatte...

Remote monitoring of environmental particulate pollution - A problem in inversion of first-kind integral equations

Science.gov (United States)

Fymat, A. L.

1975-01-01

The determination of the microstructure, chemical nature, and dynamical evolution of scattering particulates in the atmosphere is considered. A description is given of indirect sampling techniques which can circumvent most of the difficulties associated with direct sampling techniques, taking into account methods based on scattering, extinction, and diffraction of an incident light beam. Approaches for reconstructing the particulate size distribution from the direct and the scattered radiation are discussed. A new method is proposed for determining the chemical composition of the particulates and attention is given to the relevance of methods of solution involving first kind Fredholm integral equations.

Geophysical interpretation using integral equations

CERN Document Server

Eskola, L

1992-01-01

Along with the general development of numerical methods in pure and applied to apply integral equations to geophysical modelling has sciences, the ability improved considerably within the last thirty years or so. This is due to the successful derivation of integral equations that are applicable to the modelling of complex structures, and efficient numerical algorithms for their solution. A significant stimulus for this development has been the advent of fast digital computers. The purpose of this book is to give an idea of the principles by which boundary-value problems describing geophysical models can be converted into integral equations. The end results are the integral formulas and integral equations that form the theoretical framework for practical applications. The details of mathematical analysis have been kept to a minimum. Numerical algorithms are discussed only in connection with some illustrative examples involving well-documented numerical modelling results. The reader is assu­ med to have a back...

A stochastic model of multiple scattering of charged particles: process, transport equation and solutions

International Nuclear Information System (INIS)

Papiez, L.; Moskvin, V.; Tulovsky, V.

2001-01-01

The process of angular-spatial evolution of multiple scattering of charged particles can be described by a special case of Boltzmann integro-differential equation called Lewis equation. The underlying stochastic process for this evolution is the compound Poisson process on the surface of the unit sphere. The significant portion of events that constitute compound Poisson process that describes multiple scattering have diffusional character. This property allows to analyze the process of angular-spatial evolution of multiple scattering of charged particles as combination of soft and hard collision processes and compute appropriately its transition densities. These computations provide a method of the approximate solution to the Lewis equation. (orig.)

Solution of the scattering T matrix equation in discrete complex momentum space

International Nuclear Information System (INIS)

Rawitscher, G.H.; Delic, G.

1984-01-01

The scattering solution to the Lippmann-Schwinger equation is expanded into a set of spherical Bessel functions of complex wave numbers, K/sub j/, with j = 1,2 , . . . , M. The value of each K/sub j/ is determined from the condition that the spherical Bessel function smoothly matches onto an asymptotically outgoing spherical Hankel (or Coulomb) function of the correct physical wave number at a matching point R. The spherical Bessel functions thus determined are Sturmian functions, and they form a complete set in the interval 0 to R. The coefficients of the expansion of the scattering function are determined by matrix inversion of a linear set of algebraic equations, which are equivalent to the solution of the T-matrix equation in complex momentum space. In view of the presence of a matching radius, no singularities are encountered for the Green's functions, and the inclusion of Coulomb potentials offers no computational difficulties. Three numerical examples are performed in order to illustrate the convergence of the elastic scattering matrix S with M. One of these consists of a set of coupled equations which describe the breakup of a deuteron as it scatters from the nucleus on 58 Ni. A value of M of 15 or less is found sufficient to reproduce the exact S matrix element to an accuracy of four figures after the decimal point

Simulating propagation of coherent light in random media using the Fredholm type integral equation

Science.gov (United States)

Kraszewski, Maciej; Pluciński, Jerzy

2017-06-01

Studying propagation of light in random scattering materials is important for both basic and applied research. Such studies often require usage of numerical method for simulating behavior of light beams in random media. However, if such simulations require consideration of coherence properties of light, they may become a complex numerical problems. There are well established methods for simulating multiple scattering of light (e.g. Radiative Transfer Theory and Monte Carlo methods) but they do not treat coherence properties of light directly. Some variations of these methods allows to predict behavior of coherent light but only for an averaged realization of the scattering medium. This limits their application in studying many physical phenomena connected to a specific distribution of scattering particles (e.g. laser speckle). In general, numerical simulation of coherent light propagation in a specific realization of random medium is a time- and memory-consuming problem. The goal of the presented research was to develop new efficient method for solving this problem. The method, presented in our earlier works, is based on solving the Fredholm type integral equation, which describes multiple light scattering process. This equation can be discretized and solved numerically using various algorithms e.g. by direct solving the corresponding linear equations system, as well as by using iterative or Monte Carlo solvers. Here we present recent development of this method including its comparison with well-known analytical results and a finite-difference type simulations. We also present extension of the method for problems of multiple scattering of a polarized light on large spherical particles that joins presented mathematical formalism with Mie theory.

Completely integrable operator evolution equations. II

International Nuclear Information System (INIS)

Chudnovsky, D.V.

1979-01-01

The author continues the investigation of operator classical completely integrable systems. The main attention is devoted to the stationary operator non-linear Schroedinger equation. It is shown that this equation can be used for separation of variables for a large class of completely integrable equations. (Auth.)

Integral transform method for solving neutron transport problems with general anisotropic scattering in a cylinder of finite height

International Nuclear Information System (INIS)

Kumar, V.; Sahni, D.C.

1983-01-01

In this paper, the authors present the mathematical techniques that were developed for solving the integral transport equation for the criticality of a homogeneous cylinder of finite height with general anisotropic scattering. They present the integral transport equations for the Fourier transformed spherical harmonic moments of the angular flux. These moments are also represented by a series of products of spherical Bessel functions. The criticality problem is, then, posed by the matrix eigenvalue problem whose eigenvector is composed of the expansion coefficients mentioned above. An methodology of calculating the general matrix element is discussed by using the recursion relations derived in this paper. Finally, for the one-group criticality of finite cylinders, the benchmark results are generated when scattering is linearly anisotropic. Also, these benchmarks are solved and compared with the S/sub N/ method of TWOTRAN

Integral Equation Methods for Electromagnetic and Elastic Waves

CERN Document Server

Chew, Weng; Hu, Bin

2008-01-01

Integral Equation Methods for Electromagnetic and Elastic Waves is an outgrowth of several years of work. There have been no recent books on integral equation methods. There are books written on integral equations, but either they have been around for a while, or they were written by mathematicians. Much of the knowledge in integral equation methods still resides in journal papers. With this book, important relevant knowledge for integral equations are consolidated in one place and researchers need only read the pertinent chapters in this book to gain important knowledge needed for integral eq

Application and development of the Schwinger multichannel scattering theory and the partial differential equation theory of electron-molecule scattering

Science.gov (United States)

Weatherford, Charles A.

1993-01-01

One version of the multichannel theory for electron-target scattering based on the Schwinger variational principle, the SMC method, requires the introduction of a projection parameter. The role of the projection parameter a is investigated and it is shown that the principal-value operator in the SMC equation is Hermitian regardless of the value of a as long as it is real and nonzero. In a basis that is properly orthonormalizable, the matrix representation of this operator is also Hermitian. The use of such basis is consistent with the Schwinger variational principle because the Lippmann-Schwinger equation automatically builds in the correct boundary conditions. Otherwise, an auxiliary condition needs to be introduced, and Takatsuka and McKoy's original value of a is one of the three possible ways to achieve Hermiticity. In all cases but one, a can be uncoupled from the Hermiticity condition and becomes a free parameter. An equation for a based on the variational stability of the scattering amplitude is derived; its solution has an interesting property that the scattering amplitude from a converged SMC calculation is independent of the choice of a even though the SMC operator itself is a-dependent. This property provides a sensitive test of the convergence of the calculation. For a static-exchange calculation, the convergence requirement only depends on the completeness of the one-electron basis, but for a general multichannel case, the a-invariance in the scattering amplitude requires both the one-electron basis and the N plus 1-electron basis to be complete. The role of a in the SMC equation and the convergence property are illustrated using two examples: e-CO elastic scattering in the static-exchange approximation, and a two-state treatment of the e-H2 Chi(sup 1)Sigma(sub g)(+) yields b(sup 3)Sigma(sub u)(+) excitation.

Integrable dissipative nonlinear second order differential equations via factorizations and Abel equations

Energy Technology Data Exchange (ETDEWEB)

Mancas, Stefan C. [Department of Mathematics, Embry–Riddle Aeronautical University, Daytona Beach, FL 32114-3900 (United States); Rosu, Haret C., E-mail: hcr@ipicyt.edu.mx [IPICYT, Instituto Potosino de Investigacion Cientifica y Tecnologica, Apdo Postal 3-74 Tangamanga, 78231 San Luis Potosí, SLP (Mexico)

2013-09-02

We emphasize two connections, one well known and another less known, between the dissipative nonlinear second order differential equations and the Abel equations which in their first-kind form have only cubic and quadratic terms. Then, employing an old integrability criterion due to Chiellini, we introduce the corresponding integrable dissipative equations. For illustration, we present the cases of some integrable dissipative Fisher, nonlinear pendulum, and Burgers–Huxley type equations which are obtained in this way and can be of interest in applications. We also show how to obtain Abel solutions directly from the factorization of second order nonlinear equations.

An integral transform of the Salpeter equation

International Nuclear Information System (INIS)

Krolikowski, W.

1980-03-01

We find a new form of relativistic wave equation for two spin-1/2 particles, which arises by an integral transformation (in the position space) of the wave function in the Salpeter equation. The non-locality involved in this transformation is extended practically over the Compton wavelength of the lighter of two particles. In the case of equal masses the new equation assumes the form of the Breit equation with an effective integral interaction. In the one-body limit it reduces to the Dirac equation also with an effective integral interaction. (author)

B-splines and Faddeev equations

International Nuclear Information System (INIS)

Huizing, A.J.

1990-01-01

Two numerical methods for solving the three-body equations describing relativistic pion deuteron scattering have been investigated. For separable two body interactions these equations form a set of coupled one-dimensional integral equations. They are plagued by singularities which occur in the kernel of the integral equations as well as in the solution. The methods to solve these equations differ in the way they treat the singularities. First the Fuda-Stuivenberg method is discussed. The basic idea of this method is an one time iteration of the set of integral equations to treat the logarithmic singularities. In the second method, the spline method, the unknown solution is approximated by splines. Cubic splines have been used with cubic B-splines as basis. If the solution is approximated by a linear combination of basis functions, an integral equation can be transformed into a set of linear equations for the expansion coefficients. This set of linear equations is solved by standard means. Splines are determined by points called knots. A proper choice of splines to approach the solution stands for a proper choice of the knots. The solution of the three-body scattering equations has a square root behaviour at a certain point. Hence it was investigated how the knots should be chosen to approximate the square root function by cubic B-splines in an optimal way. Before applying this method to solve numerically the three-body equations describing pion-deuteron scattering, an analytically solvable example has been constructed with a singularity structure of both kernel and solution comparable to those of the three-body equations. The accuracy of the numerical solution was determined to a large extent by the accuracy of the approximation of the square root part. The results for a pion laboratory energy of 47.4 MeV agree very well with those from literature. In a complete calculation for 47.7 MeV the spline method turned out to be a factor thousand faster than the Fuda

Solving Non-Isospectral mKdV Equation and Sine-Gordon Equation Hierarchies with Self-Consistent Sources via Inverse Scattering Transform

International Nuclear Information System (INIS)

Li Qi; Zhang Dajun; Chen Dengyuan

2010-01-01

N-soliton solutions of the hierarchy of non-isospectral mKdV equation with self-consistent sources and the hierarchy of non-isospectral sine-Gordon equation with self-consistent sources are obtained via the inverse scattering transform. (general)

Semi-analytic equations to the Cox-Thompson inverse scattering method at fixed energy for special cases

International Nuclear Information System (INIS)

Palmai, T.; Apagyi, B.; Horvath, M.

2008-01-01

Solution of the Cox-Thompson inverse scattering problem at fixed energy 1-3 is reformulated resulting in semi-analytic equations. The new set of equations for the normalization constants and the nonphysical (shifted) angular momenta are free of matrix inversion operations. This simplification is a result of treating only the input phase shifts of partial waves of a given parity. Therefore, the proposed method can be applied for identical particle scattering of the bosonic type (or for certain cases of identical fermionic scattering). The new formulae are expected to be numerically more efficient than the previous ones. Based on the semi-analytic equations an approximate method is proposed for the generic inverse scattering problem, when partial waves of arbitrary parity are considered. (author)

Scattering of quantized solitary waves in the cubic Schrodinger equation

International Nuclear Information System (INIS)

Dolan, L.

1976-01-01

The quantum mechanics for N particles interacting via a delta-function potential in one space dimension and one time dimension is known. The second-quantized description of this system has for its Euler-Lagrange equations of motion the cubic Schrodinger equation. This nonlinear differential equation supports solitary wave solutions. A quantization of these solitons reproduces the weak-coupling limit to the known quantum mechanics. The phase shift for two-body scattering and the energy of the N-body bound state is derived in this approximation. The nonlinear Schrodinger equation is contrasted with the sine-Gordon theory in respect to the ideas which the classical solutions play in the description of the quantum states

Group theoretical and Hamiltonian structures of integrable evolution equations in 1x1 and 2x1 dimensions

International Nuclear Information System (INIS)

Konopel'chenko, B.G.

1983-01-01

New results in investigation of the group-theoretical and hamiltonian structure of the integrable evolution equations in 1+1 and 2+1 dimensions are briefly reviewed. Main general results, such as the form of integrable equations, Baecklund transfomations, symmetry groups, are turned out to have the same form for different spectral problems. The used generalized AKNS-method (the Ablowitz Kaup, Newell and Segur method) permits to prove that all nonlinear evolution equations considered are hamiltonians. The general condition of effective application of the ACNS mehtod to the concrete spectral problem is the possibility to calculate a recursion operator explicitly. The embedded representation is shown to be a fundamental object connected with different aspects of the inverse scattering problem

Completely integrable operator evolutionary equations

International Nuclear Information System (INIS)

Chudnovsky, D.V.

1979-01-01

The authors present natural generalizations of classical completely integrable equations where the functions are replaced by arbitrary operators. Among these equations are the non-linear Schroedinger, the Korteweg-de Vries, and the modified KdV equations. The Lax representation and the Baecklund transformations are presented. (Auth.)

Interstitial integrals in the multiple-scattering model

International Nuclear Information System (INIS)

Swanson, J.R.; Dill, D.

1982-01-01

We present an efficient method for the evaluation of integrals involving multiple-scattering wave functions over the interstitial region. Transformation of the multicenter interstitial wave functions to a single center representation followed by a geometric projection reduces the integrals to products of analytic angular integrals and numerical radial integrals. The projection function, which has the value 1 in the interstitial region and 0 elsewhere, has a closed-form partial-wave expansion. The method is tested by comparing its results with exact normalization and dipole integrals; the differences are 2% at worst and typically less than 1%. By providing an efficient means of calculating Coulomb integrals, the method allows treatment of electron correlations using a multiple scattering basis set

An explicit marching on-in-time solver for the time domain volume magnetic field integral equation

KAUST Repository

Sayed, Sadeed Bin

2014-07-01

Transient scattering from inhomogeneous dielectric objects can be modeled using time domain volume integral equations (TDVIEs). TDVIEs are oftentimes solved using marching on-in-time (MOT) techniques. Classical MOT-TDVIE solvers expand the field induced on the scatterer using local spatio-temporal basis functions. Inserting this expansion into the TDVIE and testing the resulting equation in space and time yields a system of equations that is solved by time marching. Depending on the type of the basis and testing functions and the time step, the time marching scheme can be implicit (N. T. Gres, et al., Radio Sci., 36(3), 379-386, 2001) or explicit (A. Al-Jarro, et al., IEEE Trans. Antennas Propag., 60(11), 5203-5214, 2012). Implicit MOT schemes are known to be more stable and accurate. However, under low-frequency excitation, i.e., when the time step size is large, they call for inversion of a full matrix system at very time step.

An explicit marching on-in-time solver for the time domain volume magnetic field integral equation

KAUST Repository

Sayed, Sadeed Bin; Ulku, Huseyin Arda; Bagci, Hakan

2014-01-01

Transient scattering from inhomogeneous dielectric objects can be modeled using time domain volume integral equations (TDVIEs). TDVIEs are oftentimes solved using marching on-in-time (MOT) techniques. Classical MOT-TDVIE solvers expand the field induced on the scatterer using local spatio-temporal basis functions. Inserting this expansion into the TDVIE and testing the resulting equation in space and time yields a system of equations that is solved by time marching. Depending on the type of the basis and testing functions and the time step, the time marching scheme can be implicit (N. T. Gres, et al., Radio Sci., 36(3), 379-386, 2001) or explicit (A. Al-Jarro, et al., IEEE Trans. Antennas Propag., 60(11), 5203-5214, 2012). Implicit MOT schemes are known to be more stable and accurate. However, under low-frequency excitation, i.e., when the time step size is large, they call for inversion of a full matrix system at very time step.

A time-domain finite element boundary integral approach for elastic wave scattering

Science.gov (United States)

Shi, F.; Lowe, M. J. S.; Skelton, E. A.; Craster, R. V.

2018-04-01

The response of complex scatterers, such as rough or branched cracks, to incident elastic waves is required in many areas of industrial importance such as those in non-destructive evaluation and related fields; we develop an approach to generate accurate and rapid simulations. To achieve this we develop, in the time domain, an implementation to efficiently couple the finite element (FE) method within a small local region, and the boundary integral (BI) globally. The FE explicit scheme is run in a local box to compute the surface displacement of the scatterer, by giving forcing signals to excitation nodes, which can lie on the scatterer itself. The required input forces on the excitation nodes are obtained with a reformulated FE equation, according to the incident displacement field. The surface displacements computed by the local FE are then projected, through time-domain BI formulae, to calculate the scattering signals with different modes. This new method yields huge improvements in the efficiency of FE simulations for scattering from complex scatterers. We present results using different shapes and boundary conditions, all simulated using this approach in both 2D and 3D, and then compare with full FE models and theoretical solutions to demonstrate the efficiency and accuracy of this numerical approach.

Roy-Steiner-equation analysis of pion-nucleon scattering

Science.gov (United States)

Meißner, U.-G.; Ruiz de Elvira, J.; Hoferichter, M.; Kubis, B.

2017-03-01

Low-energy pion-nucleon scattering is relevant for many areas in nuclear and hadronic physics, ranging from the scalar couplings of the nucleon to the long-range part of two-pion-exchange potentials and three-nucleon forces in Chiral Effective Field Theory. In this talk, we show how the fruitful combination of dispersion-theoretical methods, in particular in the form of Roy-Steiner equations, with modern high-precision data on hadronic atoms allows one to determine the pion-nucleon scattering amplitudes at low energies with unprecedented accuracy. Special attention will be paid to the extraction of the pion-nucleon σ-term, and we discuss in detail the current tension with recent lattice results, as well as the determination of the low-energy constants of chiral perturbation theory.

Asymptotic integration of differential and difference equations

CERN Document Server

Bodine, Sigrun

2015-01-01

This book presents the theory of asymptotic integration for both linear differential and difference equations. This type of asymptotic analysis is based on some fundamental principles by Norman Levinson. While he applied them to a special class of differential equations, subsequent work has shown that the same principles lead to asymptotic results for much wider classes of differential and also difference equations. After discussing asymptotic integration in a unified approach, this book studies how the application of these methods provides several new insights and frequent improvements to results found in earlier literature. It then continues with a brief introduction to the relatively new field of asymptotic integration for dynamic equations on time scales. Asymptotic Integration of Differential and Difference Equations is a self-contained and clearly structured presentation of some of the most important results in asymptotic integration and the techniques used in this field. It will appeal to researchers i...

Asymptotic formulae for solutions of the two-group integral neutron-transport equation

International Nuclear Information System (INIS)

Duracz, T.

1976-01-01

The steady-state, two-group integral neutron-transport equation is considered for two cases. First, for plane geometry, formulae for the asymptotic flux are obtained, under assumptions of homogeneous medium with isotropic scattering, extended to infinity (whole space and half-space), with sources vanishing at infinity as 0(esup(-IXI)). Next, for spherical geometry, the Milne problem is considered and formulae for the asymptotic flux are obtained. These formulae have the form of asymptotic expansions for small and large radii of the black sphere. (orig.) [de

Quantum-mechanical few-body scattering equations with half-on-shell energy-independent subsystem input

International Nuclear Information System (INIS)

Zeiger, E.M.

1978-01-01

New equations are presented for three- and four-body scattering, within the context of nonrelativistic quantum mechanics and a Hamiltonian scattering theory. For the three-body case Faddeev-type equations are presented which, although obtained from the rigorous Faddeev theory, only require two-body bound state wave functions and half-off-shell transition amplitudes as input. In addition, their effective potentials are independent of the three-body energy, and can easily be made real after an angular momentum decomposition. The equations are formulated in terms of physical transition amplitudes for three-body processes, except that in the breakup case the partial-wave amplitudes differ from the corresponding full amplitudes by a Watson final-state-interaction factor. Also presented are new equations for four-body scattering, obtained by generalizing our three-body formalism to the four-body case. These equations, although equivalent to those of Faddeev--Yakubovskii, are expressed in terms of singularity-free transition amplitudes, and their energy-independent effective potentials require only half-on-shell subsystem transition amplitudes (and bound state wave functions) as input. However, due to the detailed index structure of the Faddeev--Yakubovskii formalsim, the result of the generalization is considerably more complicated than in the three-body case

Roy–Steiner-equation analysis of pion–nucleon scattering

Directory of Open Access Journals (Sweden)

Meißner U.-G.

2017-01-01

Full Text Available Low-energy pion–nucleon scattering is relevant for many areas in nuclear and hadronic physics, ranging from the scalar couplings of the nucleon to the long-range part of two-pion-exchange potentials and three-nucleon forces in Chiral Effective Field Theory. In this talk, we show how the fruitful combination of dispersion-theoretical methods, in particular in the form of Roy–Steiner equations, with modern high-precision data on hadronic atoms allows one to determine the pion–nucleon scattering amplitudes at low energies with unprecedented accuracy. Special attention will be paid to the extraction of the pion–nucleon σ-term, and we discuss in detail the current tension with recent lattice results, as well as the determination of the low-energy constants of chiral perturbation theory.c

Diffusion equations and hard collisions in multiple scattering of charged particles

International Nuclear Information System (INIS)

Papiez, Lech; Tulovsky, Vladimir

1998-01-01

The processes of angular-spatial evolution of multiple scattering of charged particles are described by the Lewis (special case of Boltzmann) integro-differential equation. The underlying stochastic process for this evolution is the compound Poisson process with transition densities satisfying the Lewis equation. In this paper we derive the Lewis equation from the compound Poisson process and show that the effective method of the solution of this equation can be based on the idea of decomposition of the compound Poisson process into processes of soft and hard collisions. Formulas for transition densities of soft and hard collision processes are provided in this paper together with the formula expressing the general solution of the Lewis equation in terms of those transition densities

Diffusion equations and hard collisions in multiple scattering of charged particles

Energy Technology Data Exchange (ETDEWEB)

Papiez, Lech [Department of Radiation Oncology, Indiana University, Indianapolis, IN (United States); Tulovsky, Vladimir [Department of Mathematics, St. John' s College, Staten Island, New York, NY (United States)

1998-09-01

The processes of angular-spatial evolution of multiple scattering of charged particles are described by the Lewis (special case of Boltzmann) integro-differential equation. The underlying stochastic process for this evolution is the compound Poisson process with transition densities satisfying the Lewis equation. In this paper we derive the Lewis equation from the compound Poisson process and show that the effective method of the solution of this equation can be based on the idea of decomposition of the compound Poisson process into processes of soft and hard collisions. Formulas for transition densities of soft and hard collision processes are provided in this paper together with the formula expressing the general solution of the Lewis equation in terms of those transition densities.

Integral Parameters of the Thermal Neutron Scattering Law

International Nuclear Information System (INIS)

Purohit, S.N.

1964-09-01

Integral parameters of the thermal neutron scattering law - the thermalization binding parameter (M 2 ), the Placzek's moments of the generalized frequency spectrum of dynamical modes and the energy transfer moments of the scattering law - are theoretically discussed. A detailed study of the variation of M 2 , the thermalization time constant and the effective temperature of the vibrating atoms, with the relative weight between intra-molecular vibrations and hindered rotations for H 2 O, is presented. Theoretical results for different scattering models of H 2 O are compared with the measurements of integral experiments. A set of integral parameters for D 2 O, using Butler's model, have been obtained. Importance of the structure of hindered rotations of H 2 O and D 2 O in the study of integral parameters has also been discussed

Darboux invariants of integrable equations with variable spectral parameters

International Nuclear Information System (INIS)

Shin, H J

2008-01-01

The Darboux transformation for integrable equations with variable spectral parameters is introduced. Darboux invariant quantities are calculated, which are used in constructing the Lax pair of integrable equations. This approach serves as a systematic method for constructing inhomogeneous integrable equations and their soliton solutions. The structure functions of variable spectral parameters determine the integrability and nonlinear coupling terms. Three cases of integrable equations are treated as examples of this approach

An efficient fringe integral equation method for optimizing the antenna location on complex bodies

DEFF Research Database (Denmark)

Jørgensen, Erik; Meincke, Peter; Breinbjerg, Olav

2001-01-01

The radiation pattern of an antenna mounted nearby, or directly on, a complex three-dimensional (3D) structure can be significantly influenced by this structure. Integral equations combined with the method of moments (MoM) provide an accurate means for calculating the scattering from the structures...... in such applications. The structure is then modelled by triangular or rectangular surface patches with corresponding surface current expansion functions. A MoM matrix which is independent of the antenna location can be obtained by modelling the antenna as an impressed electric or magnetic source, e.g., a slot antenna...... can be modelled by a magnetic Hertzian dipole. For flush-mounted antennas, or antennas mounted in close vicinity of the scattering structure, the nearby impressed source induces a highly peaked surface current on the scattering structure. For the low-order basis functions usually applied...

New multidimensional partially integrable generalization of S-integrable N-wave equation

International Nuclear Information System (INIS)

Zenchuk, A. I.

2007-01-01

This paper develops a modification of the dressing method based on the inhomogeneous linear integral equation with integral operator having nonempty kernel. The method allows one to construct the systems of multidimensional partial differential equations having differential polynomial structure in any dimension n. The associated solution space is not full, although it is parametrized by certain number of arbitrary functions of (n-1) variables. We consider four-dimensional generalization of the classical (2+1)-dimensional S-integrable N-wave equation as an example

Integrable discretizations of the short pulse equation

International Nuclear Information System (INIS)

Feng Baofeng; Maruno, Ken-ichi; Ohta, Yasuhiro

2010-01-01

In this paper, we propose integrable semi-discrete and full-discrete analogues of the short pulse (SP) equation. The key construction is the bilinear form and determinant structure of solutions of the SP equation. We also give the determinant formulas of N-soliton solutions of the semi-discrete and full-discrete analogues of the SP equations, from which the multi-loop and multi-breather solutions can be generated. In the continuous limit, the full-discrete SP equation converges to the semi-discrete SP equation, and then to the continuous SP equation. Based on the semi-discrete SP equation, an integrable numerical scheme, i.e. a self-adaptive moving mesh scheme, is proposed and used for the numerical computation of the short pulse equation.

Integrable discretization s of derivative nonlinear Schroedinger equations

International Nuclear Information System (INIS)

Tsuchida, Takayuki

2002-01-01

We propose integrable discretizations of derivative nonlinear Schroedinger (DNLS) equations such as the Kaup-Newell equation, the Chen-Lee-Liu equation and the Gerdjikov-Ivanov equation by constructing Lax pairs. The discrete DNLS systems admit the reduction of complex conjugation between two dependent variables and possess bi-Hamiltonian structure. Through transformations of variables and reductions, we obtain novel integrable discretizations of the nonlinear Schroedinger (NLS), modified KdV (mKdV), mixed NLS, matrix NLS, matrix KdV, matrix mKdV, coupled NLS, coupled Hirota, coupled Sasa-Satsuma and Burgers equations. We also discuss integrable discretizations of the sine-Gordon equation, the massive Thirring model and their generalizations. (author)

Integral Parameters of the Thermal Neutron Scattering Law

Energy Technology Data Exchange (ETDEWEB)

Purohit, S N

1964-09-15

Integral parameters of the thermal neutron scattering law - the thermalization binding parameter (M{sub 2}), the Placzek's moments of the generalized frequency spectrum of dynamical modes and the energy transfer moments of the scattering law - are theoretically discussed. A detailed study of the variation of M{sub 2}, the thermalization time constant and the effective temperature of the vibrating atoms, with the relative weight between intra-molecular vibrations and hindered rotations for H{sub 2}O, is presented. Theoretical results for different scattering models of H{sub 2}O are compared with the measurements of integral experiments. A set of integral parameters for D{sub 2}O, using Butler's model, have been obtained. Importance of the structure of hindered rotations of H{sub 2}O and D{sub 2}O in the study of integral parameters has also been discussed.

Integrable peakon equations with cubic nonlinearity

International Nuclear Information System (INIS)

Hone, Andrew N W; Wang, J P

2008-01-01

We present a new integrable partial differential equation found by Vladimir Novikov. Like the Camassa-Holm and Degasperis-Procesi equations, this new equation admits peaked soliton (peakon) solutions, but it has nonlinear terms that are cubic, rather than quadratic. We give a matrix Lax pair for V Novikov's equation, and show how it is related by a reciprocal transformation to a negative flow in the Sawada-Kotera hierarchy. Infinitely many conserved quantities are found, as well as a bi-Hamiltonian structure. The latter is used to obtain the Hamiltonian form of the finite-dimensional system for the interaction of N peakons, and the two-body dynamics (N = 2) is explicitly integrated. Finally, all of this is compared with some analogous results for another cubic peakon equation derived by Zhijun Qiao. (fast track communication)

Direct Calculation of the Scattering Amplitude Without Partial Wave Analysis

Science.gov (United States)

Shertzer, J.; Temkin, A.; Fisher, Richard R. (Technical Monitor)

2001-01-01

Two new developments in scattering theory are reported. We show, in a practical way, how one can calculate the full scattering amplitude without invoking a partial wave expansion. First, the integral expression for the scattering amplitude f(theta) is simplified by an analytic integration over the azimuthal angle. Second, the full scattering wavefunction which appears in the integral expression for f(theta) is obtained by solving the Schrodinger equation with the finite element method (FEM). As an example, we calculate electron scattering from the Hartree potential. With minimal computational effort, we obtain accurate and stable results for the scattering amplitude.

On the method of inverse scattering problem and Baecklund transformations for supersymmetric equations

International Nuclear Information System (INIS)

Chaichian, M.; Kulish, P. P.

1978-04-01

Supersymmetric Liouville and sine-Gordon equations are studied. We write down for these models the system of linear equations for which the method of inverse scattering problem should be applicable. Expressions for an infinite set of conserved currents are explicitly given. Supersymmetric Baecklund transformations and generalized conservation laws are constructed. (author)

Picard-Fuchs equations of dimensionally regulated Feynman integrals

Energy Technology Data Exchange (ETDEWEB)

Zayadeh, Raphael

2013-12-15

This thesis is devoted to studying differential equations of Feynman integrals. A Feynman integral depends on a dimension D. For integer values of D it can be written as a projective integral, which is called the Feynman parameter prescription. A major complication arises from the fact that for some values of D the integral can diverge. This problem is solved within dimensional regularization by continuing the integral as a meromorphic function on the complex plane and replacing the ill-defined quantity by a Laurent series in a dimensional regularization parameter. All terms in such a Laurent expansion are periods in the sense of Kontsevich and Zagier. We describe a new method to compute differential equations of Feynman integrals. So far, the standard has been to use integration-by-parts (IBP) identities to obtain coupled systems of linear differential equations for the master integrals. Our method is based on the theory of Picard-Fuchs equations. In the case we are interested in, that of projective and quasiprojective families, a Picard-Fuchs equation can be computed by means of the Griffiths-Dwork reduction. We describe a method that is designed for fixed integer dimension. After a suitable integer shift of dimension we obtain a period of a family of hypersurfaces, hence a Picard-Fuchs equation. This equation is inhomogeneous because the domain of integration has a boundary and we only obtain a relative cycle. As a second step we shift back the dimension using Tarasov's generalized dimensional recurrence relations. Furthermore, we describe a method to directly compute the differential equation for general D without shifting the dimension. This is based on the Griffiths-Dwork reduction. The success of this method depends on the ability to solve large systems of linear equations. We give examples of two and three-loop graphs. Tarasov classifies two-loop two-point functions and we give differential equations for these. For us the most interesting example is

Picard-Fuchs equations of dimensionally regulated Feynman integrals

International Nuclear Information System (INIS)

Zayadeh, Raphael

2013-12-01

This thesis is devoted to studying differential equations of Feynman integrals. A Feynman integral depends on a dimension D. For integer values of D it can be written as a projective integral, which is called the Feynman parameter prescription. A major complication arises from the fact that for some values of D the integral can diverge. This problem is solved within dimensional regularization by continuing the integral as a meromorphic function on the complex plane and replacing the ill-defined quantity by a Laurent series in a dimensional regularization parameter. All terms in such a Laurent expansion are periods in the sense of Kontsevich and Zagier. We describe a new method to compute differential equations of Feynman integrals. So far, the standard has been to use integration-by-parts (IBP) identities to obtain coupled systems of linear differential equations for the master integrals. Our method is based on the theory of Picard-Fuchs equations. In the case we are interested in, that of projective and quasiprojective families, a Picard-Fuchs equation can be computed by means of the Griffiths-Dwork reduction. We describe a method that is designed for fixed integer dimension. After a suitable integer shift of dimension we obtain a period of a family of hypersurfaces, hence a Picard-Fuchs equation. This equation is inhomogeneous because the domain of integration has a boundary and we only obtain a relative cycle. As a second step we shift back the dimension using Tarasov's generalized dimensional recurrence relations. Furthermore, we describe a method to directly compute the differential equation for general D without shifting the dimension. This is based on the Griffiths-Dwork reduction. The success of this method depends on the ability to solve large systems of linear equations. We give examples of two and three-loop graphs. Tarasov classifies two-loop two-point functions and we give differential equations for these. For us the most interesting example is the two

Pion–nucleon scattering: from chiral perturbation theory to Roy–Steiner equations

International Nuclear Information System (INIS)

Kubis, Bastian; Hoferichter, Martin; Elvira, Jacobo Ruiz de; Meißner, Ulf-G.

2016-01-01

Ever since Weinberg’s seminal predictions of the pion–nucleon scattering amplitudes at threshold, this process has been of central interest for the study of chiral dynamics involving nucleons. The scattering lengths or the pion–nucleon σ-term are fundamental quantities characterizing the explicit breaking of chiral symmetry by means of the light quark masses. On the other hand, pion–nucleon dynamics also strongly affects the long-range part of nucleon–nucleon potentials, and hence has a far-reaching impact on nuclear physics. We discuss the fruitful combination of dispersion-theoretical methods, in the form of Roy–Steiner equations, with chiral dynamics to determine pion–nucleon scattering amplitudes at low energies with high precision.

Pion-nucleon scattering: from chiral perturbation theory to Roy-Steiner equations

Science.gov (United States)

Kubis, Bastian; Hoferichter, Martin; de Elvira, Jacobo Ruiz; Meißner, Ulf-G.

2016-11-01

Ever since Weinberg's seminal predictions of the pion-nucleon scattering amplitudes at threshold, this process has been of central interest for the study of chiral dynamics involving nucleons. The scattering lengths or the pion-nucleon σ-term are fundamental quantities characterizing the explicit breaking of chiral symmetry by means of the light quark masses. On the other hand, pion-nucleon dynamics also strongly affects the long-range part of nucleon-nucleon potentials, and hence has a far-reaching impact on nuclear physics. We discuss the fruitful combination of dispersion-theoretical methods, in the form of Roy-Steiner equations, with chiral dynamics to determine pion-nucleon scattering amplitudes at low energies with high precision.*

Scattering for wave equations with dissipative terms in layered media

Directory of Open Access Journals (Sweden)

Mitsuteru Kadowaki

2011-05-01

Full Text Available In this article, we show the existence of scattering solutions to wave equations with dissipative terms in layered media. To analyze the wave propagation in layered media, it is necessary to handle singular points called thresholds in the spectrum. Our main tools are Kato's smooth perturbation theory and some approximate operators.

Scattering state solutions of the Duffin-Kemmer-Petiau equation with the Varshni potential model

Energy Technology Data Exchange (ETDEWEB)

Oluwadare, O.J. [Federal University Oye-Ekiti, Department of Physics, Oye-Ekiti, Ekiti State (Nigeria); Oyewumi, K.J. [Federal University of Technology, Department of Physics, Minna, Niger State (Nigeria)

2017-02-15

The scattering state of the Duffin-Kemmer-Petiau equation with the Varshni potential was studied. The asymptotic wave function, the scattering phase shift and normalization constant were obtained for any J states by dealing with the centrifugal term using a suitable approximation. The analytical properties of the scattering amplitude and the bound state energy were obtained and discussed. Our numerical and graphical results indicate that the scattering phase shift depends largely on total angular momentum J, screening parameter β and potential strengths a and b. (orig.)

Some New Integrable Equations from the Self-Dual Yang-Mills Equations

International Nuclear Information System (INIS)

Ivanova, T.A.; Popov, A.D.

1994-01-01

Using the symmetry reductions of the self-dual Yang-Mills (SDYM) equations in (2+2) dimensions, we introduce new integrable equations which are 'deformations' of the chiral model in (2+1) dimensions, generalized nonlinear Schroedinger, Korteweg-de Vries, Toda lattice, Garnier, Euler-Arnold, generalized Calogero-Moser and Euler-Calogero-Moser equations. The Lax pairs for all of these equations are derived by the symmetry reductions of the Lax pair for the SDYM equations. 34 refs

Fractional Bhatnagar-Gross-Krook kinetic equation

Science.gov (United States)

Goychuk, Igor

2017-11-01

The linear Boltzmann equation (LBE) approach is generalized to describe fractional superdiffusive transport of the Lévy walk type in external force fields. The time distribution between scattering events is assumed to have a finite mean value and infinite variance. It is completely characterized by the two scattering rates, one fractional and a normal one, which defines also the mean scattering rate. We formulate a general fractional LBE approach and exemplify it with a particularly simple case of the Bohm and Gross scattering integral leading to a fractional generalization of the Bhatnagar, Gross and Krook (BGK) kinetic equation. Here, at each scattering event the particle velocity is completely randomized and takes a value from equilibrium Maxwell distribution at a given fixed temperature. We show that the retardation effects are indispensable even in the limit of infinite mean scattering rate and argue that this novel fractional kinetic equation provides a viable alternative to the fractional Kramers-Fokker-Planck (KFP) equation by Barkai and Silbey and its generalization by Friedrich et al. based on the picture of divergent mean time between scattering events. The case of divergent mean time is also discussed at length and compared with the earlier results obtained within the fractional KFP. Also a phenomenological fractional BGK equation without retardation effects is proposed in the limit of infinite scattering rates. It cannot be, however, rigorously derived from a scattering model, being rather clever postulated. It this respect, this retardationless equation is similar to the fractional KFP by Barkai and Silbey. However, it corresponds to the opposite, much more physical limit and, therefore, also presents a viable alternative.

Kaon-nucleon scattering in three-dimensional technique

International Nuclear Information System (INIS)

Salam, Agus; Fachruddin, Imam

2016-01-01

Kaon-nucleon (KN) scattering is formulated in the three-dimensional (3D) momentum space, in which the basis state is not expanded into partial waves. Based on this basis the Lippmann-Schwinger equation for the T-matrix is evaluated. We obtain as final equation for the T-matrix elements a set of two coupled integral equations in two variables, which are the momentum’s magnitude and the scattering angle. Calculations for the differential cross section and some spin observables are shown, for which we employ a hadrons exchange model with the second order contributions only.

Kaon-nucleon scattering in three-dimensional technique

Energy Technology Data Exchange (ETDEWEB)

Salam, Agus, E-mail: agus.salam@sci.ui.ac.id; Fachruddin, Imam [Departemen Fisika, FMIPA, Universitas Indonesia, Depok 16424 (Indonesia)

2016-03-11

Kaon-nucleon (KN) scattering is formulated in the three-dimensional (3D) momentum space, in which the basis state is not expanded into partial waves. Based on this basis the Lippmann-Schwinger equation for the T-matrix is evaluated. We obtain as final equation for the T-matrix elements a set of two coupled integral equations in two variables, which are the momentum’s magnitude and the scattering angle. Calculations for the differential cross section and some spin observables are shown, for which we employ a hadrons exchange model with the second order contributions only.

Calculation of the Full Scattering Amplitude without Partial Wave Decomposition II

Science.gov (United States)

Shertzer, J.; Temkin, A.

2003-01-01

As is well known, the full scattering amplitude can be expressed as an integral involving the complete scattering wave function. We have shown that the integral can be simplified and used in a practical way. Initial application to electron-hydrogen scattering without exchange was highly successful. The Schrodinger equation (SE) can be reduced to a 2d partial differential equation (pde), and was solved using the finite element method. We have now included exchange by solving the resultant SE, in the static exchange approximation. The resultant equation can be reduced to a pair of coupled pde's, to which the finite element method can still be applied. The resultant scattering amplitudes, both singlet and triplet, as a function of angle can be calculated for various energies. The results are in excellent agreement with converged partial wave results.

Construction of the Cauchy problem for solution of the integral equation of transport in a sphere with a central absorbing region

International Nuclear Information System (INIS)

Ezhov, A.A.

1978-01-01

On the basis of the integral equation for neutron transport in a homogeneous isotropically-scattering sphere with an absolutely black central part an initial value problem has been formulated which permits the construction of a numerical scheme to find the neutron flux density

Baecklund transformations for integrable lattice equations

International Nuclear Information System (INIS)

Atkinson, James

2008-01-01

We give new Baecklund transformations (BTs) for some known integrable (in the sense of being multidimensionally consistent) quadrilateral lattice equations. As opposed to the natural auto-BT inherent in every such equation, these BTs are of two other kinds. Specifically, it is found that some equations admit additional auto-BTs (with Baecklund parameter), whilst some pairs of apparently distinct equations admit a BT which connects them

PREFACE: Symmetries and Integrability of Difference Equations

Science.gov (United States)

Doliwa, Adam; Korhonen, Risto; Lafortune, Stéphane

2007-10-01

The notion of integrability was first introduced in the 19th century in the context of classical mechanics with the definition of Liouville integrability for Hamiltonian flows. Since then, several notions of integrability have been introduced for partial and ordinary differential equations. Closely related to integrability theory is the symmetry analysis of nonlinear evolution equations. Symmetry analysis takes advantage of the Lie group structure of a given equation to study its properties. Together, integrability theory and symmetry analysis provide the main method by which nonlinear evolution equations can be solved explicitly. Difference equations (DE), like differential equations, are important in numerous fields of science and have a wide variety of applications in such areas as mathematical physics, computer visualization, numerical analysis, mathematical biology, economics, combinatorics, and quantum field theory. It is thus crucial to develop tools to study and solve DEs. While the theory of symmetry and integrability for differential equations is now largely well-established, this is not yet the case for discrete equations. Although over recent years there has been significant progress in the development of a complete analytic theory of difference equations, further tools are still needed to fully understand, for instance, the symmetries, asymptotics and the singularity structure of difference equations. The series of SIDE meetings on Symmetries and Integrability of Difference Equations started in 1994. Its goal is to provide a platform for an international and interdisciplinary communication for researchers working in areas associated with integrable discrete systems, such as classical and quantum physics, computer science and numerical analysis, mathematical biology and economics, discrete geometry and combinatorics, theory of special functions, etc. The previous SIDE meetings took place in Estérel near Montréal, Canada (1994), at the University of

Coupled channel theory of pion--deuteron reaction applied to threshold scattering

International Nuclear Information System (INIS)

Mizutani, T.; Koltun, D.S.

1977-01-01

Scattering and absorption of pions by a nuclear target are treated together in a coupled channel theory. The theory is developed explicitly for the problem of pion scattering and absorption by a deuteron. The equations are presented in terms of the integral equations of three-body scattering theory. The method is then applied in an approximate from to calculate the contribution of pion absorption to the scattering length for pion--deuteron scattering. The sensitivity of the calculated results to the model assumptions and approximations is investigated

Counting master integrals. Integration by parts vs. functional equations

International Nuclear Information System (INIS)

Kniehl, Bernd A.; Tarasov, Oleg V.

2016-01-01

We illustrate the usefulness of functional equations in establishing relationships between master integrals under the integration-by-parts reduction procedure by considering a certain two-loop propagator-type diagram as an example.

Ambitwistor strings and the scattering equations at one loop

Energy Technology Data Exchange (ETDEWEB)

Adamo, Tim; Casali, Eduardo; Skinner, David [Department of Applied Mathematics & Theoretical Physics, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA (United Kingdom)

2014-04-15

Ambitwistor strings are chiral, infinite tension analogues of conventional string theory whose target space is the space of complex null geodesics and whose spectrum consists exclusively of massless states. At genus zero, these strings underpin the Cachazo-He-Yuan formulæ for tree level scattering of gravitons, gluons and scalars. In this paper we extend these formulæ in a number of directions. Firstly, we consider Ramond sector vertex operators and construct simple amplitudes involving space-time fermions. These agree with tree amplitudes in ten dimensional supergravity and super Yang-Mills. We then show that, after the usual GSO projections, the ambitwistor string partition function is modular invariant. We consider the scattering equations at genus one, and calculate one loop scattering amplitudes for NS-NS external states in the Type II ambitwistor string. We conjecture that these give new representations of (the integrand of) one loop supergravity amplitudes and we show that they have the expected behaviour under factorization of the worldsheet in both non-separating and separating degenerations.

A new treatment of nonlocality in scattering process

Science.gov (United States)

Upadhyay, N. J.; Bhagwat, A.; Jain, B. K.

2018-01-01

Nonlocality in the scattering potential leads to an integro-differential equation. In this equation nonlocality enters through an integral over the nonlocal potential kernel. The resulting Schrödinger equation is usually handled by approximating r,{r}{\\prime }-dependence of the nonlocal kernel. The present work proposes a novel method to solve the integro-differential equation. The method, using the mean value theorem of integral calculus, converts the nonhomogeneous term to a homogeneous term. The effective local potential in this equation turns out to be energy independent, but has relative angular momentum dependence. This method is accurate and valid for any form of nonlocality. As illustrative examples, the total and differential cross sections for neutron scattering off 12C, 56Fe and 100Mo nuclei are calculated with this method in the low energy region (up to 10 MeV) and are found to be in reasonable accord with the experiments.

Collinear limits beyond the leading order from the scattering equations

Energy Technology Data Exchange (ETDEWEB)

Nandan, Dhritiman; Plefka, Jan; Wormsbecher, Wadim [Institut für Physik and IRIS Adlershof, Humboldt-Universität zu Berlin,Zum Großen Windkanal 6, D-12489 Berlin (Germany)

2017-02-08

The structure of tree-level scattering amplitudes for collinear massless bosons is studied beyond their leading splitting function behavior. These near-collinear limits at sub-leading order are best studied using the Cachazo-He-Yuan (CHY) formulation of the S-matrix based on the scattering equations. We compute the collinear limits for gluons, gravitons and scalars. It is shown that the CHY integrand for an n-particle gluon scattering amplitude in the collinear limit at sub-leading order is expressed as a convolution of an (n−1)-particle gluon integrand and a collinear kernel integrand, which is universal. Our representation is shown to obey recently proposed amplitude relations in which the collinear gluons of same helicity are replaced by a single graviton. Finally, we extend our analysis to effective field theories and study the collinear limit of the non-linear sigma model, Einstein-Maxwell-Scalar and Yang-Mills-Scalar theory.

Integrable coupling system of fractional soliton equation hierarchy

Energy Technology Data Exchange (ETDEWEB)

Yu Fajun, E-mail: yfajun@163.co [College of Maths and Systematic Science, Shenyang Normal University, Shenyang 110034 (China)

2009-10-05

In this Letter, we consider the derivatives and integrals of fractional order and present a class of the integrable coupling system of the fractional order soliton equations. The fractional order coupled Boussinesq and KdV equations are the special cases of this class. Furthermore, the fractional AKNS soliton equation hierarchy is obtained.

Quantum-mechanical scattering in one dimension

International Nuclear Information System (INIS)

Boya, Luis J.

2008-01-01

The purpose of this mainly pedagogical review is to fill a lacuna in the usual treatment of scattering in quantum mechanics, by showing the essential of it in the simplest, one-dimensional setting. We define in this situation amplitudes and scattering coefficients and deal with optical and Levinson' theorems as consequences of unitarity in coordinate or momentum space. Parity waves en lieu of partial waves, integral equations and Born series, etc., are defined naturally in this frame. Several solvable examples are shown. Two topics best studied in 1d are transparent potentials and supersymmetric quantum mechanics. Elementary analytical properties and general behaviour of amplitudes give rise to study inverse problems, that is, recovering the potential from scattering data. Isospectral deformations of the wave equation give relations with some nonlinear evolution equations (Lax), solvable by the inverse scattering method (Kruskal), and we consider the KdV equation as an example. We also refer briefly to some singular potentials, where, e.g., the essence of renormalization can be read off again in the simplest setting. The whole paper emphasizes the tutorial and introductory aspects

Well Conditioned Formulations for Open Surface Scattering

National Research Council Canada - National Science Library

Ottusch, John J; Visher, John L

2008-01-01

.... This report describes an analytical preconditioner method for the EFIE on open surface PEC targets that converts the EFIE to a well conditioned, second-kind integral equation. We present theory and the results from a numerical implementation. We also discuss a 2d extension of the Poincare-Bertrand identity could be used to develop an explicitly second-kind integral equation for open surface scattering problems.

Inverse scattering solution of non-linear evolution equations in one space dimension: an introduction

International Nuclear Information System (INIS)

Alvarez-Estrada, R.F.

1979-01-01

A comprehensive review of the inverse scattering solution of certain non-linear evolution equations of physical interest in one space dimension is presented. We explain in some detail the interrelated techniques which allow to linearize exactly the following equations: (1) the Korteweg and de Vries equation; (2) the non-linear Schrodinger equation; (3) the modified Korteweg and de Vries equation; (4) the Sine-Gordon equation. We concentrate in discussing the pairs of linear operators which accomplish such an exact linearization and the solution of the associated initial value problem. The application of the method to other non-linear evolution equations is reviewed very briefly

Inequalities for differential and integral equations

CERN Document Server

Ames, William F

1997-01-01

Inequalities for Differential and Integral Equations has long been needed; it contains material which is hard to find in other books. Written by a major contributor to the field, this comprehensive resource contains many inequalities which have only recently appeared in the literature and which can be used as powerful tools in the development of applications in the theory of new classes of differential and integral equations. For researchers working in this area, it will be a valuable source of reference and inspiration. It could also be used as the text for an advanced graduate course.Key Features* Covers a variety of linear and nonlinear inequalities which find widespread applications in the theory of various classes of differential and integral equations* Contains many inequalities which have only recently appeared in literature and cannot yet be found in other books* Provides a valuable reference to engineers and graduate students

The reduced basis method for the electric field integral equation

International Nuclear Information System (INIS)

Fares, M.; Hesthaven, J.S.; Maday, Y.; Stamm, B.

2011-01-01

We introduce the reduced basis method (RBM) as an efficient tool for parametrized scattering problems in computational electromagnetics for problems where field solutions are computed using a standard Boundary Element Method (BEM) for the parametrized electric field integral equation (EFIE). This combination enables an algorithmic cooperation which results in a two step procedure. The first step consists of a computationally intense assembling of the reduced basis, that needs to be effected only once. In the second step, we compute output functionals of the solution, such as the Radar Cross Section (RCS), independently of the dimension of the discretization space, for many different parameter values in a many-query context at very little cost. Parameters include the wavenumber, the angle of the incident plane wave and its polarization.

Finite-measuring approximation of operators of scattering theory in representation of wave packets

International Nuclear Information System (INIS)

Kukulin, V.I.; Rubtsova, O.A.

2004-01-01

Several types of the packet quantization of the continuos spectrum in the scattering theory quantum problems are considered. Such a quantization leads to the convenient finite-measuring (i.e. matrix) approximation of the integral operators in the scattering theory and it makes it possible to reduce the solution of the singular integral equations, complying with the scattering theory, to the convenient purely algebraic equations on the analytical basis, whereby all the singularities are separated in the obvious form. The main attention is paid to the problems of the method practical realization [ru

Multi-parameter Analysis and Inversion for Anisotropic Media Using the Scattering Integral Method

KAUST Repository

Djebbi, Ramzi

2017-10-24

The main goal in seismic exploration is to identify locations of hydrocarbons reservoirs and give insights on where to drill new wells. Therefore, estimating an Earth model that represents the right physics of the Earth\\'s subsurface is crucial in identifying these targets. Recent seismic data, with long offsets and wide azimuth features, are more sensitive to anisotropy. Accordingly, multiple anisotropic parameters need to be extracted from the recorded data on the surface to properly describe the model. I study the prospect of applying a scattering integral approach for multi-parameter inversion for a transversely isotropic model with a vertical axis of symmetry. I mainly analyze the sensitivity kernels to understand the sensitivity of seismic data to anisotropy parameters. Then, I use a frequency domain scattering integral approach to invert for the optimal parameterization. The scattering integral approach is based on the explicit computation of the sensitivity kernels. I present a new method to compute the traveltime sensitivity kernels for wave equation tomography using the unwrapped phase. I show that the new kernels are a better alternative to conventional cross-correlation/Rytov kernels. I also derive and analyze the sensitivity kernels for a transversely isotropic model with a vertical axis of symmetry. The kernels structure, for various opening/scattering angles, highlights the trade-off regions between the parameters. For a surface recorded data, I show that the normal move-out velocity vn, ƞ and δ parameterization is suitable for a simultaneous inversion of diving waves and reflections. Moreover, when seismic data is inverted hierarchically, the horizontal velocity vh, ƞ and ϵ is the parameterization with the least trade-off. In the frequency domain, the hierarchical inversion approach is naturally implemented using frequency continuation, which makes vh, ƞ and ϵ parameterization attractive. I formulate the multi-parameter inversion using the

A conservative spectral method for the Boltzmann equation with anisotropic scattering and the grazing collisions limit

International Nuclear Information System (INIS)

Gamba, Irene M.; Haack, Jeffrey R.

2014-01-01

We present the formulation of a conservative spectral method for the Boltzmann collision operator with anisotropic scattering cross-sections. The method is an extension of the conservative spectral method of Gamba and Tharkabhushanam [17,18], which uses the weak form of the collision operator to represent the collisional term as a weighted convolution in Fourier space. The method is tested by computing the collision operator with a suitably cut-off angular cross section and comparing the results with the solution of the Landau equation. We analytically study the convergence rate of the Fourier transformed Boltzmann collision operator in the grazing collisions limit to the Fourier transformed Landau collision operator under the assumption of some regularity and decay conditions of the solution to the Boltzmann equation. Our results show that the angular singularity which corresponds to the Rutherford scattering cross section is the critical singularity for which a grazing collision limit exists for the Boltzmann operator. Additionally, we numerically study the differences between homogeneous solutions of the Boltzmann equation with the Rutherford scattering cross section and an artificial cross section, which give convergence to solutions of the Landau equation at different asymptotic rates. We numerically show the rate of the approximation as well as the consequences for the rate of entropy decay for homogeneous solutions of the Boltzmann equation and Landau equation

Reconstruction formula for a 3-d phaseless inverse scattering problem for the Schrodinger equation

OpenAIRE

Klibanov, Michael V.; Romanov, Vladimir G.

2014-01-01

The inverse scattering problem of the reconstruction of the unknown potential with compact support in the 3-d Schr\\"odinger equation is considered. Only the modulus of the scattering complex valued wave field is known, whereas the phase is unknown. It is shown that the unknown potential can be reconstructed via the inverse Radon transform. Therefore, a long standing problem posed in 1977 by K. Chadan and P.C. Sabatier in their book "Inverse Problems in Quantum Scattering Theory" is solved.

Solution of the linear transport equation, monoenergetic in multiregions with anisotopic scattering by the method F sub(N)

International Nuclear Information System (INIS)

Pontedeiro, E.M.B.D.; Maiorino, J.R.

1982-01-01

The linear equation transport, monoenergetic, with anysotropic scattering, in multiregions, by F sub(N) method, is resolved. The mathematical analysis used for this method consists in to use parcially the expansion method in singular autofunctions, or Case's method, aiming to derive a set of integral equations coupled to the angular distribution in the boundaries and interfaces, and then to approximate these distributions by polynomics of N order, aiming to derive, with the use of these boundary and continuity conditions in the interfaces, a set of algebric equations for the coef. of polynomical approximation. With the goal to obtain numerical results, a computer code (FNAM-1) with options for the number of regions, boundary conditions, F sub(N) approx order, were developed. Numerical results were then obtained for various sample problems and compared with the results published in the literature with the objective to demonstrate the precision and applicability of the F sub(N) method. (E.G.) [pt

Bounded solutions for fuzzy differential and integral equations

Energy Technology Data Exchange (ETDEWEB)

Nieto, Juan J. [Departamento de Analisis Matematico Facultad de Matematicas Universidad de Santiago de Compostela, 15782 (Spain)] e-mail: amnieto@usc.es; Rodriguez-Lopez, Rosana [Departamento de Analisis Matematico Facultad de Matematicas Universidad de Santiago de Compostela, 15782 (Spain)] e-mail: amrosana@usc.es

2006-03-01

We find sufficient conditions for the boundness of every solution of first-order fuzzy differential equations as well as certain fuzzy integral equations. Our results are based on several theorems concerning crisp differential and integral inequalities.

Gaussian-windowed frame based method of moments formulation of surface-integral-equation for extended apertures

Energy Technology Data Exchange (ETDEWEB)

Shlivinski, A., E-mail: amirshli@ee.bgu.ac.il [Department of Electrical and Computer Engineering, Ben-Gurion University of the Negev, Beer-Sheva 84105 (Israel); Lomakin, V., E-mail: vlomakin@eng.ucsd.edu [Department of Electrical and Computer Engineering, University of California, San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0407 (United States)

2016-03-01

Scattering or coupling of electromagnetic beam-field at a surface discontinuity separating two homogeneous or inhomogeneous media with different propagation characteristics is formulated using surface integral equation, which are solved by the Method of Moments with the aid of the Gabor-based Gaussian window frame set of basis and testing functions. The application of the Gaussian window frame provides (i) a mathematically exact and robust tool for spatial-spectral phase-space formulation and analysis of the problem; (ii) a system of linear equations in a transmission-line like form relating mode-like wave objects of one medium with mode-like wave objects of the second medium; (iii) furthermore, an appropriate setting of the frame parameters yields mode-like wave objects that blend plane wave properties (as if solving in the spectral domain) with Green's function properties (as if solving in the spatial domain); and (iv) a representation of the scattered field with Gaussian-beam propagators that may be used in many large (in terms of wavelengths) systems.

Algorithms For Integrating Nonlinear Differential Equations

Science.gov (United States)

Freed, A. D.; Walker, K. P.

1994-01-01

Improved algorithms developed for use in numerical integration of systems of nonhomogenous, nonlinear, first-order, ordinary differential equations. In comparison with integration algorithms, these algorithms offer greater stability and accuracy. Several asymptotically correct, thereby enabling retention of stability and accuracy when large increments of independent variable used. Accuracies attainable demonstrated by applying them to systems of nonlinear, first-order, differential equations that arise in study of viscoplastic behavior, spread of acquired immune-deficiency syndrome (AIDS) virus and predator/prey populations.

Integral equations for four identical particles in angular momentum representation

International Nuclear Information System (INIS)

Kharchenko, V.F.; Shadchin, S.A.

1975-01-01

In integral equations of motion for a system of four identical spinless particles with central pair interactions, transition is realized from the representation of relative Jacobi momenta to the representation of their moduli and relative angular moments. As a result, the variables associated with the rotation of the system as a whole are separated in the equations. The integral equations of motion for four particles are reduced to the form of an infinite system of three-demensional integral equations. The four-particle kinematic factors contained in integral kernels are expressed in terms of three-particle type kinematic factors. In the case of separable two-particle interaction, the equations of motion for four particles have the form of an infinite system of two-dimensional integral equations

Efficient Numerical Solution of Coupled Radial Differential Equations in Multichannel Scattering Problems

International Nuclear Information System (INIS)

Houfek, Karel

2008-01-01

Numerical solution of coupled radial differential equations which are encountered in multichannel scattering problems is presented. Numerical approach is based on the combination of the exterior complex scaling method and the finite-elements method with the discrete variable representation. This method can be used not only to solve multichannel scattering problem but also to find bound states and resonance positions and widths directly by diagonalization of the corresponding complex scaled Hamiltonian. Efficiency and accuracy of this method is demonstrated on an analytically solvable two-channel problem.

Feynman integrals and difference equations

International Nuclear Information System (INIS)

Moch, S.; Schneider, C.

2007-09-01

We report on the calculation of multi-loop Feynman integrals for single-scale problems by means of difference equations in Mellin space. The solution to these difference equations in terms of harmonic sums can be constructed algorithmically over difference fields, the so-called ΠΣ * -fields. We test the implementation of the Mathematica package Sigma on examples from recent higher order perturbative calculations in Quantum Chromodynamics. (orig.)

Analyses of pion-40Ca elastic scattering data using the Klein–Gordon equation

International Nuclear Information System (INIS)

Shehadeh, Z.F.

2009-01-01

The elastic scattering data for incident pion energies of 130, 163.3, 180, and 230 MeV on 40 Ca have been analyzed using the full Klein–Gordon equation (KGE), as opposed to its approximate form which renders it to the format of a Schroedinger equation with an energy-dependent potential (RSE). Calculated angular distributions, using KGE and RSE, for all four cases are nearly the same up to about 70° but differ significantly at larger angles. To fit the large-angle data of 163.3 MeV, the nature of the old potential determined by using RSE needs to be revised. The new potentials in four cases are presented and they are compatible with those determined from the inverse scattering theory at a fixed energy in the surface region. (author)

A differential operator for integrating one-loop scattering equations

Energy Technology Data Exchange (ETDEWEB)

Wang, Tianheng [Department of Physics, Nanjing University,Nanjing, Jiangsu Province (China); Chen, Gang [Department of Physics, Zhejiang Normal University,Jinhua, Zhejiang Province (China); Department of Physics and Astronomy, Uppsala University,Uppsala (Sweden); Department of Physics, Nanjing University,Nanjing, Jiangsu Province (China); Cheung, Yeuk-Kwan E. [Department of Physics, Nanjing University,Nanjing, Jiangsu Province (China); Xu, Feng [Weavi Corporation Limited, Nanjing,Jiangsu Province (China)

2017-01-09

We propose a differential operator for computing the residues associated with a class of meromorphic n-forms that frequently appear in the Cachazo-He-Yuan form of the scattering amplitudes. This differential operator is conjectured to be uniquely determined by the local duality theorem and the intersection number of the divisors in the n-form. We use the operator to evaluate the one-loop integrand of Yang-Mills theory from their generalized CHY formulae. The method can reduce the complexity of the calculation. In addition, the expression for the 1-loop four-point Yang-Mills integrand obtained in our approach has a clear correspondence with the Q-cut results.

Inverse scattering transform method and soliton solutions for Davey-Stewartson II equation

International Nuclear Information System (INIS)

Arkadiev, V.A.; Pogrebkov, A.K.; Polivanov, M.C.

1989-01-01

The inverse scattering method for Davey-Stewartson II (DS-II) equation including both soliton and continuous spectrum solutions is developed. The explicit formulae for N-soliton solutions are given. Note that our solitons decrease as |z| -2 with z tending to infinity. (author). 8 refs

Variational Integrals of a Class of Nonhomogeneous -Harmonic Equations

Directory of Open Access Journals (Sweden)

Guanfeng Li

2014-01-01

Full Text Available We introduce a class of variational integrals whose Euler equations are nonhomogeneous -harmonic equations. We investigate the relationship between the minimization problem and the Euler equation and give a simple proof of the existence of some nonhomogeneous -harmonic equations by applying direct methods of the calculus of variations. Besides, we establish some interesting results on variational integrals.

Accurate and efficient quadrature for volterra integral equations

International Nuclear Information System (INIS)

Knirk, D.L.

1976-01-01

Four quadrature schemes were tested and compared in considerable detail to determine their usefulness in the noniterative integral equation method for single-channel quantum-mechanical calculations. They are two forms of linear approximation (trapezoidal rule) and two forms of quadratic approximation (Simpson's rule). Their implementation in this method is shown, a formal discussion of error propagation is given, and tests are performed to determine actual operating characteristics on various bound and scattering problems in different potentials. The quadratic schemes are generally superior to the linear ones in terms of accuracy and efficiency. The previous implementation of Simpson's rule is shown to possess an inherent instability which requires testing on each problem for which it is used to assure its reliability. The alternative quadratic approximation does not suffer this deficiency, but still enjoys the advantages of higher order. In addition, the new scheme obeys very well an h 4 Richardson extrapolation, whereas the old one does so rather poorly. 6 figures, 11 tables

Novel approach to the Helmholtz integral equation solution by Fourier series expansion for acoustic radiation and scattering problems

CSIR Research Space (South Africa)

Shatalov, MY

2006-01-01

Full Text Available -scale structure to guarantee the numerical accuracy of solution. In the present paper the authors propose to use a novel method of solution of the Helmholtz integral equation, which is based on expansion of the integrands in double Fourier series. The main...

Seismic Imaging, One-Way Wave Equations, Pseudodifferential Operators, Path Integrals, and all that Jazz

Science.gov (United States)

Artoun, Ojenie; David-Rus, Diana; Emmett, Matthew; Fishman, Lou; Fital, Sandra; Hogan, Chad; Lim, Jisun; Lushi, Enkeleida; Marinov, Vesselin

2006-05-01

In this report we summarize an extension of Fourier analysis for the solution of the wave equation with a non-constant coefficient corresponding to an inhomogeneous medium. The underlying physics of the problem is exploited to link pseudodifferential operators and phase space path integrals to obtain a marching algorithm that incorporates the backward scattering into the evolution of the wave. This allows us to successfully apply single-sweep, one-way marching methods in inherently two-way environments, which was not achieved before through other methods for this problem.

A hierarchy of Liouville integrable discrete Hamiltonian equations

Energy Technology Data Exchange (ETDEWEB)

Xu Xixiang [College of Science, Shandong University of Science and Technology, Qingdao 266510 (China)], E-mail: xixiang_xu@yahoo.com.cn

2008-05-12

Based on a discrete four-by-four matrix spectral problem, a hierarchy of Lax integrable lattice equations with two potentials is derived. Two Hamiltonian forms are constructed for each lattice equation in the resulting hierarchy by means of the discrete variational identity. A strong symmetry operator of the resulting hierarchy is given. Finally, it is shown that the resulting lattice equations are all Liouville integrable discrete Hamiltonian systems.

Coupled-channel equations and off-shell transformations in many-body scattering

International Nuclear Information System (INIS)

Cattapan, G.; Vanzani, V.

1977-01-01

The general structure and the basic features of several many-body coupled-channel integral equations, obtained by means of the channel coupling array device, are studied in a systematic way. Particular attention is paid to the employment of symmetric transition operators. The connection between different formulations has been clarified and the role played by some off-shell transformations for many-body transition operators has been discussed. Specific choices of the coupling scheme are considered and the corresponding coupled equations are compared with similar equations previously derived. Several sets of linear relations between transition operators have also been presented and used in a three-body context to derive uncoupled integral equations with connected kernel

Complex eigenvalues for neutron transport equation with quadratically anisotropic scattering

International Nuclear Information System (INIS)

Sjoestrand, N.G.

1981-01-01

Complex eigenvalues for the monoenergetic neutron transport equation in the buckling approximation have been calculated for various combinations of linearly and quadratically anisotropic scattering. The results are discussed in terms of the time-dependent case. Tables are given of complex bucklings for real decay constants and of complex decay constants for real bucklings. The results fit nicely into the pattern of real and purely imaginary eigenvalues obtained earlier. (author)

A new path-integral representation of the T-matrix in potential scattering

International Nuclear Information System (INIS)

Carron, J.; Rosenfelder, R.

2011-01-01

We employ the method used by Barbashov and collaborators in Quantum Field Theory to derive a path-integral representation of the T-matrix in nonrelativistic potential scattering which is free of functional integration over fictitious variables as was necessary before. The resulting expression serves as a starting point for a variational approximation applied to high-energy scattering from a Gaussian potential. Good agreement with exact partial-wave calculations is found even at large scattering angles. A novel path-integral representation of the scattering length is obtained in the low-energy limit. -- Highlights: → We derive a new path-integral representation for the T-matrix in quantum scattering from a potential. → The method is based on a technique used by Barbashov and collaborators in Quantum Field Theory. → Unlike previous approaches no unphysical degrees of freedom in the path integral are needed. → The new representation is used for a variational approximation of the T-matrix at high energies. → A new expression for the scattering length at low energy is derived.

Method of mechanical quadratures for solving singular integral equations of various types

Science.gov (United States)

Sahakyan, A. V.; Amirjanyan, H. A.

2018-04-01

The method of mechanical quadratures is proposed as a common approach intended for solving the integral equations defined on finite intervals and containing Cauchy-type singular integrals. This method can be used to solve singular integral equations of the first and second kind, equations with generalized kernel, weakly singular equations, and integro-differential equations. The quadrature rules for several different integrals represented through the same coefficients are presented. This allows one to reduce the integral equations containing integrals of different types to a system of linear algebraic equations.

Feynman integrals and difference equations

Energy Technology Data Exchange (ETDEWEB)

Moch, S. [Deutsches Elektronen-Synchrotron (DESY), Zeuthen (Germany); Schneider, C. [Johannes Kepler Univ., Linz (Austria). Research Inst. for Symbolic Computation

2007-09-15

We report on the calculation of multi-loop Feynman integrals for single-scale problems by means of difference equations in Mellin space. The solution to these difference equations in terms of harmonic sums can be constructed algorithmically over difference fields, the so-called {pi}{sigma}{sup *}-fields. We test the implementation of the Mathematica package Sigma on examples from recent higher order perturbative calculations in Quantum Chromodynamics. (orig.)

Integral equation hierarchy for continuum percolation

International Nuclear Information System (INIS)

Given, J.A.

1988-01-01

In this thesis a projection operator technique is presented that yields hierarchies of integral equations satisfied exactly by the n-point connectedness functions in a continuum version of the site-bond percolation problem. The n-point connectedness functions carry the same structural information for a percolation problem as then-point correlation functions do for a thermal problem. This method extends the Potts model mapping of Fortuin and Kastelyn to the continuum by exploiting an s-state generalization of the Widom-Rowlinson model, a continuum model for phase separation. The projection operator technique is used to produce an integral equation hierarchy for percolation similar to the Born-Green heirarchy. The Kirkwood superposition approximation (SA) is extended to percolation in order to close this hierarchy and yield a nonlinear integral equation for the two-point connectedness function. The fact that this function, in the SA, is the analytic continuation to negative density of the two-point correlation function in a corresponding thermal problem is discussed. The BGY equation for percolation is solved numerically, both by an expansion in powers of the density, and by an iterative technique due to Kirkwood. It is argued both analytically and numerically, that the BYG equation for percolation, unlike its thermal counterpart, shows non-classical critical behavior, with η = 1 and γ = 0.05 ± .1. Finally a sequence of refinements to the superposition approximations based in the theory of fluids by Rice and Lekner is discussed

An analytical theory of a scattering of radio waves on meteoric ionization - II. Solution of the integro-differential equation in case of backscatter

Science.gov (United States)

Pecina, P.

2016-12-01

The integro-differential equation for the polarization vector P inside the meteor trail, representing the analytical solution of the set of Maxwell equations, is solved for the case of backscattering of radio waves on meteoric ionization. The transversal and longitudinal dimensions of a typical meteor trail are small in comparison to the distances to both transmitter and receiver and so the phase factor appearing in the kernel of the integral equation is large and rapidly changing. This allows us to use the method of stationary phase to obtain an approximate solution of the integral equation for the scattered field and for the corresponding generalized radar equation. The final solution is obtained by expanding it into the complete set of Bessel functions, which results in solving a system of linear algebraic equations for the coefficients of the expansion. The time behaviour of the meteor echoes is then obtained using the generalized radar equation. Examples are given for values of the electron density spanning a range from underdense meteor echoes to overdense meteor echoes. We show that the time behaviour of overdense meteor echoes using this method is very different from the one obtained using purely numerical solutions of the Maxwell equations. Our results are in much better agreement with the observations performed e.g. by the Ondřejov radar.

Master equations and the theory of stochastic path integrals

Science.gov (United States)

Weber, Markus F.; Frey, Erwin

2017-04-01

This review provides a pedagogic and self-contained introduction to master equations and to their representation by path integrals. Since the 1930s, master equations have served as a fundamental tool to understand the role of fluctuations in complex biological, chemical, and physical systems. Despite their simple appearance, analyses of master equations most often rely on low-noise approximations such as the Kramers-Moyal or the system size expansion, or require ad-hoc closure schemes for the derivation of low-order moment equations. We focus on numerical and analytical methods going beyond the low-noise limit and provide a unified framework for the study of master equations. After deriving the forward and backward master equations from the Chapman-Kolmogorov equation, we show how the two master equations can be cast into either of four linear partial differential equations (PDEs). Three of these PDEs are discussed in detail. The first PDE governs the time evolution of a generalized probability generating function whose basis depends on the stochastic process under consideration. Spectral methods, WKB approximations, and a variational approach have been proposed for the analysis of the PDE. The second PDE is novel and is obeyed by a distribution that is marginalized over an initial state. It proves useful for the computation of mean extinction times. The third PDE describes the time evolution of a ‘generating functional’, which generalizes the so-called Poisson representation. Subsequently, the solutions of the PDEs are expressed in terms of two path integrals: a ‘forward’ and a ‘backward’ path integral. Combined with inverse transformations, one obtains two distinct path integral representations of the conditional probability distribution solving the master equations. We exemplify both path integrals in analysing elementary chemical reactions. Moreover, we show how a well-known path integral representation of averaged observables can be recovered from

Master equations and the theory of stochastic path integrals.

Science.gov (United States)

Weber, Markus F; Frey, Erwin

2017-04-01

This review provides a pedagogic and self-contained introduction to master equations and to their representation by path integrals. Since the 1930s, master equations have served as a fundamental tool to understand the role of fluctuations in complex biological, chemical, and physical systems. Despite their simple appearance, analyses of master equations most often rely on low-noise approximations such as the Kramers-Moyal or the system size expansion, or require ad-hoc closure schemes for the derivation of low-order moment equations. We focus on numerical and analytical methods going beyond the low-noise limit and provide a unified framework for the study of master equations. After deriving the forward and backward master equations from the Chapman-Kolmogorov equation, we show how the two master equations can be cast into either of four linear partial differential equations (PDEs). Three of these PDEs are discussed in detail. The first PDE governs the time evolution of a generalized probability generating function whose basis depends on the stochastic process under consideration. Spectral methods, WKB approximations, and a variational approach have been proposed for the analysis of the PDE. The second PDE is novel and is obeyed by a distribution that is marginalized over an initial state. It proves useful for the computation of mean extinction times. The third PDE describes the time evolution of a 'generating functional', which generalizes the so-called Poisson representation. Subsequently, the solutions of the PDEs are expressed in terms of two path integrals: a 'forward' and a 'backward' path integral. Combined with inverse transformations, one obtains two distinct path integral representations of the conditional probability distribution solving the master equations. We exemplify both path integrals in analysing elementary chemical reactions. Moreover, we show how a well-known path integral representation of averaged observables can be recovered from them. Upon

Abel integral equations analysis and applications

CERN Document Server

Gorenflo, Rudolf

1991-01-01

In many fields of application of mathematics, progress is crucially dependent on the good flow of information between (i) theoretical mathematicians looking for applications, (ii) mathematicians working in applications in need of theory, and (iii) scientists and engineers applying mathematical models and methods. The intention of this book is to stimulate this flow of information. In the first three chapters (accessible to third year students of mathematics and physics and to mathematically interested engineers) applications of Abel integral equations are surveyed broadly including determination of potentials, stereology, seismic travel times, spectroscopy, optical fibres. In subsequent chapters (requiring some background in functional analysis) mapping properties of Abel integral operators and their relation to other integral transforms in various function spaces are investi- gated, questions of existence and uniqueness of solutions of linear and nonlinear Abel integral equations are treated, and for equatio...

A wavelet-based PWTD algorithm-accelerated time domain surface integral equation solver

KAUST Repository

Liu, Yang

2015-10-26

© 2015 IEEE. The multilevel plane-wave time-domain (PWTD) algorithm allows for fast and accurate analysis of transient scattering from, and radiation by, electrically large and complex structures. When used in tandem with marching-on-in-time (MOT)-based surface integral equation (SIE) solvers, it reduces the computational and memory costs of transient analysis from equation and equation to equation and equation, respectively, where Nt and Ns denote the number of temporal and spatial unknowns (Ergin et al., IEEE Trans. Antennas Mag., 41, 39-52, 1999). In the past, PWTD-accelerated MOT-SIE solvers have been applied to transient problems involving half million spatial unknowns (Shanker et al., IEEE Trans. Antennas Propag., 51, 628-641, 2003). Recently, a scalable parallel PWTD-accelerated MOT-SIE solver that leverages a hiearchical parallelization strategy has been developed and successfully applied to the transient problems involving ten million spatial unknowns (Liu et. al., in URSI Digest, 2013). We further enhanced the capabilities of this solver by implementing a compression scheme based on local cosine wavelet bases (LCBs) that exploits the sparsity in the temporal dimension (Liu et. al., in URSI Digest, 2014). Specifically, the LCB compression scheme was used to reduce the memory requirement of the PWTD ray data and computational cost of operations in the PWTD translation stage.

The Cauchy problem for the Pavlov equation

International Nuclear Information System (INIS)

Grinevich, P G; Santini, P M; Wu, D

2015-01-01

Commutation of multidimensional vector fields leads to integrable nonlinear dispersionless PDEs that arise in various problems of mathematical physics and have been intensively studied in recent literature. This report aims to solve the scattering and inverse scattering problem for integrable dispersionless PDEs, recently introduced just at a formal level, concentrating on the prototypical example of the Pavlov equation, and to justify an existence theorem for global bounded solutions of the associated Cauchy problem with small data. (paper)

Fuchs indices and the first integrals of nonlinear differential equations

International Nuclear Information System (INIS)

Kudryashov, Nikolai A.

2005-01-01

New method of finding the first integrals of nonlinear differential equations in polynomial form is presented. Basic idea of our approach is to use the scaling of solution of nonlinear differential equation and to find the dimensions of arbitrary constants in the Laurent expansion of the general solution. These dimensions allows us to obtain the scalings of members for the first integrals of nonlinear differential equations. Taking the polynomials with unknown coefficients into account we present the algorithm of finding the first integrals of nonlinear differential equations in the polynomial form. Our method is applied to look for the first integrals of eight nonlinear ordinary differential equations of the fourth order. The general solution of one of the fourth order ordinary differential equations is given

Molar mass, radius of gyration and second virial coefficient from new static light scattering equations for dilute solutions: application to 21 (macro)molecules.

Science.gov (United States)

Illien, Bertrand; Ying, Ruifeng

2009-05-11

New static light scattering (SLS) equations for dilute binary solutions are derived. Contrarily to the usual SLS equations [Carr-Zimm (CZ)], the new equations have no need for the experimental absolute Rayleigh ratio of a reference liquid and solely rely on the ratio of scattered intensities of solutions and solvent. The new equations, which are based on polarizability equations, take into account the usual refractive index increment partial differential n/partial differential rho(2) complemented by the solvent specific polarizability and a term proportional to the slope of the solution density rho versus the solute mass concentration rho(2) (density increment). Then all the equations are applied to 21 (macro)molecules with a wide range of molar mass (0.2equations clearly achieve a better agreement with supplier M values. For macromolecules (M>500 kg mol(-1)), for which the scattered intensity is no longer independent of the scattering angle, the new equations give the same value of the radius of gyration as the CZ equation and consistent values of the second virial coefficient.

Absorption line profiles in a moving atmosphere - A single scattering linear perturbation theory

Science.gov (United States)

Hays, P. B.; Abreu, V. J.

1989-01-01

An integral equation is derived which linearly relates Doppler perturbations in the spectrum of atmospheric absorption features to the wind system which creates them. The perturbation theory is developed using a single scattering model, which is validated against a multiple scattering calculation. The nature and basic properties of the kernels in the integral equation are examined. It is concluded that the kernels are well behaved and that wind velocity profiles can be recovered using standard inversion techniques.

On a Volterra Stieltjes integral equation

Directory of Open Access Journals (Sweden)

P. T. Vaz

1990-01-01

Full Text Available The paper deals with a study of linear Volterra integral equations involving Lebesgue-Stieltjes integrals in two independent variables. The authors prove an existence theorem using the Banach fixed-point principle. An explicit example is also considered.

Transmission problem for the Laplace equation and the integral equation method

Czech Academy of Sciences Publication Activity Database

Medková, Dagmar

2012-01-01

Roč. 387, č. 2 (2012), s. 837-843 ISSN 0022-247X Institutional research plan: CEZ:AV0Z10190503 Keywords : transmission problem * Laplace equation * boundary integral equation Subject RIV: BA - General Mathematics Impact factor: 1.050, year: 2012 http://www.sciencedirect.com/science/article/pii/S0022247X11008985

Numerov iteration method for second order integral-differential equation

International Nuclear Information System (INIS)

Zeng Fanan; Zhang Jiaju; Zhao Xuan

1987-01-01

In this paper, Numerov iterative method for second order integral-differential equation and system of equations are constructed. Numerical examples show that this method is better than direct method (Gauss elimination method) in CPU time and memoy requireing. Therefore, this method is an efficient method for solving integral-differential equation in nuclear physics

An integrable semi-discretization of the Boussinesq equation

International Nuclear Information System (INIS)

Zhang, Yingnan; Tian, Lixin

2016-01-01

Highlights: • A new integrable semi-discretization of the Boussinesq equation is present. • A Bäcklund transformation and a Lax pair for the differential-difference system is derived by using Hirota's bilinear method. • The soliton solutions of 'good' Boussinesq equation and numerical algorithms are investigated. - Abstract: In this paper, we present an integrable semi-discretization of the Boussinesq equation. Different from other discrete analogues, we discretize the ‘time’ variable and get an integrable differential-difference system. Under a standard limitation, the differential-difference system converges to the continuous Boussinesq equation such that the discrete system can be used to design numerical algorithms. Using Hirota's bilinear method, we find a Bäcklund transformation and a Lax pair of the differential-difference system. For the case of ‘good’ Boussinesq equation, we investigate the soliton solutions of its discrete analogue and design numerical algorithms. We find an effective way to reduce the phase shift caused by the discretization. The numerical results coincide with our analysis.

On the solution of the inverse scattering problem on a ray

International Nuclear Information System (INIS)

Egikyan, R.S.; Zhidkov, E.P.

1988-01-01

Quantum inverse scattering problem (ISP) is considered within the framework of two-particle scattering for local interaction case depending only on the scattering between particles. Constructing the solution of secondary integral equation solution of ISP is described in the clear image. Numerical calculations are conducted using a direct method

ON ASYMTOTIC APPROXIMATIONS OF FIRST INTEGRALS FOR DIFFERENTIAL AND DIFFERENCE EQUATIONS

Directory of Open Access Journals (Sweden)

W.T. van Horssen

2007-04-01

Full Text Available In this paper the concept of integrating factors for differential equations and the concept of invariance factors for difference equations to obtain first integrals or invariants will be presented. It will be shown that all integrating factors have to satisfya system of partial differential equations, and that all invariance factors have to satisfy a functional equation. In the period 1997-2001 a perturbation method based on integrating vectors was developed to approximate first integrals for systems of ordinary differential equations. This perturbation method will be reviewed shortly. Also in the paper the first results in the development of a perturbation method for difference equations based on invariance factors will be presented.

Multiple scattering of electromagnetic waves by a collection of plasma drift turbulent vortices

International Nuclear Information System (INIS)

Resendes, D.

1995-01-01

An application of the self-consistent multiple-scattering theory of electro-magnetic waves to drift turbulent vortices is presented. Using the known single-vortex solution, the integral equation describing the scattering from a finite density of drift turbulent vortices is obtained. Rather than solving this equation and then averaging, the averaging operation is taken first to obtain statistical moment equations, from which the coherent and incoherent scattering follow. These results are expressed in a Fourier basis, and the cross-section is evaluated. Limiting forms of the theory and straightforward generalizations are discussed. (Author)

Scattering of lower-hybrid waves by drift-wave density fluctuations: solutions of the radiative transfer equation

International Nuclear Information System (INIS)

Andrews, P.L.; Perkins, F.W.

1983-01-01

The investigation of the scattering of lower-hybrid waves by density fluctuations arising from drift waves in tokamaks is distinguished by the presence in the wave equation of a large, random, derivative-coupling term. The propagation of the lower-hybrid waves is well represented by a radiative transfer equation when the scale size of the density fluctuations is small compared to the overall plasma size. The radiative transfer equation is solved in two limits: first, the forward scattering limit, where the scale size of density fluctuations is large compared to the lower-hybrid perpendicular wavelength, and second, the large-angle scattering limit, where this inequality is reversed. The most important features of these solutions are well represented by analytical formulas derived by simple arguments. Based on conventional estimates for density fluctuations arising from drift waves and a parabolic density profile, the optical depth tau for scattering through a significant angle, is given by tauroughly-equal(2/N 2 /sub parallel/) (#betta#/sub p/i0/#betta#) 2 (m/sub e/c 2 /2T/sub i/)/sup 1/2/ [c/α(Ω/sub i/Ω/sub e/)/sup 1/2/ ], where #betta#/sub p/i0 is the central ion plasma frequency and T/sub i/ denotes the ion temperature near the edge of the plasma. Most of the scattering occurs near the surface. The transmission through the scattering region scales as tau - 1 and the emerging intensity has an angular spectrum proportional to cos theta, where sin theta = k/sub perpendicular/xB/sub p//(k/sub perpendicular/B/sub p/), and B/sub p/ is the poloidal field

An integral equation arising in two group neutron transport theory

International Nuclear Information System (INIS)

Cassell, J S; Williams, M M R

2003-01-01

An integral equation describing the fuel distribution necessary to maintain a flat flux in a nuclear reactor in two group transport theory is reduced to the solution of a singular integral equation. The formalism developed enables the physical aspects of the problem to be better understood and its relationship with the corresponding diffusion theory model is highlighted. The integral equation is solved by reducing it to a non-singular Fredholm equation which is then evaluated numerically

Nonlocal integrable PDEs from hierarchies of symmetry laws: The example of Pohlmeyer-Lund-Regge equation and its reflectionless potential solutions

Science.gov (United States)

Demontis, F.; Ortenzi, G.; van der Mee, C.

2018-04-01

By following the ideas presented by Fukumoto and Miyajima in Fukumoto and Miyajima (1996) we derive a generalized method for constructing integrable nonlocal equations starting from any bi-Hamiltonian hierarchy supplied with a recursion operator. This construction provides the right framework for the application of the full machinery of the inverse scattering transform. We pay attention to the Pohlmeyer-Lund-Regge equation coming from the nonlinear Schrödinger hierarchy and construct the formula for the reflectionless potential solutions which are generalizations of multi-solitons. Some explicit examples are discussed.

Calculation of the Full Scattering Amplitude without Partial Wave Decomposition II: Inclusion of Exchange

Science.gov (United States)

Shertzer, Janine; Temkin, A.

2003-01-01

As is well known, the full scattering amplitude can be expressed as an integral involving the complete scattering wave function. We have shown that the integral can be simplified and used in a practical way. Initial application to electron-hydrogen scattering without exchange was highly successful. The Schrodinger equation (SE), which can be reduced to a 2d partial differential equation (pde), was solved using the finite element method. We have now included exchange by solving the resultant SE, in the static exchange approximation, which is reducible to a pair of coupled pde's. The resultant scattering amplitudes, both singlet and triplet, calculated as a function of energy are in excellent agreement with converged partial wave results.

ON DIFFERENTIAL EQUATIONS, INTEGRABLE SYSTEMS, AND GEOMETRY

OpenAIRE

Enrique Gonzalo Reyes Garcia

2004-01-01

ON DIFFERENTIAL EQUATIONS, INTEGRABLE SYSTEMS, AND GEOMETRY Equations in partial derivatives appeared in the 18th century as essential tools for the analytic study of physical models and, later, they proved to be fundamental for the progress of mathematics. For example, fundamental results of modern differential geometry are based on deep theorems on differential equations. Reciprocally, it is possible to study differential equations through geometrical means just like it was done by o...

Harmonic R-matrices for scattering amplitudes and spectral regularization

Energy Technology Data Exchange (ETDEWEB)

Ferro, Livia; Plefka, Jan [Humboldt-Universitaet, Berlin (Germany). Inst. fuer Physik; Lukowski, Tomasz [Humboldt-Universitaet, Berlin (Germany). Inst. fuer Physik; Humboldt-Universitaet, Berlin (Germany). Inst. fuer Mathematik; Humboldt-Univ. Berlin (Germany). IRIS Adlershof; Meneghelli, Carlo [Hamburg Univ. (Germany). Fachbereich 11 - Mathematik; Deutsches Elektronen-Synchrotron (DESY), Hamburg (Germany). Theory Group; Staudacher, Matthias [Humboldt-Universitaet, Berlin (Germany). Inst. fuer Mathematik; Humboldt-Universitaet, Berlin (Germany). Inst. fuer Physik; Max-Planck-Institut fuer Gravitationsphysik (Albert-Einstein-Institut), Potsdam (Germany)

2012-12-15

Planar N=4 super Yang-Mills appears to be integrable. While this allows to find this theory's exact spectrum, integrability has hitherto been of no direct use for scattering amplitudes. To remedy this, we deform all scattering amplitudes by a spectral parameter. The deformed tree-level four-point function turns out to be essentially the one-loop R-matrix of the integrable N=4 spin chain satisfying the Yang-Baxter equation. Deformed on-shell three-point functions yield novel three-leg R-matrices satisfying bootstrap equations. Finally, we supply initial evidence that the spectral parameter might find its use as a novel symmetry-respecting regulator replacing dimensional regularization. Its physical meaning is a local deformation of particle helicity, a fact which might be useful for a much larger class of non-integrable four-dimensional field theories.

Fitting Data to Model: Structural Equation Modeling Diagnosis Using Two Scatter Plots

Science.gov (United States)

Yuan, Ke-Hai; Hayashi, Kentaro

2010-01-01

This article introduces two simple scatter plots for model diagnosis in structural equation modeling. One plot contrasts a residual-based M-distance of the structural model with the M-distance for the factor score. It contains information on outliers, good leverage observations, bad leverage observations, and normal cases. The other plot contrasts…

The unified approach to integrable relativistic equations: Soliton solutions over non-vanishing backgrounds - 1

International Nuclear Information System (INIS)

Barashenkov, I.V.; Getmanov, B.S.; Kovtun, V.E.

1992-01-01

The scheme for unified description of integrable relativistic massive systems provides an inverse scattering formalism that covers universally all (1+1)- dimensional systems of this kind. In this work we construct the N-soliton solution (over an arbitrary background) for some generic system which is associated with the sl(2,C) case of the scheme and whose reductions include the complex sine-Gordon equation, the massive Thirring model and other equations, both in the Euclidean and Minkowski spaces. Thus the N-soliton solutions for all these systems emerge in a unified form differing only in the type of constraints imposed on their parameters. In an earlier paper the case of the zero background was considered while here we concentrate on the case of the non-vanishing constant background i.e., on the N-kink solutions. (author). 18 refs

Separable expansion for realistic multichannel scattering problems

International Nuclear Information System (INIS)

Canton, L.; Cattapan, G.; Pisent, G.

1987-01-01

A new approach to the multichannel scattering problem with realistic local or nonlocal interactions is developed. By employing the negative-energy solutions of uncoupled Sturmian eigenvalue problems referring to simple auxiliary potentials, the coupling interactions appearing to the original multichannel problem are approximated by finite-rank potentials. By resorting to integral-equation tecniques the coupled-channel equations are then reduced to linear algebraic equations which can be straightforwardly solved. Compact algebraic expressions for the relevant scattering matrix elements are thus obtained. The convergence of the method is tasted in the single-channel case with realistic optical potentials. Excellent agreement is obtained with a few terms in the separable expansion for both real and absorptive interactions

Hybrid radiosity-SP3 equation based bioluminescence tomography reconstruction for turbid medium with low- and non-scattering regions

Science.gov (United States)

Chen, Xueli; Zhang, Qitan; Yang, Defu; Liang, Jimin

2014-01-01

To provide an ideal solution for a specific problem of gastric cancer detection in which low-scattering regions simultaneously existed with both the non- and high-scattering regions, a novel hybrid radiosity-SP3 equation based reconstruction algorithm for bioluminescence tomography was proposed in this paper. In the algorithm, the third-order simplified spherical harmonics approximation (SP3) was combined with the radiosity equation to describe the bioluminescent light propagation in tissues, which provided acceptable accuracy for the turbid medium with both low- and non-scattering regions. The performance of the algorithm was evaluated with digital mouse based simulations and a gastric cancer-bearing mouse based in situ experiment. Primary results demonstrated the feasibility and superiority of the proposed algorithm for the turbid medium with low- and non-scattering regions.

Hybrid radiosity-SP3 equation based bioluminescence tomography reconstruction for turbid medium with low- and non-scattering regions

International Nuclear Information System (INIS)

Chen, Xueli; Zhang, Qitan; Yang, Defu; Liang, Jimin

2014-01-01

To provide an ideal solution for a specific problem of gastric cancer detection in which low-scattering regions simultaneously existed with both the non- and high-scattering regions, a novel hybrid radiosity-SP 3 equation based reconstruction algorithm for bioluminescence tomography was proposed in this paper. In the algorithm, the third-order simplified spherical harmonics approximation (SP 3 ) was combined with the radiosity equation to describe the bioluminescent light propagation in tissues, which provided acceptable accuracy for the turbid medium with both low- and non-scattering regions. The performance of the algorithm was evaluated with digital mouse based simulations and a gastric cancer-bearing mouse based in situ experiment. Primary results demonstrated the feasibility and superiority of the proposed algorithm for the turbid medium with low- and non-scattering regions

Discrete Ordinates Approximations to the First- and Second-Order Radiation Transport Equations

International Nuclear Information System (INIS)

FAN, WESLEY C.; DRUMM, CLIFTON R.; POWELL, JENNIFER L. email wcfan@sandia.gov

2002-01-01

The conventional discrete ordinates approximation to the Boltzmann transport equation can be described in a matrix form. Specifically, the within-group scattering integral can be represented by three components: a moment-to-discrete matrix, a scattering cross-section matrix and a discrete-to-moment matrix. Using and extending these entities, we derive and summarize the matrix representations of the second-order transport equations

Discrete Ordinates Approximations to the First- and Second-Order Radiation Transport Equations

CERN Document Server

Fan, W C; Powell, J L

2002-01-01

The conventional discrete ordinates approximation to the Boltzmann transport equation can be described in a matrix form. Specifically, the within-group scattering integral can be represented by three components: a moment-to-discrete matrix, a scattering cross-section matrix and a discrete-to-moment matrix. Using and extending these entities, we derive and summarize the matrix representations of the second-order transport equations.

A Generalized Analytic Operator-Valued Function Space Integral and a Related Integral Equation

International Nuclear Information System (INIS)

Chang, K.S.; Kim, B.S.; Park, C.H.; Ryu, K.S.

2003-01-01

We introduce a generalized Wiener measure associated with a Gaussian Markov process and define a generalized analytic operator-valued function space integral as a bounded linear operator from L p into L p-ci r cumflexprime (1< p ≤ 2) by the analytic continuation of the generalized Wiener integral. We prove the existence of the integral for certain functionals which involve some Borel measures. Also we show that the generalized analytic operator-valued function space integral satisfies an integral equation related to the generalized Schroedinger equation. The resulting theorems extend the theory of operator-valued function space integrals substantially and previous theorems about these integrals are generalized by our results

An integrable semi-discretization of the Boussinesq equation

Energy Technology Data Exchange (ETDEWEB)

Zhang, Yingnan, E-mail: ynzhang@njnu.edu.cn [Jiangsu Key Laboratory for NSLSCS, School of Mathematical Sciences, Nanjing Normal University, Nanjing, Jiangsu (China); Tian, Lixin, E-mail: tianlixin@njnu.edu.cn [Jiangsu Key Laboratory for NSLSCS, School of Mathematical Sciences, Nanjing Normal University, Nanjing, Jiangsu (China); Nonlinear Scientific Research Center, Jiangsu University, Zhenjiang, Jiangsu (China)

2016-10-23

Highlights: • A new integrable semi-discretization of the Boussinesq equation is present. • A Bäcklund transformation and a Lax pair for the differential-difference system is derived by using Hirota's bilinear method. • The soliton solutions of 'good' Boussinesq equation and numerical algorithms are investigated. - Abstract: In this paper, we present an integrable semi-discretization of the Boussinesq equation. Different from other discrete analogues, we discretize the ‘time’ variable and get an integrable differential-difference system. Under a standard limitation, the differential-difference system converges to the continuous Boussinesq equation such that the discrete system can be used to design numerical algorithms. Using Hirota's bilinear method, we find a Bäcklund transformation and a Lax pair of the differential-difference system. For the case of ‘good’ Boussinesq equation, we investigate the soliton solutions of its discrete analogue and design numerical algorithms. We find an effective way to reduce the phase shift caused by the discretization. The numerical results coincide with our analysis.

Normal and adjoint integral and integrodifferential neutron transport equations. Pt. 2

International Nuclear Information System (INIS)

Velarde, G.

1976-01-01

Using the simplifying hypotheses of the integrodifferential Boltzmann equations of neutron transport, given in JEN 334 report, several integral equations, and theirs adjoint ones, are obtained. Relations between the different normal and adjoint eigenfunctions are established and, in particular, proceeding from the integrodifferential Boltzmann equation it's found out the relation between the solutions of the adjoint equation of its integral one, and the solutions of the integral equation of its adjoint one (author)

Scattering states of the Klein-Gordon equation with Coulomb-like scalar plus vector potentials in arbitrary dimension

International Nuclear Information System (INIS)

Chen Changyuan; Sun Dongsheng; Lu Falin

2004-01-01

Properties of scattering states of the Klein-Gordon equation with Coulomb-like scalar plus vector potentials are investigated in an arbitrary dimension. Exact results of normalized wave functions of scattering states in the 'k/2π scale' and formula of phase shifts are presented

On the evolution equations, solvable through the inverse scattering method

International Nuclear Information System (INIS)

Gerdjikov, V.S.; Khristov, E.Kh.

1979-01-01

The nonlinear evolution equations (NLEE), related to the one-parameter family of Dirac operators are considered in a uniform manner. The class of NLEE solvable through the inverse scatterina method and their conservation laws are described. The description of the hierarchy of Hamiltonian structures and the proof of complete integrability of the NLEE is presented. The class of Baecklund transformations for these NLEE is derived. The general formulae are illustrated by two important examples: the nonlinear Schroedinger equation and the sine-Gordon equation

On solvability of some quadratic functional-integral equation in Banach algebra

International Nuclear Information System (INIS)

Darwish, M.A.

2007-08-01

Using the technique of a suitable measure of non-compactness in Banach algebra, we prove an existence theorem for some functional-integral equations which contain, as particular cases, a lot of integral and functional-integral equations that arise in many branches of nonlinear analysis and its applications. Also, the famous Chandrasekhar's integral equation is considered as a special case. (author)

Determinantal method for complex angular momenta in potential scattering

Energy Technology Data Exchange (ETDEWEB)

Lee, B. W. [University of Pennsylvania, Philadelphia, PA (United States)

1963-01-15

In this paper I would like do describe a formulation of the complex angular momenta in potential scattering based on the Lippmann-Schwinger integral equation rather than on the Schrödinger differential equation. This is intended as a preliminary to the paper by SAWYER on the Regge poles and high energy limits in field theory (Bethe-Salpeter amplitudes), where the integral formulation is definitely more advantageous than the differential formulation.

Integrable boundary conditions and modified Lax equations

International Nuclear Information System (INIS)

Avan, Jean; Doikou, Anastasia

2008-01-01

We consider integrable boundary conditions for both discrete and continuum classical integrable models. Local integrals of motion generated by the corresponding 'transfer' matrices give rise to time evolution equations for the initial Lax operator. We systematically identify the modified Lax pairs for both discrete and continuum boundary integrable models, depending on the classical r-matrix and the boundary matrix

Integral equation for inhomogeneous condensed bosons generalizing the Gross-Pitaevskii differential equation

International Nuclear Information System (INIS)

Angilella, G.G.N.; Pucci, R.; March, N.H.

2004-01-01

We give here the derivation of a Gross-Pitaevskii-type equation for inhomogeneous condensed bosons. Instead of the original Gross-Pitaevskii differential equation, we obtain an integral equation that implies less restrictive assumptions than are made in the very recent study of Pieri and Strinati [Phys. Rev. Lett. 91, 030401 (2003)]. In particular, the Thomas-Fermi approximation and the restriction to small spatial variations of the order parameter invoked in their study are avoided

THE SIMULATION OF SCATTERING OF ELECTROMAGNETIC WAVES ON ANGULAR STRUCTURES.

Directory of Open Access Journals (Sweden)

P. A. Preobrazhensky

2017-02-01

Full Text Available The paper discusses the characteristics of scattering of electromagnetic waves on the angular diffraction structures. The solution of the problem is based on the method of integral equations. A comparative analysis of the scattering characteristics of structures with different shape is carried out.

Numerical solution of boundary-integral equations for molecular electrostatics.

Science.gov (United States)

Bardhan, Jaydeep P

2009-03-07

Numerous molecular processes, such as ion permeation through channel proteins, are governed by relatively small changes in energetics. As a result, theoretical investigations of these processes require accurate numerical methods. In the present paper, we evaluate the accuracy of two approaches to simulating boundary-integral equations for continuum models of the electrostatics of solvation. The analysis emphasizes boundary-element method simulations of the integral-equation formulation known as the apparent-surface-charge (ASC) method or polarizable-continuum model (PCM). In many numerical implementations of the ASC/PCM model, one forces the integral equation to be satisfied exactly at a set of discrete points on the boundary. We demonstrate in this paper that this approach to discretization, known as point collocation, is significantly less accurate than an alternative approach known as qualocation. Furthermore, the qualocation method offers this improvement in accuracy without increasing simulation time. Numerical examples demonstrate that electrostatic part of the solvation free energy, when calculated using the collocation and qualocation methods, can differ significantly; for a polypeptide, the answers can differ by as much as 10 kcal/mol (approximately 4% of the total electrostatic contribution to solvation). The applicability of the qualocation discretization to other integral-equation formulations is also discussed, and two equivalences between integral-equation methods are derived.

Irreducibility and co-primeness as an integrability criterion for discrete equations

International Nuclear Information System (INIS)

Kanki, Masataka; Mada, Jun; Mase, Takafumi; Tokihiro, Tetsuji

2014-01-01

We study the Laurent property, the irreducibility and co-primeness of discrete integrable and non-integrable equations. First we study a discrete integrable equation related to the Somos-4 sequence, and also a non-integrable equation as a comparison. We prove that the conditions of irreducibility and co-primeness hold only in the integrable case. Next, we generalize our previous results on the singularities of the discrete Korteweg–de Vries (dKdV) equation. In our previous paper (Kanki et al 2014 J. Phys. A: Math. Theor. 47 065201) we described the singularity confinement test (one of the integrability criteria) using the Laurent property, and the irreducibility, and co-primeness of the terms in the bilinear dKdV equation, in which we only considered simplified boundary conditions. This restriction was needed to obtain simple (monomial) relations between the bilinear form and the nonlinear form of the dKdV equation. In this paper, we prove the co-primeness of the terms in the nonlinear dKdV equation for general initial conditions and boundary conditions, by using the localization of Laurent rings and the interchange of the axes. We assert that co-primeness of the terms can be used as a new integrability criterion, which is a mathematical re-interpretation of the confinement of singularities in the case of discrete equations. (paper)

Gated integrator PXI-DAQ system for Thomson scattering diagnostics

Energy Technology Data Exchange (ETDEWEB)

Patel, Kiran, E-mail: kkpatel@ipr.res.in; Pillai, Vishal; Singh, Neha; Thomas, Jinto; Kumar, Ajai

2017-06-15

Gated Integrator (GI) PXI based data acquisition (DAQ) system has been designed and developed for the ease of acquiring fast Thomson Scattered signals (∼50 ns pulse width). The DAQ system consists of in-house designed and developed GI modules and PXI-1405 chassis with several PXI-DAQ modules. The performance of the developed system has been validated during the SST-1 campaigns. The dynamic range of the GI module depends on the integrating capacitor (C{sub i}) and the modules have been calibrated using 12 pF and 27 pF integrating capacitors. The developed GI module based data acquisition system consists of sixty four channels for simultaneous sampling using eight PXI based digitization modules having eight channels per module. The error estimation and functional tests of this unit are carried out using standard source and also with the fast detectors used for Thomson scattering diagnostics. User friendly Graphical User Interface (GUI) has been developed using LabVIEW on Windows platform to control and acquire the Thomson scattering signal. A robust, easy to operate and maintain with low power consumption, having higher dynamic range with very good sensitivity and cost effective DAQ system is developed and tested for the SST-1 Thomson scattering diagnostics.

The Cauchy problem for the Pavlov equation with large data

Science.gov (United States)

Wu, Derchyi

2017-08-01

We prove a local solvability of the Cauchy problem for the Pavlov equation with large initial data by the inverse scattering method. The Pavlov equation arises in studies Einstein-Weyl geometries and dispersionless integrable models. Our theory yields a local solvability of Cauchy problems for a quasi-linear wave equation with a characteristic initial hypersurface.

Set-Valued Stochastic Equation with Set-Valued Square Integrable Martingale

Directory of Open Access Journals (Sweden)

Li Jun-Gang

2017-01-01

Full Text Available In this paper, we shall introduce the stochastic integral of a stochastic process with respect to set-valued square integrable martingale. Then we shall give the Aumann integral measurable theorem, and give the set-valued stochastic Lebesgue integral and set-valued square integrable martingale integral equation. The existence and uniqueness of solution to set-valued stochastic integral equation are proved. The discussion will be useful in optimal control and mathematical finance in psychological factors.

Two analytic transport equation solutions for particular cases of particle history

International Nuclear Information System (INIS)

Simovic, R.

1997-01-01

For anisotropic scattering and plane geometry, the linear transport equation of particles generated by a monodirectional unit source A(x,μ) = δ(x-0)δ(μ - μ 0 ) > 0, can be stated in the form of an integral equation

Babenko’s Approach to Abel’s Integral Equations

Directory of Open Access Journals (Sweden)

Chenkuan Li

2018-03-01

Full Text Available The goal of this paper is to investigate the following Abel’s integral equation of the second kind: y ( t + λ Γ ( α ∫ 0 t ( t − τ α − 1 y ( τ d τ = f ( t , ( t > 0 and its variants by fractional calculus. Applying Babenko’s approach and fractional integrals, we provide a general method for solving Abel’s integral equation and others with a demonstration of different types of examples by showing convergence of series. In particular, we extend this equation to a distributional space for any arbitrary α ∈ R by fractional operations of generalized functions for the first time and obtain several new and interesting results that cannot be realized in the classical sense or by the Laplace transform.

Symbolic-Numeric Integration of the Dynamical Cosserat Equations

KAUST Repository

Lyakhov, Dmitry A.

2017-08-29

We devise a symbolic-numeric approach to the integration of the dynamical part of the Cosserat equations, a system of nonlinear partial differential equations describing the mechanical behavior of slender structures, like fibers and rods. This is based on our previous results on the construction of a closed form general solution to the kinematic part of the Cosserat system. Our approach combines methods of numerical exponential integration and symbolic integration of the intermediate system of nonlinear ordinary differential equations describing the dynamics of one of the arbitrary vector-functions in the general solution of the kinematic part in terms of the module of the twist vector-function. We present an experimental comparison with the well-established generalized \\\\alpha -method illustrating the computational efficiency of our approach for problems in structural mechanics.

Symbolic-Numeric Integration of the Dynamical Cosserat Equations

KAUST Repository

Lyakhov, Dmitry A.; Gerdt, Vladimir P.; Weber, Andreas G.; Michels, Dominik L.

2017-01-01

We devise a symbolic-numeric approach to the integration of the dynamical part of the Cosserat equations, a system of nonlinear partial differential equations describing the mechanical behavior of slender structures, like fibers and rods. This is based on our previous results on the construction of a closed form general solution to the kinematic part of the Cosserat system. Our approach combines methods of numerical exponential integration and symbolic integration of the intermediate system of nonlinear ordinary differential equations describing the dynamics of one of the arbitrary vector-functions in the general solution of the kinematic part in terms of the module of the twist vector-function. We present an experimental comparison with the well-established generalized \\alpha -method illustrating the computational efficiency of our approach for problems in structural mechanics.

Hybrid radiosity-SP{sub 3} equation based bioluminescence tomography reconstruction for turbid medium with low- and non-scattering regions

Energy Technology Data Exchange (ETDEWEB)

Chen, Xueli, E-mail: xlchen@xidian.edu.cn, E-mail: jimleung@mail.xidian.edu.cn; Zhang, Qitan; Yang, Defu; Liang, Jimin, E-mail: xlchen@xidian.edu.cn, E-mail: jimleung@mail.xidian.edu.cn [School of Life Science and Technology, Xidian University, Xi' an, Shaanxi 710071 (China)

2014-01-14

To provide an ideal solution for a specific problem of gastric cancer detection in which low-scattering regions simultaneously existed with both the non- and high-scattering regions, a novel hybrid radiosity-SP{sub 3} equation based reconstruction algorithm for bioluminescence tomography was proposed in this paper. In the algorithm, the third-order simplified spherical harmonics approximation (SP{sub 3}) was combined with the radiosity equation to describe the bioluminescent light propagation in tissues, which provided acceptable accuracy for the turbid medium with both low- and non-scattering regions. The performance of the algorithm was evaluated with digital mouse based simulations and a gastric cancer-bearing mouse based in situ experiment. Primary results demonstrated the feasibility and superiority of the proposed algorithm for the turbid medium with low- and non-scattering regions.

Integrability of the Gross-Pitaevskii equation with Feshbach resonance management

International Nuclear Information System (INIS)

Zhao Dun; Luo Honggang; Chai Huayue

2008-01-01

In this Letter we study the integrability of a class of Gross-Pitaevskii equations managed by Feshbach resonance in an expulsive parabolic external potential. By using WTC test, we find a condition under which the Gross-Pitaevskii equation is completely integrable. Under the present model, this integrability condition is completely consistent with that proposed by Serkin, Hasegawa, and Belyaeva [V.N. Serkin, A. Hasegawa, T.L. Belyaeva, Phys. Rev. Lett. 98 (2007) 074102]. Furthermore, this integrability can also be explicitly shown by a transformation, which can convert the Gross-Pitaevskii equation into the well-known standard nonlinear Schroedinger equation. By this transformation, each exact solution of the standard nonlinear Schroedinger equation can be converted into that of the Gross-Pitaevskii equation, which builds a systematical connection between the canonical solitons and the so-called nonautonomous ones. The finding of this transformation has a significant contribution to understanding the essential properties of the nonautonomous solitons and the dynamics of the Bose-Einstein condensates by using the Feshbach resonance technique

Model independence of scattering of three identical bosons in two dimensions

International Nuclear Information System (INIS)

Adhikari, S.K.; Tomio, L.; Delfino, A.; Frederico, T.

1992-05-01

Within the framework of scattering integral equations in momentum space we present numerical results of scattering of three-identical bosons at low energies in two dimensions for short-range separable potentials. An analysis of the present numerical results reveal the three-particle scattering observables to be independent of potential shape provided the low-energy two-particle binding energy and scattering length are held fixed throughout the investigation. (author)

First integrals of the axisymmetric shape equation of lipid membranes

Science.gov (United States)

Zhang, Yi-Heng; McDargh, Zachary; Tu, Zhan-Chun

2018-03-01

The shape equation of lipid membranes is a fourth-order partial differential equation. Under the axisymmetric condition, this equation was transformed into a second-order ordinary differential equation (ODE) by Zheng and Liu (Phys. Rev. E 48 2856 (1993)). Here we try to further reduce this second-order ODE to a first-order ODE. First, we invert the usual process of variational calculus, that is, we construct a Lagrangian for which the ODE is the corresponding Euler–Lagrange equation. Then, we seek symmetries of this Lagrangian according to the Noether theorem. Under a certain restriction on Lie groups of the shape equation, we find that the first integral only exists when the shape equation is identical to the Willmore equation, in which case the symmetry leading to the first integral is scale invariance. We also obtain the mechanical interpretation of the first integral by using the membrane stress tensor. Project supported by the National Natural Science Foundation of China (Grant No. 11274046) and the National Science Foundation of the United States (Grant No. 1515007).

Numerical treatments for solving nonlinear mixed integral equation

Directory of Open Access Journals (Sweden)

M.A. Abdou

2016-12-01

Full Text Available We consider a mixed type of nonlinear integral equation (MNLIE of the second kind in the space C[0,T]×L2(Ω,T<1. The Volterra integral terms (VITs are considered in time with continuous kernels, while the Fredholm integral term (FIT is considered in position with singular general kernel. Using the quadratic method and separation of variables method, we obtain a nonlinear system of Fredholm integral equations (NLSFIEs with singular kernel. A Toeplitz matrix method, in each case, is then used to obtain a nonlinear algebraic system. Numerical results are calculated when the kernels take a logarithmic form or Carleman function. Moreover, the error estimates, in each case, are then computed.

Quantum scattering via the discretisation of Schroedinger's equation

Energy Technology Data Exchange (ETDEWEB)

Alexopoulos, A. [Electronic Warfare and Radar Division, Defence Science and Technology Organisation (DSTO), PO Box 1500, Edinburgh 5111 (Australia)]. E-mail: aris.alexopoulos@dsto.defence.gov.au

2007-03-19

We obtain the scattering matrix for particles that encounter a quantum potential by discretising Schroedinger's time independent differential equation without the need to resort to the manipulation of the eigenfunctions directly. The singularities that arise in some solutions by other methods are handled with ease including the effects of resonances while convergence is excellent in all limits with only a small number of orders required to give accurate results. Our method compares the tunnelling probability with that of the WKB theory, exact numerical solutions and the modified Airy function method.

Neutron-triton scattering lengths for interactions reproducing low-energy trinucleon data

International Nuclear Information System (INIS)

Levashev, V.P.

1981-01-01

By solving the integral equations for four nucleons the neutron-triton scattering lengths and total cross section are calculated using different S-wave rank-one separable potentials. A number of linear correlations between the neutron-triton scattering lengths and triton binding energy are found. The scattering lengths consistent with low-energy trinucleon data. The results obtained are compared with available experimental data [ru

Magnetoconductivity of quantum wires with elastic and inelastic scattering

DEFF Research Database (Denmark)

Bruus, Henrik; Flensberg, Karsten; Smith

1993-01-01

We use a Boltzmann equation to determine the magnetoconductivity of quantum wires. The presence of a confining potential in addtion to the magnetic field removes the degeneracy of the Landau levels and allows one to associate a group velocity with each single-particle state. The distribution...... function describing the occupation of these single-particle states satisfies a Boltzmann equation, which may be solved exactly in the case of impurity scattering. In the case where the electrons scatter against both phonons and impurities we solve numerically—and in certain limits analytically—the integral...

Kwong-Wong-type integral equation on time scales

Directory of Open Access Journals (Sweden)

Baoguo Jia

2011-09-01

Full Text Available Consider the second-order nonlinear dynamic equation $$ [r(tx^Delta(ho(t]^Delta+p(tf(x(t=0, $$ where $p(t$ is the backward jump operator. We obtain a Kwong-Wong-type integral equation, that is: If $x(t$ is a nonoscillatory solution of the above equation on $[T_0,infty$, then the integral equation $$ frac{r^sigma(tx^Delta(t}{f(x^sigma(t} =P^sigma(t+int^infty_{sigma(t}frac{r^sigma(s [int^1_0f'(x_h(sdh][x^Delta(s]^2}{f(x(s f(x^sigma(s}Delta s $$ is satisfied for $tgeq T_0$, where $P^sigma(t=int^infty_{sigma(t}p(sDelta s$, and $x_h(s=x(s+hmu(sx^Delta(s$. As an application, we show that the superlinear dynamic equation $$ [r(tx^{Delta}(ho(t]^Delta+p(tf(x(t=0, $$ is oscillatory, under certain conditions.

Low equation, pion-nucleon scattering, and Castillejo-Dalitz-Dyson pole

International Nuclear Information System (INIS)

Nakano, K.; Nogami, Y.

1986-01-01

We examine the p-wave πN scattering at medium energies by means of the Low equation with a view to determining the form factor of the πN interaction. Solutions of the equation with and without a Castillejo-Dalitz-Dyson (CDD) pole are used. The solution with no CDD pole corresponds to the old Chew-Low model, whereas the one with a CDD pole corresponds to the quark version of the Chew-Low model. The πN interaction form factor is determined so that the Δ resonance is well reproduced. We find that the solution with a CDD pole leads to a softer form factor but is not as soft as those expected from the nucleon size in the quark model. Using the solutions and form factors thus determined, we also examine the pionic contributions to the nucleon magnetic moment and the nucleon mass

Parallel time domain solvers for electrically large transient scattering problems

KAUST Repository

Liu, Yang

2014-09-26

Marching on in time (MOT)-based integral equation solvers represent an increasingly appealing avenue for analyzing transient electromagnetic interactions with large and complex structures. MOT integral equation solvers for analyzing electromagnetic scattering from perfect electrically conducting objects are obtained by enforcing electric field boundary conditions and implicitly time advance electric surface current densities by iteratively solving sparse systems of equations at all time steps. Contrary to finite difference and element competitors, these solvers apply to nonlinear and multi-scale structures comprising geometrically intricate and deep sub-wavelength features residing atop electrically large platforms. Moreover, they are high-order accurate, stable in the low- and high-frequency limits, and applicable to conducting and penetrable structures represented by highly irregular meshes. This presentation reviews some recent advances in the parallel implementations of time domain integral equation solvers, specifically those that leverage multilevel plane-wave time-domain algorithm (PWTD) on modern manycore computer architectures including graphics processing units (GPUs) and distributed memory supercomputers. The GPU-based implementation achieves at least one order of magnitude speedups compared to serial implementations while the distributed parallel implementation are highly scalable to thousands of compute-nodes. A distributed parallel PWTD kernel has been adopted to solve time domain surface/volume integral equations (TDSIE/TDVIE) for analyzing transient scattering from large and complex-shaped perfectly electrically conducting (PEC)/dielectric objects involving ten million/tens of millions of spatial unknowns.

Advanced electromagnetics and scattering theory

CERN Document Server

2015-01-01

This book present the lecture notes used in two courses that the late Professor Kasra Barkeshli had offered at Sharif University of Technology, namely, Advanced Electromagnetics and Scattering Theory. The prerequisite for the sequence is vector calculus and electromagnetic fields and waves. Some familiarity with Green's functions and integral equations is desirable but not necessary. The book  provides a brief but concise introduction to classical topics in the field. It is divided into three parts including annexes. Part I covers principle of electromagnetic theory. The discussion starts with a review of the Maxwell's equations in differential and integral forms and basic boundary conditions. The solution of inhomogeneous wave equation and various field representations including Lorentz's potential functions and the Green's function method are discussed next. The solution of Helmholtz equation and wave harmonics follow. Next, the book presents plane wave propagation in dielectric and lossy media and various...

The Balescu kinetic equation with exchange interaction

International Nuclear Information System (INIS)

Belyi, V V; Kukharenko, Yu A

2009-01-01

Starting with the quantum BBGKY hierarchy for the distribution functions, we have obtained the quantum kinetic equation including the dynamical screening of the interaction potential, which exactly takes into account the exchange scattering in the plasma. The collision integral is expressed in terms of the Green function of the linearized Hartree–Fock equation. The potential energy takes into account the polarization and exchange interaction too

Initial states in integrable quantum field theory quenches from an integral equation hierarchy

Directory of Open Access Journals (Sweden)

D.X. Horváth

2016-01-01

Full Text Available We consider the problem of determining the initial state of integrable quantum field theory quenches in terms of the post-quench eigenstates. The corresponding overlaps are a fundamental input to most exact methods to treat integrable quantum quenches. We construct and examine an infinite integral equation hierarchy based on the form factor bootstrap, proposed earlier as a set of conditions determining the overlaps. Using quenches of the mass and interaction in Sinh-Gordon theory as a concrete example, we present theoretical arguments that the state has the squeezed coherent form expected for integrable quenches, and supporting an Ansatz for the solution of the hierarchy. Moreover we also develop an iterative method to solve numerically the lowest equation of the hierarchy. The iterative solution along with extensive numerical checks performed using the next equation of the hierarchy provides a strong numerical evidence that the proposed Ansatz gives a very good approximation for the solution.

Initial states in integrable quantum field theory quenches from an integral equation hierarchy

Energy Technology Data Exchange (ETDEWEB)

Horváth, D.X., E-mail: esoxluciuslinne@gmail.com [MTA-BME “Momentum” Statistical Field Theory Research Group, Budafoki út 8, 1111 Budapest (Hungary); Department of Theoretical Physics, Budapest University of Technology and Economics, Budafoki út 8, 1111 Budapest (Hungary); Sotiriadis, S., E-mail: sotiriad@sissa.it [SISSA and INFN, Via Bonomea 265, 34136 Trieste (Italy); Takács, G., E-mail: takacsg@eik.bme.hu [MTA-BME “Momentum” Statistical Field Theory Research Group, Budafoki út 8, 1111 Budapest (Hungary); Department of Theoretical Physics, Budapest University of Technology and Economics, Budafoki út 8, 1111 Budapest (Hungary)

2016-01-15

We consider the problem of determining the initial state of integrable quantum field theory quenches in terms of the post-quench eigenstates. The corresponding overlaps are a fundamental input to most exact methods to treat integrable quantum quenches. We construct and examine an infinite integral equation hierarchy based on the form factor bootstrap, proposed earlier as a set of conditions determining the overlaps. Using quenches of the mass and interaction in Sinh-Gordon theory as a concrete example, we present theoretical arguments that the state has the squeezed coherent form expected for integrable quenches, and supporting an Ansatz for the solution of the hierarchy. Moreover we also develop an iterative method to solve numerically the lowest equation of the hierarchy. The iterative solution along with extensive numerical checks performed using the next equation of the hierarchy provides a strong numerical evidence that the proposed Ansatz gives a very good approximation for the solution.

Coupling Integrable Couplings of an Equation Hierarchy

International Nuclear Information System (INIS)

Wang Hui; Xia Tie-Cheng

2013-01-01

Based on a kind of Lie algebra G proposed by Zhang, one isospectral problem is designed. Under the framework of zero curvature equation, a new kind of integrable coupling of an equation hierarchy is generated using the methods proposed by Ma and Gao. With the help of variational identity, we get the Hamiltonian structure of the hierarchy. (general)

Multi-component bi-Hamiltonian Dirac integrable equations

Energy Technology Data Exchange (ETDEWEB)

Ma Wenxiu [Department of Mathematics and Statistics, University of South Florida, Tampa, FL 33620-5700 (United States)], E-mail: mawx@math.usf.edu

2009-01-15

A specific matrix iso-spectral problem of arbitrary order is introduced and an associated hierarchy of multi-component Dirac integrable equations is constructed within the framework of zero curvature equations. The bi-Hamiltonian structure of the obtained Dirac hierarchy is presented be means of the variational trace identity. Two examples in the cases of lower order are computed.

Quadratic algebras in the noncommutative integration method of wave equation

International Nuclear Information System (INIS)

Varaksin, O.L.

1995-01-01

The paper deals with the investigation of applications of the method of noncommutative integration of linear differential equations by partial derivatives. Nontrivial example was taken for integration of three-dimensions wave equation with the use of non-Abelian quadratic algebras

Chew-Low equations as Cremoma transformations

International Nuclear Information System (INIS)

Rerikh, K.V.

1982-01-01

The Chew-Low equations for the p-wave pion-nucleon scattering with the crossing-symmetry matrix (3x3) are investigated in their well-known formulation as a system of nonlinear difference equations. These equations interpreted as geometrical transformations are shown to be a special case of the Cremona transformaions. Using the properties of the Cremona transformations we obtain the general 3-parametric functional equation on invariant algebraic and nonalgebraic curves in the space solutions of the Chew- Low equations. It is proved that there exists only one invariant algebraic curve, the parabola corresponding to the well-known solution. Analysis of the general functional equation on invariant nonalgebraic curves makes it possible to select in addition to this parabola 3 invariant forms defining implicitly 3 nonalgebraic curves and to concretize for them the general equation by means of fixing the parameters. From the transformational properties of the invariant forms with respect to the Cremona transformations, there follows an important result that the ration of these forms in proper powers is the general integral of the nonlinear system of the Chew-Low equations, which is an even antiperiodic function. The structure of the second general integral is given and the functional equations which determinne this integral are presented [ru

Integral equations for composite-particle scattering taking the Pauli principle into account

International Nuclear Information System (INIS)

Kukulin, V.I.; Neudatchin, V.G.; Pomerantsev, V.N.

1978-01-01

An approximate description of a system of three composite particles in terms of the Saito (Prog. Theor. Phys.; 41:705 (1969)) orthogonality condition model is proposed. The orthogonalising pseudopotential technique is used to derive a modified set of Fadde'ev equations where the two- and three-body exchanges due to the Pauli principle are included by orthogonalising to two-and three-body forbidden states. The scope of applicability of and the method for solving the derived equations are discussed briefly. (author)

Analytic solution of integral equations for molecular fluids

International Nuclear Information System (INIS)

Cummings, P.T.

1984-01-01

We review some recent progress in the analytic solution of integral equations for molecular fluids. The site-site Ornstein-Zernike (SSOZ) equation with approximate closures appropriate to homonuclear diatomic fluids both with and without attractive dispersion-like interactions has recently been solved in closed form analytically. In this paper, the close relationship between the SSOZ equation for homonuclear dumbells and the usual Ornstein-Zernike (OZ) equation for atomic fluids is carefully elucidated. This relationship is a key motivation for the analytic solutions of the SSOZ equation that have been obtained to date. (author)

Differential equations and integrable models: the SU(3) case

International Nuclear Information System (INIS)

Dorey, Patrick; Tateo, Roberto

2000-01-01

We exhibit a relationship between the massless a 2 (2) integrable quantum field theory and a certain third-order ordinary differential equation, thereby extending a recent result connecting the massless sine-Gordon model to the Schroedinger equation. This forms part of a more general correspondence involving A 2 -related Bethe ansatz systems and third-order differential equations. A non-linear integral equation for the generalised spectral problem is derived, and some numerical checks are performed. Duality properties are discussed, and a simple variant of the non-linear equation is suggested as a candidate to describe the finite volume ground state energies of minimal conformal field theories perturbed by the operators phi 12 , phi 21 and phi 15 . This is checked against previous results obtained using the thermodynamic Bethe ansatz

Partially integrable nonlinear equations with one higher symmetry

International Nuclear Information System (INIS)

Mikhailov, A V; Novikov, V S; Wang, J P

2005-01-01

In this letter, we present a family of second order in time nonlinear partial differential equations, which have only one higher symmetry. These equations are not integrable, but have a solution depending on one arbitrary function. (letter to the editor)

On the integrability of the generalized Fisher-type nonlinear diffusion equations

International Nuclear Information System (INIS)

Wang Dengshan; Zhang Zhifei

2009-01-01

In this paper, the geometric integrability and Lax integrability of the generalized Fisher-type nonlinear diffusion equations with modified diffusion in (1+1) and (2+1) dimensions are studied by the pseudo-spherical surface geometry method and prolongation technique. It is shown that the (1+1)-dimensional Fisher-type nonlinear diffusion equation is geometrically integrable in the sense of describing a pseudo-spherical surface of constant curvature -1 only for m = 2, and the generalized Fisher-type nonlinear diffusion equations in (1+1) and (2+1) dimensions are Lax integrable only for m = 2. This paper extends the results in Bindu et al 2001 (J. Phys. A: Math. Gen. 34 L689) and further provides the integrability information of (1+1)- and (2+1)-dimensional Fisher-type nonlinear diffusion equations for m = 2

Numerical integration of asymptotic solutions of ordinary differential equations

Science.gov (United States)

Thurston, Gaylen A.

1989-01-01

Classical asymptotic analysis of ordinary differential equations derives approximate solutions that are numerically stable. However, the analysis also leads to tedious expansions in powers of the relevant parameter for a particular problem. The expansions are replaced with integrals that can be evaluated by numerical integration. The resulting numerical solutions retain the linear independence that is the main advantage of asymptotic solutions. Examples, including the Falkner-Skan equation from laminar boundary layer theory, illustrate the method of asymptotic analysis with numerical integration.

A Fortran program for the numerical integration of the one-dimensional Schroedinger equation using exponential and Bessel fitting methods

International Nuclear Information System (INIS)

Cash, J.R.; Raptis, A.D.; Simos, T.E.

1990-01-01

An efficient algorithm is described for the accurate numerical integration of the one-dimensional Schroedinger equation. This algorithm uses a high-order, variable step Runge-Kutta like method in the region where the potential term dominates, and an exponential or Bessel fitted method in the asymptotic region. This approach can be used to compute scattering phase shifts in an efficient and reliable manner. A Fortran program which implements this algorithm is provided and some test results are given. (orig.)

Energy preserving integration of bi-Hamiltonian partial differential equations

NARCIS (Netherlands)

Karasozen, B.; Simsek, G.

2013-01-01

The energy preserving average vector field (AVF) integrator is applied to evolutionary partial differential equations (PDEs) in bi-Hamiltonian form with nonconstant Poisson structures. Numerical results for the Korteweg de Vries (KdV) equation and for the Ito type coupled KdV equation confirm the

Simplifying Differential Equations for Multiscale Feynman Integrals beyond Multiple Polylogarithms.

Science.gov (United States)

Adams, Luise; Chaubey, Ekta; Weinzierl, Stefan

2017-04-07

In this Letter we exploit factorization properties of Picard-Fuchs operators to decouple differential equations for multiscale Feynman integrals. The algorithm reduces the differential equations to blocks of the size of the order of the irreducible factors of the Picard-Fuchs operator. As a side product, our method can be used to easily convert the differential equations for Feynman integrals which evaluate to multiple polylogarithms to an ϵ form.

Discrete Painlevé equations: an integrability paradigm

International Nuclear Information System (INIS)

Grammaticos, B; Ramani, A

2014-01-01

In this paper we present a review of results on discrete Painlevé equations. We begin with an introduction which serves as a refresher on the continuous Painlevé equations. Next, in the first, main part of the paper, we introduce the discrete Painlevé equations, the various methods for their derivation, and their properties as well as their classification scheme. Along the way we present a brief summary of the two major discrete integrability detectors and of Quispel–Roberts–Thompson mapping, which plays a primordial role in the derivation of discrete Painlevé equations. The second part of the paper is more technical and focuses on the presentation of new results on what are called asymmetric discrete Painlevé equations. (comment)

Lectures on differential equations for Feynman integrals

International Nuclear Information System (INIS)

Henn, Johannes M

2015-01-01

Over the last year significant progress was made in the understanding of the computation of Feynman integrals using differential equations (DE). These lectures give a review of these developments, while not assuming any prior knowledge of the subject. After an introduction to DE for Feynman integrals, we point out how they can be simplified using algorithms available in the mathematical literature. We discuss how this is related to a recent conjecture for a canonical form of the equations. We also discuss a complementary approach that is based on properties of the space–time loop integrands, and explain how the ideas of leading singularities and d-log representations can be used to find an optimal basis for the DE. Finally, as an application of these ideas we show how single-scale integrals can be bootstrapped using the Drinfeld associator of a DE. (topical review)

Bargmann Symmetry Constraint for a Family of Liouville Integrable Differential-Difference Equations

International Nuclear Information System (INIS)

Xu Xixiang

2012-01-01

A family of integrable differential-difference equations is derived from a new matrix spectral problem. The Hamiltonian forms of obtained differential-difference equations are constructed. The Liouville integrability for the obtained integrable family is proved. Then, Bargmann symmetry constraint of the obtained integrable family is presented by binary nonliearization method of Lax pairs and adjoint Lax pairs. Under this Bargmann symmetry constraints, an integrable symplectic map and a sequences of completely integrable finite-dimensional Hamiltonian systems in Liouville sense are worked out, and every integrable differential-difference equations in the obtained family is factored by the integrable symplectic map and a completely integrable finite-dimensional Hamiltonian system. (general)

Integrable semi-discretizations of the reduced Ostrovsky equation

International Nuclear Information System (INIS)

Feng, Bao-Feng; Maruno, Ken-ichi; Ohta, Yasuhiro

2015-01-01

Based on our previous work on the reduced Ostrovsky equation (J. Phys. A: Math. Theor. 45 355203), we construct its integrable semi-discretizations. Since the reduced Ostrovsky equation admits two alternative representations, one being its original form, the other the differentiated form (the short wave limit of the Degasperis–Procesi equation) two semi-discrete analogues of the reduced Ostrovsky equation are constructed possessing the same N-loop soliton solution. The relationship between these two versions of semi-discretizations is also clarified. (paper)

Path-integral approach to resonant electron-molecule scattering

International Nuclear Information System (INIS)

Winterstetter, M.; Domcke, W.

1993-01-01

A path-integral formulation of resonant electron-molecule scattering is developed within the framework of the projection-operator formalism of scattering theory. The formation and decay of resonances is treated in real time as a quantum-mechanical electronic-tunneling process, modified by the coupling of the electronic motion with the nuclear degrees of freedom. It is shown that the electronic continuum can be summed over in the path-integral formulation, resulting formally in the path integral for an effective two-state system with coupling to vibrations. The harmonic-oscillator approximation is adopted for the vibrational motion in the present work. Approximation methods are introduced which render the numerical evaluation of the sum over paths feasible for up to ∼10 3 elementary time slices. The theory is numerically realized for simple but nontrivial models representing the 2 Π g d-wave shape resonance in e - +N 2 collisions and the 2 Σ u + p-wave shape resonance in e - +H 2 collisions, respectively. The accuracy of the path-integral results is assessed by comparison with exact numerical reference data for these models. The essential virtue of the path-integral approach is the fact that the computational effort scales at most linearly with the number of vibrational degrees of freedom. The path-integral method is thus well suited to treat electron collisions with polyatomic molecules and molecular aggregates

Block-pulse functions approach to numerical solution of Abel’s integral equation

Directory of Open Access Journals (Sweden)

Monireh Nosrati Sahlan

2015-12-01

Full Text Available This study aims to present a computational method for solving Abel’s integral equation of the second kind. The introduced method is based on the use of Block-pulse functions (BPFs via collocation method. Abel’s integral equations as singular Volterra integral equations are hard and heavy in computation, but because of the properties of BPFs, as is reported in examples, this method is more efficient and more accurate than some other methods for solving this class of integral equations. On the other hand, the benefit of this method is low cost of computing operations. The applied method transforms the singular integral equation into triangular linear algebraic system that can be solved easily. An error analysis is worked out and applications are demonstrated through illustrative examples.

Iterative numerical solution of scattering problems

International Nuclear Information System (INIS)

Tomio, L.; Adhikari, S.K.

1995-05-01

An iterative Neumann series method, employing a real auxiliary scattering integral equation, is used to calculate scattering lengths and phase shifts for the atomic Yukawa and exponential potentials. For these potentials the original Neumann series diverges. The present iterative method yields results that are far better, in convergence, stability and precision, than other momentum space methods. Accurate result is obtained in both cases with an estimated error of about 1 in 10 10 after some-8-10 iterations. (author). 31 refs, 2 tabs

Nonlinear Fredholm Integral Equation of the Second Kind with Quadrature Methods

Directory of Open Access Journals (Sweden)

M. Jafari Emamzadeh

2010-06-01

Full Text Available In this paper, a numerical method for solving the nonlinear Fredholm integral equation is presented. We intend to approximate the solution of this equation by quadrature methods and by doing so, we solve the nonlinear Fredholm integral equation more accurately. Several examples are given at the end of this paper

Integration of differential equations by the pseudo-linear (PL) approximation

International Nuclear Information System (INIS)

Bonalumi, Riccardo A.

1998-01-01

A new method of integrating differential equations was originated with the technique of approximately calculating the integrals called the pseudo-linear (PL) procedure: this method is A-stable. This article contains the following examples: 1st order ordinary differential equations (ODEs), 2nd order linear ODEs, stiff system of ODEs (neutron kinetics), one-dimensional parabolic (diffusion) partial differential equations. In this latter case, this PL method coincides with the Crank-Nicholson method

Coupled wave equations theory of surface-enhanced femtosecond stimulated Raman scattering.

Science.gov (United States)

McAnally, Michael O; McMahon, Jeffrey M; Van Duyne, Richard P; Schatz, George C

2016-09-07

We present a coupled wave semiclassical theory to describe plasmonic enhancement effects in surface-enhanced femtosecond stimulated Raman scattering (SE-FSRS). A key result is that the plasmon enhanced fields which drive the vibrational equation of motion for each normal mode results in dispersive lineshapes in the SE-FSRS spectrum. This result, which reproduces experimental lineshapes, demonstrates that plasmon-enhanced stimulated Raman methods provide unique sensitivity to a plasmonic response. Our derived SE-FSRS theory shows a plasmonic enhancement of |gpu|(2)ImχR(ω)gst (2)/ImχR(ω), where |gpu|(2) is the absolute square of the plasmonic enhancement from the Raman pump, χR(ω) is the Raman susceptibility, and gst is the plasmonic enhancement of the Stokes field in SE-FSRS. We conclude with a discussion on potential future experimental and theoretical directions for the field of plasmonically enhanced coherent Raman scattering.

On the complete integrability of the discrete Nahm equations

International Nuclear Information System (INIS)

Murray, M.K.

2000-01-01

The discrete Nahm equations, a system of matrix valued difference equations, arose in the work of Braam and Austin on half-integral mass hyperbolic monopoles. We show that the discrete Nahm equations are completely integrable in a natural sense: to any solution we can associate a spectral curve and a holomorphic line-bundle over the spectral curve, such that the discrete-time DN evolution corresponds to walking in the Jacobian of the spectral curve in a straight line through the line-bundle with steps of a fixed size. Some of the implications for hyperbolic monopoles are also discussed. (orig.)

Quantum statistics of stimulated Raman and hyper-Raman scattering by master equation approach

International Nuclear Information System (INIS)

Gupta, P.S.; Dash, J.

1991-01-01

A quantum theoretical density matrix formalism of stimulated Raman and hyper-Raman scattering using master equation approach is presented. The atomic system is described by two energy levels. The effects of upper level population and the cavity loss are incorporated. The photon statistics, coherence characteristics and the building up of the Stokes field are investigated. (author). 8 figs., 5 refs

Evaluating four-loop conformal Feynman integrals by D-dimensional differential equations

Science.gov (United States)

Eden, Burkhard; Smirnov, Vladimir A.

2016-10-01

We evaluate a four-loop conformal integral, i.e. an integral over four four-dimensional coordinates, by turning to its dimensionally regularized version and applying differential equations for the set of the corresponding 213 master integrals. To solve these linear differential equations we follow the strategy suggested by Henn and switch to a uniformly transcendental basis of master integrals. We find a solution to these equations up to weight eight in terms of multiple polylogarithms. Further, we present an analytical result for the given four-loop conformal integral considered in four-dimensional space-time in terms of single-valued harmonic polylogarithms. As a by-product, we obtain analytical results for all the other 212 master integrals within dimensional regularization, i.e. considered in D dimensions.

Evaluating four-loop conformal Feynman integrals by D-dimensional differential equations

Energy Technology Data Exchange (ETDEWEB)

Eden, Burkhard [Institut für Mathematik und Physik, Humboldt-Universität zu Berlin,Zum großen Windkanal 6, 12489 Berlin (Germany); Smirnov, Vladimir A. [Skobeltsyn Institute of Nuclear Physics, Moscow State University,119992 Moscow (Russian Federation)

2016-10-21

We evaluate a four-loop conformal integral, i.e. an integral over four four-dimensional coordinates, by turning to its dimensionally regularized version and applying differential equations for the set of the corresponding 213 master integrals. To solve these linear differential equations we follow the strategy suggested by Henn and switch to a uniformly transcendental basis of master integrals. We find a solution to these equations up to weight eight in terms of multiple polylogarithms. Further, we present an analytical result for the given four-loop conformal integral considered in four-dimensional space-time in terms of single-valued harmonic polylogarithms. As a by-product, we obtain analytical results for all the other 212 master integrals within dimensional regularization, i.e. considered in D dimensions.

Partial differential equations an introduction

CERN Document Server

Colton, David

2004-01-01

Intended for a college senior or first-year graduate-level course in partial differential equations, this text offers students in mathematics, engineering, and the applied sciences a solid foundation for advanced studies in mathematics. Classical topics presented in a modern context include coverage of integral equations and basic scattering theory. This complete and accessible treatment includes a variety of examples of inverse problems arising from improperly posed applications. Exercises at the ends of chapters, many with answers, offer a clear progression in developing an understanding of

Unconditionally stable integration of Maxwell's equations

NARCIS (Netherlands)

J.G. Verwer (Jan); M.A. Botchev

2008-01-01

htmlabstractNumerical integration of Maxwell''s equations is often based on explicit methods accepting a stability step size restriction. In literature evidence is given that there is also a need for unconditionally stable methods, as exemplified by the successful alternating direction

Is Yang-Mills equation a totally integrable system. Lecture III

International Nuclear Information System (INIS)

Chau Wang, L.L.

1981-01-01

Topics covered include: loop-space formulation of gauge theory - loop-space chiral equation; two dimensional chiral equation - conservation laws, linear system and integrability; and parallel development for the loop-space chiral equation - subtlety

Solving large sets of coupled equations iteratively by vector processing on the CYBER 205 computer

International Nuclear Information System (INIS)

Tolsma, L.D.

1985-01-01

The set of coupled linear second-order differential equations which has to be solved for the quantum-mechanical description of inelastic scattering of atomic and nuclear particles can be rewritten as an equivalent set of coupled integral equations. When some type of functions is used as piecewise analytic reference solutions, the integrals that arise in this set can be evaluated analytically. The set of integral equations can be solved iteratively. For the results mentioned an inward-outward iteration scheme has been applied. A concept of vectorization of coupled-channel Fortran programs, based on this integral method, is presented for the use on the Cyber 205 computer. It turns out that, for two heavy ion nuclear scattering test cases, this vector algorithm gives an overall speed-up of about a factor of 2 to 3 compared to a highly optimized scalar algorithm for a one vector pipeline computer

Diffraction scattering of strongly bound system

International Nuclear Information System (INIS)

Kuzmichev, V.E.

1982-04-01

The scattering of a hadron on a strongly bound system of two hadrons (dihadron) is considered in the high-energy limit for the relative hadron-dihadron motion. The dihadron scatterer motion and the internal interaction are included in our consideration. It is shown that only small values of the internal transfer momentum of dihadron particles bring the principal contribution to the three-particle propagator in eikonal approximation. On the basis of the exact analytical solution of the integral equation for the total Green function the scattering amplitude is derived. It is shown that the scattering amplitude contains only single, double, and triple scattering terms. The three new terms to the Glauber formula for the total cross section are obtained. These terms decrease both the true total hadron-hadron cross section and the screening correction. (orig.)

Integration of equations of parabolic type by the method of nets

CERN Document Server

Saul'Yev, V K; Stark, M; Ulam, S

1964-01-01

International Series of Monographs in Pure and Applied Mathematics, Volume 54: Integration of Equations of Parabolic Type by the Method of Nets deals with solving parabolic partial differential equations using the method of nets. The first part of this volume focuses on the construction of net equations, with emphasis on the stability and accuracy of the approximating net equations. The method of nets or method of finite differences (used to define the corresponding numerical method in ordinary differential equations) is one of many different approximate methods of integration of partial diff

Geometry, Heat Equation and Path Integrals on the Poincare Upper Half-Plane

OpenAIRE

Reijiro, KUBO; Research Institute for Theoretical Physics Hiroshima University

1988-01-01

Geometry, heat equation and Feynman's path integrals are studied on the Poincare upper half-plane. The fundamental solution to the heat equation ∂f/∂t=Δ_Hf is expressed in terms of a path integral defined on the upper half-plane. It is shown that Kac's statement that Feynman's path integral satisfies the Schrodinger equation is also valid for our case.

ZZ THERMOS, Multigroup P0 to P5 Thermal Scattering Kernels from ENDF/B Scattering Law Data

International Nuclear Information System (INIS)

McCrosson, F.J.; Finch, D.R.

1975-01-01

1 - Description of problem or function: Number of groups: 30-group THERMOS thermal scattering kernels. Nuclides: Molecular H 2 O, Molecular D 2 O, Graphite, Polyethylene, Benzene, Zr bound in ZrHx, H bound in ZrHx, Beryllium-9, Beryllium Oxide, Uranium Dioxide. Origin: ENDF/B library. Weighting Spectrum: yes. These data are 30-group THERMOS thermal scattering kernels for P0 to P5 Legendre orders for every temperature of every material from s(alpha,beta) data stored in the ENDF/B library. These scattering kernels were generated using the FLANGE2 computer code (NESC Abstract 368). To test the kernels, the integral properties of each set of kernels were determined by a precision integration of the diffusion length equation and compared to experimental measurements of these properties. In general, the agreement was very good. Details of the methods used and results obtained are contained in the reference. The scattering kernels are organized into a two volume magnetic tape library from which they may be retrieved easily for use in any 30-group THERMOS library. The contents of the tapes are as follows - (Material: ZA/Temperatures (degrees K)): Molecular H 2 O: 100.0/296, 350, 400, 450, 500, 600, Molecular D 2 O: 101.0/296, 350, 400, 450, 500, 600, Graphite: 6000.0/296, 400, 500, 600, 700, 800, Polyethylene: 205.0/296, 350 Benzene: 106.0/296, 350, 400, 450, 500, 600, Zr bound in ZrHx: 203.0/296, 400, 500, 600, 700, 800, H bound in ZrHx: 230.0/296, 400, 500, 600, 700, 800, Beryllium-9: 4009.0/296, 400, 500, 600, 700, 800, Beryllium Oxide: 200.0/296, 400, 500, 600, 700, 800, Uranium Dioxide: 207.0/296, 400, 500, 600, 700, 800 2 - Method of solution: Kernel generation is performed by direct integration of the thermal scattering law data to obtain the differential scattering cross sections for each Legendre order. The integral parameter calculation is done by precision integration of the diffusion length equation for several moderator absorption cross sections followed by a

Feynman path integrals - from the prodistribution definition to the calculation of glory scattering

International Nuclear Information System (INIS)

DeWitt-Morette, C.

1984-01-01

In these lectures I present a path integral calculation, starting from a global definition of Feynman path integrals and ending at a scattering cross section formula. Along the way I discuss some basic issues which had to be resolved to exploit the computational power of the proposed definition of Feynman integrals. I propose to compute the glory scattering of gravitational waves by black holes. (orig./HSI)

Fast analysis of wide-band scattering from electrically large targets with time-domain parabolic equation method

Science.gov (United States)

He, Zi; Chen, Ru-Shan

2016-03-01

An efficient three-dimensional time domain parabolic equation (TDPE) method is proposed to fast analyze the narrow-angle wideband EM scattering properties of electrically large targets. The finite difference (FD) of Crank-Nicolson (CN) scheme is used as the traditional tool to solve the time-domain parabolic equation. However, a huge computational resource is required when the meshes become dense. Therefore, the alternating direction implicit (ADI) scheme is introduced to discretize the time-domain parabolic equation. In this way, the reduced transient scattered fields can be calculated line by line in each transverse plane for any time step with unconditional stability. As a result, less computational resources are required for the proposed ADI-based TDPE method when compared with both the traditional CN-based TDPE method and the finite-different time-domain (FDTD) method. By employing the rotating TDPE method, the complete bistatic RCS can be obtained with encouraging accuracy for any observed angle. Numerical examples are given to demonstrate the accuracy and efficiency of the proposed method.

On the completeness of the natural modes for quantum mechanical potential scattering

NARCIS (Netherlands)

Hoenders, B.J.

1979-01-01

The set of natural modes, associated with quantum mechanical scattering from a central potential of finite-range is shown to be complete. The natural modes satisfy a non-Hermitian homogeneous integral equation, or alternatively, are solutions of the time independent Schrödinger equation subject to a

Uniform semiclassical approximation for absorptive scattering systems

International Nuclear Information System (INIS)

Hussein, M.S.; Pato, M.P.

1987-07-01

The uniform semiclassical approximation of the elastic scattering amplitude is generalized to absorptive systems. An integral equation is derived which connects the absorption modified amplitude to the absorption free one. Division of the amplitude into a diffractive and refractive components is then made possible. (Author) [pt

Unconditionally stable integration of Maxwell's equations

NARCIS (Netherlands)

Verwer, J.G.; Bochev, Mikhail A.

Numerical integration of Maxwell's equations is often based on explicit methods accepting a stability step size restriction. In literature evidence is given that there is also a need for unconditionally stable methods, as exemplified by the successful alternating direction implicit finite difference

Unconditionally stable integration of Maxwell's equations

NARCIS (Netherlands)

J.G. Verwer (Jan); M.A. Botchev

2009-01-01

textabstractNumerical integration of Maxwell’s equations is often based on explicit methods accepting a stability step size restriction. In literature evidence is given that there is also a need for unconditionally stable methods, as exemplified by the successful alternating direction implicit –

Stability analysis of implicit time discretizations for the Compton-scattering Fokker-Planck equation

International Nuclear Information System (INIS)

Densmore, Jeffery D.; Warsa, James S.; Lowrie, Robert B.; Morel, Jim E.

2009-01-01

The Fokker-Planck equation is a widely used approximation for modeling the Compton scattering of photons in high energy density applications. In this paper, we perform a stability analysis of three implicit time discretizations for the Compton-Scattering Fokker-Planck equation. Specifically, we examine (i) a Semi-Implicit (SI) scheme that employs backward-Euler differencing but evaluates temperature-dependent coefficients at their beginning-of-time-step values, (ii) a Fully Implicit (FI) discretization that instead evaluates temperature-dependent coefficients at their end-of-time-step values, and (iii) a Linearized Implicit (LI) scheme, which is developed by linearizing the temperature dependence of the FI discretization within each time step. Our stability analysis shows that the FI and LI schemes are unconditionally stable and cannot generate oscillatory solutions regardless of time-step size, whereas the SI discretization can suffer from instabilities and nonphysical oscillations for sufficiently large time steps. With the results of this analysis, we present time-step limits for the SI scheme that prevent undesirable behavior. We test the validity of our stability analysis and time-step limits with a set of numerical examples.

Stability analysis of implicit time discretizations for the Compton-scattering Fokker-Planck equation

Energy Technology Data Exchange (ETDEWEB)

Densmore, Jeffery D [Los Alamos National Laboratory; Warsa, James S [Los Alamos National Laboratory; Lowrie, Robert B [Los Alamos National Laboratory; Morel, Jim E [TEXAS A& M UNIV

2008-01-01

The Fokker-Planck equation is a widely used approximation for modeling the Compton scattering of photons in high energy density applications. In this paper, we perform a stability analysis of three implicit time discretizations for the Compton-Scattering Fokker-Planck equation. Specifically, we examine (i) a Semi-Implicit (SI) scheme that employs backward-Euler differencing but evaluates temperature-dependent coefficients at their beginning-of-time-step values, (ii) a Fully Implicit (FI) discretization that instead evaluates temperature-dependent coefficients at their end-of-time-step values, and (iii) a Linearized Implicit (LI) scheme, which is developed by linearizing the temperature dependence of the FI discretization within each time step. Our stability analysis shows that the FI and LI schemes are unconditionally stable and cannot generate oscillatory solutions regardless of time-step size, whereas the SI discretization can suffer from instabilities and nonphysical oscillations for sufficiently large time steps. With the results of this analysis, we present time-step limits for the SI scheme that prevent undesirable behavior. We test the validity of our stability analysis and time-step limits with a set of numerical examples.

Stability analysis of implicit time discretizations for the Compton-scattering Fokker-Planck equation

Science.gov (United States)

Densmore, Jeffery D.; Warsa, James S.; Lowrie, Robert B.; Morel, Jim E.

2009-09-01

The Fokker-Planck equation is a widely used approximation for modeling the Compton scattering of photons in high energy density applications. In this paper, we perform a stability analysis of three implicit time discretizations for the Compton-Scattering Fokker-Planck equation. Specifically, we examine (i) a Semi-Implicit (SI) scheme that employs backward-Euler differencing but evaluates temperature-dependent coefficients at their beginning-of-time-step values, (ii) a Fully Implicit (FI) discretization that instead evaluates temperature-dependent coefficients at their end-of-time-step values, and (iii) a Linearized Implicit (LI) scheme, which is developed by linearizing the temperature dependence of the FI discretization within each time step. Our stability analysis shows that the FI and LI schemes are unconditionally stable and cannot generate oscillatory solutions regardless of time-step size, whereas the SI discretization can suffer from instabilities and nonphysical oscillations for sufficiently large time steps. With the results of this analysis, we present time-step limits for the SI scheme that prevent undesirable behavior. We test the validity of our stability analysis and time-step limits with a set of numerical examples.

Numerical solution of the potential problem by integral equations without Green's functions

International Nuclear Information System (INIS)

De Mey, G.

1977-01-01

An integral equation technique will be presented to solve Laplace's equation in a two-dimensional area S. The Green's function has been replaced by a particular solution of Laplace equation in order to establish the integral equation. It is shown that accurate results can be obtained provided the pivotal elimination method is used to solve the linear algebraic set

Iterative numerical solution of scattering problems

Energy Technology Data Exchange (ETDEWEB)

Tomio, L; Adhikari, S K

1995-05-01

An iterative Neumann series method, employing a real auxiliary scattering integral equation, is used to calculate scattering lengths and phase shifts for the atomic Yukawa and exponential potentials. For these potentials the original Neumann series diverges. The present iterative method yields results that are far better, in convergence, stability and precision, than other momentum space methods. Accurate result is obtained in both cases with an estimated error of about 1 in 10{sup 10} after some-8-10 iterations. (author). 31 refs, 2 tabs.

The Cauchy problem for the Pavlov equation

Science.gov (United States)

Grinevich, P. G.; Santini, P. M.; Wu, D.

2015-10-01

Commutation of multidimensional vector fields leads to integrable nonlinear dispersionless PDEs that arise in various problems of mathematical physics and have been intensively studied in recent literature. This report aims to solve the scattering and inverse scattering problem for integrable dispersionless PDEs, recently introduced just at a formal level, concentrating on the prototypical example of the Pavlov equation, and to justify an existence theorem for global bounded solutions of the associated Cauchy problem with small data. An essential part of this work was made during the visit of the three authors to the Centro Internacional de Ciencias in Cuernavaca, Mexico in November-December 2012.

Integral propagator solvers for Vlasov-Fokker-Planck equations

International Nuclear Information System (INIS)

Donoso, J M; Rio, E del

2007-01-01

We briefly discuss the use of short-time integral propagators on solving the so-called Vlasov-Fokker-Planck equation for the dynamics of a distribution function. For this equation, the diffusion tensor is singular and the usual Gaussian representation of the short-time propagator is no longer valid. However, we prove that the path-integral approach on solving the equation is, in fact, reliable by means of our generalized propagator, which is obtained through the construction of an auxiliary solvable Fokker-Planck equation. The new representation of the grid-free advancing scheme describes the inherent cross- and self-diffusion processes, in both velocity and configuration spaces, in a natural manner, although these processes are not explicitly depicted in the differential equation. We also show that some splitting methods, as well as some finite-difference schemes, could fail in describing the aforementioned diffusion processes, governed in the whole phase space only by the velocity diffusion tensor. The short-time transition probability offers a stable and robust numerical algorithm that preserves the distribution positiveness and its norm, ensuring the smoothness of the evolving solution at any time step. (fast track communication)

The Integral Equation Method and the Neumann Problem for the Poisson Equation on NTA Domains

Czech Academy of Sciences Publication Activity Database

Medková, Dagmar

2009-01-01

Roč. 63, č. 21 (2009), s. 227-247 ISSN 0378-620X Institutional research plan: CEZ:AV0Z10190503 Keywords : Poisson equation * Neumann problem * integral equation method Subject RIV: BA - General Mathematics Impact factor: 0.477, year: 2009

Geometry, heat equation and path integrals on the Poincare upper half-plane

International Nuclear Information System (INIS)

Kubo, Reijiro.

1987-08-01

Geometry, heat equation and Feynman's path integrals are studied on the Poincare upper half-plane. The fundamental solution to the heat equation δf/δt = Δ H f is expressed in terms of a path integral defined on the upper half-plane. It is shown that Kac's proof that Feynman's path integral satisfies the Schroedinger equation is also valid for our case. (author)

On discrete 2D integrable equations of higher order

International Nuclear Information System (INIS)

Adler, V E; Postnikov, V V

2014-01-01

We study two-dimensional discrete integrable equations of order 1 with respect to one independent variable and m with respect to another one. A generalization of the multidimensional consistency property is proposed for this type of equations. The examples are related to the Bäcklund–Darboux transformations for the lattice equations of Bogoyavlensky type. (paper)

A Non-Perturbative, Finite Particle Number Approach to Relativistic Scattering Theory

Energy Technology Data Exchange (ETDEWEB)

Lindesay, James V

2001-05-11

We present integral equations for the scattering amplitudes of three scalar particles, using the Faddeev channel decomposition, which can be readily extended to any finite number of particles of any helicity. The solution of these equations, which have been demonstrated to be calculable, provide a non-perturbative way of obtaining relativistic scattering amplitudes for any finite number of particles that are Lorentz invariant, unitary, cluster decomposable and reduce unambiguously in the non-relativistic limit to the non-relativistic Faddeev equations. The aim of this program is to develop equations which explicitly depend upon physically observable input variables, and do not require ''renormalization'' or ''dressing'' of these parameters to connect them to the boundary states.

Some physical applications of fractional Schroedinger equation

International Nuclear Information System (INIS)

Guo Xiaoyi; Xu Mingyu

2006-01-01

The fractional Schroedinger equation is solved for a free particle and for an infinite square potential well. The fundamental solution of the Cauchy problem for a free particle, the energy levels and the normalized wave functions of a particle in a potential well are obtained. In the barrier penetration problem, the reflection coefficient and transmission coefficient of a particle from a rectangular potential wall is determined. In the quantum scattering problem, according to the fractional Schroedinger equation, the Green's function of the Lippmann-Schwinger integral equation is given

Exponential Convergence for Numerical Solution of Integral Equations Using Radial Basis Functions

Directory of Open Access Journals (Sweden)

Zakieh Avazzadeh

2014-01-01

Full Text Available We solve some different type of Urysohn integral equations by using the radial basis functions. These types include the linear and nonlinear Fredholm, Volterra, and mixed Volterra-Fredholm integral equations. Our main aim is to investigate the rate of convergence to solve these equations using the radial basis functions which have normic structure that utilize approximation in higher dimensions. Of course, the use of this method often leads to ill-posed systems. Thus we propose an algorithm to improve the results. Numerical results show that this method leads to the exponential convergence for solving integral equations as it was already confirmed for partial and ordinary differential equations.

Integration of the three-dimensional Vlasov equation for a magnetized plasma

International Nuclear Information System (INIS)

Cheng, C.Z.

1976-04-01

A second order splitting scheme is developed to integrate the three dimensional Vlasov equation for a plasma in a magnetic field. The integration of the Vlasov equation is divided into a series of intermediate steps and Fourier interpolation and the ASD method with a third order Taylor expansion are used to integrate the fractional equations. Numerical experiments related to cyclotron waves in 2 and 2 1 / 2 D are demonstrated with high accuracy and efficiency. The computer storage requirements are modest; for example, a typical 2D nonlinear electron plasma simulation requires only 4000 ''particles.''

Spectral parameters for scattering amplitudes in N=4 super Yang-Mills theory

International Nuclear Information System (INIS)

Ferro, Livia; Łukowski, Tomasz; Meneghelli, Carlo; Plefka, Jan; Staudacher, Matthias

2014-01-01

Planar N=4 Super Yang-Mills theory appears to be a quantum integrable four-dimensional conformal theory. This has been used to find equations believed to describe its exact spectrum of anomalous dimensions. Integrability seemingly also extends to the planar space-time scattering amplitudes of the N=4 model, which show strong signs of Yangian invariance. However, in contradistinction to the spectral problem, this has not yet led to equations determining the exact amplitudes. We propose that the missing element is the spectral parameter, ubiquitous in integrable models. We show that it may indeed be included into recent on-shell approaches to scattering amplitude integrands, providing a natural deformation of the latter. Under some constraints, Yangian symmetry is preserved. Finally we speculate that the spectral parameter might also be the regulator of choice for controlling the infrared divergences appearing when integrating the integrands in exactly four dimensions

Monograph - The Numerical Integration of Ordinary Differential Equations.

Science.gov (United States)

Hull, T. E.

The materials presented in this monograph are intended to be included in a course on ordinary differential equations at the upper division level in a college mathematics program. These materials provide an introduction to the numerical integration of ordinary differential equations, and they can be used to supplement a regular text on this…

Abecedarian School on Symmetries and Integrability of Difference Equations (ASIDE) & SIDE 12 International Conference Symmetries and Integrability of Difference Equations

CERN Document Server

Rebelo, Raphaël; Winternitz, Pavel

2017-01-01

This book shows how Lie group and integrability techniques, originally developed for differential equations, have been adapted to the case of difference equations. Difference equations are playing an increasingly important role in the natural sciences. Indeed, many phenomena are inherently discrete and thus naturally described by difference equations. More fundamentally, in subatomic physics, space-time may actually be discrete. Differential equations would then just be approximations of more basic discrete ones. Moreover, when using differential equations to analyze continuous processes, it is often necessary to resort to numerical methods. This always involves a discretization of the differential equations involved, thus replacing them by difference ones. Each of the nine peer-reviewed chapters in this volume serves as a self-contained treatment of a topic, containing introductory material as well as the latest research results and exercises. Each chapter is presented by one or more early career researchers...

On monotonic solutions of an integral equation of Abel type

International Nuclear Information System (INIS)

Darwish, Mohamed Abdalla

2007-08-01

We present an existence theorem of monotonic solutions for a quadratic integral equation of Abel type in C[0, 1]. The famous Chandrasekhar's integral equation is considered as a special case. The concept of measure of noncompactness and a fi xed point theorem due to Darbo are the main tools in carrying out our proof. (author)

Alternative integral equations and perturbation expansions for self-coupled scalar fields

International Nuclear Information System (INIS)

Ford, L.H.

1985-01-01

It is shown that the theory of a self-coupled scalar field may be expressed in terms of a class of integral equations which include the Yang-Feldman equation as a particular case. Other integral equations in this class could be used to generate alternative perturbation expansions which contain a nonanalytic dependence upon the coupling constant and are less ultraviolet divergent than the conventional perturbation expansion. (orig.)

Geometrical-integrability constraints and equations of motion in four plus extended super spaces

International Nuclear Information System (INIS)

Chau, L.L.

1987-01-01

It is pointed out that many equations of motion in physics, including gravitational and Yang-Mills equations, have a common origin: i.e. they are the results of certain geometrical integrability conditions. These integrability conditions lead to linear systems and conservation laws that are important in integrating these equations of motion

Three-body scattering problem in the fixed center approximation: The case of attraction

Energy Technology Data Exchange (ETDEWEB)

Kudryavtsev, Alexander E. [National Research Center Kurchatov Institute, Institute for Theoretical and Experimental Physics, Moscow (Russian Federation); Gani, Vakhid A. [National Research Center Kurchatov Institute, Institute for Theoretical and Experimental Physics, Moscow (Russian Federation); National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), Moscow (Russian Federation); Romanov, Alexander I. [National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), Moscow (Russian Federation)

2016-12-15

We study the scattering of a light particle on a bound pair of heavy particles (e.g., the deuteron) within the fixed center approximation in the case of light-heavy attraction, solving the integral equation for the three-body Green's function both in the coordinate and in the momentum space. The results for the three-body scattering amplitude appear to be ambiguous -they depend on a single real parameter. This parameter may be fixed by a three-body input, e.g., the three-body scattering length. We also solve the integral equation for the three-body Green function in the momentum space, introducing a finite cut-off. We show that all three approaches are equivalent. We also discuss how our approach to the problem matches with the introduction of three-body contact interaction as done by other authors. (orig.)

Differential equations for loop integrals in Baikov representation

Science.gov (United States)

Bosma, Jorrit; Larsen, Kasper J.; Zhang, Yang

2018-05-01

We present a proof that differential equations for Feynman loop integrals can always be derived in Baikov representation without involving dimension-shift identities. We moreover show that in a large class of two- and three-loop diagrams it is possible to avoid squared propagators in the intermediate steps of setting up the differential equations.

Simplified expressions of the T-matrix integrals for electromagnetic scattering.

Science.gov (United States)

Somerville, Walter R C; Auguié, Baptiste; Le Ru, Eric C

2011-09-01

The extended boundary condition method, also called the null-field method, provides a semianalytic solution to the problem of electromagnetic scattering by a particle by constructing a transition matrix (T-matrix) that links the scattered field to the incident field. This approach requires the computation of specific integrals over the particle surface, which are typically evaluated numerically. We introduce here a new set of simplified expressions for these integrals in the commonly studied case of axisymmetric particles. Simplifications are obtained using the differentiation properties of the radial functions (spherical Bessel) and angular functions (associated Legendre functions) and integrations by parts. The resulting simplified expressions not only lead to faster computations, but also reduce the risks of loss of precision and provide a simpler framework for further analytical work.

Integral equations for free-molecule ow in MEMS: recent advancements

Directory of Open Access Journals (Sweden)

Fedeli Patrick

2017-03-01

Full Text Available We address a Boundary Integral Equation (BIE approach for the analysis of gas dissipation in near-vacuum for Micro Electro Mechanical Systems (MEMS. Inspired by an analogy with the radiosity equation in computer graphics, we discuss an efficient way to compute the visible domain of integration. Moreover, we tackle the issue of near singular integrals by developing a set of analytical formulas for planar polyhedral domains. Finally a validation with experimental results taken from the literature is presented.

Nonlinear integral equations for the sausage model

Science.gov (United States)

Ahn, Changrim; Balog, Janos; Ravanini, Francesco

2017-08-01

The sausage model, first proposed by Fateev, Onofri, and Zamolodchikov, is a deformation of the O(3) sigma model preserving integrability. The target space is deformed from the sphere to ‘sausage’ shape by a deformation parameter ν. This model is defined by a factorizable S-matrix which is obtained by deforming that of the O(3) sigma model by a parameter λ. Clues for the deformed sigma model are provided by various UV and IR information through the thermodynamic Bethe ansatz (TBA) analysis based on the S-matrix. Application of TBA to the sausage model is, however, limited to the case of 1/λ integer where the coupled integral equations can be truncated to a finite number. In this paper, we propose a finite set of nonlinear integral equations (NLIEs), which are applicable to generic value of λ. Our derivation is based on T-Q relations extracted from the truncated TBA equations. For a consistency check, we compute next-leading order corrections of the vacuum energy and extract the S-matrix information in the IR limit. We also solved the NLIE both analytically and numerically in the UV limit to get the effective central charge and compared with that of the zero-mode dynamics to obtain exact relation between ν and λ. Dedicated to the memory of Petr Petrovich Kulish.

A combined finite element-boundary integral formulation for solution of two-dimensional scattering problems via CGFFT. [Conjugate Gradient Fast Fourier Transformation

Science.gov (United States)

Collins, Jeffery D.; Volakis, John L.; Jin, Jian-Ming

1990-01-01

A new technique is presented for computing the scattering by 2-D structures of arbitrary composition. The proposed solution approach combines the usual finite element method with the boundary-integral equation to formulate a discrete system. This is subsequently solved via the conjugate gradient (CG) algorithm. A particular characteristic of the method is the use of rectangular boundaries to enclose the scatterer. Several of the resulting boundary integrals are therefore convolutions and may be evaluated via the fast Fourier transform (FFT) in the implementation of the CG algorithm. The solution approach offers the principal advantage of having O(N) memory demand and employs a 1-D FFT versus a 2-D FFT as required with a traditional implementation of the CGFFT algorithm. The speed of the proposed solution method is compared with that of the traditional CGFFT algorithm, and results for rectangular bodies are given and shown to be in excellent agreement with the moment method.

Feynman path integral related to stochastic schroedinger equation

International Nuclear Information System (INIS)

Belavkin, V.P.; Smolyanov, O.G.

1998-01-01

The derivation of the Schroedinger equation describing the continuous measurement process is presented. The representation of the solution of the stochastic Schroedinger equation for continuous measurements is obtained by means of the Feynman path integral. The connection with the heuristic approach to the description of continuous measurements is considered. The connection with the Senon paradox is established [ru

On the structure of the commutative Z2 graded algebra valued integrable equations

International Nuclear Information System (INIS)

Konopelchenko, B.G.

1980-01-01

Partial differential equations integrable by the linear matrix spectral problem of arbitrary order are considered for the case that the 'potentials' take their values in the commutative infinte-dimensional Z 2 graded algebra (superalgebra). The general form of the integrable equations and their Baecklund transformations are found. The infinite sets of the integrals of the motion are constructed. The hamiltonian character of the integrable equations is proved. (orig.)

Time Domain Surface Integral Equation Solvers for Quantum Corrected Electromagnetic Analysis of Plasmonic Nanostructures

KAUST Repository

Uysal, Ismail Enes

2016-10-01

Plasmonic structures are utilized in many applications ranging from bio-medicine to solar energy generation and transfer. Numerical schemes capable of solving equations of classical electrodynamics have been the method of choice for characterizing scattering properties of such structures. However, as dimensions of these plasmonic structures reduce to nanometer scale, quantum mechanical effects start to appear. These effects cannot be accurately modeled by available classical numerical methods. One of these quantum effects is the tunneling, which is observed when two structures are located within a sub-nanometer distance of each other. At these small distances electrons “jump" from one structure to another and introduce a path for electric current to flow. Classical equations of electrodynamics and the schemes used for solving them do not account for this additional current path. This limitation can be lifted by introducing an auxiliary tunnel with material properties obtained using quantum models and applying a classical solver to the structures connected by this auxiliary tunnel. Early work on this topic focused on quantum models that are generated using a simple one-dimensional wave function to find the tunneling probability and assume a simple Drude model for the permittivity of the tunnel. These tunnel models are then used together with a classical frequency domain solver. In this thesis, a time domain surface integral equation solver for quantum corrected analysis of transient plasmonic interactions is proposed. This solver has several advantages: (i) As opposed to frequency domain solvers, it provides results at a broad band of frequencies with a single simulation. (ii) As opposed to differential equation solvers, it only discretizes surfaces (reducing number of unknowns), enforces the radiation condition implicitly (increasing the accuracy), and allows for time step selection independent of spatial discretization (increasing efficiency). The quantum model

Set-valued and fuzzy stochastic integral equations driven by semimartingales under Osgood condition

Directory of Open Access Journals (Sweden)

Malinowski Marek T.

2015-01-01

Full Text Available We analyze the set-valued stochastic integral equations driven by continuous semimartingales and prove the existence and uniqueness of solutions to such equations in the framework of the hyperspace of nonempty, bounded, convex and closed subsets of the Hilbert space L2 (consisting of square integrable random vectors. The coefficients of the equations are assumed to satisfy the Osgood type condition that is a generalization of the Lipschitz condition. Continuous dependence of solutions with respect to data of the equation is also presented. We consider equations driven by semimartingale Z and equations driven by processes A;M from decomposition of Z, where A is a process of finite variation and M is a local martingale. These equations are not equivalent. Finally, we show that the analysis of the set-valued stochastic integral equations can be extended to a case of fuzzy stochastic integral equations driven by semimartingales under Osgood type condition. To obtain our results we use the set-valued and fuzzy Maruyama type approximations and Bihari’s inequality.

Alternate Solution to Generalized Bernoulli Equations via an Integrating Factor: An Exact Differential Equation Approach

Science.gov (United States)

Tisdell, C. C.

2017-01-01

Solution methods to exact differential equations via integrating factors have a rich history dating back to Euler (1740) and the ideas enjoy applications to thermodynamics and electromagnetism. Recently, Azevedo and Valentino presented an analysis of the generalized Bernoulli equation, constructing a general solution by linearizing the problem…

ICM: an Integrated Compartment Method for numerically solving partial differential equations

Energy Technology Data Exchange (ETDEWEB)

Yeh, G.T.

1981-05-01

An integrated compartment method (ICM) is proposed to construct a set of algebraic equations from a system of partial differential equations. The ICM combines the utility of integral formulation of finite element approach, the simplicity of interpolation of finite difference approximation, and the flexibility of compartment analyses. The integral formulation eases the treatment of boundary conditions, in particular, the Neumann-type boundary conditions. The simplicity of interpolation provides great economy in computation. The flexibility of discretization with irregular compartments of various shapes and sizes offers advantages in resolving complex boundaries enclosing compound regions of interest. The basic procedures of ICM are first to discretize the region of interest into compartments, then to apply three integral theorems of vectors to transform the volume integral to the surface integral, and finally to use interpolation to relate the interfacial values in terms of compartment values to close the system. The Navier-Stokes equations are used as an example of how to derive the corresponding ICM alogrithm for a given set of partial differential equations. Because of the structure of the algorithm, the basic computer program remains the same for cases in one-, two-, or three-dimensional problems.

Multi-symplectic variational integrators for nonlinear Schrödinger equations with variable coefficients

International Nuclear Information System (INIS)

Liao Cui-Cui; Cui Jin-Chao; Liang Jiu-Zhen; Ding Xiao-Hua

2016-01-01

In this paper, we propose a variational integrator for nonlinear Schrödinger equations with variable coefficients. It is shown that our variational integrator is naturally multi-symplectic. The discrete multi-symplectic structure of the integrator is presented by a multi-symplectic form formula that can be derived from the discrete Lagrangian boundary function. As two examples of nonlinear Schrödinger equations with variable coefficients, cubic nonlinear Schrödinger equations and Gross–Pitaevskii equations are extensively studied by the proposed integrator. Our numerical simulations demonstrate that the integrator is capable of preserving the mass, momentum, and energy conservation during time evolutions. Convergence tests are presented to verify that our integrator has second-order accuracy both in time and space. (paper)

Localization of the eigenvalues of linear integral equations with applications to linear ordinary differential equations.

Science.gov (United States)

Sloss, J. M.; Kranzler, S. K.

1972-01-01

The equivalence of a considered integral equation form with an infinite system of linear equations is proved, and the localization of the eigenvalues of the infinite system is expressed. Error estimates are derived, and the problems of finding upper bounds and lower bounds for the eigenvalues are solved simultaneously.

General concepts of multichannel collision theory and their translation into the matrix formulation of few-body integral equations

International Nuclear Information System (INIS)

Sandhas, W.

1978-01-01

In the N-body problem mappings between channel states and scattering states are studied. It is shown in particular that the (2sup(N-1)-1) two-fragment MOELLER operators introduced on the whole Hilbert space are sufficient to provide all multi-fragment scattering states. Hence, each of these states is uniquely determined by (2sup(N-1)-1) Lippmann-Schwinger (LS) equations. Rewriting every set of LS equations as one matrix equation, current few-body approaches are incorporated in a rather natural way. The typical uniqueness questions of such coupled systems are discussed, and it is shown that Faddeev-type couplings lead to unique equations for arbitrary N. (author)

General concepts of multichannel collision theory and their translation into the matrix formulation of few-body integral equations

International Nuclear Information System (INIS)

Sandhas, W.

1978-04-01

In the N-body problem mappings between channel states and scattering states are studied. It is shown in particular that the (2sup(N-1)-1) two-fragment Moeller operators introduced on the whole Hilbert space are sufficient to provide all multifragment scattering states. Hence, each of these states is uniquly determined by (2sup(N-1)-1) Lippmann-Schwinger (LS) equations. Rewriting every set of LS equations as one matrix equation, current few-body approaches are incorporated in a rather natural way. The typical uniqueness questions of such coupled systems are discussed, and it si shown that Faddeev-type couplings lead to unique equations for arbitrary N. (orig.) [de

Integral equation approach to time-dependent kinematic dynamos in finite domains

International Nuclear Information System (INIS)

Xu Mingtian; Stefani, Frank; Gerbeth, Gunter

2004-01-01

The homogeneous dynamo effect is at the root of cosmic magnetic field generation. With only a very few exceptions, the numerical treatment of homogeneous dynamos is carried out in the framework of the differential equation approach. The present paper tries to facilitate the use of integral equations in dynamo research. Apart from the pedagogical value to illustrate dynamo action within the well-known picture of the Biot-Savart law, the integral equation approach has a number of practical advantages. The first advantage is its proven numerical robustness and stability. The second and perhaps most important advantage is its applicability to dynamos in arbitrary geometries. The third advantage is its intimate connection to inverse problems relevant not only for dynamos but also for technical applications of magnetohydrodynamics. The paper provides the first general formulation and application of the integral equation approach to time-dependent kinematic dynamos, with stationary dynamo sources, in finite domains. The time dependence is restricted to the magnetic field, whereas the velocity or corresponding mean-field sources of dynamo action are supposed to be stationary. For the spherically symmetric α 2 dynamo model it is shown how the general formulation is reduced to a coupled system of two radial integral equations for the defining scalars of the poloidal and toroidal field components. The integral equation formulation for spherical dynamos with general stationary velocity fields is also derived. Two numerical examples - the α 2 dynamo model with radially varying α and the Bullard-Gellman model - illustrate the equivalence of the approach with the usual differential equation method. The main advantage of the method is exemplified by the treatment of an α 2 dynamo in rectangular domains

Green function of the double-fractional Fokker-Planck equation: Path integral and stochastic differential equations

Science.gov (United States)

Kleinert, H.; Zatloukal, V.

2013-11-01

The statistics of rare events, the so-called black-swan events, is governed by non-Gaussian distributions with heavy power-like tails. We calculate the Green functions of the associated Fokker-Planck equations and solve the related stochastic differential equations. We also discuss the subject in the framework of path integration.

Integrability and structural stability of solutions to the Ginzburg-Landau equation

Science.gov (United States)

Keefe, Laurence R.

1986-01-01

The integrability of the Ginzburg-Landau equation is studied to investigate if the existence of chaotic solutions found numerically could have been predicted a priori. The equation is shown not to possess the Painleveproperty, except for a special case of the coefficients that corresponds to the integrable, nonlinear Schroedinger (NLS) equation. Regarding the Ginzburg-Landau equation as a dissipative perturbation of the NLS, numerical experiments show all but one of a family of two-tori solutions, possessed by the NLS under particular conditions, to disappear under real perturbations to the NLS coefficients of O(10 to the -6th).

Integrable discretizations of the (2+1)-dimensional sinh-Gordon equation

International Nuclear Information System (INIS)

Hu, Xing-Biao; Yu, Guo-Fu

2007-01-01

In this paper, we propose two semi-discrete equations and one fully discrete equation and study them by Hirota's bilinear method. These equations have continuum limits into a system which admits the (2+1)-dimensional generalization of the sinh-Gordon equation. As a result, two integrable semi-discrete versions and one fully discrete version for the sinh-Gordon equation are found. Baecklund transformations, nonlinear superposition formulae, determinant solution and Lax pairs for these discrete versions are presented

A New Algorithm for System of Integral Equations

Directory of Open Access Journals (Sweden)

Abdujabar Rasulov

2014-01-01

Full Text Available We develop a new algorithm to solve the system of integral equations. In this new method no need to use matrix weights. Beacause of it, we reduce computational complexity considerable. Using the new algorithm it is also possible to solve an initial boundary value problem for system of parabolic equations. To verify the efficiency, the results of computational experiments are given.

On integration of the first order differential equations in a finite terms

International Nuclear Information System (INIS)

Malykh, M D

2017-01-01

There are several approaches to the description of the concept called briefly as integration of the first order differential equations in a finite terms or symbolical integration. In the report three of them are considered: 1.) finding of a rational integral (Beaune or Poincaré problem), 2.) integration by quadratures and 3.) integration when the general solution of given differential equation is an algebraical function of a constant (Painlevé problem). Their realizations in Sage are presented. (paper)

Polynomial solutions of nonlinear integral equations

International Nuclear Information System (INIS)

Dominici, Diego

2009-01-01

We analyze the polynomial solutions of a nonlinear integral equation, generalizing the work of Bender and Ben-Naim (2007 J. Phys. A: Math. Theor. 40 F9, 2008 J. Nonlinear Math. Phys. 15 (Suppl. 3) 73). We show that, in some cases, an orthogonal solution exists and we give its general form in terms of kernel polynomials

Polynomial solutions of nonlinear integral equations

Energy Technology Data Exchange (ETDEWEB)

Dominici, Diego [Department of Mathematics, State University of New York at New Paltz, 1 Hawk Dr. Suite 9, New Paltz, NY 12561-2443 (United States)], E-mail: dominicd@newpaltz.edu

2009-05-22

We analyze the polynomial solutions of a nonlinear integral equation, generalizing the work of Bender and Ben-Naim (2007 J. Phys. A: Math. Theor. 40 F9, 2008 J. Nonlinear Math. Phys. 15 (Suppl. 3) 73). We show that, in some cases, an orthogonal solution exists and we give its general form in terms of kernel polynomials.

Thermodynamically self-consistent integral equations and the structure of liquid metals

International Nuclear Information System (INIS)

Pastore, G.; Kahl, G.

1987-01-01

We discuss the application of the new thermodynamically self-consistent integral equations for the determination of the structural properties of liquid metals. We present a detailed comparison of the structure (S(q) and g(r)) for models of liquid alkali metals as obtained from two thermodynamically self-consistent integral equations and some published exact computer simulation results; the range of states extends from the triple point to the expanded metal. The theories which only impose thermodynamic self-consistency without any fitting of external data show an excellent agreement with the simulation results, thus demonstrating that this new type of integral equation is definitely superior to the conventional ones (hypernetted chain, Percus-Yevick, mean spherical approximation, etc). (author)

THE OPTIMIZATION OF ELECTRODYNAMIC CONFIGURATION OBJECT WITH THE DESIRED CHARACTERISTICS OF SCATTERING.

Directory of Open Access Journals (Sweden)

A. P. Preobrazhensky

2017-02-01

Full Text Available This paper considers the problem of optimization of the characteristics of scattering of electromagnetic waves on periodic electrodynamic structure. The solution of the scattering problem is based on the method of integral equations, the optimization of the characteristics is based on the genetic algorithm. Recommendations on the parameters of the periodic structure under given angles are given.

Distribution theory for Schrödinger’s integral equation

NARCIS (Netherlands)

Lange, R.J.

2015-01-01

Much of the literature on point interactions in quantum mechanics has focused on the differential form of Schrödinger's equation. This paper, in contrast, investigates the integral form of Schrödinger's equation. While both forms are known to be equivalent for smooth potentials, this is not true for

Unconditionally stable integration of Maxwell’s equations

NARCIS (Netherlands)

Verwer, J.G.; Botchev, M.A.

2009-01-01

Numerical integration of Maxwell’s equations is often based on explicit methods accepting a stability step size restriction. In literature evidence is given that there is also a need for unconditionally stable methods, as exemplified by the successful alternating direction implicit - finite

Acoustic scattering by arbitrary distributions of disjoint, homogeneous cylinders or spheres.

Science.gov (United States)

Hesford, Andrew J; Astheimer, Jeffrey P; Waag, Robert C

2010-05-01

A T-matrix formulation is presented to compute acoustic scattering from arbitrary, disjoint distributions of cylinders or spheres, each with arbitrary, uniform acoustic properties. The generalized approach exploits the similarities in these scattering problems to present a single system of equations that is easily specialized to cylindrical or spherical scatterers. By employing field expansions based on orthogonal harmonic functions, continuity of pressure and normal particle velocity are directly enforced at each scatterer using diagonal, analytic expressions to eliminate the need for integral equations. The effect of a cylinder or sphere that encloses all other scatterers is simulated with an outer iterative procedure that decouples the inner-object solution from the effect of the enclosing object to improve computational efficiency when interactions among the interior objects are significant. Numerical results establish the validity and efficiency of the outer iteration procedure for nested objects. Two- and three-dimensional methods that employ this outer iteration are used to measure and characterize the accuracy of two-dimensional approximations to three-dimensional scattering of elevation-focused beams.

Quantum method of the inverse scattering problem. Pt. 1

International Nuclear Information System (INIS)

Sklyamin, E.K.; Takhtadzhyan, L.A.; Faddeev, L.D.

1978-12-01

In this work the authors use a formulation for the method of the inverse scattering problem for quantum-mechanical models of the field theory, that can be found in a quantization of these fully integrable systems. As the most important example serves the system (sinγ) 2 with the movement equation: γtt -γxx + m 2 /β sinβγ = 0 that is known under the specification Sine-Gordon-equation. (orig.) [de

Inverse scattering transform for the vector nonlinear Schroedinger equation with nonvanishing boundary conditions

International Nuclear Information System (INIS)

Prinari, Barbara; Ablowitz, Mark J.; Biondini, Gino

2006-01-01

The inverse scattering transform for the vector defocusing nonlinear Schroedinger (NLS) equation with nonvanishing boundary values at infinity is constructed. The direct scattering problem is formulated on a two-sheeted covering of the complex plane. Two out of the six Jost eigenfunctions, however, do not admit an analytic extension on either sheet of the Riemann surface. Therefore, a suitable modification of both the direct and the inverse problem formulations is necessary. On the direct side, this is accomplished by constructing two additional analytic eigenfunctions which are expressed in terms of the adjoint eigenfunctions. The discrete spectrum, bound states and symmetries of the direct problem are then discussed. In the most general situation, a discrete eigenvalue corresponds to a quartet of zeros (poles) of certain scattering data. The inverse scattering problem is formulated in terms of a generalized Riemann-Hilbert (RH) problem in the upper/lower half planes of a suitable uniformization variable. Special soliton solutions are constructed from the poles in the RH problem, and include dark-dark soliton solutions, which have dark solitonic behavior in both components, as well as dark-bright soliton solutions, which have one dark and one bright component. The linear limit is obtained from the RH problem and is shown to correspond to the Fourier transform solution obtained from the linearized vector NLS system

A new integral method for solving the point reactor neutron kinetics equations

International Nuclear Information System (INIS)

Li Haofeng; Chen Wenzhen; Luo Lei; Zhu Qian

2009-01-01

A numerical integral method that efficiently provides the solution of the point kinetics equations by using the better basis function (BBF) for the approximation of the neutron density in one time step integrations is described and investigated. The approach is based on an exact analytic integration of the neutron density equation, where the stiffness of the equations is overcome by the fully implicit formulation. The procedure is tested by using a variety of reactivity functions, including step reactivity insertion, ramp input and oscillatory reactivity changes. The solution of the better basis function method is compared to other analytical and numerical solutions of the point reactor kinetics equations. The results show that selecting a better basis function can improve the efficiency and accuracy of this integral method. The better basis function method can be used in real time forecasting for power reactors in order to prevent reactivity accidents.

Lagrangian structures, integrability and chaos for 3D dynamical equations

International Nuclear Information System (INIS)

Bustamante, Miguel D; Hojman, Sergio A

2003-01-01

In this paper, we consider the general setting for constructing action principles for three-dimensional first-order autonomous equations. We present the results for some integrable and non-integrable cases of the Lotka-Volterra equation, and show Lagrangian descriptions which are valid for systems satisfying Shil'nikov criteria on the existence of strange attractors, though chaotic behaviour has not been verified up to now. The Euler-Lagrange equations we get for these systems usually present 'time reparametrization' invariance, though other kinds of invariance may be found according to the kernel of the associated symplectic 2-form. The formulation of a Hamiltonian structure (Poisson brackets and Hamiltonians) for these systems from the Lagrangian viewpoint leads to a method of finding new constants of the motion starting from known ones, which is applied to some systems found in the literature known to possess a constant of the motion, to find the other and thus showing their integrability. In particular, we show that the so-called ABC system is completely integrable if it possesses one constant of the motion

Anisotropy of relativistic lepton coherent scattering at axial channeling

International Nuclear Information System (INIS)

Telegin, V.I.; Kanloev, A.M.; Kungurov, F.R.

1989-01-01

The contribution of the coherent and incoherent scattering of relativistic leptons passed through thin crystals in the channeling mode to their angular distribution is considered. The investigation was carried out by numerical integration of the motion equations for a great number of particles. It is shown that in the crystals with a thickness smaller than the dechanneling length the determining role in formation of distribution over the axit angles is played by the coherent scattering of particles by atomic chains. The effect of the multiple scattering on the angular distribution is negligibly small. 6 refs.; 4 figs

Universal and integrable nonlinear evolution systems of equations in 2+1 dimensions

International Nuclear Information System (INIS)

Maccari, A.

1997-01-01

Integrable systems of nonlinear partial differential equations (PDEs) are obtained from integrable equations in 2+1 dimensions, by means of a reduction method of broad applicability based on Fourier expansion and spatio endash temporal rescalings, which is asymptotically exact in the limit of weak nonlinearity. The integrability by the spectral transform is explicitly demonstrated, because the corresponding Lax pairs have been derived, applying the same reduction method to the Lax pair of the initial equation. These systems of nonlinear PDEs are likely to be of applicative relevance and have a open-quotes universalclose quotes character, inasmuch as they may be derived from a very large class of nonlinear evolution equations with a linear dispersive part. copyright 1997 American Institute of Physics

The integration of improved Monte Carlo compton scattering algorithms into the Integrated TIGER Series

International Nuclear Information System (INIS)

Quirk, Thomas J. IV

2004-01-01

The Integrated TIGER Series (ITS) is a software package that solves coupled electron-photon transport problems. ITS performs analog photon tracking for energies between 1 keV and 1 GeV. Unlike its deterministic counterpart, the Monte Carlo calculations of ITS do not require a memory-intensive meshing of phase space; however, its solutions carry statistical variations. Reducing these variations is heavily dependent on runtime. Monte Carlo simulations must therefore be both physically accurate and computationally efficient. Compton scattering is the dominant photon interaction above 100 keV and below 5-10 MeV, with higher cutoffs occurring in lighter atoms. In its current model of Compton scattering, ITS corrects the differential Klein-Nishina cross sections (which assumes a stationary, free electron) with the incoherent scattering function, a function dependent on both the momentum transfer and the atomic number of the scattering medium. While this technique accounts for binding effects on the scattering angle, it excludes the Doppler broadening the Compton line undergoes because of the momentum distribution in each bound state. To correct for these effects, Ribbefor's relativistic impulse approximation (IA) will be employed to create scattering cross section differential in both energy and angle for each element. Using the parameterizations suggested by Brusa et al., scattered photon energies and angle can be accurately sampled at a high efficiency with minimal physical data. Two-body kinematics then dictates the electron's scattered direction and energy. Finally, the atomic ionization is relaxed via Auger emission or fluorescence. Future work will extend these improvements in incoherent scattering to compounds and to adjoint calculations.

Self-consistent field theory of collisions: Orbital equations with asymptotic sources and self-averaged potentials

Energy Technology Data Exchange (ETDEWEB)

Hahn, Y.K., E-mail: ykhahn22@verizon.net

2014-12-15

The self-consistent field theory of collisions is formulated, incorporating the unique dynamics generated by the self-averaged potentials. The bound state Hartree–Fock approach is extended for the first time to scattering states, by properly resolving the principal difficulties of non-integrable continuum orbitals and imposing complex asymptotic conditions. The recently developed asymptotic source theory provides the natural theoretical basis, as the asymptotic conditions are completely transferred to the source terms and the new scattering function is made fullyintegrable. The scattering solutions can then be directly expressed in terms of bound state HF configurations, establishing the relationship between the bound and scattering state solutions. Alternatively, the integrable spin orbitals are generated by constructing the individual orbital equations that contain asymptotic sources and self-averaged potentials. However, the orbital energies are not determined by the equations, and a special channel energy fixing procedure is developed to secure the solutions. It is also shown that the variational construction of the orbital equations has intrinsic ambiguities that are generally associated with the self-consistent approach. On the other hand, when a small subset of open channels is included in the source term, the solutions are only partiallyintegrable, but the individual open channels can then be treated more simply by properly selecting the orbital energies. The configuration mixing and channel coupling are then necessary to complete the solution. The new theory improves the earlier continuum HF model. - Highlights: • First extension of HF to scattering states, with proper asymptotic conditions. • Orbital equations with asymptotic sources and integrable orbital solutions. • Construction of self-averaged potentials, and orbital energy fixing. • Channel coupling and configuration mixing, involving the new orbitals. • Critical evaluation of the

Integration Processes of Delay Differential Equation Based on Modified Laguerre Functions

Directory of Open Access Journals (Sweden)

Yeguo Sun

2012-01-01

Full Text Available We propose long-time convergent numerical integration processes for delay differential equations. We first construct an integration process based on modified Laguerre functions. Then we establish its global convergence in certain weighted Sobolev space. The proposed numerical integration processes can also be used for systems of delay differential equations. We also developed a technique for refinement of modified Laguerre-Radau interpolations. Lastly, numerical results demonstrate the spectral accuracy of the proposed method and coincide well with analysis.

An Algorithm for Computing Screened Coulomb Scattering in Geant4

OpenAIRE

Mendenhall, Marcus H.; Weller, Robert A.

2004-01-01

An algorithm has been developed for the Geant4 Monte-Carlo package for the efficient computation of screened Coulomb interatomic scattering. It explicitly integrates the classical equations of motion for scattering events, resulting in precise tracking of both the projectile and the recoil target nucleus. The algorithm permits the user to plug in an arbitrary screening function, such as Lens-Jensen screening, which is good for backscattering calculations, or Ziegler-Biersack-Littmark screenin...

Stability of negative solitary waves for an integrable modified Camassa-Holm equation

International Nuclear Information System (INIS)

Yin Jiuli; Tian Lixin; Fan Xinghua

2010-01-01

In this paper, we prove that the modified Camassa-Holm equation is Painleve integrable. We also study the orbital stability problem of negative solitary waves for this integrable equation. It is shown that the negative solitary waves are stable for arbitrary wave speed of propagation.

On a partial differential equation method for determining the free energies and coexisting phase compositions of ternary mixtures from light scattering data.

Science.gov (United States)

Ross, David S; Thurston, George M; Lutzer, Carl V

2008-08-14

In this paper we present a method for determining the free energies of ternary mixtures from light scattering data. We use an approximation that is appropriate for liquid mixtures, which we formulate as a second-order nonlinear partial differential equation. This partial differential equation (PDE) relates the Hessian of the intensive free energy to the efficiency of light scattering in the forward direction. This basic equation applies in regions of the phase diagram in which the mixtures are thermodynamically stable. In regions in which the mixtures are unstable or metastable, the appropriate PDE is the nonlinear equation for the convex hull. We formulate this equation along with continuity conditions for the transition between the two equations at cloud point loci. We show how to discretize this problem to obtain a finite-difference approximation to it, and we present an iterative method for solving the discretized problem. We present the results of calculations that were done with a computer program that implements our method. These calculations show that our method is capable of reconstructing test free energy functions from simulated light scattering data. If the cloud point loci are known, the method also finds the tie lines and tie triangles that describe thermodynamic equilibrium between two or among three liquid phases. A robust method for solving this PDE problem, such as the one presented here, can be a basis for optical, noninvasive means of characterizing the thermodynamics of multicomponent mixtures.

Fibonacci-regularization method for solving Cauchy integral equations of the first kind

Directory of Open Access Journals (Sweden)

Mohammad Ali Fariborzi Araghi

2017-09-01

Full Text Available In this paper, a novel scheme is proposed to solve the first kind Cauchy integral equation over a finite interval. For this purpose, the regularization method is considered. Then, the collocation method with Fibonacci base function is applied to solve the obtained second kind singular integral equation. Also, the error estimate of the proposed scheme is discussed. Finally, some sample Cauchy integral equations stem from the theory of airfoils in fluid mechanics are presented and solved to illustrate the importance and applicability of the given algorithm. The tables in the examples show the efficiency of the method.

Solving very large scattering problems using a parallel PWTD-enhanced surface integral equation solver

KAUST Repository

Liu, Yang; Bagci, Hakan; Michielssen, Eric

2013-01-01

numbers of temporal and spatial basis functions discretizing the current [Shanker et al., IEEE Trans. Antennas Propag., 51, 628-641, 2003]. In the past, serial versions of these solvers have been successfully applied to the analysis of scattering from

Faddeev and Glauber calculations at intermediate energies in a model for n+d scattering

International Nuclear Information System (INIS)

Elster, Ch.; Lin, T.; Gloeckle, W.; Jeschonnek, S.

2008-01-01

Obtaining cross sections for nuclear reactions at intermediate energies based on the Glauber formulation has a long tradition. Only recently the energy regime of a few hundred MeV has become accessible to ab initio Faddeev calculations of three-body scattering. In order to go to higher energies, the Faddeev equation for three-body scattering is formulated and directly solved without employing a partial wave decomposition. In the simplest form the Faddeev equation for interacting scalar particles is a three-dimensional integral equation in five variables, from which the total cross section, the cross sections for elastic scattering and breakup reactions, as well as differential cross sections are obtained. The same observables are calculated based on the Glauber formulation. The first order Glauber calculation and the Glauber rescattering corrections are compared in detail with the corresponding terms of the Faddeev multiple scattering series for projectile energies between 100 MeV and 2 GeV

Field Method for Integrating the First Order Differential Equation

Institute of Scientific and Technical Information of China (English)

JIA Li-qun; ZHENG Shi-wang; ZHANG Yao-yu

2007-01-01

An important modern method in analytical mechanics for finding the integral, which is called the field-method, is used to research the solution of a differential equation of the first order. First, by introducing an intermediate variable, a more complicated differential equation of the first order can be expressed by two simple differential equations of the first order, then the field-method in analytical mechanics is introduced for solving the two differential equations of the first order. The conclusion shows that the field-method in analytical mechanics can be fully used to find the solutions of a differential equation of the first order, thus a new method for finding the solutions of the first order is provided.

On an integrable deformed affinsphären equation. A reciprocal gasdynamic connection

International Nuclear Information System (INIS)

Rogers, C.; Huang, Yehui

2012-01-01

The integrable affinsphären equation originally arose in a geometric context but has an interesting gasdynamic connection. Here, an integrable deformed version of the affinsphären equation is derived in a novel manner via the action of reciprocal transformations on a related anisentropic gasdynamics system. A linear representation for the deformed affinsphären equation is constructed by means of the reciprocal transformations. The latter are then employed to derive a class of exact solutions in parametric form. -- Highlights: ► A deformed affinsphären equation is derived via a reciprocal transformation. ► A linear representation for the deformed affinsphären equation is constructed. ► A class of exact solutions of the deformed affinsphären equation is presented.

Analysis and regularization of the thin-wire integral equation with reduced kernel

NARCIS (Netherlands)

Beurden, van M.C.; Tijhuis, A.G.

2007-01-01

For the straight wire, modeled as a hollow tube, we establish a conditional equivalence relation between the integral equations with exact and reduced kernel. This relation allows us to examine the existence and uniqueness conditions for the integral equation with reduced kernel, based on a local

Electromagnetic wave scattering by many small particles

International Nuclear Information System (INIS)

Ramm, A.G.

2007-01-01

Scattering of electromagnetic waves by many small particles of arbitrary shapes is reduced rigorously to solving linear algebraic system of equations bypassing the usual usage of integral equations. The matrix elements of this linear algebraic system have physical meaning. They are expressed in terms of the electric and magnetic polarizability tensors. Analytical formulas are given for calculation of these tensors with any desired accuracy for homogeneous bodies of arbitrary shapes. An idea to create a 'smart' material by embedding many small particles in a given region is formulated

Derivation of the Schrodinger Equation from the Hamilton-Jacobi Equation in Feynman's Path Integral Formulation of Quantum Mechanics

Science.gov (United States)

Field, J. H.

2011-01-01

It is shown how the time-dependent Schrodinger equation may be simply derived from the dynamical postulate of Feynman's path integral formulation of quantum mechanics and the Hamilton-Jacobi equation of classical mechanics. Schrodinger's own published derivations of quantum wave equations, the first of which was also based on the Hamilton-Jacobi…

Integrable spin chains and scattering amplitudes

Energy Technology Data Exchange (ETDEWEB)

Bartels, J.; Prygarin, A. [Hamburg Univ. (Germany). 2. Inst. fuer Theoretische Physik; Lipatov, L.N. [Hamburg Univ. (Germany). 2. Inst. fuer Theoretische Physik; Petersburg Nuclear Physics Institute (Russian Federation); Sankt-Peterburgskij Univ., St. Petersburg (Russian Federation)

2011-04-15

In this review we show that the multi-particle scattering amplitudes in N=4 SYM at large N{sub c} and in the multi-Regge kinematics for some physical regions have the high energy behavior appearing from the contribution of the Mandelstam cuts in the complex angular momentum plane of the corresponding t-channel partial waves. These Mandelstam cuts or Regge cuts are resulting from gluon composite states in the adjoint representation of the gauge group SU(N{sub c}). In the leading logarithmic approximation (LLA) their contribution to the six point amplitude is in full agreement with the known two-loop result. The Hamiltonian for the Mandelstam states constructed from n gluons in LLA coincides with the local Hamiltonian of an integrable open spin chain. We construct the corresponding wave functions using the integrals of motion and the Baxter-Sklyanin approach. (orig.)

THE DECISION OF FORM FOR DIFFRACTIVE STRUCTURES IN THE PROBLEM OF SCATTERING OF RADIO WAVES.

Directory of Open Access Journals (Sweden)

A. P. Preobrazhensky

2017-02-01

Full Text Available This paper considers the problem of scattering of electromagnetic waves in different diffraction structures. The solution of the scattering problem is based on the method of integral equations. On diagrams of backscattering at various frequencies of the incident wave, the decision about the form of the object is carried out.

Numerical method for solving linear Fredholm fuzzy integral equations of the second kind

Energy Technology Data Exchange (ETDEWEB)

Abbasbandy, S. [Department of Mathematics, Imam Khomeini International University, P.O. Box 288, Ghazvin 34194 (Iran, Islamic Republic of)]. E-mail: saeid@abbasbandy.com; Babolian, E. [Faculty of Mathematical Sciences and Computer Engineering, Teacher Training University, Tehran 15618 (Iran, Islamic Republic of); Alavi, M. [Department of Mathematics, Arak Branch, Islamic Azad University, Arak 38135 (Iran, Islamic Republic of)

2007-01-15

In this paper we use parametric form of fuzzy number and convert a linear fuzzy Fredholm integral equation to two linear system of integral equation of the second kind in crisp case. We can use one of the numerical method such as Nystrom and find the approximation solution of the system and hence obtain an approximation for fuzzy solution of the linear fuzzy Fredholm integral equations of the second kind. The proposed method is illustrated by solving some numerical examples.

A new analysis of π K scattering from Roy and Steiner type equations

International Nuclear Information System (INIS)

Buettiker, P.; Descotes-Genon, S.; Moussallam, B.

2004-01-01

With the aim of generating new constraints on the OZI suppressed couplings of chiral perturbation theory a set of six equations of the Roy and Steiner type for the S- and P-waves of the πK scattering amplitudes is derived. The range of validity and the multiplicity of the solutions are discussed. Precise numerical solutions are obtained in the range E or sim 1 GeV for both πK→πK and ππ→ KK amplitudes. Our main result is the determination of a narrow allowed region for the two S-wave scattering lengths. Present experimental data below 1 GeV are found to be in generally poor agreement with our results. A set of threshold expansion parameters, as well as sub-threshold parameters are computed. For the latter, a matching with the SU(3) chiral expansion at NLO is performed. (orig.)

An integral for second-order multiple scattering perturbation theory

International Nuclear Information System (INIS)

Hoffman, G.G.

1997-01-01

This paper presents the closed form evaluation of a six-dimensional integral. The integral arises in the application to many-electron systems of a multiple scattering perturbation expansion at second order when formulated in fourier space. The resulting function can be used for the calculation of both the electron density and the effective one-electron potential in an SCF calculations. The closed form expression derived here greatly facilitates these calculations. In addition, the evaluated integral can be used for the computation of second-order corrections to the open-quotes optimized Thomas-Fermi theory.close quotes 10 refs., 2 figs

Numerical Simulation of Antennas with Improved Integral Equation Method

International Nuclear Information System (INIS)

Ma Ji; Fang Guang-You; Lu Wei

2015-01-01

Simulating antennas around a conducting object is a challenge task in computational electromagnetism, which is concerned with the behaviour of electromagnetic fields. To analyze this model efficiently, an improved integral equation-fast Fourier transform (IE-FFT) algorithm is presented in this paper. The proposed scheme employs two Cartesian grids with different size and location to enclose the antenna and the other object, respectively. On the one hand, IE-FFT technique is used to store matrix in a sparse form and accelerate the matrix-vector multiplication for each sub-domain independently. On the other hand, the mutual interaction between sub-domains is taken as the additional exciting voltage in each matrix equation. By updating integral equations several times, the whole electromagnetic system can achieve a stable status. Finally, the validity of the presented method is verified through the analysis of typical antennas in the presence of a conducting object. (paper)

Semiclassical series solution of the generalized phase shift atom--diatom scattering equations

International Nuclear Information System (INIS)

Squire, K.R.; Curtiss, C.F.

1980-01-01

A semiclassical series solution of the previously developed operator form of the generalized phase shift equations describing atom--diatom scattering is presented. This development is based on earlier work which led to a double series in powers of Planck's constant and a scaling parameter of the anisotropic portion of the intermolecular potential. The present solution is similar in that it is a double power series in Planck's constant and in the difference between the spherical radial momentum and a first order approximation. The present series solution avoids difficulties of the previous series associated with the classical turning point

Dhage Iteration Method for Generalized Quadratic Functional Integral Equations

Directory of Open Access Journals (Sweden)

Bapurao C. Dhage

2015-01-01

Full Text Available In this paper we prove the existence as well as approximations of the solutions for a certain nonlinear generalized quadratic functional integral equation. An algorithm for the solutions is developed and it is shown that the sequence of successive approximations starting at a lower or upper solution converges monotonically to the solutions of related quadratic functional integral equation under some suitable mixed hybrid conditions. We rely our main result on Dhage iteration method embodied in a recent hybrid fixed point theorem of Dhage (2014 in partially ordered normed linear spaces. An example is also provided to illustrate the abstract theory developed in the paper.

A momentum-space formulation without partial wave decomposition for scattering of two spin-half particles

Energy Technology Data Exchange (ETDEWEB)

Fachruddin, Imam, E-mail: imam.fachruddin@sci.ui.ac.id; Salam, Agus [Departemen Fisika, Universitas Indonesia, Depok 16424 (Indonesia)

2016-03-11

A new momentum-space formulation for scattering of two spin-half particles, both either identical or unidentical, is formulated. As basis states the free linear-momentum states are not expanded into the angular-momentum states, the system’s spin states are described by the product of the spin states of the two particles, and the system’s isospin states by the total isospin states of the two particles. We evaluate the Lippmann-Schwinger equations for the T-matrix elements in these basis states. The azimuthal behavior of the potential and of the T-matrix elements leads to a set of coupled integral equations for the T-matrix elements in two variables only, which are the magnitude of the relative momentum and the scattering angle. Some symmetry relations for the potential and the T-matrix elements reduce the number of the integral equations to be solved. A set of six spin operators to express any interaction of two spin-half particles is introduced. We show the spin-averaged differential cross section as being calculated in terms of the solution of the set of the integral equations.

TBA-like integral equations from quantized mirror curves

Energy Technology Data Exchange (ETDEWEB)

Okuyama, Kazumi [Department of Physics, Shinshu University,Matsumoto 390-8621 (Japan); Zakany, Szabolcs [Département de Physique Théorique, Université de Genève,Genève, CH-1211 (Switzerland)

2016-03-15

Quantizing the mirror curve of certain toric Calabi-Yau (CY) three-folds leads to a family of trace class operators. The resolvent function of these operators is known to encode topological data of the CY. In this paper, we show that in certain cases, this resolvent function satisfies a system of non-linear integral equations whose structure is very similar to the Thermodynamic Bethe Ansatz (TBA) systems. This can be used to compute spectral traces, both exactly and as a semiclassical expansion. As a main example, we consider the system related to the quantized mirror curve of local ℙ{sup 2}. According to a recent proposal, the traces of this operator are determined by the refined BPS indices of the underlying CY. We use our non-linear integral equations to test that proposal.

TBA-like integral equations from quantized mirror curves

Science.gov (United States)

Okuyama, Kazumi; Zakany, Szabolcs

2016-03-01

Quantizing the mirror curve of certain toric Calabi-Yau (CY) three-folds leads to a family of trace class operators. The resolvent function of these operators is known to encode topological data of the CY. In this paper, we show that in certain cases, this resolvent function satisfies a system of non-linear integral equations whose structure is very similar to the Thermodynamic Bethe Ansatz (TBA) systems. This can be used to compute spectral traces, both exactly and as a semiclassical expansion. As a main example, we consider the system related to the quantized mirror curve of local P2. According to a recent proposal, the traces of this operator are determined by the refined BPS indices of the underlying CY. We use our non-linear integral equations to test that proposal.

Solving Abel’s Type Integral Equation with Mikusinski’s Operator of Fractional Order

Directory of Open Access Journals (Sweden)

Ming Li

2013-01-01

Full Text Available This paper gives a novel explanation of the integral equation of Abel’s type from the point of view of Mikusinski’s operational calculus. The concept of the inverse of Mikusinski’s operator of fractional order is introduced for constructing a representation of the solution to the integral equation of Abel’s type. The proof of the existence of the inverse of the fractional Mikusinski operator is presented, providing an alternative method of treating the integral equation of Abel’s type.

Extended resolvent and inverse scattering with an application to KPI

International Nuclear Information System (INIS)

Boiti, M.; Pempinelli, F.; Pogrebkov, A.K.; Prinari, B.

2003-01-01

We present in detail an extended resolvent approach for investigating linear problems associated to 2+1 dimensional integrable equations. Our presentation is based as an example on the nonstationary Schroedinger equation with potential being a perturbation of the one-soliton potential by means of a decaying two-dimensional function. Modification of the inverse scattering theory as well as properties of the Jost solutions and spectral data as follows from the resolvent approach are given

Solving differential equations for Feynman integrals by expansions near singular points

Science.gov (United States)

Lee, Roman N.; Smirnov, Alexander V.; Smirnov, Vladimir A.

2018-03-01

We describe a strategy to solve differential equations for Feynman integrals by powers series expansions near singular points and to obtain high precision results for the corresponding master integrals. We consider Feynman integrals with two scales, i.e. non-trivially depending on one variable. The corresponding algorithm is oriented at situations where canonical form of the differential equations is impossible. We provide a computer code constructed with the help of our algorithm for a simple example of four-loop generalized sunset integrals with three equal non-zero masses and two zero masses. Our code gives values of the master integrals at any given point on the real axis with a required accuracy and a given order of expansion in the regularization parameter ɛ.

Parallel PWTD-Accelerated Explicit Solution of the Time Domain Electric Field Volume Integral Equation

KAUST Repository

Liu, Yang

2016-03-25

A parallel plane-wave time-domain (PWTD)-accelerated explicit marching-on-in-time (MOT) scheme for solving the time domain electric field volume integral equation (TD-EFVIE) is presented. The proposed scheme leverages pulse functions and Lagrange polynomials to spatially and temporally discretize the electric flux density induced throughout the scatterers, and a finite difference scheme to compute the electric fields from the Hertz electric vector potentials radiated by the flux density. The flux density is explicitly updated during time marching by a predictor-corrector (PC) scheme and the vector potentials are efficiently computed by a scalar PWTD scheme. The memory requirement and computational complexity of the resulting explicit PWTD-PC-EFVIE solver scale as ( log ) s s O N N and ( ) s t O N N , respectively. Here, s N is the number of spatial basis functions and t N is the number of time steps. A scalable parallelization of the proposed MOT scheme on distributed- memory CPU clusters is described. The efficiency, accuracy, and applicability of the resulting (parallelized) PWTD-PC-EFVIE solver are demonstrated via its application to the analysis of transient electromagnetic wave interactions on canonical and real-life scatterers represented with up to 25 million spatial discretization elements.

Parallel PWTD-Accelerated Explicit Solution of the Time Domain Electric Field Volume Integral Equation

KAUST Repository

Liu, Yang; Al-Jarro, Ahmed; Bagci, Hakan; Michielssen, Eric

2016-01-01

A parallel plane-wave time-domain (PWTD)-accelerated explicit marching-on-in-time (MOT) scheme for solving the time domain electric field volume integral equation (TD-EFVIE) is presented. The proposed scheme leverages pulse functions and Lagrange polynomials to spatially and temporally discretize the electric flux density induced throughout the scatterers, and a finite difference scheme to compute the electric fields from the Hertz electric vector potentials radiated by the flux density. The flux density is explicitly updated during time marching by a predictor-corrector (PC) scheme and the vector potentials are efficiently computed by a scalar PWTD scheme. The memory requirement and computational complexity of the resulting explicit PWTD-PC-EFVIE solver scale as ( log ) s s O N N and ( ) s t O N N , respectively. Here, s N is the number of spatial basis functions and t N is the number of time steps. A scalable parallelization of the proposed MOT scheme on distributed- memory CPU clusters is described. The efficiency, accuracy, and applicability of the resulting (parallelized) PWTD-PC-EFVIE solver are demonstrated via its application to the analysis of transient electromagnetic wave interactions on canonical and real-life scatterers represented with up to 25 million spatial discretization elements.

Constructing New Discrete Integrable Coupling System for Soliton Equation by Kronecker Product

International Nuclear Information System (INIS)

Yu Fajun; Zhang Hongqing

2008-01-01

It is shown that the Kronecker product can be applied to constructing new discrete integrable coupling system of soliton equation hierarchy in this paper. A direct application to the fractional cubic Volterra lattice spectral problem leads to a novel integrable coupling system of soliton equation hierarchy. It is also indicated that the study of discrete integrable couplings by using the Kronecker product is an efficient and straightforward method. This method can be used generally

The Magnetic Physical Optics Scattered Field in Terms of a Line Integral

DEFF Research Database (Denmark)

Meincke, Peter; Breinbjerg, Olav; Jørgensen, Erik

2000-01-01

An exact line integral representation Is derived for the magnetic physical optics field scattered by a perfectly electrically conducting planar plate illuminated by a magnetic Hertzian dipole. A numerical example is presented to illustrate the exactness of the line integral representation...

CALL FOR PAPERS: Special issue on Symmetries and Integrability of Difference Equations

Science.gov (United States)

Doliwa, Adam; Korhonen, Risto; Lafortune, Stephane

2006-10-01

This is a call for contributions to a special issue of Journal of Physics A: Mathematical and General entitled `Special issue on Symmetries and Integrability of Difference Equations' as featured at the SIDE VII meeting held during July 2006 in Melbourne (http://web.maths.unsw.edu.au/%7Eschief/side/side.html). Participants at that meeting, as well as other researchers working in the field of difference equations and discrete systems, are invited to submit a research paper to this issue. This meeting was the seventh of a series of biennial meetings devoted to the study of integrable difference equations and related topics. The notion of integrability was first introduced in the 19th century in the context of classical mechanics with the definition of Liouville integrability for Hamiltonian flows. Since then, several notions of integrability have been introduced for partial and ordinary differential equations. Closely related to integrability theory is the symmetry analysis of nonlinear evolution equations. Symmetry analysis takes advantage of the Lie group structure of a given equation to study its properties. Together, integrability theory and symmetry analysis provide the main method by which nonlinear evolution equations can be solved explicitly. Difference equations, just as differential equations, are important in numerous fields of science and have a wide variety of applications in such areas as: mathematical physics, computer visualization, numerical analysis, mathematical biology, economics, combinatorics, quantum field theory, etc. It is thus crucial to develop tools to study and solve difference equations. While the theory of symmetry and integrability for differential equations is now well-established, this is not yet the case for discrete equations. The situation has undergone impressive development in recent years and has affected a broad range of fields, including the theory of special functions, quantum integrable systems, numerical analysis, cellular

Bridging Ground Validation and Algorithms: Using Scattering and Integral Tables to Incorporate Observed DSD Correlations into Satellite Algorithms

Science.gov (United States)

Williams, C. R.

2012-12-01

The NASA Global Precipitation Mission (GPM) raindrop size distribution (DSD) Working Group is composed of NASA PMM Science Team Members and is charged to "investigate the correlations between DSD parameters using Ground Validation (GV) data sets that support, or guide, the assumptions used in satellite retrieval algorithms." Correlations between DSD parameters can be used to constrain the unknowns and reduce the degrees-of-freedom in under-constrained satellite algorithms. Over the past two years, the GPM DSD Working Group has analyzed GV data and has found correlations between the mass-weighted mean raindrop diameter (Dm) and the mass distribution standard deviation (Sm) that follows a power-law relationship. This Dm-Sm power-law relationship appears to be robust and has been observed in surface disdrometer and vertically pointing radar observations. One benefit of a Dm-Sm power-law relationship is that a three parameter DSD can be modeled with just two parameters: Dm and Nw that determines the DSD amplitude. In order to incorporate observed DSD correlations into satellite algorithms, the GPM DSD Working Group is developing scattering and integral tables that can be used by satellite algorithms. Scattering tables describe the interaction of electromagnetic waves on individual particles to generate cross sections of backscattering, extinction, and scattering. Scattering tables are independent of the distribution of particles. Integral tables combine scattering table outputs with DSD parameters and DSD correlations to generate integrated normalized reflectivity, attenuation, scattering, emission, and asymmetry coefficients. Integral tables contain both frequency dependent scattering properties and cloud microphysics. The GPM DSD Working Group has developed scattering tables for raindrops at both Dual Precipitation Radar (DPR) frequencies and at all GMI radiometer frequencies less than 100 GHz. Scattering tables include Mie and T-matrix scattering with H- and V

Integrability in three dimensions: Algebraic Bethe ansatz for anyonic models

Directory of Open Access Journals (Sweden)

Sh. Khachatryan

2015-10-01

Full Text Available We extend basic properties of two dimensional integrable models within the Algebraic Bethe Ansatz approach to 2+1 dimensions and formulate the sufficient conditions for the commutativity of transfer matrices of different spectral parameters, in analogy with Yang–Baxter or tetrahedron equations. The basic ingredient of our models is the R-matrix, which describes the scattering of a pair of particles over another pair of particles, the quark-anti-quark (meson scattering on another quark-anti-quark state. We show that the Kitaev model belongs to this class of models and its R-matrix fulfills well-defined equations for integrability.

Recursive integral equations with positive kernel for lattice calculations

International Nuclear Information System (INIS)

Illuminati, F.; Isopi, M.

1990-11-01

A Kirkwood-Salzburg integral equation, with positive defined kernel, for the states of lattice models of statistical mechanics and quantum field theory is derived. The equation is defined in the thermodynamic limit, and its iterative solution is convergent. Moreover, positivity leads to an exact a priori bound on the iteration. The equation's relevance as a reliable algorithm for lattice calculations is therefore suggested, and it is illustrated with a simple application. It should provide a viable alternative to Monte Carlo methods for models of statistical mechanics and lattice gauge theories. 10 refs

Kernel integration scatter model for parallel beam gamma camera and SPECT point source response

International Nuclear Information System (INIS)

Marinkovic, P.M.

2001-01-01

Scatter correction is a prerequisite for quantitative single photon emission computed tomography (SPECT). In this paper a kernel integration scatter Scatter correction is a prerequisite for quantitative SPECT. In this paper a kernel integration scatter model for parallel beam gamma camera and SPECT point source response based on Klein-Nishina formula is proposed. This method models primary photon distribution as well as first Compton scattering. It also includes a correction for multiple scattering by applying a point isotropic single medium buildup factor for the path segment between the point of scatter an the point of detection. Gamma ray attenuation in the object of imaging, based on known μ-map distribution, is considered too. Intrinsic spatial resolution of the camera is approximated by a simple Gaussian function. Collimator is modeled simply using acceptance angles derived from the physical dimensions of the collimator. Any gamma rays satisfying this angle were passed through the collimator to the crystal. Septal penetration and scatter in the collimator were not included in the model. The method was validated by comparison with Monte Carlo MCNP-4a numerical phantom simulation and excellent results were obtained. The physical phantom experiments, to confirm this method, are planed to be done. (author)