Approximate relativistic corrections to atomic radial wave functions

International Nuclear Information System (INIS)

Cowan, R.D.; Griffin, D.C.

1976-01-01

The mass-velocity and Darwin terms of the one-electron-atom Pauli equation have been added to the Hartree-Fock differential equations by using the HX formula to calculate a local central field potential for use in these terms. Introduction of the quantum number j is avoided by omitting the spin-orbit term of the Pauli equation. The major relativistic effects, both direct and indirect, are thereby incorporated into the wave functions, while allowing retention of the commonly used nonrelativistic formulation of energy level calculations. The improvement afforded in calculated total binding energies, excitation energies, spin-orbit parameters, and expectation values of r/sub m/ is comparable with that provided by fully relativistic Dirac-Hartree-Fock calculations

The effect of meson wave function on heavy-quark fragmentation function

Energy Technology Data Exchange (ETDEWEB)

Moosavi Nejad, S.M. [Yazd University, Faculty of Physics (Iran, Islamic Republic of); Institute for Research in Fundamental Sciences (IPM), School of Particles and Accelerators, Tehran (Iran, Islamic Republic of)

2016-05-15

We calculate the process-independent fragmentation functions (FFs) for a heavy quark to fragment into heavy mesons considering the effects of meson wave function. In all previous works, where the FFs of heavy mesons or heavy baryons were calculated, a delta function form was approximated for the wave function of hadrons. Here, for the first time, we consider a typical mesonic wave function which is different from the delta function and is the nonrelativistic limit of the solution of Bethe-Salpeter equation with the QCD kernel. We present our numerical results for the heavy FFs and show how the proposed wave function improves the previous results. As an example, we focus on the fragmentation function for c-quark to split into S-wave D{sup 0} -meson and compare our results with experimental data from BELLE and CLEO. (orig.)

Theoretical calculation of shakeup intensities using Xa--SW wave functions

International Nuclear Information System (INIS)

Tse, J.S.; Loubriel, G.

1981-01-01

The ground and 1s core hole state molecular wave functions of CH 4 , NH 3 , H 2 O, and HF obtained from Xa--SW calculations using the touching spheres (TS) and overlapping spheres (OS) approximations are used to calculate the intensity of shakeup satellites observed in their ls core level photoelectron spectra. The sudden approximation was assumed in the calculation. In case of TS Xa--SW wave functions, the one electron overlap integral inside the intersphere was calculated via Green's theorem. For OS Xa--SW wave functions, the integration over the awkwardly shaped intersphere region was circumvented by distributing the intersphere charge into the atomic spheres according to the charge partition scheme suggested by Case and Karplus. Our results show that there are no significant differences between the shakeup energies calculated from the TS and OS approximations. However, shakeup intensities calculated from TS Xa--SW wave functions are more reliable and in better numerical agreement with experiment

Approximation of wave action flux velocity in strongly sheared mean flows

Science.gov (United States)

Banihashemi, Saeideh; Kirby, James T.; Dong, Zhifei

2017-08-01

Spectral wave models based on the wave action equation typically use a theoretical framework based on depth uniform current to account for current effects on waves. In the real world, however, currents often have variations over depth. Several recent studies have made use of a depth-weighted current U˜ due to [Skop, R. A., 1987. Approximate dispersion relation for wave-current interactions. J. Waterway, Port, Coastal, and Ocean Eng. 113, 187-195.] or [Kirby, J. T., Chen, T., 1989. Surface waves on vertically sheared flows: approximate dispersion relations. J. Geophys. Res. 94, 1013-1027.] in order to account for the effect of vertical current shear. Use of the depth-weighted velocity, which is a function of wavenumber (or frequency and direction) has been further simplified in recent applications by only utilizing a weighted current based on the spectral peak wavenumber. These applications do not typically take into account the dependence of U˜ on wave number k, as well as erroneously identifying U˜ as the proper choice for current velocity in the wave action equation. Here, we derive a corrected expression for the current component of the group velocity. We demonstrate its consistency using analytic results for a current with constant vorticity, and numerical results for a measured, strongly-sheared current profile obtained in the Columbia River. The effect of choosing a single value for current velocity based on the peak wave frequency is examined, and we suggest an alternate strategy, involving a Taylor series expansion about the peak frequency, which should significantly extend the range of accuracy of current estimates available to the wave model with minimal additional programming and data transfer.

Consequences of wave function orthogonality for medium energy nuclear reactions

International Nuclear Information System (INIS)

Noble, J.V.

1978-01-01

In the usual models of high-energy bound-state to continuum transitions no account is taken of the orthogonality of the bound and continuum wave functions. This orthogonality induces considerable cancellations in the overlap integrals expressing the transition amplitudes for reactions such as (e,e'p), (γ,p), and (π,N), which are simply not included in the distorted-wave Born-approximation calculations which to date remain the only computationally feasible heirarchy of approximations. The object of this paper is to present a new formulation of the bound-state to continuum transition problem, based upon flux conservation, in which the orthogonality of wave functions is taken into account ab initio. The new formulation, while exact if exact wave functions are used, offers the possibility of using approximate wave functions for the continuum states without doing violence to the cancellations induced by orthogonality. The method is applied to single-particle states obeying the Schroedinger and Dirac equations, as well as to a coupled-channel model in which absorptive processes can be described in a fully consistent manner. Several types of absorption vertex are considered, and in the (π,N) case the equivalence of pseudoscalar and pseudovector πNN coupling is seen to follow directly from wave function orthogonality

Excited state nuclear forces from the Tamm-Dancoff approximation to time-dependent density functional theory within the plane wave basis set framework

Science.gov (United States)

Hutter, Jürg

2003-03-01

An efficient formulation of time-dependent linear response density functional theory for the use within the plane wave basis set framework is presented. The method avoids the transformation of the Kohn-Sham matrix into the canonical basis and references virtual orbitals only through a projection operator. Using a Lagrangian formulation nuclear derivatives of excited state energies within the Tamm-Dancoff approximation are derived. The algorithms were implemented into a pseudo potential/plane wave code and applied to the calculation of adiabatic excitation energies, optimized geometries and vibrational frequencies of three low lying states of formaldehyde. An overall good agreement with other time-dependent density functional calculations, multireference configuration interaction calculations and experimental data was found.

Wave equation dispersion inversion using a difference approximation to the dispersion-curve misfit gradient

KAUST Repository

Zhang, Zhendong

2016-07-26

We present a surface-wave inversion method that inverts for the S-wave velocity from the Rayleigh wave dispersion curve using a difference approximation to the gradient of the misfit function. We call this wave equation inversion of skeletonized surface waves because the skeletonized dispersion curve for the fundamental-mode Rayleigh wave is inverted using finite-difference solutions to the multi-dimensional elastic wave equation. The best match between the predicted and observed dispersion curves provides the optimal S-wave velocity model. Our method can invert for lateral velocity variations and also can mitigate the local minimum problem in full waveform inversion with a reasonable computation cost for simple models. Results with synthetic and field data illustrate the benefits and limitations of this method. © 2016 Elsevier B.V.

Parametrization of the scattering wave functions of the Paris potential

International Nuclear Information System (INIS)

Loiseau, B.; Mathelitsch, L.

1996-10-01

The neutron-proton scattering wave functions of the Paris nucleon-nucleon potential are parametrized for partial waves of total angular momenta less than 5. The inner parts of the wave functions are approximated by polynomials with a continuous transition to the outer parts, which are given by the asymptotic regime and determined by the respective phase shifts. The scattering wave functions can then be calculated at any given energy below 400 MeV. Special attention is devoted to the zero-energy limit of the low partial waves. An easy-to-use FORTRAN program, which allows the user to calculate these parametrized wave functions, is available via electronic mail. (author)

Molecular distorted-wave Born approximation for ionization of H2 by electron impact

International Nuclear Information System (INIS)

Liu, Junbo; Liu, Dejun; Zhou, Yajun

2012-01-01

The molecular distorted-wave Born approximation is proposed to study the (e, 2e) reaction for H 2 targets. The wave functions of the incoming and outgoing electrons are obtained by solving the Lippmann-Schwinger equations, and the T-matrix in the Lippmann-Schwinger equations is calculated in a momentum space static-exchange-optical model. Triple differential cross sections are computed for incident energies of 100 and 250 eV in coplanar asymmetric geometry. Comparison of the present calculated results with the available experimental data in the literature reveals that there is good agreement. (paper)

Influence of wetting layer wave functions on carrier capture in quantum dots

DEFF Research Database (Denmark)

Markussen, Troels; Kristensen, Philip; Tromborg, Bjarne

2005-01-01

This work numerically solves the effective mass Schrodinger equation and shows that the capture times are strongly influenced by details of the continuum states not accounted for by the approximate wave functions. Results show that calculations of capture time for phonon mediated carrier capture...... from a wetting layer into a quantum dot depend critically on the approximations used for the wetting layer wave functions....

Relativistic deuteron wave function on light front

International Nuclear Information System (INIS)

Karmanov, V.A.

1980-01-01

In the framework of the one boson exchange model the approximate analytical expression for the deuteron wave function (WF) at relativistic relative momenta is obtained. WF depends on extra variable having the form of a unit vector and is determined by six functions instead of two ones (S-and D-waves) in the nonrelativistic case. At moderate momenta the WF is matched with WF in the Reid model. It is emphasized the importance of indication of the qualitative observed phenomena associated with change of parametrization and spin structure of relativistic deuteron WF

Approximate Dispersion Relations for Waves on Arbitrary Shear Flows

Science.gov (United States)

Ellingsen, S. À.; Li, Y.

2017-12-01

An approximate dispersion relation is derived and presented for linear surface waves atop a shear current whose magnitude and direction can vary arbitrarily with depth. The approximation, derived to first order of deviation from potential flow, is shown to produce good approximations at all wavelengths for a wide range of naturally occuring shear flows as well as widely used model flows. The relation reduces in many cases to a 3-D generalization of the much used approximation by Skop (1987), developed further by Kirby and Chen (1989), but is shown to be more robust, succeeding in situations where the Kirby and Chen model fails. The two approximations incur the same numerical cost and difficulty. While the Kirby and Chen approximation is excellent for a wide range of currents, the exact criteria for its applicability have not been known. We explain the apparently serendipitous success of the latter and derive proper conditions of applicability for both approximate dispersion relations. Our new model has a greater range of applicability. A second order approximation is also derived. It greatly improves accuracy, which is shown to be important in difficult cases. It has an advantage over the corresponding second-order expression proposed by Kirby and Chen that its criterion of accuracy is explicitly known, which is not currently the case for the latter to our knowledge. Our second-order term is also arguably significantly simpler to implement, and more physically transparent, than its sibling due to Kirby and Chen.Plain Language SummaryIn order to answer key questions such as how the ocean surface affects the climate, erodes the coastline and transports nutrients, we must understand how waves move. This is not so easy when depth varying currents are present, as they often are in coastal waters. We have developed a modeling tool for accurately predicting wave properties in such situations, ready for use, for example, in the complex oceanographic computer models. Our

On wave functions of mesons involving the s-, c- and b-quarks

International Nuclear Information System (INIS)

Zhitnitskij, A.R.; Zhitnitskij, I.R.; Chernyak, V.L.

1983-01-01

The wave function components of pseudoscalar and vestor mesons which are antisymmertric with respect to permutation of the quark momenta are studied. The results are as follows: elt xsub(s)-xsub(u) > sub(K) approximately equal to 0.11 for the K meson, sub(K*) approximately equal to 0.15-C.20 for the K* meson, being a mean fraction of the longitudinal momentum transferred by the s(u) quark. The following estimates are obtained: / approximately equal to 0.20-0.25; / approximately equal to 0.8x10 -2 . The asymptotics of the K 0 -meson form factor and the etasub(c) → KK* decay width are found. Properties of the wave functions of mesons which contain a light and a heavy quark (D, B, ...) are considered. For the B 0 meson approximately equal to 0.10 is found. Arguments are given supporting nonenhancement of the amplitudes of the processes involving D mesons compared to similar K-meson amplitudes. A simple way is suggested to determine the asymptotic form of various wave functions

The wave equation: From eikonal to anti-eikonal approximation

Directory of Open Access Journals (Sweden)

Luis Vázquez

2016-06-01

Full Text Available When the refractive index changes very slowly compared to the wave-length we may use the eikonal approximation to the wave equation. In the opposite case, when the refractive index highly variates over the distance of one wave-length, we have what can be termed as the anti-eikonal limit. This situation is addressed in this work. The anti-eikonal limit seems to be a relevant tool in the modelling and design of new optical media. Besides, it describes a basic universal behaviour, independent of the actual values of the refractive index and, thus, of the media, for the components of a wave with wave-length much greater than the characteristic scale of the refractive index.

On the semi-classical approximation to the wave function of the universe and its stochastic interpretation

International Nuclear Information System (INIS)

Pollock, M.D.

1988-01-01

In quantum cosmology, a wave function Ψ for a given theory can be obtained by solving the Wheeler-DeWitt equation, using the semi-classical approximation to the path integral over euclidean metrics to impose the boundary condition, as described by Hawking and his collaborators. If the universe is expanding as a quasi-de Sitter space-time, then it is possible to derive a Fokker-Planck equation for the probability distribution P, as shown by Starobinsky. Arguing by analogy with quantum mechanics in flat space-time, one would expect that P ∝ ΨΨ * . We examine this assertion by reference to the scale-invariant theory L = -1/24 βR 2 , whose wave function has been calculated in mini-superspace by Horowitz, and whose classical solutions are de Sitter space-times. It appears that deviations from the relation P ∝ ΨΨ * are attributable to long-wavelength fluctuations δΦ e ≅ H/2π in the effective inflaton field Φ c =√(βR)=√(12β) H. Their existence is taken into account in the derivation of the Fokker-Planck equation, but not in the derivation of Ψ when this is restricted to mini-superspace. In the limit β → ∞, we find that δΦ e /Φ c → 0 and that P ∝ ΨΨ * . The scale-invariant theory L = (1/2εφ 2 R-1/4λΦ 4 ) can be similarly analyzed. Inclusion of a kinetic term 1/2Φ; k Φ ;k destroys this similarity, which is restored however upon addition of a term (-1/24βR 2 ). (orig.)

Optical bistability without the rotating wave approximation

Energy Technology Data Exchange (ETDEWEB)

Sharaby, Yasser A., E-mail: Yasser_Sharaby@hotmail.co [Physics Department, Faculty of Applied Sciences, Suez Canal University, Suez (Egypt); Joshi, Amitabh, E-mail: ajoshi@eiu.ed [Department of Physics, Eastern Illinois University, Charleston, IL 61920 (United States); Hassan, Shoukry S., E-mail: Shoukryhassan@hotmail.co [Mathematics Department, College of Science, University of Bahrain, P.O. Box 32038 (Bahrain)

2010-04-26

Optical bistability for two-level atomic system in a ring cavity is investigated outside the rotating wave approximation (RWA) using non-autonomous Maxwell-Bloch equations with Fourier decomposition up to first harmonic. The first harmonic output field component exhibits reversed or closed loop bistability simultaneously with the usual (anti-clockwise) bistability in the fundamental field component.

Optical bistability without the rotating wave approximation

International Nuclear Information System (INIS)

Sharaby, Yasser A.; Joshi, Amitabh; Hassan, Shoukry S.

2010-01-01

Optical bistability for two-level atomic system in a ring cavity is investigated outside the rotating wave approximation (RWA) using non-autonomous Maxwell-Bloch equations with Fourier decomposition up to first harmonic. The first harmonic output field component exhibits reversed or closed loop bistability simultaneously with the usual (anti-clockwise) bistability in the fundamental field component.

WKB wave function for many-variable systems

International Nuclear Information System (INIS)

Sakita, B.; Tzani, R.

1986-01-01

The WKB method is a non-perturbative semi-classical method in quantum mechanics. The method for a system of one degree of freedom is well known and described in standard textbooks. The method for a system with many degrees of freedom especially for quantum fields is more involved. There exist two methods: Feynman path integral and Schrodinger wave function. The Feynman path integral WKB method is essentially a stationary phase approximation for Feynman path integrals. The WKB Schrodinger wave function method is on the other hand an extension of the standard WKB to many-variable systems

On a novel iterative method to compute polynomial approximations to Bessel functions of the first kind and its connection to the solution of fractional diffusion/diffusion-wave problems

International Nuclear Information System (INIS)

Yuste, Santos Bravo; Abad, Enrique

2011-01-01

We present an iterative method to obtain approximations to Bessel functions of the first kind J p (x) (p > -1) via the repeated application of an integral operator to an initial seed function f 0 (x). The class of seed functions f 0 (x) leading to sets of increasingly accurate approximations f n (x) is considerably large and includes any polynomial. When the operator is applied once to a polynomial of degree s, it yields a polynomial of degree s + 2, and so the iteration of this operator generates sets of increasingly better polynomial approximations of increasing degree. We focus on the set of polynomial approximations generated from the seed function f 0 (x) = 1. This set of polynomials is useful not only for the computation of J p (x) but also from a physical point of view, as it describes the long-time decay modes of certain fractional diffusion and diffusion-wave problems.

Quantum scattering beyond the plane-wave approximation

Science.gov (United States)

Karlovets, Dmitry

2017-12-01

While a plane-wave approximation in high-energy physics works well in a majority of practical cases, it becomes inapplicable for scattering of the vortex particles carrying orbital angular momentum, of Airy beams, of the so-called Schrödinger cat states, and their generalizations. Such quantum states of photons, electrons and neutrons have been generated experimentally in recent years, opening up new perspectives in quantum optics, electron microscopy, particle physics, and so forth. Here we discuss the non-plane-wave effects in scattering brought about by the novel quantum numbers of these wave packets. For the well-focused electrons of intermediate energies, already available at electron microscopes, the corresponding contribution can surpass that of the radiative corrections. Moreover, collisions of the cat-like superpositions of such focused beams with atoms allow one to probe effects of the quantum interference, which have never played any role in particle scattering.

Photon wave function

OpenAIRE

Bialynicki-Birula, Iwo

2005-01-01

Photon wave function is a controversial concept. Controversies stem from the fact that photon wave functions can not have all the properties of the Schroedinger wave functions of nonrelativistic wave mechanics. Insistence on those properties that, owing to peculiarities of photon dynamics, cannot be rendered, led some physicists to the extreme opinion that the photon wave function does not exist. I reject such a fundamentalist point of view in favor of a more pragmatic approach. In my view, t...

Computing many-body wave functions with guaranteed precision: the first-order Møller-Plesset wave function for the ground state of helium atom.

Science.gov (United States)

Bischoff, Florian A; Harrison, Robert J; Valeev, Edward F

2012-09-14

We present an approach to compute accurate correlation energies for atoms and molecules using an adaptive discontinuous spectral-element multiresolution representation for the two-electron wave function. Because of the exponential storage complexity of the spectral-element representation with the number of dimensions, a brute-force computation of two-electron (six-dimensional) wave functions with high precision was not practical. To overcome the key storage bottlenecks we utilized (1) a low-rank tensor approximation (specifically, the singular value decomposition) to compress the wave function, and (2) explicitly correlated R12-type terms in the wave function to regularize the Coulomb electron-electron singularities of the Hamiltonian. All operations necessary to solve the Schrödinger equation were expressed so that the reconstruction of the full-rank form of the wave function is never necessary. Numerical performance of the method was highlighted by computing the first-order Møller-Plesset wave function of a helium atom. The computed second-order Møller-Plesset energy is precise to ~2 microhartrees, which is at the precision limit of the existing general atomic-orbital-based approaches. Our approach does not assume special geometric symmetries, hence application to molecules is straightforward.

Approximate equations at breaking for nearshore wave transformation coefficients

Digital Repository Service at National Institute of Oceanography (India)

Chandramohan, P.; Nayak, B.U.; SanilKumar, V.

Based on small amplitude wave theory approximate equations are evaluated for determining the coefficients of shoaling, refraction, bottom friction, bottom percolation and viscous dissipation at breaking. The results obtainEd. by these equations...

Parabolic approximation method for fast magnetosonic wave propagation in tokamaks

International Nuclear Information System (INIS)

Phillips, C.K.; Perkins, F.W.; Hwang, D.Q.

1985-07-01

Fast magnetosonic wave propagation in a cylindrical tokamak model is studied using a parabolic approximation method in which poloidal variations of the wave field are considered weak in comparison to the radial variations. Diffraction effects, which are ignored by ray tracing mthods, are included self-consistently using the parabolic method since continuous representations for the wave electromagnetic fields are computed directly. Numerical results are presented which illustrate the cylindrical convergence of the launched waves into a diffraction-limited focal spot on the cyclotron absorption layer near the magnetic axis for a wide range of plasma confinement parameters

Theory of magnetohydrodynamic waves: The WKB approximation revisited

International Nuclear Information System (INIS)

Barnes, A.

1992-01-01

Past treatments of the eikonal or WKB theory of the propagation of magnetohydrodynamics waves have assumed a strictly isentropic background. IF in fact there is a gradient in the background entropy, then in second order in the WKB ordering, adiabatic fluctuations (in the Lagrangian sense) are not strictly isentropic in the Eulerian sense. This means that in the second order of the WKB expansion, which determines the variation of wave amplitude along rays, the violation of isentropy must be accounted for. The present paper revisits the derivation of the WKB approximation for small-amplitude magnetohydrodynamic waves, allowing for possible spatial variation of the background entropy. The equation of variation of wave amplitude is rederived; it is a bilinear equation which, it turns out, can be recast in the action conservation form. It is shown that this action conservation equation is in fact equivalent to the action conservation law obtained from Lagrangian treatments

Neutron wave reflexions in interface media with transport equation P1 approximation

International Nuclear Information System (INIS)

Oliveira Vellozo, S. de.

1977-01-01

The propagation of neutron waves in non multiplying media is investigated employing the Telegrapher's equation obtained from the P 1 approximation of the time, space and energy dependent Boltzmann equation. Solution of the problem of propagation of sinusoidally modulated source incident on one face of the medium is obtained by analysing the Fourier component of a pulsed source introduced, for the corresponding frequency. The amplitude and the phase of the flux are computed as a function of frequency in media consisting of one, two and three regions in order to study the effects of reflection at the interfaces. The results are compared with those from the Diffusion approximation obtained by neglecting the term involving the second order time derivative. (author)

Probability function of breaking-limited surface elevation. [wind generated waves of ocean

Science.gov (United States)

Tung, C. C.; Huang, N. E.; Yuan, Y.; Long, S. R.

1989-01-01

The effect of wave breaking on the probability function of surface elevation is examined. The surface elevation limited by wave breaking zeta sub b(t) is first related to the original wave elevation zeta(t) and its second derivative. An approximate, second-order, nonlinear, non-Gaussian model for zeta(t) of arbitrary but moderate bandwidth is presented, and an expression for the probability density function zeta sub b(t) is derived. The results show clearly that the effect of wave breaking on the probability density function of surface elevation is to introduce a secondary hump on the positive side of the probability density function, a phenomenon also observed in wind wave tank experiments.

Dislocation kinetics and the acoustic-wave approximation for liquids

International Nuclear Information System (INIS)

Stout, R.B.

1983-03-01

A dislocation-dependent model for liquids describes the lattice deformation and the fluidity deformation as additive deformations. The lattice deformation represents distortions of an atom's potential energy structure and is a recoverable deformation response. The fluidity deformation represents discontinuous repositioning of atoms by dislocation kinetics in the lattice structure and is a nonrecoverable deformation response. From this model, one concludes that in liquids the acoustic-wave approximation is a description of a recoverable oscillation deformation that has dissipation because of dislocation kinetics. Other more-complex waves may exist, but such waves would rapidly disappear because of the small thermodynamic potential for dislocation kinetics in liquids

Nonlinear approximation with general wave packets

DEFF Research Database (Denmark)

Borup, Lasse; Nielsen, Morten

2005-01-01

We study nonlinear approximation in the Triebel-Lizorkin spaces with dictionaries formed by dilating and translating one single function g. A general Jackson inequality is derived for best m-term approximation with such dictionaries. In some special cases where g has a special structure, a complete...

Dynamical cluster approximation plus semiclassical approximation study for a Mott insulator and d-wave pairing

Science.gov (United States)

Kim, SungKun; Lee, Hunpyo

2017-06-01

Via a dynamical cluster approximation with N c = 4 in combination with a semiclassical approximation (DCA+SCA), we study the doped two-dimensional Hubbard model. We obtain a plaquette antiferromagnetic (AF) Mott insulator, a plaquette AF ordered metal, a pseudogap (or d-wave superconductor) and a paramagnetic metal by tuning the doping concentration. These features are similar to the behaviors observed in copper-oxide superconductors and are in qualitative agreement with the results calculated by the cluster dynamical mean field theory with the continuous-time quantum Monte Carlo (CDMFT+CTQMC) approach. The results of our DCA+SCA differ from those of the CDMFT+CTQMC approach in that the d-wave superconducting order parameters are shown even in the high doped region, unlike the results of the CDMFT+CTQMC approach. We think that the strong plaquette AF orderings in the dynamical cluster approximation (DCA) with N c = 4 suppress superconducting states with increasing doping up to strongly doped region, because frozen dynamical fluctuations in a semiclassical approximation (SCA) approach are unable to destroy those orderings. Our calculation with short-range spatial fluctuations is initial research, because the SCA can manage long-range spatial fluctuations in feasible computational times beyond the CDMFT+CTQMC tool. We believe that our future DCA+SCA calculations should supply information on the fully momentum-resolved physical properties, which could be compared with the results measured by angle-resolved photoemission spectroscopy experiments.

Discontinuous approximate molecular electronic wave-functions

International Nuclear Information System (INIS)

Stuebing, E.W.; Weare, J.H.; Parr, R.G.

1977-01-01

Following Kohn, Schlosser and Marcus and Weare and Parr an energy functional is defined for a molecular problem which is stationary in the neighborhood of the exact solution and permits the use of trial functions that are discontinuous. The functional differs from the functional of the standard Rayleigh--Ritz method in the replacement of the usual kinetic energy operators circumflex T(μ) with operators circumflex T'(μ) = circumflex T(μ) + circumflex I(μ) generates contributions from surfaces of nonsmooth behavior. If one uses the nabla PSI . nabla PSI way of writing the usual kinetic energy contributions, one must add surface integrals of the product of the average of nabla PSI and the change of PSI across surfaces of discontinuity. Various calculations are carried out for the hydrogen molecule-ion and the hydrogen molecule. It is shown that ab initio calculations on molecules can be carried out quite generally with a basis of atomic orbitals exactly obeying the zero-differential overlap (ZDO) condition, and a firm basis is thereby provided for theories of molecular electronic structure invoking the ZDO aoproximation. It is demonstrated that a valence bond theory employing orbitals exactly obeying ZDO can provide an adequate account of chemical bonding, and several suggestions are made regarding molecular orbital methods

Hierarchical wave functions revisited

International Nuclear Information System (INIS)

Li Dingping.

1997-11-01

We study the hierarchical wave functions on a sphere and on a torus. We simplify some wave functions on a sphere or a torus using the analytic properties of wave functions. The open question, the construction of the wave function for quasi electron excitation on a torus, is also solved in this paper. (author)

The choice of optimal Discrete Interaction Approximation to the kinetic integral for ocean waves

Directory of Open Access Journals (Sweden)

V. G. Polnikov

2003-01-01

Full Text Available A lot of discrete configurations for the four-wave nonlinear interaction processes have been calculated and tested by the method proposed earlier in the frame of the concept of Fast Discrete Interaction Approximation to the Hasselmann's kinetic integral (Polnikov and Farina, 2002. It was found that there are several simple configurations, which are more efficient than the one proposed originally in Hasselmann et al. (1985. Finally, the optimal multiple Discrete Interaction Approximation (DIA to the kinetic integral for deep-water waves was found. Wave spectrum features have been intercompared for a number of different configurations of DIA, applied to a long-time solution of kinetic equation. On the basis of this intercomparison the better efficiency of the configurations proposed was confirmed. Certain recommendations were given for implementation of new approximations to the wave forecast practice.

Christodoulou's nonlinear gravitational-wave memory: Evaluation in the quadrupole approximation

International Nuclear Information System (INIS)

Wiseman, A.G.; Will, C.M.

1991-01-01

Christodoulou has found a new nonlinear contribution to the net change in the wave form caused by the passage of a burst of gravity waves (''memory of the burst''). We argue that this effect is nothing but the gravitational wave form generated by the stress energy in the burst itself. We derive an explicit formula for this effect in terms of a retarded-time integral of products of time derivatives of wave-zone gravitational wave forms. The resulting effect corresponds in size to a correction 2.5 post-Newtonian orders [O((Gm/rc 2 ) 5/2 ) =(O(v/c) 5 )] beyond the quadrupole approximation, and is therefore negligible for all but the most relativistic of systems. For gravitational bremsstrahlung from two stars moving at 300 km s -1 , the effect is much less than 10 -10 of the usual linear quadrupole wave form, while for a system of coalescing binary compact objects we estimate that the effect is of order 10 -1 for two neutron stars

Diffusive Wave Approximation to the Shallow Water Equations: Computational Approach

KAUST Repository

Collier, Nathan; Radwan, Hany; Dalcin, Lisandro; Calo, Victor M.

2011-01-01

We discuss the use of time adaptivity applied to the one dimensional diffusive wave approximation to the shallow water equations. A simple and computationally economical error estimator is discussed which enables time-step size adaptivity

Nuclear Hartree-Fock approximation testing and other related approximations

International Nuclear Information System (INIS)

Cohenca, J.M.

1970-01-01

Hartree-Fock, and Tamm-Dancoff approximations are tested for angular momentum of even-even nuclei. Wave functions, energy levels and momenta are comparatively evaluated. Quadripole interactions are studied following the Elliott model. Results are applied to Ne 20 [pt

Analytic Approximations for Soliton Solutions of Short-Wave Models for Camassa-Holm and Degasperis-Procesi Equations

International Nuclear Information System (INIS)

Yang Pei; Li Zhibin; Chen Yong

2010-01-01

In this paper, the short-wave model equations are investigated, which are associated with the Camassa-Holm (CH) and Degasperis-Procesi (DP) shallow-water wave equations. Firstly, by means of the transformation of the independent variables and the travelling wave transformation, the partial differential equation is reduced to an ordinary differential equation. Secondly, the equation is solved by homotopy analysis method. Lastly, by the transformations back to the original independent variables, the solution of the original partial differential equation is obtained. The two types of solutions of the short-wave models are obtained in parametric form, one is one-cusp soliton for the CH equation while the other one is one-loop soliton for the DP equation. The approximate analytic solutions expressed by a series of exponential functions agree well with the exact solutions. It demonstrates the validity and great potential of homotopy analysis method for complicated nonlinear solitary wave problems. (general)

Analytic structure of the wave function for a hydrogen atom in an analytic potential

International Nuclear Information System (INIS)

Hill, R.N.

1984-01-01

The rate of convergence of an approximate method for solving Schroedinger's equation depends on the ability of the approximating sequence to mimic the analytic structure of the unknown exact wave function. Thus a knowledge of the analytic structure of the wave function can be of great value when approximation schemes are designed. Consider the Schroedinger equation [- 1/2 del 2 -r -1 +V(r)]Psi(r) = EPsi(r) for a hydrogen atom in a potential V(r). The general theory of elliptic partial differential equations implies that Psi is analytic at regular points, but no general theory is available at singular points. The present paper investigates the Coulomb singular point at r = 0 and shows that, if V(r) = V 1 (x, y, z)+rV 2 (x, y, z) where V 1 and V 2 are analytic functions of x, y, z at x = y = z = 0, then the wave function has the form Psi(r) = Psi 1 (x, y, z)+rPsi 2 (x, y, z) where Psi 1 and Psi 2 are analytic functions of x, y, z at x = y = z = 0

Description of the nucleon wave function as a sum of well-chosen Gaussian functions

International Nuclear Information System (INIS)

Roux, C.; Silvestre-Brac, B.

1995-01-01

We study in detail the possibility of describing the nucleon (three quark-system) wave function as a superposition of Gaussian functions. A Faddeev treatment including 8 amplitudes is performed and taken as reference for the exact values. Several approximations are proposed and compared carefully to the exact solutions. Three different potentials have been tested and several observables are considered. (author)

Energy spectra and wave function of trigonometric Rosen-Morse potential as an effective quantum chromodynamics potential in D-dimensions

Energy Technology Data Exchange (ETDEWEB)

Deta, U. A., E-mail: utamaalan@yahoo.co.id [Theoretical Physics Group, Physics Department of Post Graduate Program, Sebelas Maret University, Jl. Ir. Sutami 36A, Surakarta 57126, Indonesia and Physics Department, State University of Surabaya, Jl. Ketintang, Surabaya 60231 (Indonesia); Suparmi,; Cari,; Husein, A. S.; Yuliani, H.; Khaled, I. K. A.; Luqman, H.; Supriyanto [Theoretical Physics Group, Physics Department of Post Graduate Program, Sebelas Maret University, Jl. Ir. Sutami 36A, Surakarta 57126 (Indonesia)

2014-09-30

The Energy Spectra and Wave Function of Schrodinger equation in D-Dimensions for trigonometric Rosen-Morse potential were investigated analytically using Nikiforov-Uvarov method. This potential captures the essential traits of the quark-gluon dynamics of Quantum Chromodynamics. The approximate energy spectra are given in the close form and the corresponding approximate wave function for arbitrary l-state (l ≠ 0) in D-dimensions are formulated in the form of differential polynomials. The wave function of this potential unnormalizable for general case. The wave function of this potential unnormalizable for general case. The existence of extra dimensions (centrifugal factor) and this potential increase the energy spectra of system.

Simulation of angle-resolved photoemission spectra by approximating the final state by a plane wave: From graphene to polycyclic aromatic hydrocarbon molecules

Energy Technology Data Exchange (ETDEWEB)

Puschnig, Peter, E-mail: peter.puschnig@uni-graz.at; Lüftner, Daniel

2015-04-15

Highlights: • Computational study on angular dependent photoemission spectroscopy. • Graphene and polycyclic aromatic hydrocarbon molecules. • Plane wave final state approximation accounts for experimental findings. - Abstract: We present a computational study on the angular-resolved photoemission spectra (ARPES) from a number of polycyclic aromatic hydrocarbons and graphene. Our theoretical approach is based on ab-initio density functional theory and the one-step model where we greatly simplify the evaluation of the matrix element by assuming a plane wave for the final state. Before comparing our ARPES simulations with available experimental data, we discuss how typical approximations for the exchange-correlation energy affect orbital energies. In particular, we show that by employing a hybrid functional, considerable improvement can be obtained over semi-local functionals in terms of band widths and relative energies of π and σ states. Our ARPES simulations for graphene show that the plane wave final state approximation provides indeed an excellent description when compared to experimental band maps and constant binding energy maps. Furthermore, our ARPES simulations for a number of polycyclic aromatic molecules from the oligo-acene, oligo-phenylene, phen-anthrene families as well as for disc-shaped molecules nicely illustrate the evolution of the electronic structure from molecules with increasing size towards graphene.

Riemann zeta function from wave-packet dynamics

DEFF Research Database (Denmark)

Mack, R.; Dahl, Jens Peder; Moya-Cessa, H.

2010-01-01

We show that the time evolution of a thermal phase state of an anharmonic oscillator with logarithmic energy spectrum is intimately connected to the generalized Riemann zeta function zeta(s, a). Indeed, the autocorrelation function at a time t is determined by zeta (sigma + i tau, a), where sigma...... index of JWKB. We compare and contrast exact and approximate eigenvalues of purely logarithmic potentials. Moreover, we use a numerical method to find a potential which leads to exact logarithmic eigenvalues. We discuss possible realizations of Riemann zeta wave-packet dynamics using cold atoms...

Diagonal Born-Oppenheimer correction for coupled-cluster wave-functions

Science.gov (United States)

Shamasundar, K. R.

2018-06-01

We examine how geometry-dependent normalisation freedom of electronic wave-functions affects extraction of a meaningful diagonal Born-Oppenheimer correction (DBOC) to the ground-state Born-Oppenheimer potential energy surface (PES). By viewing this freedom as a kind of gauge-freedom, it is shown that DBOC and the resulting associated mass-dependent adiabatic PES are gauge-invariant quantities. A sum-over-states (SOS) formula for DBOC which explicitly exhibits this invariance is derived. A biorthogonal formulation suitable for DBOC computations using standard unnormalised coupled-cluster (CC) wave-functions is presented. This is shown to lead to a biorthogonal version of SOS formula with similar properties. On this basis, different computational schemes for evaluating DBOC using approximate CC wave-functions are derived. One of this agrees with the formula used in the current literature. The connection to adiabatic-to-diabatic transformations in non-adiabatic dynamics is explored and complications arising from biorthogonal nature of CC theory are identified.

Approximate optimal tracking control for near-surface AUVs with wave disturbances

Science.gov (United States)

Yang, Qing; Su, Hao; Tang, Gongyou

2016-10-01

This paper considers the optimal trajectory tracking control problem for near-surface autonomous underwater vehicles (AUVs) in the presence of wave disturbances. An approximate optimal tracking control (AOTC) approach is proposed. Firstly, a six-degrees-of-freedom (six-DOF) AUV model with its body-fixed coordinate system is decoupled and simplified and then a nonlinear control model of AUVs in the vertical plane is given. Also, an exosystem model of wave disturbances is constructed based on Hirom approximation formula. Secondly, the time-parameterized desired trajectory which is tracked by the AUV's system is represented by the exosystem. Then, the coupled two-point boundary value (TPBV) problem of optimal tracking control for AUVs is derived from the theory of quadratic optimal control. By using a recently developed successive approximation approach to construct sequences, the coupled TPBV problem is transformed into a problem of solving two decoupled linear differential sequences of state vectors and adjoint vectors. By iteratively solving the two equation sequences, the AOTC law is obtained, which consists of a nonlinear optimal feedback item, an expected output tracking item, a feedforward disturbances rejection item, and a nonlinear compensatory term. Furthermore, a wave disturbances observer model is designed in order to solve the physically realizable problem. Simulation is carried out by using the Remote Environmental Unit (REMUS) AUV model to demonstrate the effectiveness of the proposed algorithm.

Delta function excitation of waves in the earth's ionosphere

Science.gov (United States)

Vidmar, R. J.; Crawford, F. W.; Harker, K. J.

1983-01-01

Excitation of the earth's ionosphere by delta function current sheets is considered, and the temporal and spatial evolution of wave packets is analyzed for a two-component collisional F2 layer. Approximations of an inverse Fourier-Laplace transform via saddle point methods provide plots of typical wave packets. These illustrate cold plasma wave theory and may be used as a diagnostic tool since it is possible to relate specific features, e.g., the frequency of a modulation envelope, to plasma parameters such as the electron cyclotron frequency. It is also possible to deduce the propagation path length and orientation of a remote radio beacon.

Non-Hermitian wave packet approximation for coupled two-level systems in weak and intense fields

Energy Technology Data Exchange (ETDEWEB)

Puthumpally-Joseph, Raiju; Charron, Eric [Institut des Sciences Moléculaires d’Orsay (ISMO), CNRS, Univ. Paris-Sud, Université Paris-Saclay, F-91405 Orsay (France); Sukharev, Maxim [Science and Mathematics Faculty, College of Letters and Sciences, Arizona State University, Mesa, Arizona 85212 (United States)

2016-04-21

We introduce a non-Hermitian Schrödinger-type approximation of optical Bloch equations for two-level systems. This approximation provides a complete and accurate description of the coherence and decoherence dynamics in both weak and strong laser fields at the cost of losing accuracy in the description of populations. In this approach, it is sufficient to propagate the wave function of the quantum system instead of the density matrix, providing that relaxation and dephasing are taken into account via automatically adjusted time-dependent gain and decay rates. The developed formalism is applied to the problem of scattering and absorption of electromagnetic radiation by a thin layer comprised of interacting two-level emitters.

Four-body wave function of π3He-system at the threshold energy

International Nuclear Information System (INIS)

Pupyshev, V.V.; Rakityanskij, S.A.

1985-01-01

On the basis of approximate four-body equations the wave function of π 3 He-system is calculated at zero kinetic energy of the pion. In the case when distances between all four particles are comparable with the nucleus size a strong distortion of the wave function of (3N)-subsystem caused by the presence of the pion is found. The calculated four-body function is represented in a semianalytical form, which makes it possible to apply it in different calculations

Dispersion relation for Bernstein waves using a new transformation for the modified Bessel function

International Nuclear Information System (INIS)

Sato, Masumi

1985-01-01

Aitken's or Shanks' transformation of the exponent-modified Bessel function produces better approximations. Dispersion relations for the hybrid and Bernstein waves using these provide better thermal and parallel wavenumber corrections. They also predict more closely the evolution and mode-conversion of these waves. (author)

Black-hole quasinormal resonances: Wave analysis versus a geometric-optics approximation

International Nuclear Information System (INIS)

Hod, Shahar

2009-01-01

It has long been known that null unstable geodesics are related to the characteristic modes of black holes--the so-called quasinormal resonances. The basic idea is to interpret the free oscillations of a black hole in the eikonal limit in terms of null particles trapped at the unstable circular orbit and slowly leaking out. The real part of the complex quasinormal resonances is related to the angular velocity at the unstable null geodesic. The imaginary part of the resonances is related to the instability time scale (or the inverse Lyapunov exponent) of the orbit. While this geometric-optics description of the black-hole quasinormal resonances in terms of perturbed null rays is very appealing and intuitive, it is still highly important to verify the validity of this approach by directly analyzing the Teukolsky wave equation which governs the dynamics of perturbation waves in the black-hole spacetime. This is the main goal of the present paper. We first use the geometric-optics technique of perturbing a bundle of unstable null rays to calculate the resonances of near-extremal Kerr black holes in the eikonal approximation. We then directly solve the Teukolsky wave equation (supplemented by the appropriate physical boundary conditions) and show that the resultant quasinormal spectrum obtained directly from the wave analysis is in accord with the spectrum obtained from the geometric-optics approximation of perturbed null rays.

Influence of wetting-layer wave functions on phonon-mediated carrier capture into self-assembled quantum dots

DEFF Research Database (Denmark)

Markussen, Troels; Kristensen, Philip Trøst; Tromborg, Bjarne

2006-01-01

Models of carrier dynamics in quantum dots rely strongly on adequate descriptions of the carrier wave functions. In this work we numerically solve the one-band effective mass Schrodinger equation to calculate the capture times of phonon-mediated carrier capture into self-assembled quantum dots. C....... Comparing with results obtained using approximate carrier wave functions, we demonstrate that the capture times are strongly influenced by properties of the wetting layer wave functions not accounted for by earlier theoretical analyses....

Expansions for Coulomb wave functions

NARCIS (Netherlands)

Boersma, J.

1969-01-01

In this paper we derive a number of expansions for Whittaker functions, regular and irregular Coulomb wave functions. The main result consists of a new expansion for the irregular Coulomb wave functions of orders zero and one in terms of regular Coulomb wave functions. The latter expansions are

Properties of resonance wave functions.

Science.gov (United States)

More, R. M.; Gerjuoy, E.

1973-01-01

Construction and study of resonance wave functions corresponding to poles of the Green's function for several illustrative models of theoretical interest. Resonance wave functions obtained from the Siegert and Kapur-Peierls definitions of the resonance energies are compared. The comparison especially clarifies the meaning of the normalization constant of the resonance wave functions. It is shown that the wave functions may be considered renormalized in a sense analogous to that of quantum field theory. However, this renormalization is entirely automatic, and the theory has neither ad hoc procedures nor infinite quantities.

Quasi-particle excitations and dynamical structure function of trapped Bose-condensates in the WKB approximation

OpenAIRE

Csordás, András; Graham, Robert; Szépfalusy, Péter

1997-01-01

The Bogoliubov equations of the quasi-particle excitations in a weakly interacting trapped Bose-condensate are solved in the WKB approximation in an isotropic harmonic trap, determining the discrete quasi-particle energies and wave functions by torus (Bohr-Sommerfeld) quantization of the integrable classical quasi-particle dynamics. The results are used to calculate the position and strengths of the peaks in the dynamic structure function which can be observed by off-resonance inelastic light...

Classical representation of wave functions for integrable systems

International Nuclear Information System (INIS)

Kay, Kenneth G.

2004-01-01

Classical exact (CE) wave functions are certain integral representations of energy eigenfunctions that are parameterized in terms of the motion of the corresponding classical system in a semiclassically relevant way. When applied to systems for which they are not exact, such expressions serve as semiclassical approximations. Previous work identified CE wave functions for a number of specific systems and established their semiclassical usefulness. This paper explores the degree to which such representations can be found for more general systems. It is shown that CE wave functions exist, in principle, for bound states of an arbitrary integrable system that are confined to a single classically allowed region. Evidence is presented that CE representations also exist for more general states of such a system that are unbound, or that extend over more than one allowed region. The CE expressions are not unique: an innumerable variety exists for each such system. The existence proof provides a formal method for constructing CE expressions by Fourier transforming certain superpositions of energy eigenstates. The parameterization in terms of the classical motion is achieved by identifying certain quantities in these superpositions as classical action and angle variables. The semiclassical relevance of this identification is ensured by imposing some mild conditions on the coefficients in the superposition. This procedure for parameterizing exact wave functions in terms of classical variables indicates a basic relationship between the quantum and classical descriptions of states. The method of constructing CE wave functions introduced in the proof is shown to be consistent with a number of previously obtained CE formulas and is used to derive two new, closed-form, CE expressions. A simple numerical example is presented to illustrate the semiclassical application of one of these expressions and to further verify the physical significance of the classical parameterization

Wave function continuity and the diagonal Born-Oppenheimer correction at conical intersections.

Science.gov (United States)

Meek, Garrett A; Levine, Benjamin G

2016-05-14

We demonstrate that though exact in principle, the expansion of the total molecular wave function as a sum over adiabatic Born-Oppenheimer (BO) vibronic states makes inclusion of the second-derivative nonadiabatic energy term near conical intersections practically problematic. In order to construct a well-behaved molecular wave function that has density at a conical intersection, the individual BO vibronic states in the summation must be discontinuous. When the second-derivative nonadiabatic terms are added to the Hamiltonian, singularities in the diagonal BO corrections (DBOCs) of the individual BO states arise from these discontinuities. In contrast to the well-known singularities in the first-derivative couplings at conical intersections, these singularities are non-integrable, resulting in undefined DBOC matrix elements. Though these singularities suggest that the exact molecular wave function may not have density at the conical intersection point, there is no physical basis for this constraint. Instead, the singularities are artifacts of the chosen basis of discontinuous functions. We also demonstrate that continuity of the total molecular wave function does not require continuity of the individual adiabatic nuclear wave functions. We classify nonadiabatic molecular dynamics methods according to the constraints placed on wave function continuity and analyze their formal properties. Based on our analysis, it is recommended that the DBOC be neglected when employing mixed quantum-classical methods and certain approximate quantum dynamical methods in the adiabatic representation.

Micrononcasual Euclidean wave functions

International Nuclear Information System (INIS)

Enatsu, H.; Takenaka, A.; Okazaki, M.

1978-01-01

A theory which describes the internal attributes of hadrons in terms of space-time wave functions is presented. In order to develop the theory on the basis of a rather realistic model, covariant wave equations are first derived for the deuteron, in which the co-ordinates of the centre of mass of two nucleons can be defined unambiguously. Then the micro-noncasual behaviour of virtual mesons mediating between the two nucleons is expressed by means of wave functions depending only on the relative Euclidean co-ordinates with respect to the centre of mass of the two nucleons; the wave functions are assumed to obey the 0 4 and SU 2 x SU 2 groups. The properties of the wave functions under space inversion, time reversal and particle-antiparticle conjugation are investigated. It is found that the internal attributes of the mesons, such as spin, isospin, strangeness, intrinsic parity, charge parity and G-parity are explained consistently. The theory is applicable also to the case of baryons

On the asymptotic evolution of finite energy Airy wave functions.

Science.gov (United States)

Chamorro-Posada, P; Sánchez-Curto, J; Aceves, A B; McDonald, G S

2015-06-15

In general, there is an inverse relation between the degree of localization of a wave function of a certain class and its transform representation dictated by the scaling property of the Fourier transform. We report that in the case of finite energy Airy wave packets a simultaneous increase in their localization in the direct and transform domains can be obtained as the apodization parameter is varied. One consequence of this is that the far-field diffraction rate of a finite energy Airy beam decreases as the beam localization at the launch plane increases. We analyze the asymptotic properties of finite energy Airy wave functions using the stationary phase method. We obtain one dominant contribution to the long-term evolution that admits a Gaussian-like approximation, which displays the expected reduction of its broadening rate as the input localization is increased.

Probing α-particle wave functions using (rvec d,α) reactions

International Nuclear Information System (INIS)

Crosson, E.R.; Lemieux, S.K.; Ludwig, E.J.; Thompson, W.J.; Bisenberger, M.; Hertenberger, R.; Hofer, D.; Kader, H.; Schiemenz, P.; Graw, G.; Eiro, A.M.; Santos, F.D.

1993-01-01

Wave functions of the α particle corresponding to different S- and D-state deuteron-deuteron overlaps, left-angle dd|α right-angle, were investigated using exact finite-range distorted-wave Born-approximation (DWBA) analyses of (rvec d,α) reactions. Cross sections, vector, and tensor-analyzing powers were measured for (rvec d,α) reactions populating the lowest J π =7 + state in 56 Co at bombarding energies E d of 16 and 22 MeV, the lowest 7 + state in 48 Sc at E d =16 MeV, and the lowest 7 + state in 46 Sc at E d =22 MeV. We find that DWBA analyses of tensor-analyzing powers produce satisfactory agreement with the data and that A xx is especially sensitive to the D-state component of α-particle wave functions generated by different realistic nucleon-nucleon interactions

S-wave pion-nucleon phase shifts in PADE approximation

International Nuclear Information System (INIS)

Achuthan, P.; Chandramohan, T.; Venkatesan, K.

1978-01-01

The two S-wave pion nucleon pahse shifts delta( 1 ) (I = 1/2) and delta( 3 ) (I = 3/2) have been calculated in the Pade approximation using epsilon(700), rho(770), f(1260), Δ(1236) and N(938) for the energy range W = 1085 MeV - 1820 MeV in the centre of mass. Contributions from suitable resonance combinations which agree nearest with the delta( 3 ) experimental values are given. (orig.) [de

Wigner functions of s waves

International Nuclear Information System (INIS)

Dahl, J. P.; Varro, S.; Wolf, A.; Schleich, W. P.

2007-01-01

We derive explicit expressions for the Wigner function of wave functions in D dimensions which depend on the hyperradius--that is, of s waves. They are based either on the position or the momentum representation of the s wave. The corresponding Wigner function depends on three variables: the absolute value of the D-dimensional position and momentum vectors and the angle between them. We illustrate these expressions by calculating and discussing the Wigner functions of an elementary s wave and the energy eigenfunction of a free particle

Wigner functions of s waves

DEFF Research Database (Denmark)

Dahl, Jens Peder; Varro, S.; Wolf, A.

2007-01-01

We derive explicit expressions for the Wigner function of wave functions in D dimensions which depend on the hyperradius-that is, of s waves. They are based either on the position or the momentum representation of the s wave. The corresponding Wigner function depends on three variables......: the absolute value of the D-dimensional position and momentum vectors and the angle between them. We illustrate these expressions by calculating and discussing the Wigner functions of an elementary s wave and the energy eigenfunction of a free particle....

Does the instantaneous wave-free ratio approximate the fractional flow reserve?

NARCIS (Netherlands)

Johnson, Nils P.; Kirkeeide, Richard L.; Asrress, Kaleab N.; Fearon, William F.; Lockie, Timothy; Marques, Koen M. J.; Pyxaras, Stylianos A.; Rolandi, M. Cristina; van 't Veer, Marcel; de Bruyne, Bernard; Piek, Jan J.; Pijls, Nico H. J.; Redwood, Simon; Siebes, Maria; Spaan, Jos A. E.; Gould, K. Lance

2013-01-01

This study sought to examine the clinical performance of and theoretical basis for the instantaneous wave-free ratio (iFR) approximation to the fractional flow reserve (FFR). Recent work has proposed iFR as a vasodilation-free alternative to FFR for making mechanical revascularization decisions. Its

Short-distance behavior of the Bethe--Salpeter wave function in the ladder approximation

International Nuclear Information System (INIS)

Guth, A.H.; Soper, D.E.

1975-01-01

We investigate the short-distance behavior of the (Wick-rotated) Bethe--Salpeter wave function for the two spin-1/2 quarks bound by the exchange of a massive vector meson. We use the ladder-model kernel, which has the same p -4 scaling behavior as the true kernel in a theory with a fixed point of the renormalization group at g not equal to 0. For a bound state with the quantum numbers of the pion, the leading asymptotic behavior is chi (q/sup μ/) approx. cq/sup -4 + epsilon(g)/γ 5 , where epsilon (g) =1- (1-g 2 /π 2 ) 1 / 2 . Our method also provides the full asymptotic series, although it should be noted that the nonleading terms will depend on the nonleading behavior of the ladder-model kernel. A general term has the form cq - /sup a/(lnq)/sup n/phi (q/sup μ/), where c is an unknown constant, a may be integral or nonintegral, n is an integer, and phi (q/sup μ/) is a representation function of the rotation group in four dimensions

Function approximation using combined unsupervised and supervised learning.

Science.gov (United States)

Andras, Peter

2014-03-01

Function approximation is one of the core tasks that are solved using neural networks in the context of many engineering problems. However, good approximation results need good sampling of the data space, which usually requires exponentially increasing volume of data as the dimensionality of the data increases. At the same time, often the high-dimensional data is arranged around a much lower dimensional manifold. Here we propose the breaking of the function approximation task for high-dimensional data into two steps: (1) the mapping of the high-dimensional data onto a lower dimensional space corresponding to the manifold on which the data resides and (2) the approximation of the function using the mapped lower dimensional data. We use over-complete self-organizing maps (SOMs) for the mapping through unsupervised learning, and single hidden layer neural networks for the function approximation through supervised learning. We also extend the two-step procedure by considering support vector machines and Bayesian SOMs for the determination of the best parameters for the nonlinear neurons in the hidden layer of the neural networks used for the function approximation. We compare the approximation performance of the proposed neural networks using a set of functions and show that indeed the neural networks using combined unsupervised and supervised learning outperform in most cases the neural networks that learn the function approximation using the original high-dimensional data.

Legendre-tau approximations for functional differential equations

Science.gov (United States)

Ito, K.; Teglas, R.

1986-01-01

The numerical approximation of solutions to linear retarded functional differential equations are considered using the so-called Legendre-tau method. The functional differential equation is first reformulated as a partial differential equation with a nonlocal boundary condition involving time-differentiation. The approximate solution is then represented as a truncated Legendre series with time-varying coefficients which satisfy a certain system of ordinary differential equations. The method is very easy to code and yields very accurate approximations. Convergence is established, various numerical examples are presented, and comparison between the latter and cubic spline approximation is made.

Variation in Differential and Total Cross Sections Due to Different Radial Wave Functions

Science.gov (United States)

Williamson, W., Jr.; Greene, T.

1976-01-01

Three sets of analytical wave functions are used to calculate the Na (3s---3p) transition differential and total electron excitation cross sections by Born approximations. Results show expected large variations in values. (Author/CP)

Orbitally invariant internally contracted multireference unitary coupled cluster theory and its perturbative approximation: theory and test calculations of second order approximation.

Science.gov (United States)

Chen, Zhenhua; Hoffmann, Mark R

2012-07-07

A unitary wave operator, exp (G), G(+) = -G, is considered to transform a multiconfigurational reference wave function Φ to the potentially exact, within basis set limit, wave function Ψ = exp (G)Φ. To obtain a useful approximation, the Hausdorff expansion of the similarity transformed effective Hamiltonian, exp (-G)Hexp (G), is truncated at second order and the excitation manifold is limited; an additional separate perturbation approximation can also be made. In the perturbation approximation, which we refer to as multireference unitary second-order perturbation theory (MRUPT2), the Hamiltonian operator in the highest order commutator is approximated by a Mo̸ller-Plesset-type one-body zero-order Hamiltonian. If a complete active space self-consistent field wave function is used as reference, then the energy is invariant under orbital rotations within the inactive, active, and virtual orbital subspaces for both the second-order unitary coupled cluster method and its perturbative approximation. Furthermore, the redundancies of the excitation operators are addressed in a novel way, which is potentially more efficient compared to the usual full diagonalization of the metric of the excited configurations. Despite the loss of rigorous size-extensivity possibly due to the use of a variational approach rather than a projective one in the solution of the amplitudes, test calculations show that the size-extensivity errors are very small. Compared to other internally contracted multireference perturbation theories, MRUPT2 only needs reduced density matrices up to three-body even with a non-complete active space reference wave function when two-body excitations within the active orbital subspace are involved in the wave operator, exp (G). Both the coupled cluster and perturbation theory variants are amenable to large, incomplete model spaces. Applications to some widely studied model systems that can be problematic because of geometry dependent quasidegeneracy, H4, P4

Model wave functions for the deuteron

International Nuclear Information System (INIS)

Certov, A.; Mathelitsch, L.; Moravcsik, M.J.

1987-01-01

Model wave functions are constructed for the deuteron to facilitate the unambiguous exploration of dependencies on the percentage D state and on the small-, medium-, and large-distance parts of the deuteron wave function. The wave functions are constrained by those deuteron properties which are accurately known experimentally, and are in an analytic form which is easily integrable in expressions usually encountered in the use of such wave functions

Wave-function functionals for the density

International Nuclear Information System (INIS)

Slamet, Marlina; Pan Xiaoyin; Sahni, Viraht

2011-01-01

We extend the idea of the constrained-search variational method for the construction of wave-function functionals ψ[χ] of functions χ. The search is constrained to those functions χ such that ψ[χ] reproduces the density ρ(r) while simultaneously leading to an upper bound to the energy. The functionals are thereby normalized and automatically satisfy the electron-nucleus coalescence condition. The functionals ψ[χ] are also constructed to satisfy the electron-electron coalescence condition. The method is applied to the ground state of the helium atom to construct functionals ψ[χ] that reproduce the density as given by the Kinoshita correlated wave function. The expectation of single-particle operators W=Σ i r i n , n=-2,-1,1,2, W=Σ i δ(r i ) are exact, as must be the case. The expectations of the kinetic energy operator W=-(1/2)Σ i ∇ i 2 , the two-particle operators W=Σ n u n , n=-2,-1,1,2, where u=|r i -r j |, and the energy are accurate. We note that the construction of such functionals ψ[χ] is an application of the Levy-Lieb constrained-search definition of density functional theory. It is thereby possible to rigorously determine which functional ψ[χ] is closer to the true wave function.

Two-body Schrödinger wave functions in a plane-wave basis via separation of dimensions

Science.gov (United States)

Jerke, Jonathan; Poirier, Bill

2018-03-01

Using a combination of ideas, the ground and several excited electronic states of the helium atom and the hydrogen molecule are computed to chemical accuracy—i.e., to within 1-2 mhartree or better. The basic strategy is very different from the standard electronic structure approach in that the full two-electron six-dimensional (6D) problem is tackled directly, rather than starting from a single-electron Hartree-Fock approximation. Electron correlation is thus treated exactly, even though computational requirements remain modest. The method also allows for exact wave functions to be computed, as well as energy levels. From the full-dimensional 6D wave functions computed here, radial distribution functions and radial correlation functions are extracted—as well as a 2D probability density function exhibiting antisymmetry for a single Cartesian component. These calculations support a more recent interpretation of Hund's rule, which states that the lower energy of the higher spin-multiplicity states is actually due to reduced screening, rather than reduced electron-electron repulsion. Prospects for larger systems and/or electron dynamics applications appear promising.

Microscopy of electronic wave function

International Nuclear Information System (INIS)

Harb, M.

2010-01-01

This work of thesis aims to visualize, on a position sensitive detector, the spatial oscillations of slow electrons (∼ meV) emitted by a threshold photoionization in the presence of an external electric field. The interference figure obtained represents the square magnitude of electronic wavefunction. This fundamental work allows us to have access to the electronic dynamics and thus to highlight several quantum mechanisms that occur at the atomic scale (field Coulomb, electron/electron interaction..). Despite the presence an electronic core in Li atom, we have succeeded, experimentally and for the first time, in visualizing the wave function associated with the quasi-discrete Stark states coupled to the ionization continuum. Besides, using simulations of wave packet propagation, based on the 'Split-operator' method, we have conducted a comprehensive study of the H, Li and Cs atoms while revealing the significant effects of the Stark resonances. A very good agreement, on and off resonances, was obtained between simulated and experimental results. In addition, we have developed a generalized analytical model to understand deeply the function of VMI (Velocity-Map Imaging) spectrometer. This model is based on the paraxial approximation; it is based on matrix optics calculation by making an analogy between the electronic trajectory and the light beam. An excellent agreement was obtained between the model predictions and the experimental results. (author)

A partition function approximation using elementary symmetric functions.

Directory of Open Access Journals (Sweden)

Ramu Anandakrishnan

Full Text Available In statistical mechanics, the canonical partition function [Formula: see text] can be used to compute equilibrium properties of a physical system. Calculating [Formula: see text] however, is in general computationally intractable, since the computation scales exponentially with the number of particles [Formula: see text] in the system. A commonly used method for approximating equilibrium properties, is the Monte Carlo (MC method. For some problems the MC method converges slowly, requiring a very large number of MC steps. For such problems the computational cost of the Monte Carlo method can be prohibitive. Presented here is a deterministic algorithm - the direct interaction algorithm (DIA - for approximating the canonical partition function [Formula: see text] in [Formula: see text] operations. The DIA approximates the partition function as a combinatorial sum of products known as elementary symmetric functions (ESFs, which can be computed in [Formula: see text] operations. The DIA was used to compute equilibrium properties for the isotropic 2D Ising model, and the accuracy of the DIA was compared to that of the basic Metropolis Monte Carlo method. Our results show that the DIA may be a practical alternative for some problems where the Monte Carlo method converge slowly, and computational speed is a critical constraint, such as for very large systems or web-based applications.

Photo double ionization of He: C3-like wave function for the two electron continuum

Energy Technology Data Exchange (ETDEWEB)

Otranto, S.; Garibotti, C.R. [Conicet and Centro Atomico Bariloche (Argentina); Otranto, S. [Universidad Nacional del Sur, Dept. de Fisica, Bahia Blanca (Argentina)

2002-12-01

We evaluate the triply differential cross-section (TDCS) for photo double ionization (PDI) of helium. A first approximation to the final state can be obtained by neglecting the e-e interaction and the non-orthogonal kinetic energy. This leads to the C2 model which proposes as solution a product of 2 independent Coulomb wave plane waves. A better approximation is the C3 model where the C3 wave describes the e-e motion as independent of the presence of the nucleus and represents it by a Coulomb continuum wave. The C3 wave function mainly consists in the product of 3 Coulomb waves, each one representing the interaction between a pair of particles. We use a C3 final continuum wave function with an inter-electronic effective coordinate to express the nuclear screening. Comparison with the standard C3 model shows that the TDCS is enhanced in the threshold region by effect of the reduced inter-electronic repulsion introduced by the present model. A more accurate description of the intermediate energy region is also obtained. Comparison with recent experimental data shows a good overall agreement of the angular distributions. The theoretical PDI total cross-section shows a relevant improvement in the intermediate energy region relative to the C3 model, which converges to data for photon energies larger than 1 keV.

Photo double ionization of He: C3-like wave function for the two electron continuum

International Nuclear Information System (INIS)

Otranto, S.; Garibotti, C.R.; Otranto, S.

2002-01-01

We evaluate the triply differential cross-section (TDCS) for photo double ionization (PDI) of helium. A first approximation to the final state can be obtained by neglecting the e-e interaction and the non-orthogonal kinetic energy. This leads to the C2 model which proposes as solution a product of 2 independent Coulomb wave plane waves. A better approximation is the C3 model where the C3 wave describes the e-e motion as independent of the presence of the nucleus and represents it by a Coulomb continuum wave. The C3 wave function mainly consists in the product of 3 Coulomb waves, each one representing the interaction between a pair of particles. We use a C3 final continuum wave function with an inter-electronic effective coordinate to express the nuclear screening. Comparison with the standard C3 model shows that the TDCS is enhanced in the threshold region by effect of the reduced inter-electronic repulsion introduced by the present model. A more accurate description of the intermediate energy region is also obtained. Comparison with recent experimental data shows a good overall agreement of the angular distributions. The theoretical PDI total cross-section shows a relevant improvement in the intermediate energy region relative to the C3 model, which converges to data for photon energies larger than 1 keV

Twist-2 Light-Cone Pion Wave Function

OpenAIRE

Belyaev, V. M.; Johnson, Mikkel B.

1997-01-01

We present an analysis of the existing constraints for the twist-2 light-cone pion wave function. We find that existing information on the pion wave function does not exclude the possibility that the pion wave function attains its asymptotic form. New bounds on the parameters of the pion wave function are presented.

Does the Berry phase in a quantum optical system originate from the rotating wave approximation?

International Nuclear Information System (INIS)

Wang, Minghao; Wei, L.F.; Liang, J.Q.

2015-01-01

The Berry phase (BP) in a quantized light field demonstrated more than a decade ago (Fuentes-Guridi et al., 2002 [9]) has attracted considerable attention, since it plays an important role in the cavity quantum electrodynamics. However, it is argued in Larson (2012) [15] that such a BP is just due to the rotating wave approximation (RWA) and the relevant BP should vanish beyond this approximation. Based on a consistent analysis we conclude in this letter that the BP in a generic Rabi model actually exists, no matter whether the RWA is applied. The existence of BP is also generalized to a three-level atom in the quantized cavity field. - Highlights: • Non-zero Berry phases for the Rabi model (without rotating wave approximation) are verified. • A general formulation of Berry phases for both the JC model and the Rabi model is presented. • The claim of vanishing Berry phase in the Rabi model is a result of improper semiclassical approximation. • Analytic solutions for the Rabi model is presented in the semiclassical approximation

Entropy vs. Action in the (2+1)-Dimensional Hartle-Hawking Wave Function

OpenAIRE

Carlip, Steven

1992-01-01

In most attempts to compute the Hartle-Hawking ``wave function of the universe'' in Euclidean quantum gravity, two important approximations are made: the path integral is evaluated in a saddle point approximation, and only the leading (least action) extremum is taken into account. In (2+1)-dimensional gravity with a negative cosmological constant, the second assumption is shown to lead to incorrect results: although the leading extremum gives the most important single contribution to the path...

Bridging the gap between the Jaynes–Cummings and Rabi models using an intermediate rotating wave approximation

International Nuclear Information System (INIS)

Wang, Yimin; Haw, Jing Yan

2015-01-01

Highlights: • The intermediate rotating wave approximation (IRWA) is introduced. • The co-rotating and counter-rotating terms in Rabi model are expressed separately. • IRWA is applied to the near resonance case and the large detuning case. • The continuity between the Jaynes–Cummings model and the Rabi model is established. - Abstract: We present a novel approach called the intermediate rotating wave approximation (IRWA), which employs a time-averaging method to encapsulate the dynamics of light-matter interaction from strong to ultrastrong coupling regime. In contrast to the ordinary rotating wave approximation, this method addresses the co-rotating and counter-rotating terms separately to trace their physical consequences individually, and thus establishes the continuity between the Jaynes–Cummings model and the quantum Rabi model. We investigate IRWA in near resonance and large detuning cases. Our IRWA not only agrees well with both models in their respective coupling strengths, but also offers a good explanation of their differences

Digital-Control-Based Approximation of Optimal Wave Disturbances Attenuation for Nonlinear Offshore Platforms

Directory of Open Access Journals (Sweden)

Xiao-Fang Zhong

2017-12-01

Full Text Available The irregular wave disturbance attenuation problem for jacket-type offshore platforms involving the nonlinear characteristics is studied. The main contribution is that a digital-control-based approximation of optimal wave disturbances attenuation controller (AOWDAC is proposed based on iteration control theory, which consists of a feedback item of offshore state, a feedforward item of wave force and a nonlinear compensated component with iterative sequences. More specifically, by discussing the discrete model of nonlinear offshore platform subject to wave forces generated from the Joint North Sea Wave Project (JONSWAP wave spectrum and linearized wave theory, the original wave disturbances attenuation problem is formulated as the nonlinear two-point-boundary-value (TPBV problem. By introducing two vector sequences of system states and nonlinear compensated item, the solution of introduced nonlinear TPBV problem is obtained. Then, a numerical algorithm is designed to realize the feasibility of AOWDAC based on the deviation of performance index between the adjacent iteration processes. Finally, applied the proposed AOWDAC to a jacket-type offshore platform in Bohai Bay, the vibration amplitudes of the displacement and the velocity, and the required energy consumption can be reduced significantly.

Heuristic method for determining outgoing waves in many-body wave functions

International Nuclear Information System (INIS)

Redish, E.F.; Tandy, P.C.; L'Huillier, M.

1975-12-01

A new and simple method is proposed for determining the kinds of outgoing waves present in a given many-body wave function. Whether any particular wave function contains ''hidden'' rearrangement components can be determined. 1 figure

Approximate Treatment of the Dirac Equation with Hyperbolic Potential Function

Science.gov (United States)

Durmus, Aysen

2018-03-01

The time independent Dirac equation is solved analytically for equal scalar and vector hyperbolic potential function in the presence of Greene and Aldrich approximation scheme. The bound state energy equation and spinor wave functions expressed by the hypergeometric function have been obtained in detail with asymptotic iteration approach. In order to indicate the accuracy of this different approach proposed to solve second order linear differential equations, we present that in the non-relativistic limit, analytical solutions of the Dirac equation converge to those of the Schrödinger one. We introduce numerical results of the theoretical analysis for hyperbolic potential function. Bound states corresponding to arbitrary values of n and l are reported for potential parameters covering a wide range of interaction. Also, we investigate relativistic vibrational energy spectra of alkali metal diatomic molecules in the different electronic states. It is observed that theoretical vibrational energy values are consistent with experimental Rydberg-Klein-Rees (RKR) results and vibrational energies of NaK, K_2 and KRb diatomic molecules interacting with hyperbolic potential smoothly converge to the experimental dissociation limit D_e=2508cm^{-1}, 254cm^{-1} and 4221cm^{-1}, respectively.

Discovery of functional and approximate functional dependencies in relational databases

Directory of Open Access Journals (Sweden)

Ronald S. King

2003-01-01

Full Text Available This study develops the foundation for a simple, yet efficient method for uncovering functional and approximate functional dependencies in relational databases. The technique is based upon the mathematical theory of partitions defined over a relation's row identifiers. Using a levelwise algorithm the minimal non-trivial functional dependencies can be found using computations conducted on integers. Therefore, the required operations on partitions are both simple and fast. Additionally, the row identifiers provide the added advantage of nominally identifying the exceptions to approximate functional dependencies, which can be used effectively in practical data mining applications.

Quantum Chemistry on Quantum Computers: A Polynomial-Time Quantum Algorithm for Constructing the Wave Functions of Open-Shell Molecules.

Science.gov (United States)

Sugisaki, Kenji; Yamamoto, Satoru; Nakazawa, Shigeaki; Toyota, Kazuo; Sato, Kazunobu; Shiomi, Daisuke; Takui, Takeji

2016-08-18

Quantum computers are capable to efficiently perform full configuration interaction (FCI) calculations of atoms and molecules by using the quantum phase estimation (QPE) algorithm. Because the success probability of the QPE depends on the overlap between approximate and exact wave functions, efficient methods to prepare accurate initial guess wave functions enough to have sufficiently large overlap with the exact ones are highly desired. Here, we propose a quantum algorithm to construct the wave function consisting of one configuration state function, which is suitable for the initial guess wave function in QPE-based FCI calculations of open-shell molecules, based on the addition theorem of angular momentum. The proposed quantum algorithm enables us to prepare the wave function consisting of an exponential number of Slater determinants only by a polynomial number of quantum operations.

Smooth function approximation using neural networks.

Science.gov (United States)

Ferrari, Silvia; Stengel, Robert F

2005-01-01

An algebraic approach for representing multidimensional nonlinear functions by feedforward neural networks is presented. In this paper, the approach is implemented for the approximation of smooth batch data containing the function's input, output, and possibly, gradient information. The training set is associated to the network adjustable parameters by nonlinear weight equations. The cascade structure of these equations reveals that they can be treated as sets of linear systems. Hence, the training process and the network approximation properties can be investigated via linear algebra. Four algorithms are developed to achieve exact or approximate matching of input-output and/or gradient-based training sets. Their application to the design of forward and feedback neurocontrollers shows that algebraic training is characterized by faster execution speeds and better generalization properties than contemporary optimization techniques.

Quasi-fractional approximation to the Bessel functions

International Nuclear Information System (INIS)

Guerrero, P.M.L.

1989-01-01

In this paper the authors presents a simple Quasi-Fractional Approximation for Bessel Functions J ν (x), (- 1 ≤ ν < 0.5). This has been obtained by extending a method published which uses simultaneously power series and asymptotic expansions. Both functions, exact and approximated, coincide in at least two digits for positive x, and ν between - 1 and 0,4

Coordinate asymptotics of the (3→3) wave functions for a three charged particle system

International Nuclear Information System (INIS)

Merkur'ev, S.P.

1977-01-01

Coordinate asymptotics of the (3 → 3) wave functions for three particles system with Coulomb interaction in the scattering problem is plotted. (3 → 3) and (3 → 2) process cases are considered, when the particles are not connected at the initial state. For coordinate asymptotics plotting the basis functions are used which meet Schroedinger equation in the eikonal approximation. The wave functions coordinate asymptotics plotting method is described far from special directions. Wave function asymptotical form is studied in the range of special directions and (3 → 3) scattering amplitude singularities are described. All data are given in accordance with the system with 2 charged particles only. The model in question is of special interest because of the described ppn system the studying of which is of great importance in nuclear physics. Final formulae are discussed for the most general case of three charged particles. Boundary problems for Schroedinger equation are shown to give the only way of definition for the (3 → 3) wave functions. It is pointed out that in special directions wave function coordinate asymptotics is presented with accuracy that gives the possibility to set such a boundary problem

Nonlinear Ritz approximation for Fredholm functionals

Directory of Open Access Journals (Sweden)

Mudhir A. Abdul Hussain

2015-11-01

Full Text Available In this article we use the modify Lyapunov-Schmidt reduction to find nonlinear Ritz approximation for a Fredholm functional. This functional corresponds to a nonlinear Fredholm operator defined by a nonlinear fourth-order differential equation.

Symmetrized partial-wave method for density-functional cluster calculations

International Nuclear Information System (INIS)

Averill, F.W.; Painter, G.S.

1994-01-01

The computational advantage and accuracy of the Harris method is linked to the simplicity and adequacy of the reference-density model. In an earlier paper, we investigated one way the Harris functional could be extended to systems outside the limits of weakly interacting atoms by making the charge density of the interacting atoms self-consistent within the constraints of overlapping spherical atomic densities. In the present study, a method is presented for augmenting the interacting atom charge densities with symmetrized partial-wave expansions on each atomic site. The added variational freedom of the partial waves leads to a scheme capable of giving exact results within a given exchange-correlation approximation while maintaining many of the desirable convergence and stability properties of the original Harris method. Incorporation of the symmetry of the cluster in the partial-wave construction further reduces the level of computational effort. This partial-wave cluster method is illustrated by its application to the dimer C 2 , the hypothetical atomic cluster Fe 6 Al 8 , and the benzene molecule

An analytic distorted wave approximation for intermediate energy proton scattering

International Nuclear Information System (INIS)

Di Marzio, F.; Amos, K.

1982-01-01

An analytic Distorted Wave approximation has been developed for use in analyses of intermediate energy proton inelastic scattering from nuclei. Applications are made to analyse 402 and 800 MeV data from the isoscalar and isovector 1 + and 2 + states in 12 C and to the 800 MeV data from the excitation of the 2 - (8.88MeV) state in 16 O. Comparisons of predictions made using different model two-nucleon t-matrices and different models of nuclear structure are given

Correction of the near threshold behavior of electron collisional excitation cross-sections in the plane-wave Born approximation

Science.gov (United States)

Kilcrease, D. P.; Brookes, S.

2013-12-01

The modeling of NLTE plasmas requires the solution of population rate equations to determine the populations of the various atomic levels relevant to a particular problem. The equations require many cross sections for excitation, de-excitation, ionization and recombination. A simple and computational fast way to calculate electron collisional excitation cross-sections for ions is by using the plane-wave Born approximation. This is essentially a high-energy approximation and the cross section suffers from the unphysical problem of going to zero near threshold. Various remedies for this problem have been employed with varying degrees of success. We present a correction procedure for the Born cross-sections that employs the Elwert-Sommerfeld factor to correct for the use of plane waves instead of Coulomb waves in an attempt to produce a cross-section similar to that from using the more time consuming Coulomb Born approximation. We compare this new approximation with other, often employed correction procedures. We also look at some further modifications to our Born Elwert procedure and its combination with Y.K. Kim's correction of the Coulomb Born approximation for singly charged ions that more accurately approximate convergent close coupling calculations.

Capillary waves at the interface of two Bose–Einstein condensates. Long wavelengths asymptotic by trial function approach

International Nuclear Information System (INIS)

Mishonov, T.M.

2015-01-01

The dispersion relation for capillary waves at the boundary of two different Bose condensates is investigated using a trial wave-function approach applied to the Gross-Pitaevskii (GP) equations. The surface tension is expressed by the parameters of the GP equations. In the long wave-length limit the usual dispersion relation is re-derived while for wavelengths comparable to the healing length we predict significant deviations from the ω ∝ k 3/2 law which can be experimentally observed. We approximate the wave variables by a frozen order parameter, i.e. the wave function is frozen in the superfluid analogous to the magnetic field in highly conductive space plasmas. PACS codes: 67.85.Jk

The Wave Function and Quantum Reality

International Nuclear Information System (INIS)

Gao Shan

2011-01-01

We investigate the meaning of the wave function by analyzing the mass and charge density distributions of a quantum system. According to protective measurement, a charged quantum system has effective mass and charge density distributing in space, proportional to the square of the absolute value of its wave function. In a realistic interpretation, the wave function of a quantum system can be taken as a description of either a physical field or the ergodic motion of a particle. The essential difference between a field and the ergodic motion of a particle lies in the property of simultaneity; a field exists throughout space simultaneously, whereas the ergodic motion of a particle exists throughout space in a time-divided way. If the wave function is a physical field, then the mass and charge density will be distributed in space simultaneously for a charged quantum system, and thus there will exist gravitational and electrostatic self-interactions of its wave function. This not only violates the superposition principle of quantum mechanics but also contradicts experimental observations. Thus the wave function cannot be a description of a physical field but be a description of the ergodic motion of a particle. For the later there is only a localized particle with mass and charge at every instant, and thus there will not exist any self-interaction for the wave function. It is further argued that the classical ergodic models, which assume continuous motion of particles, cannot be consistent with quantum mechanics. Based on the negative result, we suggest that the wave function is a description of the quantum motion of particles, which is random and discontinuous in nature. On this interpretation, the square of the absolute value of the wave function not only gives the probability of the particle being found in certain locations, but also gives the probability of the particle being there. The suggested new interpretation of the wave function provides a natural realistic

RATIONAL APPROXIMATIONS TO GENERALIZED HYPERGEOMETRIC FUNCTIONS.

Science.gov (United States)

Under weak restrictions on the various free parameters, general theorems for rational representations of the generalized hypergeometric functions...and certain Meijer G-functions are developed. Upon specialization, these theorems yield a sequency of rational approximations which converge to the

An inductive algorithm for smooth approximation of functions

International Nuclear Information System (INIS)

Kupenova, T.N.

2011-01-01

An inductive algorithm is presented for smooth approximation of functions, based on the Tikhonov regularization method and applied to a specific kind of the Tikhonov parametric functional. The discrepancy principle is used for estimation of the regularization parameter. The principle of heuristic self-organization is applied for assessment of some parameters of the approximating function

Molecular Excitation Energies from Time-Dependent Density Functional Theory Employing Random-Phase Approximation Hessians with Exact Exchange.

Science.gov (United States)

Heßelmann, Andreas

2015-04-14

Molecular excitation energies have been calculated with time-dependent density-functional theory (TDDFT) using random-phase approximation Hessians augmented with exact exchange contributions in various orders. It has been observed that this approach yields fairly accurate local valence excitations if combined with accurate asymptotically corrected exchange-correlation potentials used in the ground-state Kohn-Sham calculations. The inclusion of long-range particle-particle with hole-hole interactions in the kernel leads to errors of 0.14 eV only for the lowest excitations of a selection of three alkene, three carbonyl, and five azabenzene molecules, thus surpassing the accuracy of a number of common TDDFT and even some wave function correlation methods. In the case of long-range charge-transfer excitations, the method typically underestimates accurate reference excitation energies by 8% on average, which is better than with standard hybrid-GGA functionals but worse compared to range-separated functional approximations.

Polynomial approximation of functions in Sobolev spaces

International Nuclear Information System (INIS)

Dupont, T.; Scott, R.

1980-01-01

Constructive proofs and several generalizations of approximation results of J. H. Bramble and S. R. Hilbert are presented. Using an averaged Taylor series, we represent a function as a polynomical plus a remainder. The remainder can be manipulated in many ways to give different types of bounds. Approximation of functions in fractional order Sobolev spaces is treated as well as the usual integer order spaces and several nonstandard Sobolev-like spaces

A local adaptive method for the numerical approximation in seismic wave modelling

Directory of Open Access Journals (Sweden)

Galuzzi Bruno G.

2017-12-01

Full Text Available We propose a new numerical approach for the solution of the 2D acoustic wave equation to model the predicted data in the field of active-source seismic inverse problems. This method consists in using an explicit finite difference technique with an adaptive order of approximation of the spatial derivatives that takes into account the local velocity at the grid nodes. Testing our method to simulate the recorded seismograms in a marine seismic acquisition, we found that the low computational time and the low approximation error of the proposed approach make it suitable in the context of seismic inversion problems.

Risk analysis of breakwater caisson under wave attack using load surface approximation

Science.gov (United States)

Kim, Dong Hyawn

2014-12-01

A new load surface based approach to the reliability analysis of caisson-type breakwater is proposed. Uncertainties of the horizontal and vertical wave loads acting on breakwater are considered by using the so-called load surfaces, which can be estimated as functions of wave height, water level, and so on. Then, the first-order reliability method (FORM) can be applied to determine the probability of failure under the wave action. In this way, the reliability analysis of breakwaters with uncertainties both in wave height and in water level is possible. Moreover, the uncertainty in wave breaking can be taken into account by considering a random variable for wave height ratio which relates the significant wave height to the maximum wave height. The proposed approach is applied numerically to the reliability analysis of caisson breakwater under wave attack that may undergo partial or full wave breaking.

ONETEP: linear-scaling density-functional theory with plane-waves

International Nuclear Information System (INIS)

Haynes, P D; Mostof, A A; Skylaris, C-K; Payne, M C

2006-01-01

This paper provides a general overview of the methodology implemented in onetep (Order-N Electronic Total Energy Package), a parallel density-functional theory code for largescale first-principles quantum-mechanical calculations. The distinctive features of onetep are linear-scaling in both computational effort and resources, obtained by making well-controlled approximations which enable simulations to be performed with plane-wave accuracy. Titanium dioxide clusters of increasing size designed to mimic surfaces are studied to demonstrate the accuracy and scaling of onetep

On root mean square approximation by exponential functions

OpenAIRE

Sharipov, Ruslan

2014-01-01

The problem of root mean square approximation of a square integrable function by finite linear combinations of exponential functions is considered. It is subdivided into linear and nonlinear parts. The linear approximation problem is solved. Then the nonlinear problem is studied in some particular example.

Green function iterative solution of ground state wave function for Yukawa potential

International Nuclear Information System (INIS)

Zhang Zhao

2003-01-01

The newly developed single trajectory quadrature method is applied to solve central potentials. First, based on the series expansion method an exact analytic solution of the ground state for Hulthen potential and an approximate solution for Yukawa potential are obtained respectively. Second, the newly developed iterative method based on Green function defined by quadratures along the single trajectory is applied to solve Yukawa potential using the Coulomb solution and Hulthen solution as the trial functions respectively. The results show that a more proper choice of the trial function will give a better convergence. To further improve the convergence the iterative method is combined with the variational method to solve the ground state wave function for Yukawa potential, using variational solutions of the Coulomb and Hulthen potentials as the trial functions. The results give much better convergence. Finally, the obtained critical screen coefficient is applied to discuss the dissociate temperature of J/ψ in high temperature QGP

Connection of relativistic and nonrelativistic wave functions in the calculation of leptonic widths

International Nuclear Information System (INIS)

Durand, B.; Durand, L.

1984-01-01

We generalize our previous JWKB relations between the relativistic qq-bar wave function at the origin and (a) the inverse density of states of the qq-bar system and (b) the nonrelativistic qq-bar wave function at the origin, to the case of potentials with a Coulomb singularity. We show that the square of the Bethe-Salpeter wave function at the the origin is given approximately for 1 - states by for M/sub n/>2m/sub q/, where F(v) = (4πα/sub s//3v)[1-exp(-4πα /sub s//3v)] -1 is the usual Coulomb factor and g(v)approx. =1 is associated with the lowest-order gluonic radiative corrections. We present numerical evidence for the remarkable accuracy of these relations, which have important implications for the use of nonrelativistic potential models to describe quarkonium systems. We also discuss some subtleties in the v and α/sub s/ dependence of corrections to leptonic widths

Modeling shock waves in an ideal gas: combining the Burnett approximation and Holian's conjecture.

Science.gov (United States)

He, Yi-Guang; Tang, Xiu-Zhang; Pu, Yi-Kang

2008-07-01

We model a shock wave in an ideal gas by combining the Burnett approximation and Holian's conjecture. We use the temperature in the direction of shock propagation rather than the average temperature in the Burnett transport coefficients. The shock wave profiles and shock thickness are compared with other theories. The results are found to agree better with the nonequilibrium molecular dynamics (NEMD) and direct simulation Monte Carlo (DSMC) data than the Burnett equations and the modified Navier-Stokes theory.

Using function approximation to determine neural network accuracy

International Nuclear Information System (INIS)

Wichman, R.F.; Alexander, J.

2013-01-01

Many, if not most, control processes demonstrate nonlinear behavior in some portion of their operating range and the ability of neural networks to model non-linear dynamics makes them very appealing for control. Control of high reliability safety systems, and autonomous control in process or robotic applications, however, require accurate and consistent control and neural networks are only approximators of various functions so their degree of approximation becomes important. In this paper, the factors affecting the ability of a feed-forward back-propagation neural network to accurately approximate a non-linear function are explored. Compared to pattern recognition using a neural network for function approximation provides an easy and accurate method for determining the network's accuracy. In contrast to other techniques, we show that errors arising in function approximation or curve fitting are caused by the neural network itself rather than scatter in the data. A method is proposed that provides improvements in the accuracy achieved during training and resulting ability of the network to generalize after training. Binary input vectors provided a more accurate model than with scalar inputs and retraining using a small number of the outlier x,y pairs improved generalization. (author)

Modeling of pseudoacoustic P-waves in orthorhombic media with a low-rank approximation

KAUST Repository

Song, Xiaolei

2013-06-04

Wavefield extrapolation in pseudoacoustic orthorhombic anisotropic media suffers from wave-mode coupling and stability limitations in the parameter range. We use the dispersion relation for scalar wave propagation in pseudoacoustic orthorhombic media to model acoustic wavefields. The wavenumber-domain application of the Laplacian operator allows us to propagate the P-waves exclusively, without imposing any conditions on the parameter range of stability. It also allows us to avoid dispersion artifacts commonly associated with evaluating the Laplacian operator in space domain using practical finite-difference stencils. To handle the corresponding space-wavenumber mixed-domain operator, we apply the low-rank approximation approach. Considering the number of parameters necessary to describe orthorhombic anisotropy, the low-rank approach yields space-wavenumber decomposition of the extrapolator operator that is dependent on space location regardless of the parameters, a feature necessary for orthorhombic anisotropy. Numerical experiments that the proposed wavefield extrapolator is accurate and practically free of dispersion. Furthermore, there is no coupling of qSv and qP waves because we use the analytical dispersion solution corresponding to the P-wave.

Antiferromagnetism and d-wave superconductivity in (doped) Mott insulators: A wave function approach

OpenAIRE

Weng, Z. Y.; Zhou, Y.; Muthukumar, V. N.

2003-01-01

We propose a class of wave functions that provide a unified description of antiferromagnetism and d-wave superconductivity in (doped) Mott insulators. The wave function has a Jastrow form and prohibits double occupancies. In the absence of holes, the wave function describes antiferromagnetism accurately. Off diagonal long range order develops at finite doping and the superconducting order parameter has d-wave symmetry. We also show how nodal quasiparticles and neutral spin excitations can be ...

Diffusive Wave Approximation to the Shallow Water Equations: Computational Approach

KAUST Repository

Collier, Nathan

2011-05-14

We discuss the use of time adaptivity applied to the one dimensional diffusive wave approximation to the shallow water equations. A simple and computationally economical error estimator is discussed which enables time-step size adaptivity. This robust adaptive time discretization corrects the initial time step size to achieve a user specified bound on the discretization error and allows time step size variations of several orders of magnitude. In particular, in the one dimensional results presented in this work feature a change of four orders of magnitudes for the time step over the entire simulation.

A simple and realistic triton wave function

International Nuclear Information System (INIS)

Lomnitz-Adler, J.; Pandharipande, V.R.

1980-01-01

We propose a simple triton wave function that consists of a product of three correlation operators operating on a three-body spin-isospin state. This wave function is formally similar to that used in the recent variational theories of nuclear matter, the main difference being in the long-range behavior of the correlation operators. Variational calculations are carried out with the Reid potential, using this wave function in the so-called 'symmetrized product' and 'independent pair' forms. The triton energy and density distributions obtained with the symmetrized product wave function agree with those obtained in Faddeev and other variational calculations using harmonic oscillator states. The proposed wave function and calculational methods can be easily generalized to treat the four-nucleon α-particle. (orig.)

Approximating Exponential and Logarithmic Functions Using Polynomial Interpolation

Science.gov (United States)

Gordon, Sheldon P.; Yang, Yajun

2017-01-01

This article takes a closer look at the problem of approximating the exponential and logarithmic functions using polynomials. Either as an alternative to or a precursor to Taylor polynomial approximations at the precalculus level, interpolating polynomials are considered. A measure of error is given and the behaviour of the error function is…

Subsystem density functional theory with meta-generalized gradient approximation exchange-correlation functionals.

Science.gov (United States)

Śmiga, Szymon; Fabiano, Eduardo; Laricchia, Savio; Constantin, Lucian A; Della Sala, Fabio

2015-04-21

We analyze the methodology and the performance of subsystem density functional theory (DFT) with meta-generalized gradient approximation (meta-GGA) exchange-correlation functionals for non-bonded molecular systems. Meta-GGA functionals depend on the Kohn-Sham kinetic energy density (KED), which is not known as an explicit functional of the density. Therefore, they cannot be directly applied in subsystem DFT calculations. We propose a Laplacian-level approximation to the KED which overcomes this limitation and provides a simple and accurate way to apply meta-GGA exchange-correlation functionals in subsystem DFT calculations. The so obtained density and energy errors, with respect to the corresponding supermolecular calculations, are comparable with conventional approaches, depending almost exclusively on the approximations in the non-additive kinetic embedding term. An embedding energy error decomposition explains the accuracy of our method.

Deep inelastic scattering in formalism with wave functions of rest compound system

International Nuclear Information System (INIS)

Sisakyan, A.N.; Kvinikhidze, A.N.; Khvedelidze, A.M.

1987-01-01

One of the most simple examples of interaction of compound systems: deep inelastic scattering of the point particle on hadron is considered. By choosing the compound particle (hadron) rest system the corresponding cross section is expressed in terms of more usual from the view point of nonrelativistic quantum mechanics wave functions of the rest bound state. A new variant of structure functions expansion into a series in terms of the coupling constant is suggested. Each therm of a series due to correct account of the energy conservation law in any order of the perturbation theory possess spectral property. Analysis in QCD shows that in the bound state rest system (P-vector=0) the pulse approximation though satisfies the requirements of scale invariance is insufficient for correct description of elastic limit x Bj →1 by contrast to P Z →∞ system. It means that parton model is equivalent to pulse approximation only in P Z →∞ system. To obtain the leading in asymptotic region x Bj →1 terms account of component interaction in the finite state is necessary. The simplicity and physical evidence of the wave functions are attained due to the seeming complication of calculations according to the perturbation theory

Electron Correlation from the Adiabatic Connection for Multireference Wave Functions

Science.gov (United States)

Pernal, Katarzyna

2018-01-01

An adiabatic connection (AC) formula for the electron correlation energy is derived for a broad class of multireference wave functions. The AC expression recovers dynamic correlation energy and assures a balanced treatment of the correlation energy. Coupling the AC formalism with the extended random phase approximation allows one to find the correlation energy only from reference one- and two-electron reduced density matrices. If the generalized valence bond perfect pairing model is employed a simple closed-form expression for the approximate AC formula is obtained. This results in the overall M5 scaling of the computation cost making the method one of the most efficient multireference approaches accounting for dynamic electron correlation also for the strongly correlated systems.

On propagation of axisymmetric waves in pressurized functionally graded elastomeric hollow cylinders

Science.gov (United States)

Wu, Bin; Su, Yipin; Liu, Dongying; Chen, Weiqiu; Zhang, Chuanzeng

2018-05-01

Soft materials can be designed with a functionally graded (FG) property for specific applications. Such material inhomogeneity can also be found in many soft biological tissues whose functionality is only partly understood to date. In this paper, we analyze the axisymmetric guided wave propagation in a pressurized FG elastomeric hollow cylinder. The cylinder is subjected to a combined action of axial pre-stretch and pressure difference applied to the inner and outer cylindrical surfaces. We consider both torsional waves and longitudinal waves propagating in the FG cylinder made of incompressible isotropic elastomer, which is characterized by the Mooney-Rivlin strain energy function but with the material parameters varying with the radial coordinate in an affine way. The pressure difference generates an inhomogeneous deformation field in the FG cylinder, which dramatically complicates the superimposed wave problem described by the small-on-large theory. A particularly efficient approach is hence employed which combines the state-space formalism for the incremental wave motion with the approximate laminate or multi-layer technique. Dispersion relations for the two types of axisymmetric guided waves are then derived analytically. The accuracy and convergence of the proposed approach is validated numerically. The effects of the pressure difference, material gradient, and axial pre-stretch on both the torsional and the longitudinal wave propagation characteristics are discussed in detail through numerical examples. It is found that the frequency of axisymmetric waves depends nonlinearly on the pressure difference and the material gradient, and an increase in the material gradient enhances the capability of the pressure difference to adjust the wave behavior in the FG cylinder. This work provides a theoretical guidance for characterizing FG soft materials by in-situ ultrasonic nondestructive evaluation and for designing tunable waveguides via material tailoring along

Simulating propagation of decoupled elastic waves using low-rank approximate mixed-domain integral operators for anisotropic media

KAUST Repository

Cheng, Jiubing; Alkhalifah, Tariq Ali; Wu, Zedong; Zou, Peng; Wang, Chenlong

2016-01-01

In elastic imaging, the extrapolated vector fields are decoupled into pure wave modes, such that the imaging condition produces interpretable images. Conventionally, mode decoupling in anisotropic media is costly because the operators involved are dependent on the velocity, and thus they are not stationary. We have developed an efficient pseudospectral approach to directly extrapolate the decoupled elastic waves using low-rank approximate mixed-domain integral operators on the basis of the elastic displacement wave equation. We have applied k-space adjustment to the pseudospectral solution to allow for a relatively large extrapolation time step. The low-rank approximation was, thus, applied to the spectral operators that simultaneously extrapolate and decompose the elastic wavefields. Synthetic examples on transversely isotropic and orthorhombic models showed that our approach has the potential to efficiently and accurately simulate the propagations of the decoupled quasi-P and quasi-S modes as well as the total wavefields for elastic wave modeling, imaging, and inversion.

Simulating propagation of decoupled elastic waves using low-rank approximate mixed-domain integral operators for anisotropic media

KAUST Repository

Cheng, Jiubing

2016-03-15

In elastic imaging, the extrapolated vector fields are decoupled into pure wave modes, such that the imaging condition produces interpretable images. Conventionally, mode decoupling in anisotropic media is costly because the operators involved are dependent on the velocity, and thus they are not stationary. We have developed an efficient pseudospectral approach to directly extrapolate the decoupled elastic waves using low-rank approximate mixed-domain integral operators on the basis of the elastic displacement wave equation. We have applied k-space adjustment to the pseudospectral solution to allow for a relatively large extrapolation time step. The low-rank approximation was, thus, applied to the spectral operators that simultaneously extrapolate and decompose the elastic wavefields. Synthetic examples on transversely isotropic and orthorhombic models showed that our approach has the potential to efficiently and accurately simulate the propagations of the decoupled quasi-P and quasi-S modes as well as the total wavefields for elastic wave modeling, imaging, and inversion.

A phase space approach to wave propagation with dispersion.

Science.gov (United States)

Ben-Benjamin, Jonathan S; Cohen, Leon; Loughlin, Patrick J

2015-08-01

A phase space approximation method for linear dispersive wave propagation with arbitrary initial conditions is developed. The results expand on a previous approximation in terms of the Wigner distribution of a single mode. In contrast to this previously considered single-mode case, the approximation presented here is for the full wave and is obtained by a different approach. This solution requires one to obtain (i) the initial modal functions from the given initial wave, and (ii) the initial cross-Wigner distribution between different modal functions. The full wave is the sum of modal functions. The approximation is obtained for general linear wave equations by transforming the equations to phase space, and then solving in the new domain. It is shown that each modal function of the wave satisfies a Schrödinger-type equation where the equivalent "Hamiltonian" operator is the dispersion relation corresponding to the mode and where the wavenumber is replaced by the wavenumber operator. Application to the beam equation is considered to illustrate the approach.

Prestack wavefield approximations

KAUST Repository

Alkhalifah, Tariq

2013-01-01

The double-square-root (DSR) relation offers a platform to perform prestack imaging using an extended single wavefield that honors the geometrical configuration between sources, receivers, and the image point, or in other words, prestack wavefields. Extrapolating such wavefields, nevertheless, suffers from limitations. Chief among them is the singularity associated with horizontally propagating waves. I have devised highly accurate approximations free of such singularities which are highly accurate. Specifically, I use Padé expansions with denominators given by a power series that is an order lower than that of the numerator, and thus, introduce a free variable to balance the series order and normalize the singularity. For the higher-order Padé approximation, the errors are negligible. Additional simplifications, like recasting the DSR formula as a function of scattering angle, allow for a singularity free form that is useful for constant-angle-gather imaging. A dynamic form of this DSR formula can be supported by kinematic evaluations of the scattering angle to provide efficient prestack wavefield construction. Applying a similar approximation to the dip angle yields an efficient 1D wave equation with the scattering and dip angles extracted from, for example, DSR ray tracing. Application to the complex Marmousi data set demonstrates that these approximations, although they may provide less than optimal results, allow for efficient and flexible implementations. © 2013 Society of Exploration Geophysicists.

Prestack wavefield approximations

KAUST Repository

Alkhalifah, Tariq

2013-09-01

The double-square-root (DSR) relation offers a platform to perform prestack imaging using an extended single wavefield that honors the geometrical configuration between sources, receivers, and the image point, or in other words, prestack wavefields. Extrapolating such wavefields, nevertheless, suffers from limitations. Chief among them is the singularity associated with horizontally propagating waves. I have devised highly accurate approximations free of such singularities which are highly accurate. Specifically, I use Padé expansions with denominators given by a power series that is an order lower than that of the numerator, and thus, introduce a free variable to balance the series order and normalize the singularity. For the higher-order Padé approximation, the errors are negligible. Additional simplifications, like recasting the DSR formula as a function of scattering angle, allow for a singularity free form that is useful for constant-angle-gather imaging. A dynamic form of this DSR formula can be supported by kinematic evaluations of the scattering angle to provide efficient prestack wavefield construction. Applying a similar approximation to the dip angle yields an efficient 1D wave equation with the scattering and dip angles extracted from, for example, DSR ray tracing. Application to the complex Marmousi data set demonstrates that these approximations, although they may provide less than optimal results, allow for efficient and flexible implementations. © 2013 Society of Exploration Geophysicists.

Rapidity resummation for B-meson wave functions

Directory of Open Access Journals (Sweden)

Shen Yue-Long

2014-01-01

Full Text Available Transverse-momentum dependent (TMD hadronic wave functions develop light-cone divergences under QCD corrections, which are commonly regularized by the rapidity ζ of gauge vector defining the non-light-like Wilson lines. The yielding rapidity logarithms from infrared enhancement need to be resummed for both hadronic wave functions and short-distance functions, to achieve scheme-independent calculations of physical quantities. We briefly review the recent progress on the rapidity resummation for B-meson wave functions which are the key ingredients of TMD factorization formulae for radiative-leptonic, semi-leptonic and non-leptonic B-meson decays. The crucial observation is that rapidity resummation induces a strong suppression of B-meson wave functions at small light-quark momentum, strengthening the applicability of TMD factorization in exclusive B-meson decays. The phenomenological consequence of rapidity-resummation improved B-meson wave functions is further discussed in the context of B → π transition form factors at large hadronic recoil.

Dynamic equations for gauge-invariant wave functions

International Nuclear Information System (INIS)

Kapshaj, V.N.; Skachkov, N.B.; Solovtsov, I.L.

1984-01-01

The Bethe-Salpeter and quasipotential dynamic equations for wave functions of relative quark motion, have been derived. Wave functions are determined by the gauge invariant method. The V.A. Fock gauge condition is used in the construction. Despite the transl tional noninvariance of the gauge condition the standard separation of variables has been obtained and wave function doesn't contain gauge exponents

Uniform analytic approximation of Wigner rotation matrices

Science.gov (United States)

Hoffmann, Scott E.

2018-02-01

We derive the leading asymptotic approximation, for low angle θ, of the Wigner rotation matrix elements, dm1m2 j(θ ) , uniform in j, m1, and m2. The result is in terms of a Bessel function of integer order. We numerically investigate the error for a variety of cases and find that the approximation can be useful over a significant range of angles. This approximation has application in the partial wave analysis of wavepacket scattering.

Angular momentum projection on a mesh of cranked Hartree-Fock wave functions

International Nuclear Information System (INIS)

Baye, D.; Heenen, P.

1984-01-01

A method for projecting on angular momentum wave functions discretized on a three-dimensional Cartesian mesh is presented. The method is based on a matrix representation of the rotation operator. It is applied to cranked Hartree-Fock wave functions calculated for 24 Mg with a simple interaction. In this case, the accuracy of the projected matrix elements is estimated to be of the order of 0.1%. An extensive comparison of the projected and cranking energies is made. The validity of the cranking method as an approximation to a variation-after-projection calculation seems to be wider than usually expected. The study of the fission barrier of 24 Mg for the channel 4 He- 16 O- 4 He shows that the cranking predictions for these very deformed states are quite reliable

Wave equation of a nonlinear triatomic molecule and the adiabatic correction to the Born--Oppenheimer approximation

International Nuclear Information System (INIS)

Bardo, R.D.; Wolfsberg, M.

1977-01-01

The wave equation for a nonlinear polyatomic molecule is formulated in molecule-fixed coordinates by a method originally due to Hirschfelder and Wigner. Application is made to a triatomic molecule, and the wave equation is explicitly presented in a useful molecule-fixed coordinate system. The formula for the adiabatic correction to the Born--Oppenheimer approximation for a triatomic molecule is obtained. The extension of the present formulation to larger polyatomic molecules is pointed out. Some terms in the triatomic molecule wave equation are discussed in detail

Statistical wave function

International Nuclear Information System (INIS)

Levine, R.D.

1988-01-01

Statistical considerations are applied to quantum mechanical amplitudes. The physical motivation is the progress in the spectroscopy of highly excited states. The corresponding wave functions are strongly mixed. In terms of a basis set of eigenfunctions of a zeroth-order Hamiltonian with good quantum numbers, such wave functions have contributions from many basis states. The vector x is considered whose components are the expansion coefficients in that basis. Any amplitude can be written as a dagger x x. It is argued that the components of x and hence other amplitudes can be regarded as random variables. The maximum entropy formalism is applied to determine the corresponding distribution function. Two amplitudes a dagger x x and b dagger x x are independently distributed if b dagger x a = 0. It is suggested that the theory of quantal measurements implies that, in general, one can one determine the distribution of amplitudes and not the amplitudes themselves

Approximated solutions to Born-Infeld dynamics

International Nuclear Information System (INIS)

Ferraro, Rafael; Nigro, Mauro

2016-01-01

The Born-Infeld equation in the plane is usefully captured in complex language. The general exact solution can be written as a combination of holomorphic and anti-holomorphic functions. However, this solution only expresses the potential in an implicit way. We rework the formulation to obtain the complex potential in an explicit way, by means of a perturbative procedure. We take care of the secular behavior common to this kind of approach, by resorting to a symmetry the equation has at the considered order of approximation. We apply the method to build approximated solutions to Born-Infeld electrodynamics. We solve for BI electromagnetic waves traveling in opposite directions. We study the propagation at interfaces, with the aim of searching for effects susceptible to experimental detection. In particular, we show that a reflected wave is produced when a wave is incident on a semi-space containing a magnetostatic field.

Approximated solutions to Born-Infeld dynamics

Energy Technology Data Exchange (ETDEWEB)

Ferraro, Rafael [Instituto de Astronomía y Física del Espacio (IAFE, CONICET-UBA),Casilla de Correo 67, Sucursal 28, 1428 Buenos Aires (Argentina); Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires,Ciudad Universitaria, Pabellón I, 1428 Buenos Aires (Argentina); Nigro, Mauro [Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires,Ciudad Universitaria, Pabellón I, 1428 Buenos Aires (Argentina)

2016-02-01

The Born-Infeld equation in the plane is usefully captured in complex language. The general exact solution can be written as a combination of holomorphic and anti-holomorphic functions. However, this solution only expresses the potential in an implicit way. We rework the formulation to obtain the complex potential in an explicit way, by means of a perturbative procedure. We take care of the secular behavior common to this kind of approach, by resorting to a symmetry the equation has at the considered order of approximation. We apply the method to build approximated solutions to Born-Infeld electrodynamics. We solve for BI electromagnetic waves traveling in opposite directions. We study the propagation at interfaces, with the aim of searching for effects susceptible to experimental detection. In particular, we show that a reflected wave is produced when a wave is incident on a semi-space containing a magnetostatic field.

Comparison of four support-vector based function approximators

NARCIS (Netherlands)

de Kruif, B.J.; de Vries, Theodorus J.A.

2004-01-01

One of the uses of the support vector machine (SVM), as introduced in V.N. Vapnik (2000), is as a function approximator. The SVM and approximators based on it, approximate a relation in data by applying interpolation between so-called support vectors, being a limited number of samples that have been

Three-Dimensional Visualization of Wave Functions for Rotating Molecule: Plot of Spherical Harmonics

Science.gov (United States)

Nagaoka, Shin-ichi; Teramae, Hiroyuki; Nagashima, Umpei

2013-01-01

At an early stage of learning quantum chemistry, undergraduate students usually encounter the concepts of the particle in a box, the harmonic oscillator, and then the particle on a sphere. Rotational levels of a diatomic molecule can be well approximated by the energy levels of the particle on a sphere. Wave functions for the particle in a…

Boundary conditions of the exact impulse wave function

International Nuclear Information System (INIS)

Gravielle, M.; Miraglia, J.E.

1997-01-01

The behavior of the exact impulse wave function is investigated at intermediate and high impact energies. Numerical details of the wave function and its perturbative potential are reported. We conclude that the impulse wave function does not tend to the proper Coulomb asymptotic limit. For electron capture, however, it is shown that the impulse wave function produces reliable probabilities even for intermediate velocities and symmetric collision systems. copyright 1997 The American Physical Society

CDW (continuum distorted wave) type approximation for electron capture at large angles

International Nuclear Information System (INIS)

Fojon, O.A.; Maidagan, J.M.

1990-01-01

A calculation is made for the probability of electron capture in shell K at great angles using a second order symmetrical model used related to the continuum distorted wave (CDW) approximation. The influence of Coulomb distortion of nuclei is studied and compared with OBK and CIS type calculations. Numerical results are compared with experimental results of the collision of H + on C at intermediate energies. (Author). 19 refs., 2 figs

Function approximation of tasks by neural networks

International Nuclear Information System (INIS)

Gougam, L.A.; Chikhi, A.; Mekideche-Chafa, F.

2008-01-01

For several years now, neural network models have enjoyed wide popularity, being applied to problems of regression, classification and time series analysis. Neural networks have been recently seen as attractive tools for developing efficient solutions for many real world problems in function approximation. The latter is a very important task in environments where computation has to be based on extracting information from data samples in real world processes. In a previous contribution, we have used a well known simplified architecture to show that it provides a reasonably efficient, practical and robust, multi-frequency analysis. We have investigated the universal approximation theory of neural networks whose transfer functions are: sigmoid (because of biological relevance), Gaussian and two specified families of wavelets. The latter have been found to be more appropriate to use. The aim of the present contribution is therefore to use a m exican hat wavelet a s transfer function to approximate different tasks relevant and inherent to various applications in physics. The results complement and provide new insights into previously published results on this problem

Soft dipole mode of neutron-rich light nuclei in asymptotic potential approximation

International Nuclear Information System (INIS)

Filippov, G.F.; Lashko, Yu.A.; Shvedov, L.P.

2000-01-01

Completely antisymmetrized 1''-continuum wave functions as well as the ground state wave function for ''6He have been constructed in asymptotic potential approximation. The behaviour of two-channel S-matrix elements shows on the existence of 1''- resonant state just above the three-body decay threshold of ''6He

Deep inelastic scattering and light-cone wave functions

International Nuclear Information System (INIS)

Belyaev, V.M.; Johnson, M.B.

1996-01-01

In the framework of light-cone QCD rules, we study the valence quark distribution function q(x B ) of a pion for moderate x B . The sum rule with the leading twist-2 wave function gives q(x B ) = φ π (x B ). Twist-4 wave functions give about 30% for x B ∼0.5. It is shown that QCD sum rule predictions, with the asymptotic pion wave function, are in good agreement with experimental data. We found that a two-hump profile for the twist-2 wave function leads to a valence quark distribution function that contradicts experimental data

Nonlinear Dispersive Elastic Waves in Solids: Exact, Approximate, and Numerical Solutions

Science.gov (United States)

Khajehtourian, Romik

extended method, which is based on a standard transfer-matrix formulation augmented with a nonlinear enrichment at the constitutive material level, yields an approximate band structure that is accurate to an amplitude that is roughly one eighth of the unit cell length. This approach represents a new paradigm for examining the balance between periodicity and nonlinearity in shaping the nature of wave motion.

On single nucleon wave functions in nuclei

International Nuclear Information System (INIS)

Talmi, Igal

2011-01-01

The strong and singular interaction between nucleons, makes the nuclear many body theory very complicated. Still, nuclei exhibit simple and regular features which are simply described by the shell model. Wave functions of individual nucleons may be considered just as model wave functions which bear little resemblance to the real ones. There is, however, experimental evidence for the reality of single nucleon wave functions. There is a simple method of constructing such wave functions for valence nucleons. It is shown that this method can be improved by considering the polarization of the core by the valence nucleon. This gives rise to some rearrangement energy which affects the single valence nucleon energy within the nucleus.

Relation between equal-time and light-front wave functions

International Nuclear Information System (INIS)

Miller, Gerald A.; Tiburzi, Brian C.

2010-01-01

The relation between equal-time and light-front wave functions is studied using models for which the four-dimensional solution of the Bethe-Salpeter wave function can be obtained. The popular prescription of defining the longitudinal momentum fraction using the instant-form free kinetic energy and third component of momentum is found to be incorrect except in the nonrelativistic limit. One may obtain light-front wave functions from rest-frame, instant-form wave functions by boosting the latter wave functions to the infinite momentum frame. Despite this difficulty, we prove a relation between certain integrals of the equal-time and light-front wave functions.

Low-energy P-wave phaseshifts for positron-hydrogen elastic scattering using an adiabatic approximation

International Nuclear Information System (INIS)

Armour, E.A.G.; Beker, C.A.; Farina, J.E.G.

1981-01-01

P-wave phaseshifts for positron-hydrogen elastic scattering are calculated using a new adiabatic approximation in which the length of the radius vector from the proton to the positron is fixed but its direction is allowed to vary. This adiabatic approximation makes possible the full inclusion in the calculation of virtual states in which angular momentum is transferred to the target H atom. The results obtained agree qualitatively with the highly accurate results of Bhatia and co-workers (Phys. Rev.; A9:219 (1974)) and are much closer to them than the results obtained using the usual adiabatic approximation in which the radius vector from the proton to the positron is fixed. (author)

Taylor-series method for four-nucleon wave functions

International Nuclear Information System (INIS)

Sandulescu, A.; Tarnoveanu, I.; Rizea, M.

1977-09-01

Taylor-series method for transforming the infinite or finite well two-nucleon wave functions from individual coordinates to relative and c.m. coordinates, by expanding the single particle shell model wave functions around c.m. of the system, is generalized to four-nucleon wave functions. Also the connections with the Talmi-Moshinsky method for two and four harmonic oscillator wave functions are deduced. For both methods Fortran IV programs for the expansion coefficients have been written and the equivalence of corresponding expressions numerically proved. (author)

Covariance Function for Nearshore Wave Assimilation Systems

Science.gov (United States)

2018-01-30

which is applicable for any spectral wave model. The four dimensional variational (4DVar) assimilation methods are based on the mathematical ...covariance can be modeled by a parameterized Gaussian function, for nearshore wave assimilation applications , the covariance function depends primarily on...SPECTRAL ACTION DENSITY, RESPECTIVELY. ............................ 5 FIGURE 2. TOP ROW: STATISTICAL ANALYSIS OF THE WAVE-FIELD PROPERTIES AT THE

Green function for three-wave coupling problems

International Nuclear Information System (INIS)

Molevich, N E

2001-01-01

The Green function is found for three-wave coupling problems. The function was used for analysis of parametric amplification in dissipative and active media. It is shown that the parametric increment in active media can become exponential. As an example, the nonstationary stimulated scattering of electromagnetic waves by sound and temperatures waves is considered. (nonlinear optical phenomena)

Test of distorted wave kinematic coupling approximation calculations for knockout reactions

International Nuclear Information System (INIS)

Jain, A.K.

1990-01-01

A test has been devised to check the validity of conventional distorted-wave impulse approximation (DWIA) treatment of knockout reactions. The conventional DWIA formalism separates the three-body final state Schroedinger equation for a knockout reaction into two two-body Schroedinger equations by assuming an asymptotic constant value for the three-body coupling term commonly known as the kinematic coupling approximation (KCA). In the test case, which consists of an extreme asymmetric situation where one of the distorting optical potentials is assumed to vanish, the three-body final state Schroedinger equation can be solved exactly as a product of two two-body solutions using one particular set of relative coordinates. Large influence of the three-body coupling term is seen in the comparison of the exact and KCA results for (α,2α) and (p,pα) knockout reactions when the distorting optical potentials are weakly absorbing

Mathieu functions describing particles evolving in electromagnetic waves

Science.gov (United States)

Mihu, Denisa-Andreea; Dariescu, Marina-Aura

2017-12-01

Solutions of Klein-Gordon equation for particles moving in a standing wave configuration bring into attention an intricate and complicated category of special functions, namely the Mathieu functions. The stability of the solutions governed by the intercorrelation between Mathieu equation' parameters is discussed. For specific intervals of the wave number, the instability regime installs, pointing out the tendency of exponential growth for the oscillatory wave functions, as a consequence of parametric resonance phenomenon. The expression of the wave function allows the computation of the four-dimensional conserved current density components.

Six Impossible Things: Fractional Charge From Laughlin's Wave Function

International Nuclear Information System (INIS)

Shrivastava, Keshav N.

2010-01-01

The Laughlin's wave function is found to be the zero-energy ground state of a δ-function Hamiltonian. The finite negative value of the ground state energy which is 91 per cent of Wigner value, can be obtained only when Coulomb correlations are introduced. The Laughlin's wave function is of short range and it overlaps with that of the exact wave functions of small (number of electrons 2 or 5) systems. (i) It is impossible to obtain fractional charge from Laughlin's wave function. (ii) It is impossible to prove that the Laughlin's wave function gives the ground state of the Coulomb Hamiltonian. (iii) It is impossible to have particle-hole symmetry in the Laughlin's wave function. (iv) It is impossible to derive the value of m in the Laughlin's wave function. The value of m in ψ m can not be proved to be 3 or 5. (v) It is impossible to prove that the Laughlin's state is incompressible because the compressible states are also likely. (vi) It is impossible for the Laughlin's wave function to have spin. This effort is directed to explain the experimental data of quantum Hall effect in GaAs/AlGaAs.

Nonlinear Schroedinger Approximations for Partial Differential Equations with Quadratic and Quasilinear Terms

Science.gov (United States)

Cummings, Patrick

We consider the approximation of solutions of two complicated, physical systems via the nonlinear Schrodinger equation (NLS). In particular, we discuss the evolution of wave packets and long waves in two physical models. Due to the complicated nature of the equations governing many physical systems and the in-depth knowledge we have for solutions of the nonlinear Schrodinger equation, it is advantageous to use approximation results of this kind to model these physical systems. The approximations are simple enough that we can use them to understand the qualitative and quantitative behavior of the solutions, and by justifying them we can show that the behavior of the approximation captures the behavior of solutions to the original equation, at least for long, but finite time. We first consider a model of the water wave equations which can be approximated by wave packets using the NLS equation. We discuss a new proof that both simplifies and strengthens previous justification results of Schneider and Wayne. Rather than using analytic norms, as was done by Schneider and Wayne, we construct a modified energy functional so that the approximation holds for the full interval of existence of the approximate NLS solution as opposed to a subinterval (as is seen in the analytic case). Furthermore, the proof avoids problems associated with inverting the normal form transform by working with a modified energy functional motivated by Craig and Hunter et al. We then consider the Klein-Gordon-Zakharov system and prove a long wave approximation result. In this case there is a non-trivial resonance that cannot be eliminated via a normal form transform. By combining the normal form transform for small Fourier modes and using analytic norms elsewhere, we can get a justification result on the order 1 over epsilon squared time scale.

High energy approximations for nuclear knockout form factors at small momentum transfer

International Nuclear Information System (INIS)

Amado, R.D.; Cannata, F.; Dedonder, J.P.

1985-01-01

We obtain an explicit approximate expression for the nucleon knockout form factor at small momentum transfer induced by a scalar probe in a single particle model in terms of the momentum space bound state wave function. Our form preserves the orthogonality constraint without using explicitly the final state scattering wave function. We examine the leading large momentum behavior of the momentum space wave function and of correction terms to our expression for the form factor in the case where the bound state is an s state

Bound-state wave functions at rest in describing deep inelastic scattering

International Nuclear Information System (INIS)

Khvedelidze, A.M.; Kvinikhidze, A.N.

1991-01-01

The deep inelastic process of the lepton-hadron scattering is studied in the bound-state rest frame. A new version of expanding structure functions in interaction constant powers is proposed, each term in it having spectral properties. This expansion makes it possible to consider contributions of composites in the final state to the cross section. It is shown that, as compared with the system P z →∞, the impulse approximation is insufficient for describing correctly the elastic limit in the composite particle rest frame. The leading asymptotics of structure functions as χ Bj →1 can be obtained by taking into account the interaction of contituents in the final state. It is shown that in contrast to the 'light-cone' formalism the ratio F 2 en (χ)/F 2 ep (χ) as χ Bj →1 depends on the explicit form of the spatial part of the nucleon wave function and, in particular, assuming the relativistic character of internal motion, it may be lower than the well-known prediction (i.e. 3/7). This is due to the correct consideration of spin degrees of freedom of the wave function of the nucleon at rest. (orig.)

NN wave function at small distances and hard bremsstrahlung in the process pp→ppγ

International Nuclear Information System (INIS)

Neudatchin, V. G.; Khokhlov, N. A.; Shirokov, A. M.; Knyr, V. A.

1997-01-01

Various possibilities of studying the NN wave function at small distances--and in particular, quark degrees of freedom in the NN system--are discussed. It is shown that there is such a possibility at moderate energies--namely, hard bremsstrahlung in the process pp→ppγ at proton-beam energies in the range 350-450 MeV permits distinguishing between the pp wave function with nodes in S and P waves that corresponds to the Moscow potential of NN interaction from functions obtained with repulsive-core mesonic potentials. In the regions where photon energies in the c.m.s. are maximal (forward and backward photon emission angles in the laboratory frame), the pp→ppγ cross section calculated with the Moscow potential has maxima at which it is approximately five times larger than the analogous cross section calculated with repulsive-core mesonic potentials. The coordinate-representation formalism of the theory of bremsstrahlung is expounded

Computing wave functions in multichannel collisions with non-local potentials using the R-matrix method

Science.gov (United States)

Bonitati, Joey; Slimmer, Ben; Li, Weichuan; Potel, Gregory; Nunes, Filomena

2017-09-01

The calculable form of the R-matrix method has been previously shown to be a useful tool in approximately solving the Schrodinger equation in nuclear scattering problems. We use this technique combined with the Gauss quadrature for the Lagrange-mesh method to efficiently solve for the wave functions of projectile nuclei in low energy collisions (1-100 MeV) involving an arbitrary number of channels. We include the local Woods-Saxon potential, the non-local potential of Perey and Buck, a Coulomb potential, and a coupling potential to computationally solve for the wave function of two nuclei at short distances. Object oriented programming is used to increase modularity, and parallel programming techniques are introduced to reduce computation time. We conclude that the R-matrix method is an effective method to predict the wave functions of nuclei in scattering problems involving both multiple channels and non-local potentials. Michigan State University iCER ACRES REU.

Noncommuting limits of oscillator wave functions

International Nuclear Information System (INIS)

Daboul, J.; Pogosyan, G. S.; Wolf, K. B.

2007-01-01

Quantum harmonic oscillators with spring constants k > 0 plus constant forces f exhibit rescaled and displaced Hermite-Gaussian wave functions, and discrete, lower bound spectra. We examine their limits when (k, f) → (0, 0) along two different paths. When f → 0 and then k → 0, the contraction is standard: the system becomes free with a double continuous, positive spectrum, and the wave functions limit to plane waves of definite parity. On the other hand, when k → 0 first, the contraction path passes through the free-fall system, with a continuous, nondegenerate, unbounded spectrum and displaced Airy wave functions, while parity is lost. The subsequent f → 0 limit of the nonstandard path shows the dc hysteresis phenomenon of noncommuting contractions: the lost parity reappears as an infinitely oscillating superposition of the two limiting solutions that are related by the symmetry

Calculation of deuteron wave functions with relativistic interactions

International Nuclear Information System (INIS)

Buck, W.W. III.

1976-01-01

Deuteron wave functions with a repulsive core are obtained numerically from a fully relativistic wave equation introduced by Gross. The numerical technique enables analytic solutions for classes of interactions composed of the relativistic exchanges of a single pion and a single phenomenological meson, sigma. The pion is chosen to interact as a mixture of pseudoscalar and pseudovector. The amount of mixture is determined by a free mixing parameter, lambda, ranging between 1 (pure pseudoscalar) and (pure pseudovector). Each value of lambda corresponds, then, to a different interaction. Solutions are found for lambda = 1, .9, .8, .6, and 0. The wave functions for each interaction come in a group of four. Of the four wave functions, two are the usual S and D state wave functions, while the remaining two, arising out of the relativistic prescription, are identified as 3 P 1 and 1 P 1 wave functions (P state wave functions). For the interactions solved for, the D state probabilities ranged between 5.1 percent and 6.3 percent, while the total P state probabilities ranged between 0.7 percent and 2.7 percent. The method of obtaining solutions was to adjust the sigma meson parameters to give the correct binding energy and a good quadrupole moment. All wave functions obtained are applied to relativistic N-d scattering in the backward direction where the effect of the P states is quite measurable

Density functional theory of electron transfer beyond the Born-Oppenheimer approximation: Case study of LiF

Science.gov (United States)

Li, Chen; Requist, Ryan; Gross, E. K. U.

2018-02-01

We perform model calculations for a stretched LiF molecule, demonstrating that nonadiabatic charge transfer effects can be accurately and seamlessly described within a density functional framework. In alkali halides like LiF, there is an abrupt change in the ground state electronic distribution due to an electron transfer at a critical bond length R = Rc, where an avoided crossing of the lowest adiabatic potential energy surfaces calls the validity of the Born-Oppenheimer approximation into doubt. Modeling the R-dependent electronic structure of LiF within a two-site Hubbard model, we find that nonadiabatic electron-nuclear coupling produces a sizable elongation of the critical Rc by 0.5 bohr. This effect is very accurately captured by a simple and rigorously derived correction, with an M-1 prefactor, to the exchange-correlation potential in density functional theory, M = reduced nuclear mass. Since this nonadiabatic term depends on gradients of the nuclear wave function and conditional electronic density, ∇Rχ(R) and ∇Rn(r, R), it couples the Kohn-Sham equations at neighboring R points. Motivated by an observed localization of nonadiabatic effects in nuclear configuration space, we propose a local conditional density approximation—an approximation that reduces the search for nonadiabatic density functionals to the search for a single function y(n).

9Be scattering with microscopic wave functions and the continuum-discretized coupled-channel method

Science.gov (United States)

Descouvemont, P.; Itagaki, N.

2018-01-01

We use microscopic 9Be wave functions defined in a α +α +n multicluster model to compute 9Be+target scattering cross sections. The parameter sets describing 9Be are generated in the spirit of the stochastic variational method, and the optimal solution is obtained by superposing Slater determinants and by diagonalizing the Hamiltonian. The 9Be three-body continuum is approximated by square-integral wave functions. The 9Be microscopic wave functions are then used in a continuum-discretized coupled-channel (CDCC) calculation of 9Be+208Pb and of 9Be+27Al elastic scattering. Without any parameter fitting, we obtain a fair agreement with experiment. For a heavy target, the influence of 9Be breakup is important, while it is weaker for light targets. This result confirms previous nonmicroscopic CDCC calculations. One of the main advantages of the microscopic CDCC is that it is based on nucleon-target interactions only; there is no adjustable parameter. The present work represents a first step towards more ambitious calculations involving heavier Be isotopes.

Approximate self-consistent potentials for density-functional-theory exchange-correlation functionals

International Nuclear Information System (INIS)

Cafiero, Mauricio; Gonzalez, Carlos

2005-01-01

We show that potentials for exchange-correlation functionals within the Kohn-Sham density-functional-theory framework may be written as potentials for simpler functionals multiplied by a factor close to unity, and in a self-consistent field calculation, these effective potentials find the correct self-consistent solutions. This simple theory is demonstrated with self-consistent exchange-only calculations of the atomization energies of some small molecules using the Perdew-Kurth-Zupan-Blaha (PKZB) meta-generalized-gradient-approximation (meta-GGA) exchange functional. The atomization energies obtained with our method agree with or surpass previous meta-GGA calculations performed in a non-self-consistent manner. The results of this work suggest the utility of this simple theory to approximate exchange-correlation potentials corresponding to energy functionals too complicated to generate closed forms for their potentials. We hope that this method will encourage the development of complex functionals which have correct boundary conditions and are free of self-interaction errors without the worry that the functionals are too complex to differentiate to obtain potentials

Improved wave functions for large-N expansions

International Nuclear Information System (INIS)

Imbo, T.; Sukhatme, U.

1985-01-01

Existing large-N expansions of radial wave functions phi/sub n/,l(r) are only accurate near the minimum of the effective potential. Within the framework of the shifted 1/N expansion, we use known analytic results to motivate a simple modification so that the improved wave functions are accurate over a wide range of r and any choice of quantum numbers n and l. It is shown that these wave functions yield simple and accurate analytic expressions for certain quantities of interest in quarkonium physics

Approximation of Analytic Functions by Bessel's Functions of Fractional Order

Directory of Open Access Journals (Sweden)

Soon-Mo Jung

2011-01-01

Full Text Available We will solve the inhomogeneous Bessel's differential equation x2y″(x+xy′(x+(x2-ν2y(x=∑m=0∞amxm, where ν is a positive nonintegral number and apply this result for approximating analytic functions of a special type by the Bessel functions of fractional order.

Sequential function approximation on arbitrarily distributed point sets

Science.gov (United States)

Wu, Kailiang; Xiu, Dongbin

2018-02-01

We present a randomized iterative method for approximating unknown function sequentially on arbitrary point set. The method is based on a recently developed sequential approximation (SA) method, which approximates a target function using one data point at each step and avoids matrix operations. The focus of this paper is on data sets with highly irregular distribution of the points. We present a nearest neighbor replacement (NNR) algorithm, which allows one to sample the irregular data sets in a near optimal manner. We provide mathematical justification and error estimates for the NNR algorithm. Extensive numerical examples are also presented to demonstrate that the NNR algorithm can deliver satisfactory convergence for the SA method on data sets with high irregularity in their point distributions.

Nonstandard jump functions for radically symmetric shock waves

International Nuclear Information System (INIS)

Baty, Roy S.; Tucker, Don H.; Stanescu, Dan

2008-01-01

Nonstandard analysis is applied to derive generalized jump functions for radially symmetric, one-dimensional, magnetogasdynamic shock waves. It is assumed that the shock wave jumps occur on infinitesimal intervals and the jump functions for the physical parameters occur smoothly across these intervals. Locally integrable predistributions of the Heaviside function are used to model the flow variables across a shock wave. The equations of motion expressed in nonconservative form are then applied to derive unambiguous relationships between the jump functions for the physical parameters for two families of self-similar flows. It is shown that the microstructures for these families of radially symmetric, magnetogasdynamic shock waves coincide in a nonstandard sense for a specified density jump function.

Integral transform technique for meson wave functions

International Nuclear Information System (INIS)

Bakulev, A.P.; Mikhajlov, S.V.

1996-01-01

In a recent paper [1] we proposed a new approach for extracting the wave function of the π-meson φ π (x) and the masses and wave functions of its first resonances from the new QCD sum rules for nondiagonal correlators obtained in [2]. Here, we test our approach using an exactly solvable toy model as an illustrating example. We demonstrate the validity of the method and suggest a pure algebraic procedure for extracting the masses and wave functions relating to the case under investigation. We also explore the stability of the procedure under perturbations of the theoretical part of the sum rule. In application to the pion case, this results not only in the mass and wave function of the first resonance (π'), but also in the estimation of π''-mass. 17 refs., 11 figs

Expansion of continuum functions on resonance wave functions and amplitudes

International Nuclear Information System (INIS)

Bang, J.; Gareev, F.A.; Gizzatkulov, M.H.; Goncharov, S.A.

1978-01-01

To overcome difficulties encountered with wave functions of continuum spectrum (for example, in a shell model with continuum) the pole expansion (by the Mittag-Leffler theorem) of wave functions, scattering amplitudes and the Green functions with positive energies are considered. It is shown that resonance functions (the Gamov functions) form a complete set over which the continuum functions could be expanded. The general view of these expansions for final potentials and for the Coulomb repulsion potential are obtained and discussed. It is shown that the application of the method to nuclear structure calculations leads to simple algebraic equations

Approximated treatment of the Pauli principle effects in elastic collisons

International Nuclear Information System (INIS)

Schechter, H.

1984-08-01

Exact microscopic methods like the RGM (Resonanting Group Method) and the GCM (Generator Coordinate Method) and approximate methods like the OCM (Orthogonality Condition Model) are used to study the effects of Pauli Principle in α- 16 O elastic scattering. Using V2 and BL nucleon-nucleon interactions, nucleus-nucleus effective potentials are obtained from RGM 'exact' wave functions and also from an approximate method developed previoulsy. Using these potentials in the OCM Saito Equation phase-shifts are calculated for partial waves Λ = 0, 1, ... 11, in the energy range 0 [pt

Cheap contouring of costly functions: the Pilot Approximation Trajectory algorithm

International Nuclear Information System (INIS)

Huttunen, Janne M J; Stark, Philip B

2012-01-01

The Pilot Approximation Trajectory (PAT) contour algorithm can find the contour of a function accurately when it is not practical to evaluate the function on a grid dense enough to use a standard contour algorithm, for instance, when evaluating the function involves conducting a physical experiment or a computationally intensive simulation. PAT relies on an inexpensive pilot approximation to the function, such as interpolating from a sparse grid of inexact values, or solving a partial differential equation (PDE) numerically using a coarse discretization. For each level of interest, the location and ‘trajectory’ of an approximate contour of this pilot function are used to decide where to evaluate the original function to find points on its contour. Those points are joined by line segments to form the PAT approximation of the contour of the original function. Approximating a contour numerically amounts to estimating a lower level set of the function, the set of points on which the function does not exceed the contour level. The area of the symmetric difference between the true lower level set and the estimated lower level set measures the accuracy of the contour. PAT measures its own accuracy by finding an upper confidence bound for this area. In examples, PAT can estimate a contour more accurately than standard algorithms, using far fewer function evaluations than standard algorithms require. We illustrate PAT by constructing a confidence set for viscosity and thermal conductivity of a flowing gas from simulated noisy temperature measurements, a problem in which each evaluation of the function to be contoured requires solving a different set of coupled nonlinear PDEs. (paper)

Approximate formulas for elasticity of the Tornquist functions and some their advantages

Science.gov (United States)

Issin, Meyram

2017-09-01

In this article functions of demand for prime necessity, second necessity and luxury goods depending on the income are considered. These functions are called Tornquist functions. By means of the return model the demand for prime necessity goods and second necessity goods are approximately described. Then on the basis of a method of the smallest squares approximate formulas for elasticity of these Tornquist functions are received. To receive an approximate formula for elasticity of function of demand for luxury goods, the linear asymptotic formula is constructed for this function. Some benefits of approximate formulas for elasticity of Tornquist functions are specified.

High-Dimensional Function Approximation With Neural Networks for Large Volumes of Data.

Science.gov (United States)

Andras, Peter

2018-02-01

Approximation of high-dimensional functions is a challenge for neural networks due to the curse of dimensionality. Often the data for which the approximated function is defined resides on a low-dimensional manifold and in principle the approximation of the function over this manifold should improve the approximation performance. It has been show that projecting the data manifold into a lower dimensional space, followed by the neural network approximation of the function over this space, provides a more precise approximation of the function than the approximation of the function with neural networks in the original data space. However, if the data volume is very large, the projection into the low-dimensional space has to be based on a limited sample of the data. Here, we investigate the nature of the approximation error of neural networks trained over the projection space. We show that such neural networks should have better approximation performance than neural networks trained on high-dimensional data even if the projection is based on a relatively sparse sample of the data manifold. We also find that it is preferable to use a uniformly distributed sparse sample of the data for the purpose of the generation of the low-dimensional projection. We illustrate these results considering the practical neural network approximation of a set of functions defined on high-dimensional data including real world data as well.

Improved Wave-vessel Transfer Functions by Uncertainty Modelling

DEFF Research Database (Denmark)

Nielsen, Ulrik Dam; Fønss Bach, Kasper; Iseki, Toshio

2016-01-01

This paper deals with uncertainty modelling of wave-vessel transfer functions used to calculate or predict wave-induced responses of a ship in a seaway. Although transfer functions, in theory, can be calculated to exactly reflect the behaviour of the ship when exposed to waves, uncertainty in inp...

Application of the N-quantum approximation method to bound state problems

International Nuclear Information System (INIS)

Raychaudhuri, A.

1977-01-01

The N-quantum approximation (NQA) method is examined in the light of its application to bound state problems. Bound state wave functions are obtained as expansion coefficients in a truncated Haag expansion. From the equations of motion for the Heisenberg field and the NQA expansion, an equation satisfied by the wave function is derived. Two different bound state systems are considered. In one case, the bound state problem of two identical scalars by scalar exchange is analyzed using the NQA. An integral equation satisfied by the wave function is derived. In the nonrelativistic limit, the equation is shown to reduce to the Schroedinger equation. The equation is solved numerically, and the results compared with those obtained for this system by other methods. The NQA method is also applied to the bound state of two spin 1/2 particles with electromagnetic interaction. The integral equation for the wave function is shown to agree with the corresponding Bethe Salpeter equation in the nonrelativistic limit. Using the Dirac (4 x 4) matrices the wave function is expanded in terms of structure functions and the equation for the wave function is reduced to two disjoint sets of coupled equation for the structure functions

A finite range coupled channel Born approximation code

International Nuclear Information System (INIS)

Nagel, P.; Koshel, R.D.

1978-01-01

The computer code OUKID calculates differential cross sections for direct transfer nuclear reactions in which multistep processes, arising from strongly coupled inelastic states in both the target and residual nuclei, are possible. The code is designed for heavy ion reactions where full finite range and recoil effects are important. Distorted wave functions for the elastic and inelastic scattering are calculated by solving sets of coupled differential equations using a Matrix Numerov integration procedure. These wave functions are then expanded into bases of spherical Bessel functions by the plane-wave expansion method. This approach allows the six-dimensional integrals for the transition amplitude to be reduced to products of two one-dimensional integrals. Thus, the inelastic scattering is treated in a coupled channel formalism while the transfer process is treated in a finite range born approximation formalism. (Auth.)

Calculation of the nucleon structure function from the nucleon wave function

Science.gov (United States)

Hussar, Paul E.

1993-01-01

Harmonic oscillator wave functions have played an historically important role in our understanding of the structure of the nucleon, most notably by providing insight into the mass spectra of the low-lying states. High energy scattering experiments are known to give us a picture of the nucleon wave function at high-momentum transfer and in a frame in which the nucleon is traveling fast. A simple model that crosses the twin bridges of momentum scale and Lorentz frame that separate the pictures of the nucleon wave function provided by the deep inelastic scattering data and by the oscillator model is presented.

Beyond the relativistic mean-field approximation. II. Configuration mixing of mean-field wave functions projected on angular momentum and particle number

International Nuclear Information System (INIS)

Niksic, T.; Vretenar, D.; Ring, P.

2006-01-01

The framework of relativistic self-consistent mean-field models is extended to include correlations related to the restoration of broken symmetries and to fluctuations of collective variables. The generator coordinate method is used to perform configuration mixing of angular-momentum and particle-number projected relativistic wave functions. The geometry is restricted to axially symmetric shapes, and the intrinsic wave functions are generated from the solutions of the relativistic mean-field+Lipkin-Nogami BCS equations, with a constraint on the mass quadrupole moment. The model employs a relativistic point-coupling (contact) nucleon-nucleon effective interaction in the particle-hole channel, and a density-independent δ-interaction in the pairing channel. Illustrative calculations are performed for 24 Mg, 32 S, and 36 Ar, and compared with results obtained employing the model developed in the first part of this work, i.e., without particle-number projection, as well as with the corresponding nonrelativistic models based on Skyrme and Gogny effective interactions

A pair density functional theory utilizing the correlated wave function

International Nuclear Information System (INIS)

Higuchi, M; Higuchi, K

2009-01-01

We propose a practical scheme for calculating the ground-state pair density (PD) by utilizing the correlated wave function. As the correlated wave function, we adopt a linear combination of the single Slater determinants that are constructed from the solutions of the initial scheme [Higuchi M and Higuchi K 2007 Physica B 387, 117]. The single-particle equation is derived by performing the variational principle within the set of PDs that are constructed from such correlated wave functions. Since the search region of the PD is substantially extended as compared with the initial scheme, it is expected that the present scheme can cover more correlation effects. The single-particle equation is practical, and may be easily applied to actual calculations.

Charge symmetry of electron wave functions in a quantized electromagnetic wave field

Energy Technology Data Exchange (ETDEWEB)

Fedorov, M V [AN SSSR, Moscow. Fizicheskij Inst.

1975-01-01

An attempt to clear up the reasons of the electron charge symmetry violation in the quantum wave field was made in this article. For this purpose the connection between the Dirac equation and the electron wave functions in the external field with the exact equation of quantum electrodynamics is established. Attention is paid to the fact that a number of equations for single-electron wave functions can be used in the framework of the same assumptions. It permits the construction of the charge-symmetric solutions in particular.

Faddeev wave function decomposition using bipolar harmonics

International Nuclear Information System (INIS)

Friar, J.L.; Tomusiak, E.L.; Gibson, B.F.; Payne, G.L.

1981-01-01

The standard partial wave (channel) representation for the Faddeev solution to the Schroedinger equation for the ground state of 3 nucleons is written in terms of functions which couple the interacting pair and spectator angular momenta to give S, P, and D waves. For each such coupling there are three terms, one for each of the three cyclic permutations of the nucleon coordinates. A series of spherical harmonic identities is developed which allows writing the Faddeev solution in terms of a basis set of 5 bipolar harmonics: 1 for S waves; 1 for P waves; and 3 for D waves. The choice of a D-wave basis is largely arbitrary, and specific choices correspond to the decomposition schemes of Derrick and Blatt, Sachs, Gibson and Schiff, and Bolsterli and Jezak. The bipolar harmonic form greatly simplifies applications which utilize the wave function, and we specifically discuss the isoscalar charge (or mass) density and the 3 He Coulomb energy

Optimization of nonlinear wave function parameters

International Nuclear Information System (INIS)

Shepard, R.; Minkoff, M.; Chemistry

2006-01-01

An energy-based optimization method is presented for our recently developed nonlinear wave function expansion form for electronic wave functions. This expansion form is based on spin eigenfunctions, using the graphical unitary group approach (GUGA). The wave function is expanded in a basis of product functions, allowing application to closed-shell and open-shell systems and to ground and excited electronic states. Each product basis function is itself a multiconfigurational function that depends on a relatively small number of nonlinear parameters called arc factors. The energy-based optimization is formulated in terms of analytic arc factor gradients and orbital-level Hamiltonian matrices that correspond to a specific kind of uncontraction of each of the product basis functions. These orbital-level Hamiltonian matrices give an intuitive representation of the energy in terms of disjoint subsets of the arc factors, they provide for an efficient computation of gradients of the energy with respect to the arc factors, and they allow optimal arc factors to be determined in closed form for subspaces of the full variation problem. Timings for energy and arc factor gradient computations involving expansion spaces of > 10 24 configuration state functions are reported. Preliminary convergence studies and molecular dissociation curves are presented for some small molecules

Tur\\'an type inequalities for regular Coulomb wave functions

OpenAIRE

Baricz, Árpád

2015-01-01

Tur\\'an, Mitrinovi\\'c-Adamovi\\'c and Wilker type inequalities are deduced for regular Coulomb wave functions. The proofs are based on a Mittag-Leffler expansion for the regular Coulomb wave function, which may be of independent interest. Moreover, some complete monotonicity results concerning the Coulomb zeta functions and some interlacing properties of the zeros of Coulomb wave functions are given.

Reducing Approximation Error in the Fourier Flexible Functional Form

Directory of Open Access Journals (Sweden)

Tristan D. Skolrud

2017-12-01

Full Text Available The Fourier Flexible form provides a global approximation to an unknown data generating process. In terms of limiting function specification error, this form is preferable to functional forms based on second-order Taylor series expansions. The Fourier Flexible form is a truncated Fourier series expansion appended to a second-order expansion in logarithms. By replacing the logarithmic expansion with a Box-Cox transformation, we show that the Fourier Flexible form can reduce approximation error by 25% on average in the tails of the data distribution. The new functional form allows for nested testing of a larger set of commonly implemented functional forms.

Relativistic bound state wave functions

International Nuclear Information System (INIS)

Micu, L.

2005-01-01

A particular method of writing the bound state wave functions in relativistic form is applied to the solutions of the Dirac equation with confining potentials in order to obtain a relativistic description of a quark antiquark bound system representing a given meson. Concerning the role of the effective constituent in the present approach we first observe that without this additional constituent we couldn't expand the bound state wave function in terms of products of free states. Indeed, we notice that if the wave function depends on the relative coordinates only, all the expansion coefficients would be infinite. Secondly we remark that the effective constituent enabled us to give a Lorentz covariant meaning to the potential energy of the bound system which is now seen as the 4th component of a 4-momentum. On the other side, by relating the effective constituent to the quantum fluctuations of the background field which generate the binding, we provided a justification for the existence of some spatial degrees of freedom accompanying the interaction potential. These ones, which are quite unusual in quantum mechanics, in our model are the natural consequence of the the independence of the quarks and can be seen as the effect of the imperfect cancellation of the vector momenta during the quantum fluctuations. Related with all these we remark that the adequate representation for the relativistic description of a bound system is the momentum representation, because of the transparent and easy way of writing the conservation laws and the transformation properties of the wave functions. The only condition to be fulfilled is to find a suitable way to take into account the potential energy of the bound system. A particular feature of the present approach is that the confining forces are due to a kind of glue where both quarks are embedded. This recalls other bound state models where the wave function is factorized in terms of constituent wave functions and the confinement is

Study of Ion Acoustic Wave Damping through Green's Functions

DEFF Research Database (Denmark)

Hsuan, H.C.S.; Jensen, Vagn Orla

1973-01-01

Green's function analyses of ion acoustic waves in streaming plasmas show that, in general, the waves damp algebraically rather than exponentially with distance from exciter.......Green's function analyses of ion acoustic waves in streaming plasmas show that, in general, the waves damp algebraically rather than exponentially with distance from exciter....

Cavity losses for the dissipative Jaynes-Cummings Hamiltonian beyond rotating wave approximation

International Nuclear Information System (INIS)

Scala, M; Militello, B; Messina, A; Maniscalco, S; Piilo, J; Suominen, K-A

2007-01-01

A microscopic derivation of the master equation for the Jaynes-Cummings model with cavity losses is given, taking into account the terms in the dissipator which vary with frequencies of the order of the vacuum Rabi frequency. Our approach allows us to single out physical contexts wherein the usual phenomenological dissipator turns out to be fully justified and constitutes an extension of our previous analysis (Scala et al 2007 Phys. Rev. A 75 013811), where a microscopic derivation was given in the framework of the rotating wave approximation

Wave function of free electron in a strong laser plasma

International Nuclear Information System (INIS)

Zhu Shitong; Shen Wenda; Guo Qizhi

1993-01-01

The wave function of free electron in a strong laser plasma is obtained by solving exactly the Dirac equation in a curved space-time with optical metric for the laser plasma. When the laser field is diminished to zero, the wave function is naturally reduced to relativistic wave function of free electron. The possible application of the wave function is discussed

Precise analytic approximations for the Bessel function J1 (x)

Science.gov (United States)

Maass, Fernando; Martin, Pablo

2018-03-01

Precise and straightforward analytic approximations for the Bessel function J1 (x) have been found. Power series and asymptotic expansions have been used to determine the parameters of the approximation, which is as a bridge between both expansions, and it is a combination of rational and trigonometric functions multiplied with fractional powers of x. Here, several improvements with respect to the so called Multipoint Quasirational Approximation technique have been performed. Two procedures have been used to determine the parameters of the approximations. The maximum absolute errors are in both cases smaller than 0.01. The zeros of the approximation are also very precise with less than 0.04 per cent for the first one. A second approximation has been also determined using two more parameters, and in this way the accuracy has been increased to less than 0.001.

Analytic approximation for the modified Bessel function I -2/3(x)

Science.gov (United States)

Martin, Pablo; Olivares, Jorge; Maass, Fernando

2017-12-01

In the present work an analytic approximation to modified Bessel function of negative fractional order I -2/3(x) is presented. The validity of the approximation is for every positive value of the independent variable. The accuracy is high in spite of the small number (4) of parameters used. The approximation is a combination of elementary functions with rational ones. Power series and assymptotic expansions are simultaneously used to obtain the approximation.

On quantum mechanical phase-space wave functions

DEFF Research Database (Denmark)

Wlodarz, Joachim J.

1994-01-01

An approach to quantum mechanics based on the notion of a phase-space wave function is proposed within the Weyl-Wigner-Moyal representation. It is shown that the Schrodinger equation for the phase-space wave function is equivalent to the quantum Liouville equation for the Wigner distribution...... function. The relationship to the recent results by Torres-Vega and Frederick [J. Chem. Phys. 98, 3103 (1993)] is also discussed....

Computational resources to filter gravitational wave data with P-approximant templates

International Nuclear Information System (INIS)

Porter, Edward K

2002-01-01

The prior knowledge of the gravitational waveform from compact binary systems makes matched filtering an attractive detection strategy. This detection method involves the filtering of the detector output with a set of theoretical waveforms or templates. One of the most important factors in this strategy is knowing how many templates are needed in order to reduce the loss of possible signals. In this study, we calculate the number of templates and computational power needed for a one-step search for gravitational waves from inspiralling binary systems. We build on previous works by first expanding the post-Newtonian waveforms to 2.5-PN order and second, for the first time, calculating the number of templates needed when using P-approximant waveforms. The analysis is carried out for the four main first-generation interferometers, LIGO, GEO600, VIRGO and TAMA. As well as template number, we also calculate the computational cost of generating banks of templates for filtering GW data. We carry out the calculations for two initial conditions. In the first case we assume a minimum individual mass of 1 M o-dot and in the second, we assume a minimum individual mass of 5 M o-dot . We find that, in general, we need more P-approximant templates to carry out a search than if we use standard PN templates. This increase varies according to the order of PN-approximation, but can be as high as a factor of 3 and is explained by the smaller span of the P-approximant templates as we go to higher masses. The promising outcome is that for 2-PN templates, the increase is small and is outweighed by the known robustness of the 2-PN P-approximant templates

Integral approximants for functions of higher monodromic dimension

Energy Technology Data Exchange (ETDEWEB)

Baker, G.A. Jr.

1987-01-01

In addition to the description of multiform, locally analytic functions as covering a many sheeted version of the complex plane, Riemann also introduced the notion of considering them as describing a space whose ''monodromic'' dimension is the number of linearly independent coverings by the monogenic analytic function at each point of the complex plane. I suggest that this latter concept is natural for integral approximants (sub-class of Hermite-Pade approximants) and discuss results for both ''horizontal'' and ''diagonal'' sequences of approximants. Some theorems are now available in both cases and make clear the natural domain of convergence of the horizontal sequences is a disk centered on the origin and that of the diagonal sequences is a suitably cut complex-plane together with its identically cut pendant Riemann sheets. Some numerical examples have also been computed.

A constrained approximation for nuclear barrier penetration and fission

International Nuclear Information System (INIS)

Tang, H.H.K.; Negele, J.W.; Massachusetts Inst. of Tech., Cambridge; Massachusetts Inst. of Tech., Cambridge

1983-01-01

An approximation to the time-dependent mean-field theory for barrier penetration by a nucleus is obtained in terms of constrained Hartree-Fock wave functions and a coherent velocity field. A discrete approximation to the continuum theory suitable for practical numerical calculations is presented and applied to three illustrative models. Potential application of the theory to the study of nuclear fission is discussed. (orig.)

Approximate convex hull of affine iterated function system attractors

International Nuclear Information System (INIS)

Mishkinis, Anton; Gentil, Christian; Lanquetin, Sandrine; Sokolov, Dmitry

2012-01-01

Highlights: ► We present an iterative algorithm to approximate affine IFS attractor convex hull. ► Elimination of the interior points significantly reduces the complexity. ► To optimize calculations, we merge the convex hull images at each iteration. ► Approximation by ellipses increases speed of convergence to the exact convex hull. ► We present a method of the output convex hull simplification. - Abstract: In this paper, we present an algorithm to construct an approximate convex hull of the attractors of an affine iterated function system (IFS). We construct a sequence of convex hull approximations for any required precision using the self-similarity property of the attractor in order to optimize calculations. Due to the affine properties of IFS transformations, the number of points considered in the construction is reduced. The time complexity of our algorithm is a linear function of the number of iterations and the number of points in the output approximate convex hull. The number of iterations and the execution time increases logarithmically with increasing accuracy. In addition, we introduce a method to simplify the approximate convex hull without loss of accuracy.

HIGHLY PRECISE APPROXIMATION OF FREE SURFACE GREEN FUNCTION AND ITS HIGH ORDER DERIVATIVES BASED ON REFINED SUBDOMAINS

Directory of Open Access Journals (Sweden)

Jiameng Wu

2018-01-01

Full Text Available The infinite depth free surface Green function (GF and its high order derivatives for diffraction and radiation of water waves are considered. Especially second order derivatives are essential requirements in high-order panel method. In this paper, concerning the classical representation, composed of a semi-infinite integral involving a Bessel function and a Cauchy singularity, not only the GF and its first order derivatives but also second order derivatives are derived from four kinds of analytical series expansion and refined division of whole calculation domain. The approximations of special functions, particularly the hypergeometric function and the algorithmic applicability with different subdomains are implemented. As a result, the computation accuracy can reach 10-9 in whole domain compared with conventional methods based on direct numerical integration. Furthermore, numerical efficiency is almost equivalent to that with the classical method.

Approximation Of Multi-Valued Inverse Functions Using Clustering And Sugeno Fuzzy Inference

Science.gov (United States)

Walden, Maria A.; Bikdash, Marwan; Homaifar, Abdollah

1998-01-01

Finding the inverse of a continuous function can be challenging and computationally expensive when the inverse function is multi-valued. Difficulties may be compounded when the function itself is difficult to evaluate. We show that we can use fuzzy-logic approximators such as Sugeno inference systems to compute the inverse on-line. To do so, a fuzzy clustering algorithm can be used in conjunction with a discriminating function to split the function data into branches for the different values of the forward function. These data sets are then fed into a recursive least-squares learning algorithm that finds the proper coefficients of the Sugeno approximators; each Sugeno approximator finds one value of the inverse function. Discussions about the accuracy of the approximation will be included.

Neutrino wave function and oscillation suppression

International Nuclear Information System (INIS)

Dolgov, A.D.; Lychkovskiy, O.V.; Mamonov, A.A.; Okun, L.B.; Schepkin, M.G.

2005-01-01

We consider a thought experiment, in which a neutrino is produced by an electron on a nucleus in a crystal. The wave function of the oscillating neutrino is calculated assuming that the electron is described by a wave packet. If the electron is relativistic and the spatial size of its wave packet is much larger than the size of the crystal cell, then the wave packet of the produced neutrino has essentially the same size as the wave packet of the electron. We investigate the suppression of neutrino oscillations at large distances caused by two mechanisms: (1) spatial separation of wave packets corresponding to different neutrino masses; (2) neutrino energy dispersion for given neutrino mass eigenstates. We resolve the contributions of these two mechanisms. (orig.)

Comparative study on spreading function for directional wave spectra

Digital Repository Service at National Institute of Oceanography (India)

Bhat, S.S.; Anand, N.M.; Nayak, B.U.

-dimensional wave energy S(f) and the directional spreading function D(f, theta). This paper reviews various spreading functions proposed in the past for estimating the directional wave energy and presents their application to the Indian wave condition. It is found...

Simultaneous projection of particle-number and angular momentum BCS wave-functions in the rare-earth nuclei

International Nuclear Information System (INIS)

Oudih, M.R.; Benhamouda, N.; Fellah, M.; Allal, N.H.

2000-01-01

A method of simultaneous particle-number and angular-momentum projection of the BCS wave-function is presented. The particle number projection method is of FBCS type. In the frame work of the adiabatic approximation, the rotational energies of the axially symmetric even-even nuclei are established and numerically calculated for the rare-earth region. (author)

Second-order Born approximation for the ionization of molecules by electron and positron impact

Energy Technology Data Exchange (ETDEWEB)

Dal Cappello, C. [Universite Paul Verlaine-Metz, Laboratoire de Physique Moleculaire et des Collisions, Institut Jean Barriol (FR2843), 1 Boulevard Arago, F-57078 Metz Cedex 3 (France); Rezkallah, Z.; Houamer, S. [Laboratoire de Physique Quantique et Systemes Dynamiques, Departement de Physique, Faculte des Sciences Universite Ferhat Abbas, Setif 19000 (Algeria); Charpentier, I. [Universite Paul Verlaine-Metz, Laboratoire de Physique et Mecanique des Materiaux UMR 7554, Ile du Saulcy, F-57045 Metz Cedex 1 (France); Hervieux, P. A. [Institut de Physique et Chimie des Materiaux de Strasbourg, 23 Rue du Loess, BP 43, F-67034 Strasbourg Cedex 2 (France); Ruiz-Lopez, M. F. [Nancy-University, Equipe de Chimie et Biochimie Theoriques, UMR CNRS-UHP 7565, BP 239, F-54506 Vandoeuvre-les-Nancy (France); Dey, R. [Max-Planck Institut fuer Plasmaphysik, Boltzmannstr. 2, D-85748 Garching (Germany); Roy, A. C. [School of Mathematical Sciences, Ramakrishna Mission Vivekananda University, Belur Math 711202, West Bengal (India)

2011-09-15

Second-order Born approximation is applied to study the ionization of molecules. The initial and final states are described by single-center wave functions. For the initial state a Gaussian wave function is used while for the ejected electron it is a distorted wave. Results of the present model are compared with recent (e,2e) experiments on the water molecule. Preliminary results are also presented for the ionization of the thymine molecule by electrons and positrons.

Variational random phase approximation for the anharmonic oscillator

International Nuclear Information System (INIS)

Dukelsky, J.; Schuck, P.

1990-04-01

The recently derived Variational Random Phase Approximation is examined using the anharmonic oscillator model. Special attention is paid to the ground state RPA wave function and the convergence of the proposed truncation scheme to obtain the diagonal density matrix. Comparison with the standard Coupled Cluster method is made

Time adaptivity in the diffusive wave approximation to the shallow water equations

KAUST Repository

Collier, Nathan; Radwan, Hany; Dalcí n, Lisandro D.; Calo, Victor M.

2013-01-01

We discuss the use of time adaptivity applied to the one dimensional diffusive wave approximation to the shallow water equations. A simple and computationally economical error estimator is discussed which enables time-step size adaptivity. This robust adaptive time discretization corrects the initial time step size to achieve a user specified bound on the discretization error and allows time step size variations of several orders of magnitude. In particular, the one dimensional results presented in this work feature a change of four orders of magnitudes for the time step over the entire simulation. © 2011 Elsevier B.V.

Time adaptivity in the diffusive wave approximation to the shallow water equations

KAUST Repository

Collier, Nathan

2013-05-01

We discuss the use of time adaptivity applied to the one dimensional diffusive wave approximation to the shallow water equations. A simple and computationally economical error estimator is discussed which enables time-step size adaptivity. This robust adaptive time discretization corrects the initial time step size to achieve a user specified bound on the discretization error and allows time step size variations of several orders of magnitude. In particular, the one dimensional results presented in this work feature a change of four orders of magnitudes for the time step over the entire simulation. © 2011 Elsevier B.V.

Efficient approximation of black-box functions and Pareto sets

NARCIS (Netherlands)

Rennen, G.

2009-01-01

In the case of time-consuming simulation models or other so-called black-box functions, we determine a metamodel which approximates the relation between the input- and output-variables of the simulation model. To solve multi-objective optimization problems, we approximate the Pareto set, i.e. the

Semiclassical multicomponent wave function

NARCIS (Netherlands)

Mostovoy, M.V.

A consistent method for obtaining the semiclassical multicomponent wave function for any value of adiabatic parameter is discussed and illustrated by examining the motion of a neutral particle in a nonuniform magnetic field. The method generalizes the Bohr-Sommerfeld quantization rule to

A new way of obtaining analytic approximations of Chandrasekhar's H function

International Nuclear Information System (INIS)

Vukanic, J.; Arsenovic, D.; Davidovic, D.

2007-01-01

Applying the mean value theorem for definite integrals in the non-linear integral equation for Chandrasekhar's H function describing conservative isotropic scattering, we have derived a new, simple analytic approximation for it, with a maximal relative error below 2.5%. With this new function as a starting-point, after a single iteration in the corresponding integral equation, we have obtained a new, highly accurate analytic approximation for the H function. As its maximal relative error is below 0.07%, it significantly surpasses the accuracy of other analytic approximations

Square well approximation to the optical potential

International Nuclear Information System (INIS)

Jain, A.K.; Gupta, M.C.; Marwadi, P.R.

1976-01-01

Approximations for obtaining T-matrix elements for a sum of several potentials in terms of T-matrices for individual potentials are studied. Based on model calculations for S-wave for a sum of two separable non-local potentials of Yukawa type form factors and a sum of two delta function potentials, it is shown that the T-matrix for a sum of several potentials can be approximated satisfactorily over all the energy regions by the sum of T-matrices for individual potentials. Based on this, an approximate method for finding T-matrix for any local potential by approximating it by a sum of suitable number of square wells is presented. This provides an interesting way to calculate the T-matrix for any arbitary potential in terms of Bessel functions to a good degree of accuracy. The method is applied to the Saxon-Wood potentials and good agreement with exact results is found. (author)

On Approximate Solutions of Functional Equations in Vector Lattices

Directory of Open Access Journals (Sweden)

Bogdan Batko

2014-01-01

Full Text Available We provide a method of approximation of approximate solutions of functional equations in the class of functions acting into a Riesz space (algebra. The main aim of the paper is to provide a general theorem that can act as a tool applicable to a possibly wide class of functional equations. The idea is based on the use of the Spectral Representation Theory for Riesz spaces. The main result will be applied to prove the stability of an alternative Cauchy functional equation F(x+y+F(x+F(y≠0⇒F(x+y=F(x+F(y in Riesz spaces, the Cauchy equation with squares F(x+y2=(F(x+F(y2 in f-algebras, and the quadratic functional equation F(x+y+F(x-y=2F(x+2F(y in Riesz spaces.

Simulating propagation of decomposed elastic waves using low-rank approximate mixed-domain integral operators for heterogeneous transversely isotropic media

KAUST Repository

Cheng, Jiubing

2014-08-05

In elastic imaging, the extrapolated vector fields are decomposed into pure wave modes, such that the imaging condition produces interpretable images, which characterize reflectivity of different reflection types. Conventionally, wavefield decomposition in anisotropic media is costly as the operators involved is dependent on the velocity, and thus not stationary. In this abstract, we propose an efficient approach to directly extrapolate the decomposed elastic waves using lowrank approximate mixed space/wavenumber domain integral operators for heterogeneous transverse isotropic (TI) media. The low-rank approximation is, thus, applied to the pseudospectral extrapolation and decomposition at the same time. The pseudo-spectral implementation also allows for relatively large time steps in which the low-rank approximation is applied. Synthetic examples show that it can yield dispersionfree extrapolation of the decomposed quasi-P (qP) and quasi- SV (qSV) modes, which can be used for imaging, as well as the total elastic wavefields.

Compact two-electron wave function for bond dissociation and Van der Waals interactions: a natural amplitude assessment.

Science.gov (United States)

Giesbertz, Klaas J H; van Leeuwen, Robert

2014-05-14

Electron correlations in molecules can be divided in short range dynamical correlations, long range Van der Waals type interactions, and near degeneracy static correlations. In this work, we analyze for a one-dimensional model of a two-electron system how these three types of correlations can be incorporated in a simple wave function of restricted functional form consisting of an orbital product multiplied by a single correlation function f (r12) depending on the interelectronic distance r12. Since the three types of correlations mentioned lead to different signatures in terms of the natural orbital (NO) amplitudes in two-electron systems, we make an analysis of the wave function in terms of the NO amplitudes for a model system of a diatomic molecule. In our numerical implementation, we fully optimize the orbitals and the correlation function on a spatial grid without restrictions on their functional form. Due to this particular form of the wave function, we can prove that none of the amplitudes vanishes and moreover that it displays a distinct sign pattern and a series of avoided crossings as a function of the bond distance in agreement with the exact solution. This shows that the wave function ansatz correctly incorporates the long range Van der Waals interactions. We further show that the approximate wave function gives an excellent binding curve and is able to describe static correlations. We show that in order to do this the correlation function f (r12) needs to diverge for large r12 at large internuclear distances while for shorter bond distances it increases as a function of r12 to a maximum value after which it decays exponentially. We further give a physical interpretation of this behavior.

On the construction of translationally invariant deformed wave functions

International Nuclear Information System (INIS)

Guardiola, R.

1975-01-01

Translationally invariant nuclear wave functions are constructed from deformed harmonic oscillator shell-model wave functions, with an exact projection of angular momentum quantum numbers. It is shown that the computation of matrix elements with the translationally invariant wave functions is as simple as the standard calculation, and formulae are obtained for (i) the potential energy, (ii) the kinetic energy and rms radius, and (iii) the charge form factor. (Auth.)

Application of a simplified calculation for full-wave microtremor H/ V spectral ratio based on the diffuse field approximation to identify underground velocity structures

Science.gov (United States)

Wu, Hao; Masaki, Kazuaki; Irikura, Kojiro; Sánchez-Sesma, Francisco José

2017-12-01

Under the diffuse field approximation, the full-wave (FW) microtremor H/ V spectral ratio ( H/ V) is modeled as the square root of the ratio of the sum of imaginary parts of the Green's function of the horizontal components to that of the vertical one. For a given layered medium, the FW H/ V can be well approximated with only surface waves (SW) H/ V of the "cap-layered" medium which consists of the given layered medium and a new larger velocity half-space (cap layer) at large depth. Because the contribution of surface waves can be simply obtained by the residue theorem, the computation of SW H/ V of cap-layered medium is faster than that of FW H/ V evaluated by discrete wavenumber method and contour integration method. The simplified computation of SW H/ V was then applied to identify the underground velocity structures at six KiK-net strong-motion stations. The inverted underground velocity structures were used to evaluate FW H/ Vs which were consistent with the SW H/ Vs of corresponding cap-layered media. The previous study on surface waves H/ Vs proposed with the distributed surface sources assumption and a fixed Rayleigh-to-Love waves amplitude ratio for horizontal motions showed a good agreement with the SW H/ Vs of our study. The consistency between observed and theoretical spectral ratios, such as the earthquake motions of H/ V spectral ratio and spectral ratio of horizontal motions between surface and bottom of borehole, indicated that the underground velocity structures identified from SW H/ V of cap-layered medium were well resolved by the new method.[Figure not available: see fulltext.

Structure and representation of correlation functions and the density matrix for a statistical wave field in optics

International Nuclear Information System (INIS)

Sudarshan, E.C.G.; Mukunda, N.

1978-03-01

A systematic structure analysis of the correlation functions of statistical quantum optics is carried out. From a suitably defined auxiliary two-point function identification of the excited modes in the wave field is found. The relative simplicity of the higher order correlation functions emerges as a by-product and the conditions under which these are made pure are derived. These results depend in a crucial manner on the notion of coherence indices aand of unimodular coherence indices. A new class of approximate expressions for the density operator of a statistical wave field is worked out based on discrete characteristic sets. These are even more economical than the diagonal coherent state representations. An appreciation of the subtleties of quantum theory obtains. Certain implications for the physics of light beams are cited. 28 references

Wave-function reconstruction in a graded semiconductor superlattice

DEFF Research Database (Denmark)

Lyssenko, V. G.; Hvam, Jørn Märcher; Meinhold, D.

2004-01-01

We reconstruct a test wave function in a strongly coupled, graded well-width superlattice by resolving the spatial extension of the interband polarisation and deducing the wave function employing non-linear optical spectroscopy. The graded gap superlattice allows us to precisely control the dista...

Optimized implementations of rational approximations for the Voigt and complex error function

International Nuclear Information System (INIS)

Schreier, Franz

2011-01-01

Rational functions are frequently used as efficient yet accurate numerical approximations for real and complex valued functions. For the complex error function w(x+iy), whose real part is the Voigt function K(x,y), code optimizations of rational approximations are investigated. An assessment of requirements for atmospheric radiative transfer modeling indicates a y range over many orders of magnitude and accuracy better than 10 -4 . Following a brief survey of complex error function algorithms in general and rational function approximations in particular the problems associated with subdivisions of the x, y plane (i.e., conditional branches in the code) are discussed and practical aspects of Fortran and Python implementations are considered. Benchmark tests of a variety of algorithms demonstrate that programming language, compiler choice, and implementation details influence computational speed and there is no unique ranking of algorithms. A new implementation, based on subdivision of the upper half-plane in only two regions, combining Weideman's rational approximation for small |x|+y<15 and Humlicek's rational approximation otherwise is shown to be efficient and accurate for all x, y.

Light-front wave function of composite system with spin

International Nuclear Information System (INIS)

Karmanov, V.A.

1979-01-01

The method to construct the relativistic wave function with spin on the light front is developed. The spin structure of the deuteron wave function in relativistic range is found. The calculation methods are illustrated by the calculation of elastic pd-scattering cross section. The consideration carried out is equivalent to the solution of the problem of taking into account the spins and angular momenta in the parton wave functions in the infinite momentum frame

Multidimensional stochastic approximation using locally contractive functions

Science.gov (United States)

Lawton, W. M.

1975-01-01

A Robbins-Monro type multidimensional stochastic approximation algorithm which converges in mean square and with probability one to the fixed point of a locally contractive regression function is developed. The algorithm is applied to obtain maximum likelihood estimates of the parameters for a mixture of multivariate normal distributions.

Hadronic wave functions and high momentum transfer interactions in quantum chromodynamics

International Nuclear Information System (INIS)

Brodsky, S.J.; Huang, T.; Lepage, G.P.

1983-01-01

This chapter emphasizes the utility of a Fock state representation of the meson and baryon wave functions as a means not only to parametrize the effects of bound state dynamics in QCD phenomena, but also to interrelate exclusive, inclusive, and higher twist processes. Discusses hadronic wave functions in QCD, measures of hadronic wave functions (form factors of composite systems, form factors of mesons, the meson distribution amplitude); large momentum transfer exclusive processes (two-photon processes); deep inelastic lepton scattering; and the phenomenology of hadronic wave functions (measures of hadron wave functions, constraints on the pion and proton valence wave function, quark jet diffraction excitation, the ''unveiling'' of the hadronic wave function and intrinsic charm). Finds that the testing ground of perturbative QCD where rigorous, definitive tests of the theory can be made can now be extended throughout a large domain of large momentum transfer exclusive and inclusive lepton, photon, and hadron reactions

Long-range-corrected Rung 3.5 density functional approximations

Science.gov (United States)

Janesko, Benjamin G.; Proynov, Emil; Scalmani, Giovanni; Frisch, Michael J.

2018-03-01

Rung 3.5 functionals are a new class of approximations for density functional theory. They provide a flexible intermediate between exact (Hartree-Fock, HF) exchange and semilocal approximations for exchange. Existing Rung 3.5 functionals inherit semilocal functionals' limitations in atomic cores and density tails. Here we address those limitations using range-separated admixture of HF exchange. We present three new functionals. LRC-ωΠLDA combines long-range HF exchange with short-range Rung 3.5 ΠLDA exchange. SLC-ΠLDA combines short- and long-range HF exchange with middle-range ΠLDA exchange. LRC-ωΠLDA-AC incorporates a combination of HF, semilocal, and Rung 3.5 exchange in the short range, based on an adiabatic connection. We test these in a new Rung 3.5 implementation including up to analytic fourth derivatives. LRC-ωΠLDA and SLC-ΠLDA improve atomization energies and reaction barriers by a factor of 8 compared to the full-range ΠLDA. LRC-ωΠLDA-AC brings further improvement approaching the accuracy of standard long-range corrected schemes LC-ωPBE and SLC-PBE. The new functionals yield highest occupied orbital energies closer to experimental ionization potentials and describe correctly the weak charge-transfer complex of ethylene and dichlorine and the hole-spin distribution created by an Al defect in quartz. This study provides a framework for more flexible range-separated Rung 3.5 approximations.

Analysis of a finite PML approximation to the three dimensional elastic wave scattering problem

KAUST Repository

Bramble, James H.

2010-01-01

We consider the application of a perfectly matched layer (PML) technique to approximate solutions to the elastic wave scattering problem in the frequency domain. The PML is viewed as a complex coordinate shift in spherical coordinates which leads to a variable complex coefficient equation for the displacement vector posed on an infinite domain (the complement of the scatterer). The rapid decay of the PML solution suggests truncation to a bounded domain with a convenient outer boundary condition and subsequent finite element approximation (for the truncated problem). We prove existence and uniqueness of the solutions to the infinite domain and truncated domain PML equations (provided that the truncated domain is sufficiently large). We also show exponential convergence of the solution of the truncated PML problem to the solution of the original scattering problem in the region of interest. We then analyze a Galerkin numerical approximation to the truncated PML problem and prove that it is well posed provided that the PML damping parameter and mesh size are small enough. Finally, computational results illustrating the efficiency of the finite element PML approximation are presented. © 2010 American Mathematical Society.

An Approximate Proximal Bundle Method to Minimize a Class of Maximum Eigenvalue Functions

Directory of Open Access Journals (Sweden)

Wei Wang

2014-01-01

Full Text Available We present an approximate nonsmooth algorithm to solve a minimization problem, in which the objective function is the sum of a maximum eigenvalue function of matrices and a convex function. The essential idea to solve the optimization problem in this paper is similar to the thought of proximal bundle method, but the difference is that we choose approximate subgradient and function value to construct approximate cutting-plane model to solve the above mentioned problem. An important advantage of the approximate cutting-plane model for objective function is that it is more stable than cutting-plane model. In addition, the approximate proximal bundle method algorithm can be given. Furthermore, the sequences generated by the algorithm converge to the optimal solution of the original problem.

Approximate analytical solutions of Klein-Gordon equation with Hulthen potentials for nonzero angular momentum

International Nuclear Information System (INIS)

Chen Changyuan; Sun Dongsheng; Lu Falin

2007-01-01

Using the exponential function transformation approach along with an approximation for the centrifugal potential, the radial Klein-Gordon equation with the vector and scalar Hulthen potential is transformed to a hypergeometric differential equation. The approximate analytical solutions of bound states are attained for different l. The analytical energy equation and the unnormalized radial wave functions expressed in terms of hypergeometric polynomials are given

Efficient and Flexible Computation of Many-Electron Wave Function Overlaps.

Science.gov (United States)

Plasser, Felix; Ruckenbauer, Matthias; Mai, Sebastian; Oppel, Markus; Marquetand, Philipp; González, Leticia

2016-03-08

A new algorithm for the computation of the overlap between many-electron wave functions is described. This algorithm allows for the extensive use of recurring intermediates and thus provides high computational efficiency. Because of the general formalism employed, overlaps can be computed for varying wave function types, molecular orbitals, basis sets, and molecular geometries. This paves the way for efficiently computing nonadiabatic interaction terms for dynamics simulations. In addition, other application areas can be envisaged, such as the comparison of wave functions constructed at different levels of theory. Aside from explaining the algorithm and evaluating the performance, a detailed analysis of the numerical stability of wave function overlaps is carried out, and strategies for overcoming potential severe pitfalls due to displaced atoms and truncated wave functions are presented.

Pade approximants for entire functions with regularly decreasing Taylor coefficients

International Nuclear Information System (INIS)

Rusak, V N; Starovoitov, A P

2002-01-01

For a class of entire functions the asymptotic behaviour of the Hadamard determinants D n,m as 0≤m≤m(n)→∞ and n→∞ is described. This enables one to study the behaviour of parabolic sequences from Pade and Chebyshev tables for many individual entire functions. The central result of the paper is as follows: for some sequences {(n,m(n))} in certain classes of entire functions (with regular Taylor coefficients) the Pade approximants {π n,m(n) }, which provide the locally best possible rational approximations, converge to the given function uniformly on the compact set D={z:|z|≤1} with asymptotically best rate

Bohmian Conditional Wave Functions (and the status of the quantum state)

International Nuclear Information System (INIS)

Norsen, Travis

2016-01-01

The de Broglie - Bohm pilot-wave theory - uniquely among realistic candidate quantum theories - allows a straightforward and simple definition of the wave function of a subsystem of some larger system (such as the entire universe). Such sub-system wave functions are called “Conditional Wave Functions” (CWFs). Here we explain this concept and indicate the CWF's role in the Bohmian explanation of the usual quantum formalism, and then develop (and motivate) the more speculative idea that something like single-particle wave functions could replace the (ontologically problematical) universal wave function in some future, empirically adequate, pilot-wave-type theory. Throughout the presentation is pedagogical and points are illustrated with simple toy models. (paper)

Quantifying Diastolic Function: From E-Waves as Triangles to Physiologic Contours via the 'Geometric Method'.

Science.gov (United States)

Golman, Mikhail; Padovano, William; Shmuylovich, Leonid; Kovács, Sándor J

2018-03-01

Conventional echocardiographic diastolic function (DF) assessment approximates transmitral flow velocity contours (Doppler E-waves) as triangles, with peak (E peak ), acceleration time (AT), and deceleration time (DT) as indexes. These metrics have limited value because they are unable to characterize the underlying physiology. The parametrized diastolic filling (PDF) formalism provides a physiologic, kinematic mechanism based characterization of DF by extracting chamber stiffness (k), relaxation (c), and load (x o ) from E-wave contours. We derive the mathematical relationship between the PDF parameters and E peak , AT, DT and thereby introduce the geometric method (GM) that computes the PDF parameters using E peak , AT, and DT as input. Numerical experiments validated GM by analysis of 208 E-waves from 31 datasets spanning the full range of clinical diastolic function. GM yielded indistinguishable average parameter values per subject vs. the gold-standard PDF method (k: R 2  = 0.94, c: R 2  = 0.95, x o : R 2  = 0.95, p PDF method to quantify DF in terms of physiologic chamber properties.

Bessel harmonic analysis and approximation of functions on the half-line

International Nuclear Information System (INIS)

Platonov, Sergei S

2007-01-01

We study problems of approximation of functions on [0,+∞) in the metric of L p with power weight using generalized Bessel shifts. We prove analogues of direct Jackson theorems for the modulus of smoothness of arbitrary order defined in terms of generalized Bessel shifts. We establish the equivalence of the modulus of smoothness and the K-functional. We define function spaces of Nikol'skii-Besov type and describe them in terms of best approximations. As a tool for approximation, we use a certain class of entire functions of exponential type. In this class, we prove analogues of Bernstein's inequality and others for the Bessel differential operator and its fractional powers. The main tool we use to solve these problems is Bessel harmonic analysis

Approximation methods for the partition functions of anharmonic systems

International Nuclear Information System (INIS)

Lew, P.; Ishida, T.

1979-07-01

The analytical approximations for the classical, quantum mechanical and reduced partition functions of the diatomic molecule oscillating internally under the influence of the Morse potential have been derived and their convergences have been tested numerically. This successful analytical method is used in the treatment of anharmonic systems. Using Schwinger perturbation method in the framework of second quantization formulism, the reduced partition function of polyatomic systems can be put into an expression which consists separately of contributions from the harmonic terms, Morse potential correction terms and interaction terms due to the off-diagonal potential coefficients. The calculated results of the reduced partition function from the approximation method on the 2-D and 3-D model systems agree well with the numerical exact calculations

Multiquark masses and wave functions through modified Green's function Monte Carlo method

International Nuclear Information System (INIS)

Kerbikov, B.O.; Polikarpov, M.I.; Shevchenko, L.V.

1987-01-01

The Modified Green's function Monte Carlo method (MGFMC) is used to calculate the masses and ground-state wave functions of multiquark systems in the potential model. The previously developed MGFMC is generalized in order to treat systems containing quarks with inequal masses. The obtained results are presented with the Cornell potential for the masses and the wave functions of light and heavy flavoured baryons and multiquark states (N=6, 9, 12) made of light quarks

Approximation of the Doppler broadening function by Frobenius method

International Nuclear Information System (INIS)

Palma, Daniel A.P.; Martinez, Aquilino S.; Silva, Fernando C.

2005-01-01

An analytical approximation of the Doppler broadening function ψ(x,ξ) is proposed. This approximation is based on the solution of the differential equation for ψ(x,ξ) using the methods of Frobenius and the parameters variation. The analytical form derived for ψ(x,ξ) in terms of elementary functions is very simple and precise. It can be useful for applications related to the treatment of nuclear resonances mainly for the calculations of multigroup parameters and self-protection factors of the resonances, being the last used to correct microscopic cross-sections measurements by the activation technique. (author)

Wave function of the Universe as a leaking system

International Nuclear Information System (INIS)

Suen, W.; Young, K.

1989-01-01

We propose a path-integral formulation for the wave function of the Universe which requires neither the Euclidean nor the conformal rotation. The boundary condition is taken to be that ''all possible boundaries are included.'' The resulting wave function in a simple model is shown to have the following properties: (i) the wave function tends to zero as the scale factor of the Universe tends to zero; (ii) in the semiclassical regime, it contains only the expanding component; (iii) it favors inflation

Wave function for harmonically confined electrons in time-dependent electric and magnetostatic fields.

Science.gov (United States)

Zhu, Hong-Ming; Chen, Jin-Wang; Pan, Xiao-Yin; Sahni, Viraht

2014-01-14

We derive via the interaction "representation" the many-body wave function for harmonically confined electrons in the presence of a magnetostatic field and perturbed by a spatially homogeneous time-dependent electric field-the Generalized Kohn Theorem (GKT) wave function. In the absence of the harmonic confinement - the uniform electron gas - the GKT wave function reduces to the Kohn Theorem wave function. Without the magnetostatic field, the GKT wave function is the Harmonic Potential Theorem wave function. We further prove the validity of the connection between the GKT wave function derived and the system in an accelerated frame of reference. Finally, we provide examples of the application of the GKT wave function.

Effect of Forcing Function on Nonlinear Acoustic Standing Waves

Science.gov (United States)

Finkheiner, Joshua R.; Li, Xiao-Fan; Raman, Ganesh; Daniels, Chris; Steinetz, Bruce

2003-01-01

Nonlinear acoustic standing waves of high amplitude have been demonstrated by utilizing the effects of resonator shape to prevent the pressure waves from entering saturation. Experimentally, nonlinear acoustic standing waves have been generated by shaking an entire resonating cavity. While this promotes more efficient energy transfer than a piston-driven resonator, it also introduces complicated structural dynamics into the system. Experiments have shown that these dynamics result in resonator forcing functions comprised of a sum of several Fourier modes. However, previous numerical studies of the acoustics generated within the resonator assumed simple sinusoidal waves as the driving force. Using a previously developed numerical code, this paper demonstrates the effects of using a forcing function constructed with a series of harmonic sinusoidal waves on resonating cavities. From these results, a method will be demonstrated which allows the direct numerical analysis of experimentally generated nonlinear acoustic waves in resonators driven by harmonic forcing functions.

Corrected Fourier series and its application to function approximation

Directory of Open Access Journals (Sweden)

Qing-Hua Zhang

2005-01-01

Full Text Available Any quasismooth function f(x in a finite interval [0,x0], which has only a finite number of finite discontinuities and has only a finite number of extremes, can be approximated by a uniformly convergent Fourier series and a correction function. The correction function consists of algebraic polynomials and Heaviside step functions and is required by the aperiodicity at the endpoints (i.e., f(0≠f(x0 and the finite discontinuities in between. The uniformly convergent Fourier series and the correction function are collectively referred to as the corrected Fourier series. We prove that in order for the mth derivative of the Fourier series to be uniformly convergent, the order of the polynomial need not exceed (m+1. In other words, including the no-more-than-(m+1 polynomial has eliminated the Gibbs phenomenon of the Fourier series until its mth derivative. The corrected Fourier series is then applied to function approximation; the procedures to determine the coefficients of the corrected Fourier series are illustrated in detail using examples.

Approximate angular momentum projection from cranked intrinsic states

International Nuclear Information System (INIS)

Goodman, A.L.

1979-01-01

High-spin spectra are determined by approximately projecting states of good angular momentum from cranked Hartree-Fock-Bogoliubov (CHFB) wave functions. For each J the projected energy is E/sub PROJ/ approx. = E/sub CHFB/ - (ΔJ) 2 /2 J/sub CHFB/, where the moment of inertia J and the fluctuation ΔJ are spin dependent. For /sup 168,170/Yb and 174 Hf the projected J is less than the CHFB value for all J. Consequently approximate projection increases all yrast excitation energies for these nuclei

Projected potential profiles across interfaces obtained by reconstructing the exit face wave function from through focal series

International Nuclear Information System (INIS)

Bhattacharyya, Somnath; Koch, Christoph T.; Ruehle, Manfred

2006-01-01

An iterative method for reconstructing the exit face wave function from a through focal series of transmission electron microscopy image line profiles across an interface is presented. Apart from high-resolution images recorded with small changes in defocus, this method works also well for a large defocus range as used for Fresnel imaging. Using the phase-object approximation the projected electrostatic as well as the absorptive potential profiles across an interface are determined from this exit face wave function. A new experimental image alignment procedure was developed in order to align images with large relative defocus shift. The performance of this procedure is shown to be superior to other image alignment procedures existing in the literature. The reconstruction method is applied to both simulated and experimental images

Measurement of light-cone wave functions by diffractive dissociation

Energy Technology Data Exchange (ETDEWEB)

Asheri, D. [Tel Aviv Univ., School of Physics and Astronomy, Sackler Faculty of Exact Science (Israel)

2005-07-01

The measurement of the pion light-cone wave function is revisited and results for the Gegenbauer coefficients are presented. Measurements of the photon electromagnetic and hadronic wave functions are described and results are presented. (authors)

An Approximate Redistributed Proximal Bundle Method with Inexact Data for Minimizing Nonsmooth Nonconvex Functions

Directory of Open Access Journals (Sweden)

Jie Shen

2015-01-01

Full Text Available We describe an extension of the redistributed technique form classical proximal bundle method to the inexact situation for minimizing nonsmooth nonconvex functions. The cutting-planes model we construct is not the approximation to the whole nonconvex function, but to the local convexification of the approximate objective function, and this kind of local convexification is modified dynamically in order to always yield nonnegative linearization errors. Since we only employ the approximate function values and approximate subgradients, theoretical convergence analysis shows that an approximate stationary point or some double approximate stationary point can be obtained under some mild conditions.

Finite rank separable approximation for Skyrme interactions: spin-isospin excitations

International Nuclear Information System (INIS)

Severyukhin, A.P.; Voronov, V.V.; Borzov, I.N.; Nguyen Van Giai

2012-01-01

A finite rank separable approximation for the quasiparticle random phase approximation with the Skyrme interactions is applied for the case of charge-exchange nuclear modes. The coupling between one- and two-phonon terms in the wave functions is taken into account. It has been shown that the approximation reproduces reasonably well the full charge-exchange RPA results for the spin-dipole resonances in 132 Sn. As an illustration of the method, the phonon-phonon coupling effect on the β-decay half-life of 78 Ni is considered

Complexity of Gaussian-Radial-Basis Networks Approximating Smooth Functions

Czech Academy of Sciences Publication Activity Database

Kainen, P.C.; Kůrková, Věra; Sanguineti, M.

2009-01-01

Roč. 25, č. 1 (2009), s. 63-74 ISSN 0885-064X R&D Projects: GA ČR GA201/08/1744 Institutional research plan: CEZ:AV0Z10300504 Keywords : Gaussian-radial-basis-function networks * rates of approximation * model complexity * variation norms * Bessel and Sobolev norms * tractability of approximation Subject RIV: IN - Informatics, Computer Science Impact factor: 1.227, year: 2009

Counterstreaming magnetized plasmas. II. Perpendicular wave propagation

International Nuclear Information System (INIS)

Tautz, R.C.; Schlickeiser, R.

2006-01-01

The properties of longitudinal and transverse oscillations in magnetized symmetric counterstreaming Maxwellian plasmas with equal thermal velocities for waves propagating perpendicular to the stream direction are investigated on the basis of Maxwell equations and the nonrelativistic Vlasov equation. With the constraint of vanishing particle flux in the stream direction, three distinct dispersion relations are known, which are the ordinary-wave mode, the Bernstein wave mode, and the extraordinary electromagnetic wave mode, where the latter two are only approximations. In this article, all three dispersion relations are evaluated for a counterstreaming Maxwellian distribution function in terms of the hypergeometric function 2 F 2 . The growth rates for the ordinary-wave mode are compared to earlier results by Bornatici and Lee [Phys. Fluids 13, 3007 (1970)], who derived approximate results, whereas in this article the exact dispersion relation is solved numerically. The original results are therefore improved and show differences of up to 21% to the results obtained in this article

Quantal density functional theory II. Approximation methods and applications

International Nuclear Information System (INIS)

Sahni, Viraht

2010-01-01

This book is on approximation methods and applications of Quantal Density Functional Theory (QDFT), a new local effective-potential-energy theory of electronic structure. What distinguishes the theory from traditional density functional theory is that the electron correlations due to the Pauli exclusion principle, Coulomb repulsion, and the correlation contribution to the kinetic energy -- the Correlation-Kinetic effects -- are separately and explicitly defined. As such it is possible to study each property of interest as a function of the different electron correlations. Approximations methods based on the incorporation of different electron correlations, as well as a many-body perturbation theory within the context of QDFT, are developed. The applications are to the few-electron inhomogeneous electron gas systems in atoms and molecules, as well as to the many-electron inhomogeneity at metallic surfaces. (orig.)

Wind wave source functions in opposing seas

KAUST Repository

Langodan, Sabique

2015-08-26

The Red Sea is a challenge for wave modeling because of its unique two opposed wave systems, forced by opposite winds and converging at its center. We investigate the different physical aspects of wave evolution and propagation in the convergence zone. The two opposing wave systems have similar amplitude and frequency, each driven by the action of its own wind. Wave patterns at the centre of the Red Sea, as derived from extensive tests and intercomparison between model and measured data, suggest that the currently available wave model source functions may not properly represent the evolution of the local fields that appear to be characterized by a less effective wind input and an enhanced white-capping. We propose and test a possible simple solution to improve the wave-model simulation under opposing winds and waves condition. This article is protected by copyright. All rights reserved.

Reduced density matrix functional theory via a wave function based approach

Energy Technology Data Exchange (ETDEWEB)

Schade, Robert; Bloechl, Peter [Institute for Theoretical Physics, Clausthal University of Technology, Clausthal (Germany); Pruschke, Thomas [Institute for Theoretical Physics, University of Goettingen, Goettingen (Germany)

2016-07-01

We propose a new method for the calculation of the electronic and atomic structure of correlated electron systems based on reduced density matrix functional theory (rDMFT). The density-matrix functional is evaluated on the fly using Levy's constrained search formalism. The present implementation rests on a local approximation of the interaction reminiscent to that of dynamical mean field theory (DMFT). We focus here on additional approximations to the exact density-matrix functional in the local approximation and evaluate their performance.

Optimized Perturbation Theory for Wave Functions of Quantum Systems

International Nuclear Information System (INIS)

Hatsuda, T.; Tanaka, T.; Kunihiro, T.

1997-01-01

The notion of the optimized perturbation, which has been successfully applied to energy eigenvalues, is generalized to treat wave functions of quantum systems. The key ingredient is to construct an envelope of a set of perturbative wave functions. This leads to a condition similar to that obtained from the principle of minimal sensitivity. Applications of the method to the quantum anharmonic oscillator and the double well potential show that uniformly valid wave functions with correct asymptotic behavior are obtained in the first-order optimized perturbation even for strong couplings. copyright 1997 The American Physical Society

Conformal invariance and pion wave functions of nonleading twist

International Nuclear Information System (INIS)

Braun, V.M.; Filyanov, I.E.

1989-01-01

The restrictions are studied for the general structure of pion wave functions of twist 3 and twist 4 imposed by the conformal symmetry and the equations of motion. A systematic expansion of wave functions in the conformal spin is built and the first order corrections to asymptotic formulae are calculated by the QCD sum rule method. In particular, we have found a multiplicatively renormalizable contribution into the two-particle wave function of twist 4 which cannot be expanded in a finite set of Gegenbauer polynomials. 19 refs.; 5 figs

General Forms of Wave Functions for Dipositronium, Ps2

Science.gov (United States)

Schrader, D.M.

2007-01-01

The consequences of particle interchange symmetry for the structure of wave functions of the states of dipositronium was recently discussed by the author [I]. In the present work, the methodology is simply explained, and the wave functions are explicitly given.

Coherent molecular transistor: control through variation of the gate wave function.

Science.gov (United States)

Ernzerhof, Matthias

2014-03-21

In quantum interference transistors (QUITs), the current through the device is controlled by variation of the gate component of the wave function that interferes with the wave function component joining the source and the sink. Initially, mesoscopic QUITs have been studied and more recently, QUITs at the molecular scale have been proposed and implemented. Typically, in these devices the gate lead is subjected to externally adjustable physical parameters that permit interference control through modifications of the gate wave function. Here, we present an alternative model of a molecular QUIT in which the gate wave function is directly considered as a variable and the transistor operation is discussed in terms of this variable. This implies that we specify the gate current as well as the phase of the gate wave function component and calculate the resulting current through the source-sink channel. Thus, we extend on prior works that focus on the phase of the gate wave function component as a control parameter while having zero or certain discrete values of the current. We address a large class of systems, including finite graphene flakes, and obtain analytic solutions for how the gate wave function controls the transistor.

Coherent molecular transistor: Control through variation of the gate wave function

International Nuclear Information System (INIS)

Ernzerhof, Matthias

2014-01-01

In quantum interference transistors (QUITs), the current through the device is controlled by variation of the gate component of the wave function that interferes with the wave function component joining the source and the sink. Initially, mesoscopic QUITs have been studied and more recently, QUITs at the molecular scale have been proposed and implemented. Typically, in these devices the gate lead is subjected to externally adjustable physical parameters that permit interference control through modifications of the gate wave function. Here, we present an alternative model of a molecular QUIT in which the gate wave function is directly considered as a variable and the transistor operation is discussed in terms of this variable. This implies that we specify the gate current as well as the phase of the gate wave function component and calculate the resulting current through the source-sink channel. Thus, we extend on prior works that focus on the phase of the gate wave function component as a control parameter while having zero or certain discrete values of the current. We address a large class of systems, including finite graphene flakes, and obtain analytic solutions for how the gate wave function controls the transistor

Collapse of the wave function models, ontology, origin, and implications

CERN Document Server

2018-01-01

This is the first single volume about the collapse theories of quantum mechanics, which is becoming a very active field of research in both physics and philosophy. In standard quantum mechanics, it is postulated that when the wave function of a quantum system is measured, it no longer follows the Schrödinger equation, but instantaneously and randomly collapses to one of the wave functions that correspond to definite measurement results. However, why and how a definite measurement result appears is unknown. A promising solution to this problem are collapse theories in which the collapse of the wave function is spontaneous and dynamical. Chapters written by distinguished physicists and philosophers of physics discuss the origin and implications of wave-function collapse, the controversies around collapse models and their ontologies, and new arguments for the reality of wave function collapse. This is an invaluable resource for students and researchers interested in the philosophy of physics and foundations of ...

Punctual Pade approximants as a regularization procedure for divergent and oscillatory partial wave expansions of the scattering amplitude

International Nuclear Information System (INIS)

Garibotti, C.R.; Grinstein, F.F.

1978-01-01

Previous theorems on the convergence of the [n,n+m] punctual Pade approximants to the scattering amplitude are extended. The new proofs include the cases of nonforward and backward scattering corresponding to potentials having 1/r and 1/r 2 long-range behaviors, for which the partial wave expansions are divergent and oscillatory, respectively. In this way, the ability of the approximation scheme as a summation method is established for all of the long-range potentials of interest in potential scattering

Punctual Pade Approximants as a regularization procedure for divergent and oscillatory partial wave expansions of the scattering amplitude

International Nuclear Information System (INIS)

Garibotti, C.R.; Grinstein, F.F.

1978-01-01

Previous theorems on the convergence of the [n, n+m] Punctual Pade Approximants to the scattering amplitude are extended. The new proofs include the cases of non-forward and backward scattering corresponding to potentials having 1/r and 1/r 2 long range behaviours, for which the partial wave expansions are divergent and oscillatory, respectively. In this way, the ability of the approximation scheme as a summation method is established for all of the long range potentials of interest in potential scattering [pt

Monotone Approximations of Minimum and Maximum Functions and Multi-objective Problems

International Nuclear Information System (INIS)

Stipanović, Dušan M.; Tomlin, Claire J.; Leitmann, George

2012-01-01

In this paper the problem of accomplishing multiple objectives by a number of agents represented as dynamic systems is considered. Each agent is assumed to have a goal which is to accomplish one or more objectives where each objective is mathematically formulated using an appropriate objective function. Sufficient conditions for accomplishing objectives are derived using particular convergent approximations of minimum and maximum functions depending on the formulation of the goals and objectives. These approximations are differentiable functions and they monotonically converge to the corresponding minimum or maximum function. Finally, an illustrative pursuit-evasion game example with two evaders and two pursuers is provided.

Monotone Approximations of Minimum and Maximum Functions and Multi-objective Problems

Energy Technology Data Exchange (ETDEWEB)

Stipanovic, Dusan M., E-mail: dusan@illinois.edu [University of Illinois at Urbana-Champaign, Coordinated Science Laboratory, Department of Industrial and Enterprise Systems Engineering (United States); Tomlin, Claire J., E-mail: tomlin@eecs.berkeley.edu [University of California at Berkeley, Department of Electrical Engineering and Computer Science (United States); Leitmann, George, E-mail: gleit@berkeley.edu [University of California at Berkeley, College of Engineering (United States)

2012-12-15

In this paper the problem of accomplishing multiple objectives by a number of agents represented as dynamic systems is considered. Each agent is assumed to have a goal which is to accomplish one or more objectives where each objective is mathematically formulated using an appropriate objective function. Sufficient conditions for accomplishing objectives are derived using particular convergent approximations of minimum and maximum functions depending on the formulation of the goals and objectives. These approximations are differentiable functions and they monotonically converge to the corresponding minimum or maximum function. Finally, an illustrative pursuit-evasion game example with two evaders and two pursuers is provided.

Proposal for determining the energy content of gravitational waves by using approximate symmetries of differential equations

International Nuclear Information System (INIS)

Hussain, Ibrar; Qadir, Asghar; Mahomed, F. M.

2009-01-01

Since gravitational wave spacetimes are time-varying vacuum solutions of Einstein's field equations, there is no unambiguous means to define their energy content. However, Weber and Wheeler had demonstrated that they do impart energy to test particles. There have been various proposals to define the energy content, but they have not met with great success. Here we propose a definition using 'slightly broken' Noether symmetries. We check whether this definition is physically acceptable. The procedure adopted is to appeal to 'approximate symmetries' as defined in Lie analysis and use them in the limit of the exact symmetry holding. A problem is noted with the use of the proposal for plane-fronted gravitational waves. To attain a better understanding of the implications of this proposal we also use an artificially constructed time-varying nonvacuum metric and evaluate its Weyl and stress-energy tensors so as to obtain the gravitational and matter components separately and compare them with the energy content obtained by our proposal. The procedure is also used for cylindrical gravitational wave solutions. The usefulness of the definition is demonstrated by the fact that it leads to a result on whether gravitational waves suffer self-damping.

The potential-free approach to the construction of the NN-wave functions

International Nuclear Information System (INIS)

Troitsky, V.E.

1984-01-01

The traditional approaches to the nonrelativistic NN-interaction use local and nonlocal potentials of the kind defined by different dynamical speculations. The wave functions are obtained then from the Schroedinger equation with the chosen potential. Here the author obtains the wave functions (scattering wave function and bound state wave function) directly from the scattering phases in the frame of a dispersion approach without use of potential. (Auth.)

Analytical approximations to seawater optical phase functions of scattering

Science.gov (United States)

Haltrin, Vladimir I.

2004-11-01

This paper proposes a number of analytical approximations to the classic and recently measured seawater light scattering phase functions. The three types of analytical phase functions are derived: individual representations for 15 Petzold, 41 Mankovsky, and 91 Gulf of Mexico phase functions; collective fits to Petzold phase functions; and analytical representations that take into account dependencies between inherent optical properties of seawater. The proposed phase functions may be used for problems of radiative transfer, remote sensing, visibility and image propagation in natural waters of various turbidity.

The Green-function transform and wave propagation

Directory of Open Access Journals (Sweden)

Colin eSheppard

2014-11-01

Full Text Available Fourier methods well known in signal processing are applied to three-dimensional wave propagation problems. The Fourier transform of the Green function, when written explicitly in terms of a real-valued spatial frequency, consists of homogeneous and inhomogeneous components. Both parts are necessary to result in a pure out-going wave that satisfies causality. The homogeneous component consists only of propagating waves, but the inhomogeneous component contains both evanescent and propagating terms. Thus we make a distinction between inhomogeneous waves and evanescent waves. The evanescent component is completely contained in the region of the inhomogeneous component outside the k-space sphere. Further, propagating waves in the Weyl expansion contain both homogeneous and inhomogeneous components. The connection between the Whittaker and Weyl expansions is discussed. A list of relevant spherically symmetric Fourier transforms is given.

Approximation by rational functions as processing method, analysis and transformation of neutron data

International Nuclear Information System (INIS)

Gaj, E.V.; Badikov, S.A.; Gusejnov, M.A.; Rabotnov, N.S.

1988-01-01

Possible applications of rational functions in the analysis of neutron cross sections, angular distributions and neutron constants generation are described. Results of investigations made in this direction, which have been obtained after the preceding conference in Kiev, are presented: the method of simultaneous treatment of several cross sections for one compound nucleus in the resonance range; the use of the Pade approximation for elastically scattered neutron angular distribution approximation; obtaining of subgroup constants on the basis of rational approximation of cross section functional dependence on dilution cross section; the first experience in function approximation by two variables

Bridge density functional approximation for non-uniform hard core repulsive Yukawa fluid

International Nuclear Information System (INIS)

Zhou Shiqi

2008-01-01

In this work, a bridge density functional approximation (BDFA) (J. Chem. Phys. 112, 8079 (2000)) for a non-uniform hard-sphere fluid is extended to a non-uniform hard-core repulsive Yukawa (HCRY) fluid. It is found that the choice of a bulk bridge functional approximation is crucial for both a uniform HCRY fluid and a non-uniform HCRY fluid. A new bridge functional approximation is proposed, which can accurately predict the radial distribution function of the bulk HCRY fluid. With the new bridge functional approximation and its associated bulk second order direct correlation function as input, the BDFA can be used to well calculate the density profile of the HCRY fluid subjected to the influence of varying external fields, and the theoretical predictions are in good agreement with the corresponding simulation data. The calculated results indicate that the present BDFA captures quantitatively the phenomena such as the coexistence of solid-like high density phase and low density gas phase, and the adsorption properties of the HCRY fluid, which qualitatively differ from those of the fluids combining both hard-core repulsion and an attractive tail. (condensed matter: structure, thermal and mechanical properties)

Simulating propagation of decomposed elastic waves using low-rank approximate mixed-domain integral operators for heterogeneous transversely isotropic media

KAUST Repository

Cheng, Jiubing; Wu, Zedong; Alkhalifah, Tariq Ali

2014-01-01

decomposition in anisotropic media is costly as the operators involved is dependent on the velocity, and thus not stationary. In this abstract, we propose an efficient approach to directly extrapolate the decomposed elastic waves using lowrank approximate mixed

On Love's approximation for fluid-filled elastic tubes

International Nuclear Information System (INIS)

Caroli, E.; Mainardi, F.

1980-01-01

A simple procedure is set up to introduce Love's approximation for wave propagation in thin-walled fluid-filled elastic tubes. The dispersion relation for linear waves and the radial profile for fluid pressure are determined in this approximation. It is shown that the Love approximation is valid in the low-frequency regime. (author)

Simulation of wind wave growth with reference source functions

Science.gov (United States)

Badulin, Sergei I.; Zakharov, Vladimir E.; Pushkarev, Andrei N.

2013-04-01

We present results of extensive simulations of wind wave growth with the so-called reference source function in the right-hand side of the Hasselmann equation written as follows First, we use Webb's algorithm [8] for calculating the exact nonlinear transfer function Snl. Second, we consider a family of wind input functions in accordance with recent consideration [9] ( )s S = ?(k)N , ?(k) = ? ? ?- f (?). in k 0 ?0 in (2) Function fin(?) describes dependence on angle ?. Parameters in (2) are tunable and determine magnitude (parameters ?0, ?0) and wave growth rate s [9]. Exponent s plays a key role in this study being responsible for reference scenarios of wave growth: s = 4-3 gives linear growth of wave momentum, s = 2 - linear growth of wave energy and s = 8-3 - constant rate of wave action growth. Note, the values are close to ones of conventional parameterizations of wave growth rates (e.g. s = 1 for [7] and s = 2 for [5]). Dissipation function Sdiss is chosen as one providing the Phillips spectrum E(?) ~ ?5 at high frequency range [3] (parameter ?diss fixes a dissipation scale of wind waves) Sdiss = Cdissμ4w?N (k)θ(? - ?diss) (3) Here frequency-dependent wave steepness μ2w = E(?,?)?5-g2 makes this function to be heavily nonlinear and provides a remarkable property of stationary solutions at high frequencies: the dissipation coefficient Cdiss should keep certain value to provide the observed power-law tails close to the Phillips spectrum E(?) ~ ?-5. Our recent estimates [3] give Cdiss ? 2.0. The Hasselmann equation (1) with the new functions Sin, Sdiss (2,3) has a family of self-similar solutions of the same form as previously studied models [1,3,9] and proposes a solid basis for further theoretical and numerical study of wave evolution under action of all the physical mechanisms: wind input, wave dissipation and nonlinear transfer. Simulations of duration- and fetch-limited wind wave growth have been carried out within the above model setup to check its

Balancing Exchange Mixing in Density-Functional Approximations for Iron Porphyrin.

Science.gov (United States)

Berryman, Victoria E J; Boyd, Russell J; Johnson, Erin R

2015-07-14

Predicting the correct ground-state multiplicity for iron(II) porphyrin, a high-spin quintet, remains a significant challenge for electronic-structure methods, including commonly employed density functionals. An even greater challenge for these methods is correctly predicting favorable binding of O2 to iron(II) porphyrin, due to the open-shell singlet character of the adduct. In this work, the performance of a modest set of contemporary density-functional approximations is assessed and the results interpreted using Bader delocalization indices. It is found that inclusion of greater proportions of Hartree-Fock exchange, in hybrid or range-separated hybrid functionals, has opposing effects; it improves the ability of the functional to identify the ground state but is detrimental to predicting favorable dioxygen binding. Because of the uncomplementary nature of these properties, accurate prediction of both the relative spin-state energies and the O2 binding enthalpy eludes conventional density-functional approximations.

On approximation of functions by product operators

Directory of Open Access Journals (Sweden)

Hare Krishna Nigam

2013-12-01

Full Text Available In the present paper, two quite new reults on the degree of approximation of a function f belonging to the class Lip(α,r, 1≤ r <∞ and the weighted class W(Lr,ξ(t, 1≤ r <∞ by (C,2(E,1 product operators have been obtained. The results obtained in the present paper generalize various known results on single operators.

Simple functional-differential equations for the bound-state wave-function components

International Nuclear Information System (INIS)

Kamuntavicius, G.P.

1986-01-01

The author presents a new method of a direct derivation of differential equations for the wave-function components of identical-particles systems. The method generates in a simple manner all the possible variants of these equations. In some cases they are the differential equations of Faddeev or Yakubovskii. It is shown that the case of the bound states allows to formulate very simple equations for the components which are equivalent to the Schroedinger equation for the complete wave function. The components with a minimal antisymmetry are defined and the corresponding equations are derived. (Auth.)

Hartree-Fock implementation using a Laguerre-based wave function for the ground state and correlation energies of two-electron atoms.

Science.gov (United States)

King, Andrew W; Baskerville, Adam L; Cox, Hazel

2018-03-13

An implementation of the Hartree-Fock (HF) method using a Laguerre-based wave function is described and used to accurately study the ground state of two-electron atoms in the fixed nucleus approximation, and by comparison with fully correlated (FC) energies, used to determine accurate electron correlation energies. A variational parameter A is included in the wave function and is shown to rapidly increase the convergence of the energy. The one-electron integrals are solved by series solution and an analytical form is found for the two-electron integrals. This methodology is used to produce accurate wave functions, energies and expectation values for the helium isoelectronic sequence, including at low nuclear charge just prior to electron detachment. Additionally, the critical nuclear charge for binding two electrons within the HF approach is calculated and determined to be Z HF C =1.031 177 528.This article is part of the theme issue 'Modern theoretical chemistry'. © 2018 The Author(s).

Adaptive Linear and Normalized Combination of Radial Basis Function Networks for Function Approximation and Regression

Directory of Open Access Journals (Sweden)

Yunfeng Wu

2014-01-01

Full Text Available This paper presents a novel adaptive linear and normalized combination (ALNC method that can be used to combine the component radial basis function networks (RBFNs to implement better function approximation and regression tasks. The optimization of the fusion weights is obtained by solving a constrained quadratic programming problem. According to the instantaneous errors generated by the component RBFNs, the ALNC is able to perform the selective ensemble of multiple leaners by adaptively adjusting the fusion weights from one instance to another. The results of the experiments on eight synthetic function approximation and six benchmark regression data sets show that the ALNC method can effectively help the ensemble system achieve a higher accuracy (measured in terms of mean-squared error and the better fidelity (characterized by normalized correlation coefficient of approximation, in relation to the popular simple average, weighted average, and the Bagging methods.

Distorted-wave Born approximation in the case of an optical scattering potential

International Nuclear Information System (INIS)

Mytnichenko, Sergey V.

2005-01-01

Application of the distorted-wave Born approximation in the conventional form developed for the case of a real scattering potential is shown to cause significant errors in calculating X-ray diffuse scattering from non-ideal crystals, superlattices, multilayers and other objects if energy dissipation (photoabsorption, inelastic scattering, and so on) is not negligible, or in other words, in the case of an optical (complex) scattering potential. We show how a correct expression for the X-ray diffuse-scattering cross-section can be obtained in this case. Generally, the diffuse-scattering cross-section from an optical potential is not T-invariant, i.e. the reciprocity principle is violated. Violations of T-invariance are more evident when the dynamical nature of the diffraction is more critical

The Hartree-Fock seniority approximation

International Nuclear Information System (INIS)

Gomez, J.M.G.; Prieto, C.

1986-01-01

A new self-consistent method is used to take into account the mean-field and the pairing correlations in nuclei at the same time. We call it the Hartree-Fock seniority approximation, because the long-range and short-range correlations are treated in the frameworks of Hartree-Fock theory and the seniority scheme. The method is developed in detail for a minimum-seniority variational wave function in the coordinate representation for an effective interaction of the Skyrme type. An advantage of the present approach over the Hartree-Fock-Bogoliubov theory is the exact conservation of angular momentum and particle number. Furthermore, the computational effort required in the Hartree-Fock seniority approximation is similar to that ofthe pure Hartree-Fock picture. Some numerical calculations for Ca isotopes are presented. (orig.)

Electron capture by alpha particles from helium atoms in a Coulomb-Born distorted-wave approximation

International Nuclear Information System (INIS)

Ghanbari-Adivi, E; Ghavaminia, H

2012-01-01

A three-body Coulomb-Born continuum distorted-wave approximation is applied to calculate the differential and total cross sections for single-electron exchange in the collision of fast alpha particles with helium atoms in their ground states. The applied first-order distorted wave theory satisfies correct Coulomb boundary conditions. Both post and prior forms of the transition amplitude are calculated. The nuclear-screening effect of the passive electron on the differential and total cross sections is investigated. The results are compared with those of other theories and with the available experimental data. For differential cross sections, the comparisons show a reasonable agreement with empirical measurements at higher impact energies. The agreement between experimental data and the present calculations for total cross sections with the average of the post and prior forms of the transition amplitude is reasonable at all the specified energies.

Rational function approximation method for discrete ordinates problems in slab geometry

International Nuclear Information System (INIS)

Leal, Andre Luiz do C.; Barros, Ricardo C.

2009-01-01

In this work we use rational function approaches to obtain the transfer functions that appear in the spectral Green's function (SGF) auxiliary equations for one-speed isotropic scattering SN equations in one-dimensional Cartesian geometry. For this task we use the computation of the Pade approximants to compare the results with the standard SGF method's applied to deep penetration problems in homogeneous domains. This work is a preliminary investigation of a new proposal for handling leakage terms that appear in the two transverse integrated one-dimensional SN equations in the exponential SGF method (SGF-ExpN). Numerical results are presented to illustrate the rational function approximation accuracy. (author)

Numerical approximations of difference functional equations and applications

Directory of Open Access Journals (Sweden)

Zdzisław Kamont

2005-01-01

Full Text Available We give a theorem on the error estimate of approximate solutions for difference functional equations of the Volterra type. We apply this general result in the investigation of the stability of difference schemes generated by nonlinear first order partial differential functional equations and by parabolic problems. We show that all known results on difference methods for initial or initial boundary value problems can be obtained as particular cases of this general and simple result. We assume that the right hand sides of equations satisfy nonlinear estimates of the Perron type with respect to functional variables.

A convenient analytical form for the triton wave function

International Nuclear Information System (INIS)

Hajduk, C.; Green, A.M.; Sainio, M.E.

1979-01-01

The triton wave function obtained by solving the Faddeev equations with the Reid soft core potential is parametrized in a symmetrized cluster form. As a test the 3 He charge form factor is calculated for the exact and the parametrized wave functions and reasonable agreement between the two is found. (author)

Application of modified analytical function for approximation and computer simulation of diffraction profile

International Nuclear Information System (INIS)

Marrero, S. I.; Turibus, S. N.; Assis, J. T. De; Monin, V. I.

2011-01-01

Data processing of the most of diffraction experiments is based on determination of diffraction line position and measurement of broadening of diffraction profile. High precision and digitalisation of these procedures can be resolved by approximation of experimental diffraction profiles by analytical functions. There are various functions for these purposes both simples, like Gauss function, but no suitable for wild range of experimental profiles and good approximating functions but complicated for practice using, like Vougt or PersonVII functions. Proposed analytical function is modified Cauchy function which uses two variable parameters allowing describing any experimental diffraction profile. In the presented paper modified function was applied for approximation of diffraction lines of steels after various physical and mechanical treatments and simulation of diffraction profiles applied for study of stress gradients and distortions of crystal structure. (Author)

BAYESZ, S-Wave, P-Wave Resonance Level Spacing and Strength Functions

International Nuclear Information System (INIS)

Moore, M.S.

1982-01-01

A - Description of problem or function: BAYESZ calculates average s- and p-wave level spacings, strength functions, and average radiation widths of a mixed sequence of s- and p-wave resonances whose parameters are supplied as input. The code is based on two physical assumptions: 1) The neutron reduced width distribution for each open channel is a chi-squared distribution with one degree of freedom, i.e. Porter-Thomas. 2) The spacing distribution follows the Gaussian Orthogonal Ensemble. This property is used, however, only to fix the s- to p-wave level density ratio as proportional to (2J+1) with a spin cut-off correction. B - Method of solution: The method used is an extension of that described by Moore et al. in reference (1), and is based on the method of moments of a truncated Porter-Thomas distribution. C - Restrictions on the complexity of the problem: Parameters for a maximum of 500 individual resonances can be specified. This restriction can be relaxed by increasing array dimensions

ESTIMA, Neutron Width Level Spacing, Neutron Strength Function of S- Wave, P-Wave Resonances

International Nuclear Information System (INIS)

Fort, E.

1982-01-01

1 - Description of problem or function: ESTIMA calculates level spacing and neutron strength function of a mixed sequence of s- and p-wave resonances given a set of neutron widths as input parameters. Three algorithms are used, two of which calculate s-wave average parameters and assume that the reduced widths obey a Porter-Thomas distribution truncated by a minimum detection threshold. The third performs a maximum likelihood fit to a truncated chi-squared distribution of any specified number of degrees of freedom, i.e. it can be used for calculating s-wave or p-wave average parameters. Resonances of undeclared angular orbital momentum are divided into groups of probable s-wave and probable p-wave by a simple application of Bayes' Theorem. 2 - Method of solution: Three algorithms are used: i) GAMN method, based on simple moments properties of a Porter-Thomas distribution. ii) Missing Level Estimator, a simplified version of the algorithm used by the program BAYESZ. iii) ESTIMA, a maximum likelihood fit. 3 - Restrictions on the complexity of the problem: A maximum of 400 resonances is allowed in the version available from NEADB, however this restriction can be relaxed by increasing array dimensions

Four-body correlation embedded in antisymmetrized geminal power wave function.

Science.gov (United States)

Kawasaki, Airi; Sugino, Osamu

2016-12-28

We extend the Coleman's antisymmetrized geminal power (AGP) to develop a wave function theory that can incorporate up to four-body correlation in a region of strong correlation. To facilitate the variational determination of the wave function, the total energy is rewritten in terms of the traces of geminals. This novel trace formula is applied to a simple model system consisting of one dimensional Hubbard ring with a site of strong correlation. Our scheme significantly improves the result obtained by the AGP-configuration interaction scheme of Uemura et al. and also achieves more efficient compression of the degrees of freedom of the wave function. We regard the result as a step toward a first-principles wave function theory for a strongly correlated point defect or adsorbate embedded in an AGP-based mean-field medium.

On- and off-resonance radiation-atom-coupling matrix elements involving extended atomic wave functions

Science.gov (United States)

Komninos, Yannis; Mercouris, Theodoros; Nicolaides, Cleanthes A.

2014-01-01

In continuation of our earlier works, we present results concerning the computation of matrix elements of the multipolar Hamiltonian (MPH) between extended wave functions that are obtained numerically. The choice of the MPH is discussed in connection with the broader issue of the form of radiation-atom (or -molecule) interaction that is appropriate for the systematic solution of various problems of matter-radiation interaction. We derive analytic formulas, in terms of the sine-integral function and spherical Bessel functions of various orders, for the cumulative radial integrals that were obtained and calculated by Komninos, Mercouris, and Nicolaides [Phys. Rev. A 71, 023410 (2005), 10.1103/PhysRevA.71.023410]. This development allows the much faster and more accurate computation of such matrix elements, a fact that enhances the efficiency with which the time-dependent Schrödinger equation is solved nonperturbatively, in the framework of the state-specific expansion approach. The formulas are applicable to the general case where a pair of orbitals with angular parts |ℓ1,m1> and |ℓ2,m2> are coupled radiatively. As a test case, we calculate the matrix elements of the electric field and of the paramagnetic operators for on- and off-resonance transitions, between hydrogenic circular states of high angular momentum, whose quantum numbers are chosen so as to satisfy electric dipole and electric quadrupole selection rules. Because of the nature of their wave function (they are nodeless and the large centrifugal barrier keeps their overwhelming part at large distances from the nucleus), the validity of the electric dipole approximation in various applications where the off-resonance couplings must be considered becomes precarious. For example, for the transition from the circular state with n = 20 to that with n = 21, for which ≈400 a.u., the dipole approximation starts to fail already at XUV wavelengths (λ <125nm).

A Proton-Cyclotron Wave Storm Generated by Unstable Proton Distribution Functions in the Solar Wind

Science.gov (United States)

Wicks, R. T.; Alexander, R. L.; Stevens, M.; Wilson, L. B., III; Moya, P. S.; Vinas, A.; Jian, L. K.; Roberts, D. A.; O’Modhrain, S.; Gilbert, J. A.;

Accuracy of three-body wave functions obtained with the correlation-function hyperspherical-harmonic method

International Nuclear Information System (INIS)

Haftel, M.I.; Mandelzweig, V.B.

1990-01-01

The local convergence and accuracy of wave functions obtained by direct solution of the Schroedinger equation with the help of the correlation-function hyperspherical-harmonic method are analyzed for ground and excited states of the helium atom and for the ground state of the positronium negative ion. The inclusion of the cusp conditions into the correlation function is shown to be of crucial importance, not only near the coalescence points, but also away from them. The proper inclusion of all cusps yields for the ground state of the helium atom the local wave-function accuracy of about 10 -7 for different interparticle distances. The omission of one of the cusps in the excited helium atom reduces the wave-function precision to 10 -2 near the corresponding coalescence point and to 10 -4 --10 -5 away from it

QCD Phenomenology and Light-Front Wave Functions

International Nuclear Information System (INIS)

Brodsky, St.J.

2001-01-01

A natural calculus for describing the bound-state structure of relativistic composite systems in quantum field theory is the light-front Fock expansion which encodes the properties of a hadrons in terms of a set of frame-independent n-particle wave functions. Light-front quantization in the doubly-transverse light-cone gauge has a number of remarkable advantages, including explicit unitarity, a physical Fock expansion, the absence of ghost degrees of freedom, and the decoupling properties needed to prove factorization theorems in high momentum transfer inclusive and exclusive reactions. A number of applications are discussed in these lectures, including semileptonic B decays, two-photon exclusive reactions, diffractive dissociation into jets, and deeply virtual Compton scattering. The relation of the intrinsic sea to the light-front wave functions is discussed. Light-front quantization can also be used in the Hamiltonian form to construct an event generator for high energy physics reactions at the amplitude level. The light-cone partition function, summed over exponentially-weighted light-cone energies, has simple boost properties which may be useful for studies in heavy ion collisions. I also review recent work which shows that the structure functions measured in deep inelastic lepton scattering are affected by final-state rescattering, thus modifying their connection to light-front probability distributions. In particular, the shadowing of nuclear structure functions is due to destructive interference effects from leading-twist diffraction of the virtual photon, physics not included in the nuclear light-cone wave functions. (author)

On approximation and energy estimates for delta 6-convex functions.

Science.gov (United States)

Saleem, Muhammad Shoaib; Pečarić, Josip; Rehman, Nasir; Khan, Muhammad Wahab; Zahoor, Muhammad Sajid

2018-01-01

The smooth approximation and weighted energy estimates for delta 6-convex functions are derived in this research. Moreover, we conclude that if 6-convex functions are closed in uniform norm, then their third derivatives are closed in weighted [Formula: see text]-norm.

Cerebral functional connectivity and Mayer waves in mice: Phenomena and separability.

Science.gov (United States)

Bumstead, Jonathan R; Bauer, Adam Q; Wright, Patrick W; Culver, Joseph P

2017-02-01

Resting-state functional connectivity is a growing neuroimaging approach that analyses the spatiotemporal structure of spontaneous brain activity, often using low-frequency (Mayer waves. Despite how close in frequency these phenomena exist, there is little research on how vasomotion and Mayer waves are related to or affect resting-state functional connectivity. In this study, we analyze spontaneous hemodynamic fluctuations over the mouse cortex using optical intrinsic signal imaging. We found spontaneous occurrence of oscillatory hemodynamics ∼0.2 Hz consistent with the properties of Mayer waves reported in the literature. Across a group of mice (n = 19), there was a large variability in the magnitude of Mayer waves. However, regardless of the magnitude of Mayer waves, functional connectivity patterns could be recovered from hemodynamic signals when filtered to the lower frequency band, 0.01-0.08 Hz. Our results demonstrate that both Mayer waves and resting-state functional connectivity patterns can co-exist simultaneously, and that they can be separated by applying bandpass filters.

Degenerate RS perturbation theory. [Rayleigh-Schroedinger energies and wave functions

Science.gov (United States)

Hirschfelder, J. O.; Certain, P. R.

1974-01-01

A concise, systematic procedure is given for determining the Rayleigh-Schroedinger energies and wave functions of degenerate states to arbitrarily high orders even when the degeneracies of the various states are resolved in arbitrary orders. The procedure is expressed in terms of an iterative cycle in which the energy through the (2n + 1)-th order is expressed in terms of the partially determined wave function through the n-th order. Both a direct and an operator derivation are given. The two approaches are equivalent and can be transcribed into each other. The direct approach deals with the wave functions (without the use of formal operators) and has the advantage that it resembles the usual treatment of nondegenerate perturbations and maintains close contact with the basic physics. In the operator approach, the wave functions are expressed in terms of infinite-order operators which are determined by the successive resolution of the space of the zeroth-order functions.

Deep-inelastic structure functions in an approximation to the bag theory

International Nuclear Information System (INIS)

Jaffe, R.L.

1975-01-01

A cavity approximation to the bag theory developed earlier is extended to the treatment of forward virtual Compton scattering. In the Bjorken limit and for small values of ω (ω = vertical-bar2p center-dot q/q 2 vertical-bar) it is argued that the operator nature of the bag boundaries might be ignored. Structure functions are calculated in one and three dimensions. Bjorken scaling is obtained. The model provides a realization of light-cone current algebra and possesses a parton interpretation. The structure functions show a quasielastic peak. The spreading of the structure functions about the peak is associated with confinement. As expected, Regge behavior is not obtained for large ω. The ''momentum sum rule'' is saturated, indicating that the hadron's charged constituents carry all the momentum in this model. νW/subL/ is found to scale and is calculable. Application of the model to the calculation of spin-dependent and chiral-symmetry--violating structure functions is proposed. The nature of the intermediate states in this approximation is discussed. Problems associated with the cavity approximation are also discussed

On the approximation of the limit cycles function

Directory of Open Access Journals (Sweden)

L. Cherkas

2007-11-01

Full Text Available We consider planar vector fields depending on a real parameter. It is assumed that this vector field has a family of limit cycles which can be described by means of the limit cycles function $l$. We prove a relationship between the multiplicity of a limit cycle of this family and the order of a zero of the limit cycles function. Moreover, we present a procedure to approximate $l(x$, which is based on the Newton scheme applied to the Poincaré function and represents a continuation method. Finally, we demonstrate the effectiveness of the proposed procedure by means of a Liénard system.

General Rytov approximation.

Science.gov (United States)

Potvin, Guy

2015-10-01

We examine how the Rytov approximation describing log-amplitude and phase fluctuations of a wave propagating through weak uniform turbulence can be generalized to the case of turbulence with a large-scale nonuniform component. We show how the large-scale refractive index field creates Fermat rays using the path integral formulation for paraxial propagation. We then show how the second-order derivatives of the Fermat ray action affect the Rytov approximation, and we discuss how a numerical algorithm would model the general Rytov approximation.

Reduction of the Breit Coulomb equation to an equivalent Schroedinger equation, and investigation of the behavior of the wave function near the origin

International Nuclear Information System (INIS)

Malenfant, J.

1988-01-01

The Breit equation for two equal-mass spin-1/2 particles interacting through an attractive Coulomb potential is separated into its angular and radial parts, obtaining coupled sets of first-order differential equations for the radial wave functions. The radial equations for the 1 J/sub J/, 3 J/sub J/, and 3 P 0 states are further reduced to a single, one-dimensional Schroedinger equation with a relatively simple effective potential. No approximations, other than the initial one of an instantaneous Coulomb interaction, are made in deriving this equation; it accounts for all relativistic effects, as well as for mixing between different components of the wave function. Approximate solutions are derived for this Schroedinger equation, which gives the correct O(α 4 ) term for the 1 1 S 0 energy and for the n 1 J/sub J/ energies, for J>0. The radial equations for the 3 (J +- 1)/sub J/ states are reduced to two second-order coupled equations. At small r, the Breit Coulomb wave functions behave as r/sup ν//sup -1/, where ν is either √J(J+1)+1-α 2 /4 or √J(J+1)-α 2 /4 . The 1 S 0 and 3 P 0 wave functions therefore diverge at the origin as r/sup //sup √//sup 1-//sup α//sup <2//4 -1$. This divergence of the J = 0 states, however, does not occur when the spin-spin interaction, -(α/r)αxα, is added to the Coulomb potential

Strong semiclassical approximation of Wigner functions for the Hartree dynamics

KAUST Repository

Athanassoulis, Agissilaos; Paul, Thierry; Pezzotti, Federica; Pulvirenti, Mario

2011-01-01

We consider the Wigner equation corresponding to a nonlinear Schrödinger evolution of the Hartree type in the semiclassical limit h → 0. Under appropriate assumptions on the initial data and the interaction potential, we show that the Wigner function is close in L 2 to its weak limit, the solution of the corresponding Vlasov equation. The strong approximation allows the construction of semiclassical operator-valued observables, approximating their quantum counterparts in Hilbert-Schmidt topology. The proof makes use of a pointwise-positivity manipulation, which seems necessary in working with the L 2 norm and the precise form of the nonlinearity. We employ the Husimi function as a pivot between the classical probability density and the Wigner function, which - as it is well known - is not pointwise positive in general.

Determination of the S-wave scattering shape parameter P from the zero-energy wave function

International Nuclear Information System (INIS)

Kermode, M.W.; van Dijk, W.

1990-01-01

We show that for S-wave scattering at an energy k 2 by a local potential which supports no more than one bound state, the shape parameter P and coefficients of higher powers of k 2 in the effective range expansion function cotδ=-1/a+1/2 r 0 k 2 -Pr 0 3 k 3 +Qr 0 5 k 6 +..., where δ is the phase shift, may be obtained from the zero-energy wave function, u 0 (r). Thus δ itself may be determined from u 0 . We show that Pr 0 3 =∫ 0 R [β(r)u 0 2 (r)-bar β(r)bar u 0 2 (r)]dr, where r 0 is the effective range, β(r) is determined from an integral involving the wave function, and bar β(r) is a simple function of r which involves the scattering length and effective range

Mean-field approximation for spacing distribution functions in classical systems

Science.gov (United States)

González, Diego Luis; Pimpinelli, Alberto; Einstein, T. L.

2012-01-01

We propose a mean-field method to calculate approximately the spacing distribution functions p(n)(s) in one-dimensional classical many-particle systems. We compare our method with two other commonly used methods, the independent interval approximation and the extended Wigner surmise. In our mean-field approach, p(n)(s) is calculated from a set of Langevin equations, which are decoupled by using a mean-field approximation. We find that in spite of its simplicity, the mean-field approximation provides good results in several systems. We offer many examples illustrating that the three previously mentioned methods give a reasonable description of the statistical behavior of the system. The physical interpretation of each method is also discussed.

Mathieu functions and its useful approximation for elliptical waveguides

Science.gov (United States)

Pillay, Shamini; Kumar, Deepak

2017-11-01

The standard form of the Mathieu differential equation is where a and q are real parameters and q > 0. In this paper we obtain closed formula for the generic term of expansions of modified Mathieu functions in terms of Bessel and modified Bessel functions in the following cases: Let ξ0 = ξ0, where i can take the values 1 and 2 corresponding to the first and the second boundary. These approximations also provide alternative methods for numerical evaluation of Mathieu functions.

Generalized finite polynomial approximation (WINIMAX) to the reduced partition function of isotopic molecules

International Nuclear Information System (INIS)

Lee, M.W.; Bigeleisen, J.

1978-01-01

The MINIMAX finite polynomial approximation to an arbitrary function has been generalized to include a weighting function (WINIMAX). It is suggested that an exponential is a reasonable weighting function for the logarithm of the reduced partition function of a harmonic oscillator. Comparison of the error function for finite orthogonal polynomial (FOP), MINIMAX, and WINIMAX expansions of the logarithm of the reduced vibrational partition function show WINIMAX to be the best of the three approximations. A condensed table of WINIMAX coefficients is presented. The FOP, MINIMAX, and WINIMAX approximations are compared with exact calculations of the logarithm of the reduced partition function ratios for isotopic substitution in H 2 O, CH 4 , CH 2 O, C 2 H 4 , and C 2 H 6 at 300 0 K. Both deuterium and heavy atom isotope substitution are studied. Except for a third order expansion involving deuterium substitution, the WINIMAX method is superior to FOP and MINIMAX. At the level of a second order expansion WINIMAX approximations to ln(s/s')f are good to 2.5% and 6.5% for deuterium and heavy atom substitution, respectively

Regression analysis and transfer function in estimating the parameters of central pulse waves from brachial pulse wave.

Science.gov (United States)

Chai Rui; Li Si-Man; Xu Li-Sheng; Yao Yang; Hao Li-Ling

2017-07-01

This study mainly analyzed the parameters such as ascending branch slope (A_slope), dicrotic notch height (Hn), diastolic area (Ad) and systolic area (As) diastolic blood pressure (DBP), systolic blood pressure (SBP), pulse pressure (PP), subendocardial viability ratio (SEVR), waveform parameter (k), stroke volume (SV), cardiac output (CO) and peripheral resistance (RS) of central pulse wave invasively and non-invasively measured. These parameters extracted from the central pulse wave invasively measured were compared with the parameters measured from the brachial pulse waves by a regression model and a transfer function model. The accuracy of the parameters which were estimated by the regression model and the transfer function model was compared too. Our findings showed that in addition to the k value, the above parameters of the central pulse wave and the brachial pulse wave invasively measured had positive correlation. Both the regression model parameters including A_slope, DBP, SEVR and the transfer function model parameters had good consistency with the parameters invasively measured, and they had the same effect of consistency. The regression equations of the three parameters were expressed by Y'=a+bx. The SBP, PP, SV, CO of central pulse wave could be calculated through the regression model, but their accuracies were worse than that of transfer function model.

Approximation of the exponential integral (well function) using sampling methods

Science.gov (United States)

Baalousha, Husam Musa

2015-04-01

Exponential integral (also known as well function) is often used in hydrogeology to solve Theis and Hantush equations. Many methods have been developed to approximate the exponential integral. Most of these methods are based on numerical approximations and are valid for a certain range of the argument value. This paper presents a new approach to approximate the exponential integral. The new approach is based on sampling methods. Three different sampling methods; Latin Hypercube Sampling (LHS), Orthogonal Array (OA), and Orthogonal Array-based Latin Hypercube (OA-LH) have been used to approximate the function. Different argument values, covering a wide range, have been used. The results of sampling methods were compared with results obtained by Mathematica software, which was used as a benchmark. All three sampling methods converge to the result obtained by Mathematica, at different rates. It was found that the orthogonal array (OA) method has the fastest convergence rate compared with LHS and OA-LH. The root mean square error RMSE of OA was in the order of 1E-08. This method can be used with any argument value, and can be used to solve other integrals in hydrogeology such as the leaky aquifer integral.

Animating Nested Taylor Polynomials to Approximate a Function

Science.gov (United States)

Mazzone, Eric F.; Piper, Bruce R.

2010-01-01

The way that Taylor polynomials approximate functions can be demonstrated by moving the center point while keeping the degree fixed. These animations are particularly nice when the Taylor polynomials do not intersect and form a nested family. We prove a result that shows when this nesting occurs. The animations can be shown in class or…

P wave dispersion and maximum P wave duration are independently associated with rapid renal function decline.

Science.gov (United States)

Su, Ho-Ming; Tsai, Wei-Chung; Lin, Tsung-Hsien; Hsu, Po-Chao; Lee, Wen-Hsien; Lin, Ming-Yen; Chen, Szu-Chia; Lee, Chee-Siong; Voon, Wen-Chol; Lai, Wen-Ter; Sheu, Sheng-Hsiung

2012-01-01

The P wave parameters measured by 12-lead electrocardiogram (ECG) are commonly used as noninvasive tools to assess for left atrial enlargement. There are limited studies to evaluate whether P wave parameters are independently associated with decline in renal function. Accordingly, the aim of this study is to assess whether P wave parameters are independently associated with progression to renal end point of ≥25% decline in estimated glomerular filtration rate (eGFR). This longitudinal study included 166 patients. The renal end point was defined as ≥25% decline in eGFR. We measured two ECG P wave parameters corrected by heart rate, i.e. corrected P wave dispersion (PWdisperC) and corrected P wave maximum duration (PWdurMaxC). Heart function and structure were measured from echocardiography. Clinical data, P wave parameters, and echocardiographic measurements were compared and analyzed. Forty-three patients (25.9%) reached renal end point. Kaplan-Meier curves for renal end point-free survival showed PWdisperC > median (63.0 ms) (log-rank P = 0.004) and PWdurMaxC > median (117.9 ms) (log-rank Pfunction decline.

On approximation and energy estimates for delta 6-convex functions

Directory of Open Access Journals (Sweden)

Muhammad Shoaib Saleem

2018-02-01

Full Text Available Abstract The smooth approximation and weighted energy estimates for delta 6-convex functions are derived in this research. Moreover, we conclude that if 6-convex functions are closed in uniform norm, then their third derivatives are closed in weighted L2 $L^{2}$-norm.

BCS wave function, matrix product states, and the Ising conformal field theory

Science.gov (United States)

Montes, Sebastián; Rodríguez-Laguna, Javier; Sierra, Germán

2017-11-01

We present a characterization of the many-body lattice wave functions obtained from the conformal blocks (CBs) of the Ising conformal field theory (CFT). The formalism is interpreted as a matrix product state using continuous ancillary degrees of freedom. We provide analytic and numerical evidence that the resulting states can be written as BCS states. We give a complete proof that the translationally invariant 1D configurations have a BCS form and we find suitable parent Hamiltonians. In particular, we prove that the ground state of the finite-size critical Ising transverse field (ITF) Hamiltonian can be obtained with this construction. Finally, we study 2D configurations using an operator product expansion (OPE) approximation. We associate these states to the weak pairing phase of the p +i p superconductor via the scaling of the pairing function and the entanglement spectrum.

Correlation energy functional within the GW -RPA: Exact forms, approximate forms, and challenges

Science.gov (United States)

Ismail-Beigi, Sohrab

2010-05-01

In principle, the Luttinger-Ward Green’s-function formalism allows one to compute simultaneously the total energy and the quasiparticle band structure of a many-body electronic system from first principles. We present approximate and exact expressions for the correlation energy within the GW -random-phase approximation that are more amenable to computation and allow for developing efficient approximations to the self-energy operator and correlation energy. The exact form is a sum over differences between plasmon and interband energies. The approximate forms are based on summing over screened interband transitions. We also demonstrate that blind extremization of such functionals leads to unphysical results: imposing physical constraints on the allowed solutions (Green’s functions) is necessary. Finally, we present some relevant numerical results for atomic systems.

Systematic theoretical investigation of the zero-field splitting in Gd(III) complexes: Wave function and density functional approaches

Energy Technology Data Exchange (ETDEWEB)

Khan, Shehryar, E-mail: sherkhan@fysik.su.se; Odelius, Michael, E-mail: odelius@fysik.su.se [Department of Physics, Stockholm University, AlbaNova University Center, S-106 91 Stockholm (Sweden); Kubica-Misztal, Aleksandra [Institute of Physics, Jagiellonian University, ul. Reymonta 4, PL-30-059 Krakow (Poland); Kruk, Danuta [Faculty of Mathematics and Computer Science, University of Warmia and Mazury in Olsztyn, Sloneczna 54, Olsztyn PL-10710 (Poland); Kowalewski, Jozef [Department of Materials and Environmental Chemistry, Arrhenius Laboratory, Stockholm University, S-106 91 Stockholm (Sweden)

2015-01-21

The zero-field splitting (ZFS) of the electronic ground state in paramagnetic ions is a sensitive probe of the variations in the electronic and molecular structure with an impact on fields ranging from fundamental physical chemistry to medical applications. A detailed analysis of the ZFS in a series of symmetric Gd(III) complexes is presented in order to establish the applicability and accuracy of computational methods using multiconfigurational complete-active-space self-consistent field wave functions and of density functional theory calculations. The various computational schemes are then applied to larger complexes Gd(III)DOTA(H{sub 2}O){sup −}, Gd(III)DTPA(H{sub 2}O){sup 2−}, and Gd(III)(H{sub 2}O){sub 8}{sup 3+} in order to analyze how the theoretical results compare to experimentally derived parameters. In contrast to approximations based on density functional theory, the multiconfigurational methods produce results for the ZFS of Gd(III) complexes on the correct order of magnitude.

Efficient approximation of the Struve functions Hn occurring in the calculation of sound radiation quantities.

Science.gov (United States)

Aarts, Ronald M; Janssen, Augustus J E M

2016-12-01

The Struve functions H n (z), n=0, 1, ...  are approximated in a simple, accurate form that is valid for all z≥0. The authors previously treated the case n = 1 that arises in impedance calculations for the rigid-piston circular radiator mounted in an infinite planar baffle [Aarts and Janssen, J. Acoust. Soc. Am. 113, 2635-2637 (2003)]. The more general Struve functions occur when other acoustical quantities and/or non-rigid pistons are considered. The key step in the paper just cited is to express H 1 (z) as (2/π)-J 0 (z)+(2/π) I(z), where J 0 is the Bessel function of order zero and the first kind and I(z) is the Fourier cosine transform of [(1-t)/(1+t)] 1/2 , 0≤t≤1. The square-root function is optimally approximated by a linear function ĉt+d̂, 0≤t≤1, and the resulting approximated Fourier integral is readily computed explicitly in terms of sin z/z and (1-cos z)/z 2 . The same approach has been used by Maurel, Pagneux, Barra, and Lund [Phys. Rev. B 75, 224112 (2007)] to approximate H 0 (z) for all z≥0. In the present paper, the square-root function is optimally approximated by a piecewise linear function consisting of two linear functions supported by [0,t̂ 0 ] and [t̂ 0 ,1] with t̂ 0 the optimal take-over point. It is shown that the optimal two-piece linear function is actually continuous at the take-over point, causing a reduction of the additional complexity in the resulting approximations of H 0 and H 1 . Furthermore, this allows analytic computation of the optimal two-piece linear function. By using the two-piece instead of the one-piece linear approximation, the root mean square approximation error is reduced by roughly a factor of 3 while the maximum approximation error is reduced by a factor of 4.5 for H 0 and of 2.6 for H 1 . Recursion relations satisfied by Struve functions, initialized with the approximations of H 0 and H 1 , yield approximations for higher order Struve functions.

Wave function collapse implies divergence of average displacement

OpenAIRE

Marchewka, A.; Schuss, Z.

2005-01-01

We show that propagating a truncated discontinuous wave function by Schr\\"odinger's equation, as asserted by the collapse axiom, gives rise to non-existence of the average displacement of the particle on the line. It also implies that there is no Zeno effect. On the other hand, if the truncation is done so that the reduced wave function is continuous, the average coordinate is finite and there is a Zeno effect. Therefore the collapse axiom of measurement needs to be revised.

Horizon wave-function and the quantum cosmic censorship

OpenAIRE

Casadio, RobertoDipartimento di Fisica e Astronomia, Alma Mater Università di Bologna, via Irnerio 46, Bologna, 40126, Italy; Micu, Octavian(Institute of Space Science, Bucharest, P.O. Box MG-23, Bucharest-Magurele, RO-077125, Romania); Stojkovic, Dejan(HEPCOS, Department of Physics, SUNY at Buffalo, Buffalo, NY, 14260-1500, United States)

2015-01-01

We investigate the Cosmic Censorship Conjecture by means of the horizon wave-function (HWF) formalism. We consider a charged massive particle whose quantum mechanical state is represented by a spherically symmetric Gaussian wave-function, and restrict our attention to the superxtremal case (with charge-to-mass ratio $\\alpha>1$), which is the prototype of a naked singularity in the classical theory. We find that one can still obtain a normalisable HWF for $\\alpha^2 2$, and the uncertainty in t...

Wave function of the Universe in the early stage of its evolution

International Nuclear Information System (INIS)

Maydanyuk, Sergei P.

2008-01-01

In quantum cosmological models, constructed in the framework of Friedmann-Robertson-Walker metrics, a nucleation of the Universe with its further expansion is described as a tunneling transition through an effective barrier between regions with small and large values of the scale factor a at non-zero (or zero) energy. The approach for describing this tunneling consists of constructing a wave function satisfying an appropriate boundary condition. There are various ways for defining the boundary condition that lead to different estimates of the barrier penetrability and the tunneling time. In order to describe the escape from the tunneling region as accurately as possible and to construct the total wave function on the basis of its two partial solutions unambiguously, we use the tunneling boundary condition that the total wave function must represent only the outgoing wave at the point of escape from the barrier, where the following definition for the wave is introduced: the wave is represented by the wave function whose modulus changes minimally under a variation of the scale factor a. We construct a new method for a direct non-semiclassical calculation of the total stationary wave function of the Universe, analyze the behavior of this wave function in the tunneling region, near the escape point and in the asymptotic region, and estimate the barrier penetrability. We observe oscillations of the modulus of the wave function in the external region starting from the turning point which decrease with increasing of a and which are not shown in semiclassical calculations. The period of such an oscillation decreases uniformly with increasing a and can be used as a fully quantum dynamical characteristic of the expansion of the Universe. (orig.)

Relativistic amplitudes in terms of wave functions

International Nuclear Information System (INIS)

Karmanov, V.A.

1978-01-01

In the framework of the invariant diagram technique which arises at the formulation of the fueld theory on the light front the question about conditions at which the relativistic amplitudes may be expressed through the wave functions is investigated. The amplitudes obtained depend on four-vector ω, determining the light front surface. The way is shown to find such values of the four-vector ω, at which the contribution of diagrams not expressed through wave functions is minimal. The investigation carried out is equivalent to the study of the dependence of amplitudes of the old-fashioned perturbation theory in the in the infinite momentum frame on direction of the infinite momentum

Antisymmetrized four-body wave function and coexistence of single particle and cluster structures

International Nuclear Information System (INIS)

Sasakawa, T.

1979-01-01

It is shown that each Yakubovski component of the totally antisymmetric four-body wave function satisfies the same equation as the unantisymmetric wave function. In the antisymmetric total wave function, the wave functions belonging to the same kind of partition are totally antisymmetric among themselves. This leads to the coexistence of cluster models, including the single particle model as a special case of the cluster model, as a sum

Approximation solutions for indifference pricing under general utility functions

NARCIS (Netherlands)

Chen, An; Pelsser, Antoon; Vellekoop, M.H.

2008-01-01

With the aid of Taylor-based approximations, this paper presents results for pricing insurance contracts by using indifference pricing under general utility functions. We discuss the connection between the resulting "theoretical" indifference prices and the pricing rule-of-thumb that practitioners

Approximate Solutions for Indifference Pricing under General Utility Functions

NARCIS (Netherlands)

Chen, A.; Pelsser, A.; Vellekoop, M.

2007-01-01

With the aid of Taylor-based approximations, this paper presents results for pricing insurance contracts by using indifference pricing under general utility functions. We discuss the connection between the resulting "theoretical" indifference prices and the pricing rule-of-thumb that practitioners

A wave propagation model of blood flow in large vessels using an approximate velocity profile function

NARCIS (Netherlands)

Bessems, D.; Rutten, M.C.M.; Vosse, van de F.N.

2007-01-01

Lumped-parameter models (zero-dimensional) and wave-propagation models (one-dimensional) for pressure and flow in large vessels, as well as fully three-dimensional fluid–structure interaction models for pressure and velocity, can contribute valuably to answering physiological and patho-physiological

Relativistic plasma dispersion functions

International Nuclear Information System (INIS)

Robinson, P.A.

1986-01-01

The known properties of plasma dispersion functions (PDF's) for waves in weakly relativistic, magnetized, thermal plasmas are reviewed and a large number of new results are presented. The PDF's required for the description of waves with small wave number perpendicular to the magnetic field (Dnestrovskii and Shkarofsky functions) are considered in detail; these functions also arise in certain quantum electrodynamical calculations involving strongly magnetized plasmas. Series, asymptotic series, recursion relations, integral forms, derivatives, differential equations, and approximations for these functions are discussed as are their analytic properties and connections with standard transcendental functions. In addition a more general class of PDF's relevant to waves of arbitrary perpendicular wave number is introduced and a range of properties of these functions are derived

The distribution of waves in the inner magnetosphere as a function of solar wind parameters

Science.gov (United States)

Aryan, Homayon; Balikhin, Michael A.; Agapitov, Oleksiy; Krasnoselskikh, Vladimir; Yearby, Keith

Energetic electrons within the Earth’s radiation belts represent a serious hazard to geostationary satellites. The interactions of electrons with chorus waves play an important role in both the acceleration and loss of radiation belt electrons. Studies of the evolution of energetic electron fluxes rely heavily on numerical codes in order to model energy and pitch angle diffusion due to electron interaction with plasma waves in the frame of quasilinear approximation. Application of these codes requires knowledge of statistical wave models to present wave distributions in the magnetosphere. A number of such models are based on CRESS, Cluster, THEMIS and other mission data. These models present wave distributions as a function of L-shell, magnetic local time, magnetic latitude and geomagnetic activity expressed by geomagnetic indices (Kp or Ae). However, it has been shown by G. Reeves and co-authors that only 50% of geomagnetic storms increase flux of relativistic electrons at GEO while 20% cause a decrease. This emphasizes the importance of including solar wind parameters in addition to geomagnetic indices. The present study examines almost four years (01, January, 2004 to 29, September, 2007) of STAFF (Spatio-Temporal Analysis of Field Fluctuation) data from Double Star TC1 combined with geomagnetic indices and solar wind parameters from OMNI database in order to present a comprehensive model of chorus wave intensities as a function of L-shell, magnetic local time, magnetic latitude, geomagnetic indices and solar wind parameters. The results show that chorus emission is not only sub-storm dependent but also dependent upon solar wind parameters with solar wind velocity evidently the most influential solar wind parameter. The largest peak intensities are observed for lower band chorus during active conditions, high solar wind velocity, low density and high pressure.

Wave functions constructed from an invariant sum over histories satisfy constraints

International Nuclear Information System (INIS)

Halliwell, J.J.; Hartle, J.B.

1991-01-01

Invariance of classical equations of motion under a group parametrized by functions of time implies constraints between canonical coordinates and momenta. In the Dirac formulation of quantum mechanics, invariance is normally imposed by demanding that physical wave functions are annihilated by the operator versions of these constraints. In the sum-over-histories quantum mechanics, however, wave functions are specified, directly, by appropriate functional integrals. It therefore becomes an interesting question whether the wave functions so specified obey the operator constraints of the Dirac theory. In this paper, we show for a wide class of theories, including gauge theories, general relativity, and first-quantized string theories, that wave functions constructed from a sum over histories are, in fact, annihilated by the constraints provided that the sum over histories is constructed in a manner which respects the invariance generated by the constraints. By this we mean a sum over histories defined with an invariant action, invariant measure, and an invariant class of paths summed over

Tropical Cyclogenesis in a Tropical Wave Critical Layer: Easterly Waves

Science.gov (United States)

Dunkerton, T. J.; Montgomery, M. T.; Wang, Z.

2009-01-01

The development of tropical depressions within tropical waves over the Atlantic and eastern Pacific is usually preceded by a "surface low along the wave" as if to suggest a hybrid wave-vortex structure in which flow streamlines not only undulate with the waves, but form a closed circulation in the lower troposphere surrounding the low. This structure, equatorward of the easterly jet axis, is identified herein as the familiar critical layer of waves in shear flow, a flow configuration which arguably provides the simplest conceptual framework for tropical cyclogenesis resulting from tropical waves, their interaction with the mean flow, and with diabatic processes associated with deep moist convection. The recirculating Kelvin cat's eye within the critical layer represents a sweet spot for tropical cyclogenesis in which a proto-vortex may form and grow within its parent wave. A common location for storm development is given by the intersection of the wave's critical latitude and trough axis at the center of the cat's eye, with analyzed vorticity centroid nearby. The wave and vortex live together for a time, and initially propagate at approximately the same speed. In most cases this coupled propagation continues for a few days after a tropical depression is identified. For easterly waves, as the name suggests, the propagation is westward. It is shown that in order to visualize optimally the associated Lagrangian motions, one should view the flow streamlines, or stream function, in a frame of reference translating horizontally with the phase propagation of the parent wave. In this co-moving frame, streamlines are approximately equivalent to particle trajectories. The closed circulation is quasi-stationary, and a dividing streamline separates air within the cat's eye from air outside.

General Relativistic Radiant Shock Waves in the Post-Quasistatic Approximation

Energy Technology Data Exchange (ETDEWEB)

H, Jorge A Rueda [Centro de Fisica Fundamental, Universidad de Los Andes, Merida 5101, Venezuela Escuela de Fisica, Universidad Industrial de Santander, A.A. 678, Bucaramanga (Colombia); Nunez, L A [Centro de Fisica Fundamental, Universidad de Los Andes, Merida 5101, Venezuela Centro Nacional de Calculo Cientifico, Universidad de Los Andes, CeCalCULA, Corporacion Parque Tecnologico de Merida, Merida 5101, Venezuela (Venezuela)

2007-05-15

An evolution of radiant shock wave front is considered in the framework of a recently presented method to study self-gravitating relativistic spheres, whose rationale becomes intelligible and finds full justification within the context of a suitable definition of the post-quasistatic approximation. The spherical matter configuration is divided into two regions by the shock and each side of the interface having a different equation of state and anisotropic phase. In order to simulate dissipation effects due to the transfer of photons and/or neutrinos within the matter configuration, we introduce the flux factor, the variable Eddington factor and a closure relation between them. As we expected the strong of the shock increases the speed of the fluid to relativistic ones and for some critical values is larger than light speed. In addition, we find that energy conditions are very sensible to the anisotropy, specially the strong energy condition. As a special feature of the model, we find that the contribution of the matter and radiation to the radial pressure are the same order of magnitude as in the mant as in the core, moreover, in the core radiation pressure is larger than matter pressure.

General Relativistic Radiant Shock Waves in the Post-Quasistatic Approximation

International Nuclear Information System (INIS)

H, Jorge A Rueda; Nunez, L A

2007-01-01

An evolution of radiant shock wave front is considered in the framework of a recently presented method to study self-gravitating relativistic spheres, whose rationale becomes intelligible and finds full justification within the context of a suitable definition of the post-quasistatic approximation. The spherical matter configuration is divided into two regions by the shock and each side of the interface having a different equation of state and anisotropic phase. In order to simulate dissipation effects due to the transfer of photons and/or neutrinos within the matter configuration, we introduce the flux factor, the variable Eddington factor and a closure relation between them. As we expected the strong of the shock increases the speed of the fluid to relativistic ones and for some critical values is larger than light speed. In addition, we find that energy conditions are very sensible to the anisotropy, specially the strong energy condition. As a special feature of the model, we find that the contribution of the matter and radiation to the radial pressure are the same order of magnitude as in the mant as in the core, moreover, in the core radiation pressure is larger than matter pressure

Irregular wave functions of a hydrogen atom in a uniform magnetic field

Science.gov (United States)

Wintgen, D.; Hoenig, A.

1989-01-01

The highly excited irregular wave functions of a hydrogen atom in a uniform magnetic field are investigated analytically, with wave function scarring by periodic orbits considered quantitatively. The results obtained confirm that the contributions of closed classical orbits to the spatial wave functions vanish in the semiclassical limit. Their disappearance, however, is slow. This discussion is illustrated by numerical examples.

Applications exponential approximation by integer shifts of Gaussian functions

Directory of Open Access Journals (Sweden)

S. M. Sitnik

2013-01-01

Full Text Available In this paper we consider approximations of functions using integer shifts of Gaussians – quadratic exponentials. A method is proposed to find coefficients of node functions by solving linear systems of equations. The explicit formula for the determinant of the system is found, based on it solvability of linear system under consideration is proved and uniqueness of its solution. We compare results with known ones and briefly indicate applications to signal theory.

Spin-Density Functionals from Current-Density Functional Theory and Vice Versa: A Road towards New Approximations

International Nuclear Information System (INIS)

Capelle, K.; Gross, E.

1997-01-01

It is shown that the exchange-correlation functional of spin-density functional theory is identical, on a certain set of densities, with the exchange-correlation functional of current-density functional theory. This rigorous connection is used to construct new approximations of the exchange-correlation functionals. These include a conceptually new generalized-gradient spin-density functional and a nonlocal current-density functional. copyright 1997 The American Physical Society

Prestack traveltime approximations

KAUST Repository

Alkhalifah, Tariq Ali

2011-01-01

Most prestack traveltime relations we tend work with are based on homogeneous (or semi-homogenous, possibly effective) media approximations. This includes the multi-focusing or double square-root (DSR) and the common reflection stack (CRS) equations. Using the DSR equation, I analyze the associated eikonal form in the general source-receiver domain. Like its wave-equation counterpart, it suffers from a critical singularity for horizontally traveling waves. As a result, I derive expansion based solutions of this eikonal based on polynomial expansions in terms of the reflection and dip angles in a generally inhomogenous background medium. These approximate solutions are free of singularities and can be used to estimate travetimes for small to moderate offsets (or reflection angles) in a generally inhomogeneous medium. A Marmousi example demonstrates the usefulness of the approach. © 2011 Society of Exploration Geophysicists.

Response functions of free mass gravitational wave antennas

Science.gov (United States)

Estabrook, F. B.

1985-01-01

The work of Gursel, Linsay, Spero, Saulson, Whitcomb and Weiss (1984) on the response of a free-mass interferometric antenna is extended. Starting from first principles, the earlier work derived the response of a 2-arm gravitational wave antenna to plane polarized gravitational waves. Equivalent formulas (generalized slightly to allow for arbitrary elliptical polarization) are obtained by a simple differencing of the '3-pulse' Doppler response functions of two 1-arm antennas. A '4-pulse' response function is found, with quite complicated angular dependences for arbitrary incident polarization. The differencing method can as readily be used to write exact response functions ('3n+1 pulse') for antennas having multiple passes or more arms.

A point-value enhanced finite volume method based on approximate delta functions

Science.gov (United States)

Xuan, Li-Jun; Majdalani, Joseph

2018-02-01

We revisit the concept of an approximate delta function (ADF), introduced by Huynh (2011) [1], in the form of a finite-order polynomial that holds identical integral properties to the Dirac delta function when used in conjunction with a finite-order polynomial integrand over a finite domain. We show that the use of generic ADF polynomials can be effective at recovering and generalizing several high-order methods, including Taylor-based and nodal-based Discontinuous Galerkin methods, as well as the Correction Procedure via Reconstruction. Based on the ADF concept, we then proceed to formulate a Point-value enhanced Finite Volume (PFV) method, which stores and updates the cell-averaged values inside each element as well as the unknown quantities and, if needed, their derivatives on nodal points. The sharing of nodal information with surrounding elements saves the number of degrees of freedom compared to other compact methods at the same order. To ensure conservation, cell-averaged values are updated using an identical approach to that adopted in the finite volume method. Here, the updating of nodal values and their derivatives is achieved through an ADF concept that leverages all of the elements within the domain of integration that share the same nodal point. The resulting scheme is shown to be very stable at successively increasing orders. Both accuracy and stability of the PFV method are verified using a Fourier analysis and through applications to the linear wave and nonlinear Burgers' equations in one-dimensional space.

Analytic perturbation theory for screened Coulomb potential: full continuum wave function

International Nuclear Information System (INIS)

Bechler, A.; Ennan, Mc J.; Pratt, R.H.

1979-01-01

An analytic perturbation theory developed previously is used to find a continuum screened-Coulomb wave function characterized by definite asymptotic momentum. This wave function satisfies an inhomogeneous partial differential equation which is solved in parabolic coordinates; the solution depends on both parabolic variables. We calculate partial wave projections of this solution and show that we can choose to add a solution of the homogeneous equation such that the partial wave projections become equal to the normalized continuum radial function found previously. However, finding the unique solution with given asymptotic linear momentum will require either using boundary conditions to determine the unique needed solution of the homogeneous equation or equivalently specifying the screened-Coulomb phase-shifts. (author)

Methods and Algorithms for Approximating the Gamma Function and Related Functions. A survey. Part I: Asymptotic Series

Directory of Open Access Journals (Sweden)

Cristinel Mortici

2015-01-01

Full Text Available In this survey we present our recent results on analysis of gamma function and related functions. The results obtained are in the theory of asymptotic analysis, approximation of gamma and polygamma functions, or in the theory of completely monotonic functions. The motivation of this first part is the work of C. Mortici [Product Approximations via Asymptotic Integration Amer. Math. Monthly 117 (2010 434-441] where a simple strategy for constructing asymptotic series is presented. The classical asymptotic series associated to Stirling, Wallis, Glaisher-Kinkelin are rediscovered. In the second section we discuss some new inequalities related to Landau constants and we establish some asymptotic formulas.

Local density approximation for exchange in excited-state density functional theory

OpenAIRE

Harbola, Manoj K.; Samal, Prasanjit

2004-01-01

Local density approximation for the exchange energy is made for treatment of excited-states in density-functional theory. It is shown that taking care of the state-dependence of the LDA exchange energy functional leads to accurate excitation energies.

On Born approximation in black hole scattering

Science.gov (United States)

Batic, D.; Kelkar, N. G.; Nowakowski, M.

2011-12-01

A massless field propagating on spherically symmetric black hole metrics such as the Schwarzschild, Reissner-Nordström and Reissner-Nordström-de Sitter backgrounds is considered. In particular, explicit formulae in terms of transcendental functions for the scattering of massless scalar particles off black holes are derived within a Born approximation. It is shown that the conditions on the existence of the Born integral forbid a straightforward extraction of the quasi normal modes using the Born approximation for the scattering amplitude. Such a method has been used in literature. We suggest a novel, well defined method, to extract the large imaginary part of quasinormal modes via the Coulomb-like phase shift. Furthermore, we compare the numerically evaluated exact scattering amplitude with the Born one to find that the approximation is not very useful for the scattering of massless scalar, electromagnetic as well as gravitational waves from black holes.

Short time propagation of a singular wave function: Some surprising results

Science.gov (United States)

Marchewka, A.; Granot, E.; Schuss, Z.

2007-08-01

The Schrödinger evolution of an initially singular wave function was investigated. First it was shown that a wide range of physical problems can be described by initially singular wave function. Then it was demonstrated that outside the support of the initial wave function the time evolution is governed to leading order by the values of the wave function and its derivatives at the singular points. Short-time universality appears where it depends only on a single parameter—the value at the singular point (not even on its derivatives). It was also demonstrated that the short-time evolution in the presence of an absorptive potential is different than in the presence of a nonabsorptive one. Therefore, this dynamics can be harnessed to the determination whether a potential is absorptive or not simply by measuring only the transmitted particles density.

New approximation for Glauber theory on stripping of relativistic deuterons

International Nuclear Information System (INIS)

Nissen-Meyer, S.A.

1978-03-01

The momentum distribution of forward protons from relativistic collisions of deuterons with nuclei is computed from a Glauber theoretical Ansatz of Bertocchi and Tekou. The outgoing proton-neutron scattering state (disintegrated deuteron) with a plane wave minus the components of this plane wave along the deuteron bound state vector is approximated. With no fitted parameters good agreement is found with data from the reaction d + C 12 → p + X in the region corresponding to nonrelativistic Fermi momenta in the forward direction. At more relativistic Fermi momenta, the model deviates more from the data, which can be due to incorrect choice of the short distance part of the deuteron wave function as well as off-shell effects in the deuteron

Slow Growth and Optimal Approximation of Pseudoanalytic Functions on the Disk

Directory of Open Access Journals (Sweden)

Devendra Kumar

2013-07-01

Full Text Available Pseudoanalytic functions (PAF are constructed as complex combination of real-valued analytic solutions to the Stokes-Betrami System. These solutions include the generalized biaxisymmetric potentials. McCoy [10] considered the approximation of pseudoanalytic functions on the disk. Kumar et al. [9] studied the generalized order and generalized type of PAF in terms of the Fourier coefficients occurring in its local expansion and optimal approximation errors in Bernstein sense on the disk. The aim of this paper is to improve the results of McCoy [10] and Kumar et al. [9]. Our results apply satisfactorily for slow growth.

Linear Water Waves

Science.gov (United States)

Kuznetsov, N.; Maz'ya, V.; Vainberg, B.

2002-08-01

This book gives a self-contained and up-to-date account of mathematical results in the linear theory of water waves. The study of waves has many applications, including the prediction of behavior of floating bodies (ships, submarines, tension-leg platforms etc.), the calculation of wave-making resistance in naval architecture, and the description of wave patterns over bottom topography in geophysical hydrodynamics. The first section deals with time-harmonic waves. Three linear boundary value problems serve as the approximate mathematical models for these types of water waves. The next section uses a plethora of mathematical techniques in the investigation of these three problems. The techniques used in the book include integral equations based on Green's functions, various inequalities between the kinetic and potential energy and integral identities which are indispensable for proving the uniqueness theorems. The so-called inverse procedure is applied to constructing examples of non-uniqueness, usually referred to as 'trapped nodes.'

Deorbitalization strategies for meta-generalized-gradient-approximation exchange-correlation functionals

Science.gov (United States)

Mejia-Rodriguez, Daniel; Trickey, S. B.

2017-11-01

We explore the simplification of widely used meta-generalized-gradient approximation (mGGA) exchange-correlation functionals to the Laplacian level of refinement by use of approximate kinetic-energy density functionals (KEDFs). Such deorbitalization is motivated by the prospect of reducing computational cost while recovering a strictly Kohn-Sham local potential framework (rather than the usual generalized Kohn-Sham treatment of mGGAs). A KEDF that has been rather successful in solid simulations proves to be inadequate for deorbitalization, but we produce other forms which, with parametrization to Kohn-Sham results (not experimental data) on a small training set, yield rather good results on standard molecular test sets when used to deorbitalize the meta-GGA made very simple, Tao-Perdew-Staroverov-Scuseria, and strongly constrained and appropriately normed functionals. We also study the difference between high-fidelity and best-performing deorbitalizations and discuss possible implications for use in ab initio molecular dynamics simulations of complicated condensed phase systems.

A full scale approximation of covariance functions for large spatial data sets

KAUST Repository

Sang, Huiyan

2011-10-10

Gaussian process models have been widely used in spatial statistics but face tremendous computational challenges for very large data sets. The model fitting and spatial prediction of such models typically require O(n 3) operations for a data set of size n. Various approximations of the covariance functions have been introduced to reduce the computational cost. However, most existing approximations cannot simultaneously capture both the large- and the small-scale spatial dependence. A new approximation scheme is developed to provide a high quality approximation to the covariance function at both the large and the small spatial scales. The new approximation is the summation of two parts: a reduced rank covariance and a compactly supported covariance obtained by tapering the covariance of the residual of the reduced rank approximation. Whereas the former part mainly captures the large-scale spatial variation, the latter part captures the small-scale, local variation that is unexplained by the former part. By combining the reduced rank representation and sparse matrix techniques, our approach allows for efficient computation for maximum likelihood estimation, spatial prediction and Bayesian inference. We illustrate the new approach with simulated and real data sets. © 2011 Royal Statistical Society.

A full scale approximation of covariance functions for large spatial data sets

KAUST Repository

Sang, Huiyan; Huang, Jianhua Z.

2011-01-01

Gaussian process models have been widely used in spatial statistics but face tremendous computational challenges for very large data sets. The model fitting and spatial prediction of such models typically require O(n 3) operations for a data set of size n. Various approximations of the covariance functions have been introduced to reduce the computational cost. However, most existing approximations cannot simultaneously capture both the large- and the small-scale spatial dependence. A new approximation scheme is developed to provide a high quality approximation to the covariance function at both the large and the small spatial scales. The new approximation is the summation of two parts: a reduced rank covariance and a compactly supported covariance obtained by tapering the covariance of the residual of the reduced rank approximation. Whereas the former part mainly captures the large-scale spatial variation, the latter part captures the small-scale, local variation that is unexplained by the former part. By combining the reduced rank representation and sparse matrix techniques, our approach allows for efficient computation for maximum likelihood estimation, spatial prediction and Bayesian inference. We illustrate the new approach with simulated and real data sets. © 2011 Royal Statistical Society.

The Positron-Electron Correlation Energy In ZnO Calculated With The Modified Single Wave Function Of Positron

International Nuclear Information System (INIS)

Chau Van Tao; Trinh Hoa Lang; Le Hoang Chien; Nguyen Huu Loc; Nguyen Anh Tuan

2011-01-01

Positron-electron correlation energy of the ZnO - positron system is studied on assumption that positron binds with the outer shell electrons of Zinc and Oxygen to form the pseudo ZnO - positron molecule before it annihilates with one of these electrons. In this work, the single wave function for positron is form by LCAO approximation and is modified according to the principle of linear superposition, and by using Variational Quantum Monte Carlo method (VQMC) [7] the correlation energy of this system is estimated with the value E c e-p = - 9.3 ± 1.1 eV. In the theoretical aspect it turns out that this result is more reasonable and closer to those of other methods [3] than the one which is done without modifying the wave function of positron [1]. To confirm this legitimate approach, however, the further calculations of positron annihilation rate in ZnO have to be carried out in our next work. (author)

Anomalies of the free loop wave equation in the WKB approximation

International Nuclear Information System (INIS)

Weisz, P.; Luescher, M.; Symanzik, K.

1980-04-01

We derive a well-defined, reparametrization invariant expression for the next to leading term in the small h/2π expansion of the Euclidean loop Green's functional PSI(C). To this order in h/2π, we then verify that PSI(C) satisfies a renormalized loop wave equation, which involves a number of local, but non-harmonic anomalous terms. Also, we find that the quantum fluctuations of the string give rise, in 3 + 1 dimensions, to a correction of the static quark potential by an attractive Coulomb potential of universal strength αsub(string) = π/12. (orig.)

Complex-valued derivative propagation method with approximate Bohmian trajectories: Application to electronic nonadiabatic dynamics

Science.gov (United States)

Wang, Yu; Chou, Chia-Chun

2018-05-01

The coupled complex quantum Hamilton-Jacobi equations for electronic nonadiabatic transitions are approximately solved by propagating individual quantum trajectories in real space. Equations of motion are derived through use of the derivative propagation method for the complex actions and their spatial derivatives for wave packets moving on each of the coupled electronic potential surfaces. These equations for two surfaces are converted into the moving frame with the same grid point velocities. Excellent wave functions can be obtained by making use of the superposition principle even when nodes develop in wave packet scattering.

Finite-rank potential that reproduces the Pade approximant

International Nuclear Information System (INIS)

Tani, S.

1979-01-01

If a scattering potential is of a finite rank, say N, the exact solution of the problem can be obtained from the Born series, if the potential strength is within the radius of convergence; the exact solution can be obtained from the analytical continuation of the formal Born series outside the radius of convergence. Beyond the first 2N terms of the Born series, an individual term of the Born series depends on the first 2N terms, and the [N/N] Pade approximant and the exact solution agree with each other. The above-mentioned features of a finite-rank problem are relevant to scattering theory in general, because most scattering problems may be handled as an extension of the rank-N problem, in which the rank N tends to infinity. The foregoing aspects of scattering theory will be studied in depth in the present paper, and in so doing we proceed in the opposite direction. Namely, given a potential, we calculate the first 2N terms of the Born series for the K matrix and the first N terms of the Born series for the wave function. Using these data, a special rank-N potential is constructed in such a way that it reproduces the [N/N] Pade approximant of the K matrix of the original scattering problem. One great advantage of obtaining such a rank-N potential is that the wave function of the system may be approximated in the same spirit as done for the K matrix; hence, we can introduce a new approximation method for dealing with an off-shell T matrix. A part of the mathematical work is incomplete, but the physical aspects are thoroughly discussed

Near-Shore Floating Wave Energy Converters

DEFF Research Database (Denmark)

Ruol, Piero; Zanuttigh, Barbara; Martinelli, Luca

2011-01-01

and transmission characteristics are approximated to functions of wave height, period and obliquity. Their order of magnitude are 20% and 80%, respectively. It is imagined that an array of DEXA is deployed in front of Marina di Ravenna beach (IT), a highly touristic site of the Adriatic Coast. Based on the CERC......Aim of this note is to analyse the possible application of a Wave Energy Converter (WEC) as a combined tool to protect the coast and harvest energy. Physical model tests are used to evaluate wave transmission past a near-shore floating WEC of the wave activated body type, named DEXA. Efficiency...

Validity of the broken-pair approximation for N = 50, even-A nuclei

International Nuclear Information System (INIS)

Haq, S.; Gambhir, Y.K.

1977-01-01

The validity of the broken-pair approximation as an approximation to the seniority shell model is investigated. The results of the broken-pair approximation and the seniority shell model, obtained by employing identical input information (single-particle levels and their energies, effective two-body matrix elements, 88 Sr inert core) for N = 50, even-A nuclei are compared. A close agreement obtained between the calculated broken-pair approximation and the seniority shell model energies for 90 Zr, 92 Mo, 94 Ru, and 96 Pd nuclei and large (95--100 %) overlaps between the broken-pair approximation and the senority shell model wave functions for 92 Mo, demonstrates the validity of the broken-pair approximation in this region and in general its usefulness as a good approximation to the seniority shell model

Period functions for Maass wave forms and cohomology

CERN Document Server

Bruggeman, R; Zagier, D; Bruggeman, R W; Zagier, D

2015-01-01

The authors construct explicit isomorphisms between spaces of Maass wave forms and cohomology groups for discrete cofinite groups \\Gamma\\subset\\mathrm{PSL}_2({\\mathbb{R}}). In the case that \\Gamma is the modular group \\mathrm{PSL}_2({\\mathbb{Z}}) this gives a cohomological framework for the results in Period functions for Maass wave forms. I, of J. Lewis and D. Zagier in Ann. Math. 153 (2001), 191-258, where a bijection was given between cuspidal Maass forms and period functions. The authors introduce the concepts of mixed parabolic cohomology group and semi-analytic vectors in principal serie

Evolution of wave function in a dissipative system

Science.gov (United States)

Yu, Li-Hua; Sun, Chang-Pu

1994-01-01

For a dissipative system with Ohmic friction, we obtain a simple and exact solution for the wave function of the system plus the bath. It is described by the direct product in two independent Hilbert space. One of them is described by an effective Hamiltonian, the other represents the effect of the bath, i.e., the Brownian motion, thus clarifying the structure of the wave function of the system whose energy is dissipated by its interaction with the bath. No path integral technology is needed in this treatment. The derivation of the Weisskopf-Wigner line width theory follows easily.

Cut-off frequencies of circumferential horizontal shear waves in various functionally graded cylinder shells.

Science.gov (United States)

Shen, Xiaoqin; Ren, Dawei; Cao, Xiaoshan; Wang, Ji

2018-03-01

In this study, cut-off frequencies of the circumferential SH waves in functionally graded piezoelectric-piezomagnetic material (FGPPM) cylinder shells with traction free, electrical and magnetic open boundary conditions are investigated analytically. The Wentzel-Kramers-Brillouin (WKB) method is employed for solving differential equations with variable coefficients for general cases. For comparison, Bessel functions and Kummer functions are used for solving cut-off frequency problems in homogenous and ideal FGPPM cylinder shells. It is shown that the WKB solution for the cut-off frequencies has good precise. The set of cut-off frequencies is a series of approximate arithmetic progressions, for which the difference is a function of the density and the effective elastic parameter. The relationship between the difference and the gradient coefficient is described. These results provide theoretical guidance for the non-destructive evaluation of curved shells based on the cut-off frequencies. Copyright © 2017 Elsevier B.V. All rights reserved.

Correlated wave functions for three-particle systems with Coulomb interaction - The muonic helium atom

Science.gov (United States)

Huang, K.-N.

1977-01-01

A computational procedure for calculating correlated wave functions is proposed for three-particle systems interacting through Coulomb forces. Calculations are carried out for the muonic helium atom. Variational wave functions which explicitly contain interparticle coordinates are presented for the ground and excited states. General Hylleraas-type trial functions are used as the basis for the correlated wave functions. Excited-state energies of the muonic helium atom computed from 1- and 35-term wave functions are listed for four states.

Longitudinal wave function control in single quantum dots with an applied magnetic field

Science.gov (United States)

Cao, Shuo; Tang, Jing; Gao, Yunan; Sun, Yue; Qiu, Kangsheng; Zhao, Yanhui; He, Min; Shi, Jin-An; Gu, Lin; Williams, David A.; Sheng, Weidong; Jin, Kuijuan; Xu, Xiulai

2015-01-01

Controlling single-particle wave functions in single semiconductor quantum dots is in demand to implement solid-state quantum information processing and spintronics. Normally, particle wave functions can be tuned transversely by an perpendicular magnetic field. We report a longitudinal wave function control in single quantum dots with a magnetic field. For a pure InAs quantum dot with a shape of pyramid or truncated pyramid, the hole wave function always occupies the base because of the less confinement at base, which induces a permanent dipole oriented from base to apex. With applying magnetic field along the base-apex direction, the hole wave function shrinks in the base plane. Because of the linear changing of the confinement for hole wave function from base to apex, the center of effective mass moves up during shrinking process. Due to the uniform confine potential for electrons, the center of effective mass of electrons does not move much, which results in a permanent dipole moment change and an inverted electron-hole alignment along the magnetic field direction. Manipulating the wave function longitudinally not only provides an alternative way to control the charge distribution with magnetic field but also a new method to tune electron-hole interaction in single quantum dots. PMID:25624018

Longitudinal wave function control in single quantum dots with an applied magnetic field.

Science.gov (United States)

Cao, Shuo; Tang, Jing; Gao, Yunan; Sun, Yue; Qiu, Kangsheng; Zhao, Yanhui; He, Min; Shi, Jin-An; Gu, Lin; Williams, David A; Sheng, Weidong; Jin, Kuijuan; Xu, Xiulai

2015-01-27

Controlling single-particle wave functions in single semiconductor quantum dots is in demand to implement solid-state quantum information processing and spintronics. Normally, particle wave functions can be tuned transversely by an perpendicular magnetic field. We report a longitudinal wave function control in single quantum dots with a magnetic field. For a pure InAs quantum dot with a shape of pyramid or truncated pyramid, the hole wave function always occupies the base because of the less confinement at base, which induces a permanent dipole oriented from base to apex. With applying magnetic field along the base-apex direction, the hole wave function shrinks in the base plane. Because of the linear changing of the confinement for hole wave function from base to apex, the center of effective mass moves up during shrinking process. Due to the uniform confine potential for electrons, the center of effective mass of electrons does not move much, which results in a permanent dipole moment change and an inverted electron-hole alignment along the magnetic field direction. Manipulating the wave function longitudinally not only provides an alternative way to control the charge distribution with magnetic field but also a new method to tune electron-hole interaction in single quantum dots.

Wave drag as the objective function in transonic fighter wing optimization

Science.gov (United States)

Phillips, P. S.

1984-01-01

The original computational method for determining wave drag in a three dimensional transonic analysis method was replaced by a wave drag formula based on the loss in momentum across an isentropic shock. This formula was used as the objective function in a numerical optimization procedure to reduce the wave drag of a fighter wing at transonic maneuver conditions. The optimization procedure minimized wave drag through modifications to the wing section contours defined by a wing profile shape function. A significant reduction in wave drag was achieved while maintaining a high lift coefficient. Comparisons of the pressure distributions for the initial and optimized wing geometries showed significant reductions in the leading-edge peaks and shock strength across the span.

Are there approximate relations among transverse momentum dependent distribution functions?

Energy Technology Data Exchange (ETDEWEB)

Harutyun AVAKIAN; Anatoli Efremov; Klaus Goeke; Andreas Metz; Peter Schweitzer; Tobias Teckentrup

2007-10-11

Certain {\\sl exact} relations among transverse momentum dependent parton distribution functions due to QCD equations of motion turn into {\\sl approximate} ones upon the neglect of pure twist-3 terms. On the basis of available data from HERMES we test the practical usefulness of one such ``Wandzura-Wilczek-type approximation'', namely of that connecting $h_{1L}^{\\perp(1)a}(x)$ to $h_L^a(x)$, and discuss how it can be further tested by future CLAS and COMPASS data.

Discrete expansions of continuum wave functions

International Nuclear Information System (INIS)

Bang, J.; Ershov, S.N.; Gareev, F.A.; Kazacha, G.S.

1980-01-01

Different methods of expanding continuum wave functions in terms of discrete basis sets are discussed. The convergence properties of these expansions are investigated, both from a mathematical and a numerical point of view, for the case of potentials of Woods-Saxon and square well type. (orig.)

Nonlinear instability and chaos in plasma wave-wave interactions, I., Introduction

International Nuclear Information System (INIS)

Kueny, C.S.; Morrison, P.J.

1994-11-01

Conventional linear stability analyses may fail for fluid systems with an indefinite free energy functional. When such a system is linearly stable, it is said to possess negative energy modes. Instability may then occur either via dissipation of the negative energy modes, or nonlinearly via resonant wave-wave coupling, leading to explosive growth. In the dissipationless case, it is conjectured that intrinsic chaotic behavior may allow initially nonresonant systems to reach resonance by diffusion in phase space. In this and a companion paper [submitted to Physics of Plasmas], this phenomenon is demonstrated for a simple equilibrium involving cold counterstreaming ions. The system is described in the fluid approximation by a Hamiltonian functional and associated noncanonical Poisson bracket. By Fourier decomposition and appropriate coordinate transformations, the Hamiltonian for the perturbed energy is expressed in action-angle form. The normal modes correspond to Doppler-shifted ion-acoustic waves of positive and negative energy. Nonlinear coupling leads to decay instability via two-wave interactions, and to either decay or explosive instability via three-wave interactions. These instabilities are described for various (integrable) systems of waves interacting via single nonlinear terms. This discussion provides the foundation for the treatment of nonintegrable systems in the companion paper

Nonlinear instability and chaos in plasma wave--wave interactions. I. Introduction

International Nuclear Information System (INIS)

Kueny, C.S.; Morrison, P.J.

1995-01-01

Conventional linear stability analyses may fail for fluid systems with an indefinite free-energy functional. When such a system is linearly stable, it is said to possess negative energy modes. Instability may then occur either via dissipation of the negative energy modes, or nonlinearly via resonant wave--wave coupling, leading to explosive growth. In the dissipationless case, it is conjectured that intrinsic chaotic behavior may allow initially nonresonant systems to reach resonance by diffusion in phase space. In this and a companion paper (submitted to Phys. Plasmas), this phenomenon is demonstrated for a simple equilibrium involving cold counterstreaming ions. The system is described in the fluid approximation by a Hamiltonian functional and associated noncanonical Poisson bracket. By Fourier decomposition and appropriate coordinate transformations, the Hamiltonian for the perturbed energy is expressed in action-angle form. The normal modes correspond to Doppler-shifted ion-acoustic waves of positive and negative energy. Nonlinear coupling leads to decay instability via two-wave interactions, and to either decay or explosive instability via three-wave interactions. These instabilities are described for various integrable systems of waves interacting via single nonlinear terms. This discussion provides the foundation for the treatment of nonintegrable systems in the companion paper. copyright 1995 American Institute of Physics

Piecewise quadratic Lyapunov functions for stability verification of approximate explicit MPC

Directory of Open Access Journals (Sweden)

Morten Hovd

2010-04-01

Full Text Available Explicit MPC of constrained linear systems is known to result in a piecewise affine controller and therefore also piecewise affine closed loop dynamics. The complexity of such analytic formulations of the control law can grow exponentially with the prediction horizon. The suboptimal solutions offer a trade-off in terms of complexity and several approaches can be found in the literature for the construction of approximate MPC laws. In the present paper a piecewise quadratic (PWQ Lyapunov function is used for the stability verification of an of approximate explicit Model Predictive Control (MPC. A novel relaxation method is proposed for the LMI criteria on the Lyapunov function design. This relaxation is applicable to the design of PWQ Lyapunov functions for discrete-time piecewise affine systems in general.

Linear density response function in the projector augmented wave method

DEFF Research Database (Denmark)

Yan, Jun; Mortensen, Jens Jørgen; Jacobsen, Karsten Wedel

2011-01-01

We present an implementation of the linear density response function within the projector-augmented wave method with applications to the linear optical and dielectric properties of both solids, surfaces, and interfaces. The response function is represented in plane waves while the single...... functions of Si, C, SiC, AlP, and GaAs compare well with previous calculations. While optical properties of semiconductors, in particular excitonic effects, are generally not well described by ALDA, we obtain excellent agreement with experiments for the surface loss function of graphene and the Mg(0001...

Hierarchical low-rank approximation for high dimensional approximation

KAUST Repository

Nouy, Anthony

2016-01-07

Tensor methods are among the most prominent tools for the numerical solution of high-dimensional problems where functions of multiple variables have to be approximated. Such high-dimensional approximation problems naturally arise in stochastic analysis and uncertainty quantification. In many practical situations, the approximation of high-dimensional functions is made computationally tractable by using rank-structured approximations. In this talk, we present algorithms for the approximation in hierarchical tensor format using statistical methods. Sparse representations in a given tensor format are obtained with adaptive or convex relaxation methods, with a selection of parameters using crossvalidation methods.

Hierarchical low-rank approximation for high dimensional approximation

KAUST Repository

Nouy, Anthony

2016-01-01

Tensor methods are among the most prominent tools for the numerical solution of high-dimensional problems where functions of multiple variables have to be approximated. Such high-dimensional approximation problems naturally arise in stochastic analysis and uncertainty quantification. In many practical situations, the approximation of high-dimensional functions is made computationally tractable by using rank-structured approximations. In this talk, we present algorithms for the approximation in hierarchical tensor format using statistical methods. Sparse representations in a given tensor format are obtained with adaptive or convex relaxation methods, with a selection of parameters using crossvalidation methods.

Order in large and chaos in small components of nuclear wave functions

International Nuclear Information System (INIS)

Soloviev, V.G.

1992-06-01

An investigation of the order and chaos of the nuclear excited states has shown that there is order in the large and chaos in the small quasiparticle or phonon components of the nuclear wave functions. The order-to-chaos transition is treated as a transition from the large to the small components of the nuclear wave function. The analysis has shown that relatively large many-quasiparticle components of the wave function at an excitation energy (4-8)MeV may exist. The large many-quasiparticle components of the wave functions of the neutron resonances are responsible for enhanced E1-, M1- and E2-transition probabilities from neutron resonance to levels lying (1-2)MeV below them. (author)

Performance of wave function and density functional methods for water hydrogen bond spin-spin coupling constants.

Science.gov (United States)

García de la Vega, J M; Omar, S; San Fabián, J

2017-04-01

Spin-spin coupling constants in water monomer and dimer have been calculated using several wave function and density functional-based methods. CCSD, MCSCF, and SOPPA wave functions methods yield similar results, specially when an additive approach is used with the MCSCF. Several functionals have been used to analyze their performance with the Jacob's ladder and a set of functionals with different HF exchange were tested. Functionals with large HF exchange appropriately predict 1 J O H , 2 J H H and 2h J O O couplings, while 1h J O H is better calculated with functionals that include a reduced fraction of HF exchange. Accurate functionals for 1 J O H and 2 J H H have been tested in a tetramer water model. The hydrogen bond effects on these intramolecular couplings are additive when they are calculated by SOPPA(CCSD) wave function and DFT methods. Graphical Abstract Evaluation of the additive effect of the hydrogen bond on spin-spin coupling constants of water using WF and DFT methods.

Generalized squeezing rotating-wave approximation to the isotropic and anisotropic Rabi model in the ultrastrong-coupling regime

OpenAIRE

Zhang, Yu-Yu

2016-01-01

Generalized squeezing rotating-wave approximation (GSRWA) is proposed by employing both the displacement and the squeezing transformations. A solvable Hamiltonian is reformulated in the same form as the ordinary RWA ones. For a qubit coupled to oscillators experiment, a well-defined Schr\\"{o}dinger-cat-like entangled state is given by the displaced-squeezed oscillator state instead of the original displaced state. For the isotropic Rabi case, the mean photon number and the ground-state energy...

Evolution of the orbitals Dy-4f in the DyB2 compound using the LDA, PBE approximations, and the PBE0 hybrid functional

Science.gov (United States)

Rasero Causil, Diego; Ortega López, César; Espitia Rico, Miguel

2018-04-01

Computational calculations of total energy based on density functional theory were used to investigate the structural, electronic, and magnetic properties of the DyB2 compounds in the hexagonal structure. The calculations were carried out by means of the full-potential linearized augmented plane wave (FP-LAPW) method, employing the computational Wien2k package. The local density approximation (LDA) and the generalized gradient approximation (GGA) were used for the electron-electron interactions. Additionally, we used the functional hybrid PBE0 for a better description the electronic and magnetic properties, because the DyB2 compound is a strongly-correlated system. We found that the calculated lattice constant agrees well with the values reported theoretically and experimentally. The density of states (DOS) calculation shows that the compound exhibits a metallic behavior and has magnetic properties, with a total magnetic moment of 5.47 μ0/cell determined mainly by the 4f states of the rare earth elements. The functional PBE0 shows a strong localization of the Dy-4f orbitals.

Approximation of functions in two variables by some linear positive operators

Directory of Open Access Journals (Sweden)

Mariola Skorupka

1995-12-01

Full Text Available We introduce some linear positive operators of the Szasz-Mirakjan type in the weighted spaces of continuous functions in two variables. We study the degree of the approximation of functions by these operators. The similar results for functions in one variable are given in [5]. Some operators of the Szasz-Mirakjan type are examined also in [3], [4].

Special software for computing the special functions of wave catastrophes

Directory of Open Access Journals (Sweden)

Andrey S. Kryukovsky

2015-01-01

Full Text Available The method of ordinary differential equations in the context of calculating the special functions of wave catastrophes is considered. Complementary numerical methods and algorithms are described. The paper shows approaches to accelerate such calculations using capabilities of modern computing systems. Methods for calculating the special functions of wave catastrophes are considered in the framework of parallel computing and distributed systems. The paper covers the development process of special software for calculating of special functions, questions of portability, extensibility and interoperability.

Propagating wave correlations in complex systems

International Nuclear Information System (INIS)

Creagh, Stephen C; Gradoni, Gabriele; Hartmann, Timo; Tanner, Gregor

2017-01-01

We describe a novel approach for computing wave correlation functions inside finite spatial domains driven by complex and statistical sources. By exploiting semiclassical approximations, we provide explicit algorithms to calculate the local mean of these correlation functions in terms of the underlying classical dynamics. By defining appropriate ensemble averages, we show that fluctuations about the mean can be characterised in terms of classical correlations. We give in particular an explicit expression relating fluctuations of diagonal contributions to those of the full wave correlation function. The methods have a wide range of applications both in quantum mechanics and for classical wave problems such as in vibro-acoustics and electromagnetism. We apply the methods here to simple quantum systems, so-called quantum maps, which model the behaviour of generic problems on Poincaré sections. Although low-dimensional, these models exhibit a chaotic classical limit and share common characteristics with wave propagation in complex structures. (paper)

Inverse Schroedinger equation and the exact wave function

International Nuclear Information System (INIS)

Nakatsuji, Hiroshi

2002-01-01

Using the inverse of the Hamiltonian, we introduce the inverse Schroedinger equation (ISE) that is equivalent to the ordinary Schroedinger equation (SE). The ISE has the variational principle and the H-square group of equations as the SE has. When we use a positive Hamiltonian, shifting the energy origin, the inverse energy becomes monotonic and we further have the inverse Ritz variational principle and cross-H-square equations. The concepts of the SE and the ISE are combined to generalize the theory for calculating the exact wave function that is a common eigenfunction of the SE and ISE. The Krylov sequence is extended to include the inverse Hamiltonian, and the complete Krylov sequence is introduced. The iterative configuration interaction (ICI) theory is generalized to cover both the SE and ISE concepts and four different computational methods of calculating the exact wave function are presented in both analytical and matrix representations. The exact wave-function theory based on the inverse Hamiltonian can be applied to systems that have singularities in the Hamiltonian. The generalized ICI theory is applied to the hydrogen atom, giving the exact solution without any singularity problem

A surprise in the first Born approximation for electron scattering

International Nuclear Information System (INIS)

Treacy, M.M.J.; Van Dyck, D.

2012-01-01

A standard textbook derivation for the scattering of electrons by a weak potential under the first Born approximation suggests that the far-field scattered wave should be in phase with the incident wave. However, it is well known that waves scattered from a weak phase object should be phase-shifted by π/2 relative to the incident wave. A disturbing consequence of this missing phase is that, according to the Optical Theorem, the total scattering cross section would be zero in the first Born approximation. We resolve this mystery pedagogically by showing that the first Born approximation fails to conserve electrons even to first order. Modifying the derivation to conserve electrons introduces the correct phase without changing the scattering amplitude. We also show that the far-field expansion for the scattered waves used in many texts is inappropriate for computing an exit wave from a sample, and that the near-field expansion also give the appropriately phase-shifted result. -- Highlights: ► The first Born approximation is usually invoked as the theoretical physical basis for kinematical electron scattering theory. ► Although it predicts the correct scattering amplitude, it predicts the wrong phase; the scattered wave is missing a prefactor of i. ► We show that this arises because the standard textbook version of the first Born approximation does not conserve electrons. ► We show how this can be fixed.

Convergence of repeated quantum nondemolition measurements and wave-function collapse

International Nuclear Information System (INIS)

Bauer, Michel; Bernard, Denis

2011-01-01

Motivated by recent experiments on quantum trapped fields, we give a rigorous proof that repeated indirect quantum nondemolition (QND) measurements converge to the collapse of the wave function as predicted by the postulates of quantum mechanics for direct measurements. We also relate the rate of convergence toward the collapsed wave function to the relative entropy of each indirect measurement, a result which makes contact with information theory.

Construction of Bethe Salpeter wave functions and applications in QCD

International Nuclear Information System (INIS)

Gromes, D.

1993-01-01

We suggest an ansatz for the Bethe Salpeter wave function which is strictly covariant, obeys the spectrum conditions, and has the correct non relativistic limit. As a first simple application we present a wave function for the pion. It contains two parameters, one of them being the quark mass. The decay constant and the form factor derived from this are in excellent agreement with the data. (orig.)

Application of Coupled-Wave Wentzel-Kramers-Brillouin Approximation to Ground Penetrating Radar

Directory of Open Access Journals (Sweden)

Igor Prokopovich

2017-12-01

Full Text Available This paper deals with bistatic subsurface probing of a horizontally layered dielectric half-space by means of ultra-wideband electromagnetic waves. In particular, the main objective of this work is to present a new method for the solution of the two-dimensional back-scattering problem arising when a pulsed electromagnetic signal impinges on a non-uniform dielectric half-space; this scenario is of interest for ground penetrating radar (GPR applications. For the analytical description of the signal generated by the interaction of the emitted pulse with the environment, we developed and implemented a novel time-domain version of the coupled-wave Wentzel-Kramers-Brillouin approximation. We compared our solution with finite-difference time-domain (FDTD results, achieving a very good agreement. We then applied the proposed technique to two case studies: in particular, our method was employed for the post-processing of experimental radargrams collected on Lake Chebarkul, in Russia, and for the simulation of GPR probing of the Moon surface, to detect smooth gradients of the dielectric permittivity in lunar regolith. The main conclusions resulting from our study are that our semi-analytical method is accurate, radically accelerates calculations compared to simpler mathematical formulations with a mostly numerical nature (such as the FDTD technique, and can be effectively used to aid the interpretation of GPR data. The method is capable to correctly predict the protracted return signals originated by smooth transition layers of the subsurface dielectric medium. The accuracy and numerical efficiency of our computational approach make promising its further development.

Quantum chaos induced by nonadiabatic coupling in wave-packet dynamics

International Nuclear Information System (INIS)

Higuchi, Hisashi; Takatsuka, Kazuo

2002-01-01

The effect of nonadiabatic coupling due to breakdown of the Born-Oppenheimer approximation on chaos is investigated. A couple of measures (indicators) that detect the extent of chaos in wave-packet dynamics on coupled potential functions are devised. Using them, we show that chaos is indeed induced by a nonadiabatic coupling in individual time-dependent wave-packet dynamics. This chaos is genuinely of quantum nature, since it arises from bifurcation and merging of a wave packet at the quasicrossing region of two coupled potential functions

The deuteron bound state wave function with tensor forces

International Nuclear Information System (INIS)

Takemasa, Tadashi

1991-01-01

A FORTRAN program named DEUTERON is developed to calculate the binding energy and wave function of a deuteron, when the interaction between two nucleons is described in terms of central, tensor, spin-orbit, and quadratic LS potentials with or without a hard core. An important use of the program is to provide the deuteron wave function required in nuclear reaction calculations involving a deuteron. Also, this program may be employed in nuclear Hartree-Fock calculations using an effective nucleon-nucleon interaction with a tensor component. (author)

Extracting a shape function for a signal with intra-wave frequency modulation.

Science.gov (United States)

Hou, Thomas Y; Shi, Zuoqiang

2016-04-13

In this paper, we develop an effective and robust adaptive time-frequency analysis method for signals with intra-wave frequency modulation. To handle this kind of signals effectively, we generalize our data-driven time-frequency analysis by using a shape function to describe the intra-wave frequency modulation. The idea of using a shape function in time-frequency analysis was first proposed by Wu (Wu 2013 Appl. Comput. Harmon. Anal. 35, 181-199. (doi:10.1016/j.acha.2012.08.008)). A shape function could be any smooth 2π-periodic function. Based on this model, we propose to solve an optimization problem to extract the shape function. By exploring the fact that the shape function is a periodic function with respect to its phase function, we can identify certain low-rank structure of the signal. This low-rank structure enables us to extract the shape function from the signal. Once the shape function is obtained, the instantaneous frequency with intra-wave modulation can be recovered from the shape function. We demonstrate the robustness and efficiency of our method by applying it to several synthetic and real signals. One important observation is that this approach is very stable to noise perturbation. By using the shape function approach, we can capture the intra-wave frequency modulation very well even for noise-polluted signals. In comparison, existing methods such as empirical mode decomposition/ensemble empirical mode decomposition seem to have difficulty in capturing the intra-wave modulation when the signal is polluted by noise. © 2016 The Author(s).

Entanglement of two atoms interacting with a dissipative coherent cavity field without rotating wave approximation

International Nuclear Information System (INIS)

Kang Guo-Dong; Fang Mao-Fa; Ouyang Xi-Cheng; Deng Xiao-Juan

2010-01-01

Considering two identical two-level atoms interacting with a single-model dissipative coherent cavity field without rotating wave approximation, we explore the entanglement dynamics of the two atoms prepared in different states using concurrence. Interestingly, our results show that the entanglement between the two atoms that initially disentangled will come up to a large constant rapidly, and then keeps steady in the following time or always has its maximum when prepared in some special Bell states. The model considered in this study is a good candidate for quantum information processing especially for quantum computation as steady high-degree atomic entanglement resource obtained in dissipative cavity

Numerical analysis of different neural transfer functions used for best approximation

International Nuclear Information System (INIS)

Gougam, L.A.; Chikhi, A.; Biskri, S.; Chafa, F.

2006-01-01

It is widely recognised that the choice of transfer functions in neural networks is of en importance to their performance. In this paper, different neural transfer functions usec approximation are discussed. We begin with sigmoi'dal functions used most often by diffi authors . At a second step, we use Gaussian functions as previously suggested in refere Finally, we deal with a specified wavelet family. A comparison between the three cases < above is made exhibiting therefore the advantages of each transfer function. The approa< function improves as the dimension N of the elementary task basis increases

Validity of the Glauber approximation with account of the deuteron quark structure

International Nuclear Information System (INIS)

Zakharov, B.G.; Kopeliovich, B.Z.

1985-01-01

Corrections to the hadron-deuteron elastic scattering amplitude due to the deuteron quark structore are calculated. The two-gluon pomeron model which provides a good description of the hadron-hardron scattering data has been used. It is shown that the Glauber approximation with elastic corrections provides with a good description of the data on the total and differential pd-scattering cross sections, if the admixture of the six-quark bag state in the deuteron S 6 wave function is small ( 6 approximately 0.8 fm

Comparison of Regression Analysis and Transfer Function in Estimating the Parameters of Central Pulse Waves from Brachial Pulse Wave.

Science.gov (United States)

Chai, Rui; Xu, Li-Sheng; Yao, Yang; Hao, Li-Ling; Qi, Lin

2017-01-01

This study analyzed ascending branch slope (A_slope), dicrotic notch height (Hn), diastolic area (Ad) and systolic area (As) diastolic blood pressure (DBP), systolic blood pressure (SBP), pulse pressure (PP), subendocardial viability ratio (SEVR), waveform parameter (k), stroke volume (SV), cardiac output (CO), and peripheral resistance (RS) of central pulse wave invasively and non-invasively measured. Invasively measured parameters were compared with parameters measured from brachial pulse waves by regression model and transfer function model. Accuracy of parameters estimated by regression and transfer function model, was compared too. Findings showed that k value, central pulse wave and brachial pulse wave parameters invasively measured, correlated positively. Regression model parameters including A_slope, DBP, SEVR, and transfer function model parameters had good consistency with parameters invasively measured. They had same effect of consistency. SBP, PP, SV, and CO could be calculated through the regression model, but their accuracies were worse than that of transfer function model.

Distorted wave method in reactions with composite particles

International Nuclear Information System (INIS)

Zelenskaya, N.S.; Teplov, I.B.

1980-01-01

The work deals with the distorbed wave method with a finite radius of interaction (DWBAFR) as applied to quantitative analysis of direct nuclear reactions with composite particles (including heavy ions) considering the reaction mechanisms other than the cluster stripping mechanism, in particular the exchange processes. The accurate equations of the distorbed-wave method in the three-body problem and the general formula dor calculating differential cross-sections of arbitrary binary reactions by DWBAFR are presented. Accurate and approximate methods allowing for finite interaction radius are discussed. Two main versions of exact account of recoil effects: separation of variables in wave functions of relative motion of particles and in interaction potentials and separation of variables in distorted waves are analysed. Given is a characteristic of the known calculated programs approximately and exactly taking account of recoil effects for direct and exchange processes [ru

Approximately analytical solutions of the Manning-Rosen potential with the spin-orbit coupling term and spin symmetry

International Nuclear Information System (INIS)

Wei Gaofeng; Dong Shihai

2008-01-01

In this Letter the approximately analytical bound state solutions of the Dirac equation with the Manning-Rosen potential for arbitrary spin-orbit coupling quantum number k are carried out by taking a properly approximate expansion for the spin-orbit coupling term. In the case of exact spin symmetry, the associated two-component spinor wave functions of the Dirac equation for arbitrary spin-orbit quantum number k are presented and the corresponding bound state energy equation is derived. We study briefly two special cases; the general s-wave problem and the equal scalar and vector Manning-Rosen potential

Semi-analytical Karhunen-Loeve representation of irregular waves based on the prolate spheroidal wave functions

Science.gov (United States)

Lee, Gibbeum; Cho, Yeunwoo

2018-01-01

A new semi-analytical approach is presented to solving the matrix eigenvalue problem or the integral equation in Karhunen-Loeve (K-L) representation of random data such as irregular ocean waves. Instead of direct numerical approach to this matrix eigenvalue problem, which may suffer from the computational inaccuracy for big data, a pair of integral and differential equations are considered, which are related to the so-called prolate spheroidal wave functions (PSWF). First, the PSWF is expressed as a summation of a small number of the analytical Legendre functions. After substituting them into the PSWF differential equation, a much smaller size matrix eigenvalue problem is obtained than the direct numerical K-L matrix eigenvalue problem. By solving this with a minimal numerical effort, the PSWF and the associated eigenvalue of the PSWF differential equation are obtained. Then, the eigenvalue of the PSWF integral equation is analytically expressed by the functional values of the PSWF and the eigenvalues obtained in the PSWF differential equation. Finally, the analytically expressed PSWFs and the eigenvalues in the PWSF integral equation are used to form the kernel matrix in the K-L integral equation for the representation of exemplary wave data such as ordinary irregular waves. It is found that, with the same accuracy, the required memory size of the present method is smaller than that of the direct numerical K-L representation and the computation time of the present method is shorter than that of the semi-analytical method based on the sinusoidal functions.

Approximate Coulomb effects in the three-body scattering problem

International Nuclear Information System (INIS)

Haftel, M.I.; Zankel, H.

1981-01-01

From the momentum space Faddeev equations we derive approximate expressions which describe the Coulomb-nuclear interference in the three-body elastic scattering, rearrangement, and breakup problems and apply the formalism to p-d elastic scattering. The approximations treat the Coulomb interference as mainly a two-body effect, but we allow for the charge distribution of the deuteron in the p-d calculations. Real and imaginary parts of the Coulomb correction to the elastic scattering phase shifts are described in terms of on-shell quantities only. In the case of pure Coulomb breakup we recover the distorted-wave Born approximation result. Comparing the derived approximation with the full Faddeev p-d elastic scattering calculation, which includes the Coulomb force, we obtain good qualitative agreement in S and P waves, but disagreement in repulsive higher partial waves. The on-shell approximation investigated is found to be superior to other current approximations. The calculated differential cross sections at 10 MeV raise the question of whether there is a significant Coulomb-nuclear interference at backward angles

Quantum mean-field approximations for nuclear bound states and tunneling

International Nuclear Information System (INIS)

Negele, J.W.; Levit, S.; Paltiel, Z.; Massachusetts Inst. of Tech., Cambridge

1979-01-01

A conceptual framework has been presented in which observables are approximated in terms of a self-consistent quantum mean-field theory. Since the SPA (Stationary Phase Approximation) determines the optimal mean field to approximate a given observable, it is natural that when one changes the observable, the best mean field to describe it changes as well. Although the theory superficially appears applicable to any observable expressible in terms of an evolution operator, for example an S-matrix element, one would have to go far beyond the SPA to adequately approximate the overlap of two many-body wave functions. The most salient open problems thus concern quantitative assessment of the accuracy of the SPA, reformulation of the theory to accomodate hard cores, and selection of sensible expectation values of few-body operators to address in scattering problems

Double photoionization of helium: A new correlated double continuum wave function

Energy Technology Data Exchange (ETDEWEB)

Macri, P.A.; Kornberg, M.A.; Miraglia, J.E. [Consejo Nacional de Investigaciones Cientificas y Tecnicas, Buenos Aires (Argentina). Inst. de Astron. y Fisica del Espacio; Garibotti, C.R.; Gasaneo, G.; Colavecchia, F.D. [Centro Atomico Bariloche, Comision Nacional de Energia Atomica, 8400 S.C. de Bariloche, Rio Negro (Argentina)

1997-10-01

In this work we discuss the failures and goodness of using the product of two and three Coulomb waves to represent the double-continuum wave function of two electrons in the field of an ion. Furthermore, we present a new wave function for the double continuum, which takes into account the non-diagonal part of the kinetic energy. It satisfies the correct boundary conditions for large particle separations, and treats the electronic interaction in a more realistic way than the previously enunciated models. (orig.). 14 refs.

Wave-equation reflection traveltime inversion

KAUST Repository

Zhang, Sanzong

2011-01-01

The main difficulty with iterative waveform inversion using a gradient optimization method is that it tends to get stuck in local minima associated within the waveform misfit function. This is because the waveform misfit function is highly nonlinear with respect to changes in the velocity model. To reduce this nonlinearity, we present a reflection traveltime tomography method based on the wave equation which enjoys a more quasi-linear relationship between the model and the data. A local crosscorrelation of the windowed downgoing direct wave and the upgoing reflection wave at the image point yields the lag time that maximizes the correlation. This lag time represents the reflection traveltime residual that is back-projected into the earth model to update the velocity in the same way as wave-equation transmission traveltime inversion. No travel-time picking is needed and no high-frequency approximation is assumed. The mathematical derivation and the numerical examples are presented to partly demonstrate its efficiency and robustness. © 2011 Society of Exploration Geophysicists.

When Density Functional Approximations Meet Iron Oxides.

Science.gov (United States)

Meng, Yu; Liu, Xing-Wu; Huo, Chun-Fang; Guo, Wen-Ping; Cao, Dong-Bo; Peng, Qing; Dearden, Albert; Gonze, Xavier; Yang, Yong; Wang, Jianguo; Jiao, Haijun; Li, Yongwang; Wen, Xiao-Dong

2016-10-11

Three density functional approximations (DFAs), PBE, PBE+U, and Heyd-Scuseria-Ernzerhof screened hybrid functional (HSE), were employed to investigate the geometric, electronic, magnetic, and thermodynamic properties of four iron oxides, namely, α-FeOOH, α-Fe 2 O 3 , Fe 3 O 4 , and FeO. Comparing our calculated results with available experimental data, we found that HSE (a = 0.15) (containing 15% "screened" Hartree-Fock exchange) can provide reliable values of lattice constants, Fe magnetic moments, band gaps, and formation energies of all four iron oxides, while standard HSE (a = 0.25) seriously overestimates the band gaps and formation energies. For PBE+U, a suitable U value can give quite good results for the electronic properties of each iron oxide, but it is challenging to accurately get other properties of the four iron oxides using the same U value. Subsequently, we calculated the Gibbs free energies of transformation reactions among iron oxides using the HSE (a = 0.15) functional and plotted the equilibrium phase diagrams of the iron oxide system under various conditions, which provide reliable theoretical insight into the phase transformations of iron oxides.

Approximate Dynamic Programming: Combining Regional and Local State Following Approximations.

Science.gov (United States)

Deptula, Patryk; Rosenfeld, Joel A; Kamalapurkar, Rushikesh; Dixon, Warren E

2018-06-01

An infinite-horizon optimal regulation problem for a control-affine deterministic system is solved online using a local state following (StaF) kernel and a regional model-based reinforcement learning (R-MBRL) method to approximate the value function. Unlike traditional methods such as R-MBRL that aim to approximate the value function over a large compact set, the StaF kernel approach aims to approximate the value function in a local neighborhood of the state that travels within a compact set. In this paper, the value function is approximated using a state-dependent convex combination of the StaF-based and the R-MBRL-based approximations. As the state enters a neighborhood containing the origin, the value function transitions from being approximated by the StaF approach to the R-MBRL approach. Semiglobal uniformly ultimately bounded (SGUUB) convergence of the system states to the origin is established using a Lyapunov-based analysis. Simulation results are provided for two, three, six, and ten-state dynamical systems to demonstrate the scalability and performance of the developed method.

Wave Function Engineering in CdSe/PbS Core/Shell Quantum Dots.

Science.gov (United States)

Wieliczka, Brian M; Kaledin, Alexey L; Buhro, William E; Loomis, Richard A

2018-05-25

The synthesis of epitaxial CdSe/PbS core/shell quantum dots (QDs) is reported. The PbS shell grows in a rock salt structure on the zinc blende CdSe core, thereby creating a crystal structure mismatch through additive growth. Absorption and photoluminescence (PL) band edge features shift to lower energies with increasing shell thickness, but remain above the CdSe bulk band gap. Nevertheless, the profiles of the absorption spectra vary with shell growth, indicating that the overlap of the electron and hole wave functions is changing significantly. This leads to over an order of magnitude reduction of absorption near the band gap and a large, tunable energy shift, of up to 550 meV, between the onset of strong absorption and the band edge PL. While the bulk valence and conduction bands adopt an inverse type-I alignment, the observed spectroscopic behavior is consistent with a transition between quasi-type-I and quasi-type-II behavior depending on shell thickness. Three effective mass approximation models support this hypothesis and suggest that the large difference in effective masses between the core and shell results in hole localization in the CdSe core and a delocalization of the electron across the entire QD. These results show the tuning of wave functions and transition energies in CdSe/PbS nanoheterostructures with prospects for use in optoelectronic devices for luminescent solar concentration or multiexciton generation.

On the interpretation of wave function overlaps in quantum dots

DEFF Research Database (Denmark)

Stobbe, Søren; Hvam, Jørn Märcher; Lodahl, Peter

2011-01-01

The spontaneous emission rate of excitons strongly confined in quantum dots (QDs) is proportional to the overlap integral of electron and hole envelope wave functions. A common and intuitive interpretation of this result is that the spontaneous emission rate is proportional to the probability...... that the electron and the hole are located at the same point or region in space, i.e., they must coincide spatially to recombine. Here, we show that this interpretation is not correct even loosely speaking. By general mathematical considerations we compare the envelope wave function overlap, the exchange overlap...... integral, and the probability of electrons and holes coinciding, and find that the frequency dependence of the envelope wave function overlap integral is very different from that expected from the common interpretation. We show that these theoretical considerations lead to predictions for measurements. We...

Photon wave function formalism for analysis of Mach–Zehnder interferometer and sum-frequency generation

Energy Technology Data Exchange (ETDEWEB)

Ritboon, Atirach, E-mail: atirach.3.14@gmail.com [School of Physics and Astronomy, University of Glasgow, Glasgow, G12 8QQ (United Kingdom); Department of Physics, Faculty of Science, Prince of Songkla University, Hat Yai 90112 (Thailand); Daengngam, Chalongrat, E-mail: chalongrat.d@psu.ac.th [Department of Physics, Faculty of Science, Prince of Songkla University, Hat Yai 90112 (Thailand); Pengpan, Teparksorn, E-mail: teparksorn.p@psu.ac.th [Department of Physics, Faculty of Science, Prince of Songkla University, Hat Yai 90112 (Thailand)

2016-08-15

Biakynicki-Birula introduced a photon wave function similar to the matter wave function that satisfies the Schrödinger equation. Its second quantization form can be applied to investigate nonlinear optics at nearly full quantum level. In this paper, we applied the photon wave function formalism to analyze both linear optical processes in the well-known Mach–Zehnder interferometer and nonlinear optical processes for sum-frequency generation in dispersive and lossless medium. Results by photon wave function formalism agree with the well-established Maxwell treatments and existing experimental verifications.

Photon wave function formalism for analysis of Mach–Zehnder interferometer and sum-frequency generation

International Nuclear Information System (INIS)

Ritboon, Atirach; Daengngam, Chalongrat; Pengpan, Teparksorn

2016-01-01

Biakynicki-Birula introduced a photon wave function similar to the matter wave function that satisfies the Schrödinger equation. Its second quantization form can be applied to investigate nonlinear optics at nearly full quantum level. In this paper, we applied the photon wave function formalism to analyze both linear optical processes in the well-known Mach–Zehnder interferometer and nonlinear optical processes for sum-frequency generation in dispersive and lossless medium. Results by photon wave function formalism agree with the well-established Maxwell treatments and existing experimental verifications.

Sum rules for baryonic vertex functions and the proton wave function in QCD

International Nuclear Information System (INIS)

Lavelle, M.J.

1985-01-01

We consider light-cone sum rules for vertex functions involving baryon-meson couplings. These sum rules relate the non-perturbative, and experimentally known, coupling constants to the moments of the wave function of the proton state. Our results for these moments are consistent with those obtained from QCD sum rules for two-point functions. (orig.)

Many-body perturbation theory using the density-functional concept: beyond the GW approximation.

Science.gov (United States)

Bruneval, Fabien; Sottile, Francesco; Olevano, Valerio; Del Sole, Rodolfo; Reining, Lucia

2005-05-13

We propose an alternative formulation of many-body perturbation theory that uses the density-functional concept. Instead of the usual four-point integral equation for the polarizability, we obtain a two-point one, which leads to excellent optical absorption and energy-loss spectra. The corresponding three-point vertex function and self-energy are then simply calculated via an integration, for any level of approximation. Moreover, we show the direct impact of this formulation on the time-dependent density-functional theory. Numerical results for the band gap of bulk silicon and solid argon illustrate corrections beyond the GW approximation for the self-energy.

The exact solutions and approximate analytic solutions of the (2 + 1)-dimensional KP equation based on symmetry method.

Science.gov (United States)

Gai, Litao; Bilige, Sudao; Jie, Yingmo

2016-01-01

In this paper, we successfully obtained the exact solutions and the approximate analytic solutions of the (2 + 1)-dimensional KP equation based on the Lie symmetry, the extended tanh method and the homotopy perturbation method. In first part, we obtained the symmetries of the (2 + 1)-dimensional KP equation based on the Wu-differential characteristic set algorithm and reduced it. In the second part, we constructed the abundant exact travelling wave solutions by using the extended tanh method. These solutions are expressed by the hyperbolic functions, the trigonometric functions and the rational functions respectively. It should be noted that when the parameters are taken as special values, some solitary wave solutions are derived from the hyperbolic function solutions. Finally, we apply the homotopy perturbation method to obtain the approximate analytic solutions based on four kinds of initial conditions.

Site-occupation embedding theory using Bethe ansatz local density approximations

Science.gov (United States)

Senjean, Bruno; Nakatani, Naoki; Tsuchiizu, Masahisa; Fromager, Emmanuel

2018-06-01

Site-occupation embedding theory (SOET) is an alternative formulation of density functional theory (DFT) for model Hamiltonians where the fully interacting Hubbard problem is mapped, in principle exactly, onto an impurity-interacting (rather than a noninteracting) one. It provides a rigorous framework for combining wave-function (or Green function)-based methods with DFT. In this work, exact expressions for the per-site energy and double occupation of the uniform Hubbard model are derived in the context of SOET. As readily seen from these derivations, the so-called bath contribution to the per-site correlation energy is, in addition to the latter, the key density functional quantity to model in SOET. Various approximations based on Bethe ansatz and perturbative solutions to the Hubbard and single-impurity Anderson models are constructed and tested on a one-dimensional ring. The self-consistent calculation of the embedded impurity wave function has been performed with the density-matrix renormalization group method. It has been shown that promising results are obtained in specific regimes of correlation and density. Possible further developments have been proposed in order to provide reliable embedding functionals and potentials.

The effects of coulomb distortion on the first, second, and third sturcture functions for (e, e'p) reactions

International Nuclear Information System (INIS)

Kim, K. S.; Cheoun, Myung Ki; Cheon, Il Tong; Chung, Yeun Gun

1998-01-01

In this paper, we study the electron Coulomb distortion effects on the first, second, and third structure functions for the exclusive reaction (e, e'p) in the quasielastic region. For a heavy target ( 208 Pb) or a light nucleus ( 16 O), these structure functions calculated using the distorted wave Born approximation for the electron Coulomb distortion have shapes similar to those calculated using the plane wave Born approximation, but the effects are changed in magnitude. We use the approximate Moeller potential which has a 'plane-wave-like' form and hence permits the separation of the cross section into five structure functions. We investigate the dependence of the azimuthal angle for the outgoing proton on each structure functions. In this calculation, we use the Dirac-Hartree single particle wave functions for the ground state and the relativistic optical wave functions for the continuum proton

On Approximation of Hyper-geometric Function Values of a Special Class

Directory of Open Access Journals (Sweden)

P. L. Ivankov

2017-01-01

Full Text Available Investigations of arithmetic properties of the hyper-geometric function values make it possible to single out two trends, namely, Siegel’s method and methods based on the effective construction of a linear approximating form. There are also methods combining both approaches mentioned.  The Siegel’s method allows obtaining the most general results concerning the abovementioned problems. In many cases it was used to establish the algebraic independence of the values of corresponding functions. Although the effective methods do not allow obtaining propositions of such generality they have nevertheless some advantages. Among these advantages one can distinguish at least two: a higher precision of the quantitative results obtained by effective methods and a possibility to study the hyper-geometric functions with irrational parameters.In this paper we apply the effective construction to estimate a measure of the linear independence of the hyper-geometric function values over the imaginary quadratic field. The functions themselves were chosen by a special way so that it could be possible to demonstrate a new approach to the effective construction of a linear approximating form. This approach makes it possible also to extend the well-known effective construction methods of the linear approximating forms for poly-logarithms to the functions of more general type.To obtain the arithmetic result we had to establish a linear independence of the functions under consideration over the field of rational functions. It is apparently impossible to apply directly known theorems containing sufficient (and in some cases needful and sufficient conditions for the system of functions appearing in the theorems mentioned. For this reason, a special technique has been developed to solve this problem.The paper presents the obtained arithmetic results concerning the values of integral functions, but, with appropriate alterations, the theorems proved can be adapted to