Derivation and Numerical analysis of an Attenuation Operator for non-relativistic waves.

Quantum mechanical models for particles are strictly dependent on the Schrödinger equation, where the solutions and the Hermitian polynomials form a mathematical foundation to derive expectation values for observables. As for all quantum systems, the solutions are derived in discrete energy levels, and yield probability density, the kinetic energy and average momentum. In this study however, an attenuation Hamiltonian is derived by the algebraic relation of the momentum and position operators, and the derived equation, where the attenuation of kinetic energy is the eigenvalue, is studied numerically. The numerical solutions suggest that the change in kinetic energy from one transition to the next proceeds in an undular fashion, and not in a definite manner. This suggests that any sub-atomic particle which experiences a transition from one level to the next, does so by both gaining and losing energy in an undular manner before reaching an equilibrium with a new and stabilized kinetic energy. The results show also that the phase of the change in kinetic energy between transitions differs between high and low momenta and that higher levels of momentum attenuate more smoothly than transitions between lower energy levels. The investigated attenuation operator may be important for future pinning and quasipinning approaches and play a role in future quantum information processing. Future research is required on the spectrum of the operator and on its potential analytical solutions.

Introduction

Wavefunction studies are important for the development of various technologies in science, particularly quantum computers, quantum holography systems and laser technologies. Particularly, assigning the values and properties of wavefunctions to algorithms and quantum computational structures is of paramount importance, and several groups have developed various models in the last years[1-3]. Weinacht et al. reported a laser-based system where the ionization signals of excited atoms was assigned to specific wavefunction integrals and thus formed the rational structure on the wavefunction-based computational quantum algorithm[1]. Sau et al. (2010) recently reported a two-dimensional model and the wavefunction solutions for the v = 5/2 fractional quantum Hall system under a high magnetic field as a topological qubit for quantum computational algorithms[2]. Paskaukas and You (2001) showed that the probability density of the wavefunction is a good property to discern between identical particles[3] and hence can be used to construct wavefunction-sensitive algorithms for quantum computation. These various properties are ultimately applied in the quantum computer which is engineered to form a system of coherent adaptation of the quantum mechanical wavefunction to execute a task[4,5]. For a computational model to make a practical application in a quantum computing system, it must be scalable and hence its solutions (wavefunction, spins, etc.) which are individual qbits, must be represented in a Hilbert space composed of infinitely many dimensions[4]. The wavefunctions can describe electrons, anyons, photons, ions or atoms[5-9] and can be pinned[10-12]. In order to discern between wavefunctions and particles, operators are paramount. Several groups have developed and studied new quantum operators in recent years, such as non-Hermitian Hamiltonians with special symmetric properties[13], and complex-scaled Hamiltonians[14]. In this study however, an attenuation Hamiltonian for quantum waves is investigated, where the behaviour of waves between transitions is mapped. Transitions between energy levels can for instance be caused by the release or absorption of a photon, and are given by the formwhich describes the wavefunction as a mixture of two states, m and n. The wave-particle entity which is subjected to the transition between the energy levels, m and n, experiences a change in kinetic energy, which can be described by the form:

The one-dimensional Schrödinger equation forms the foundation for the derivation of the expectation values of the position, momentum and kinetic energy, which are derive by deploying the algebraic norm on the operators , and , which are all included as levels of integration in the Hamiltonian kinetic operator:from the Schrödinger equation Hψ = Eψ.

The Schrödinger Hamiltonian:derives from the general wave equation, which in turn originates from the algebraic properties of the harmonic representation, where the representation is constructed in terms of the two operators P and Q, which respectively define the momentum and the position. The representation of the Schrödinger model, also known as the oscillator representation with an attractive potential is the norm of the operators defined in the association:which is diagonalized by the Hermite functions on [15]. Other combinations of the norms are given[15], however such representation describe repulsive potentials (i.e. −P2).

In order to develop a model for the attenuation in kinetic energy, P3, during transitions between states m and n, we consider electrons entirely as waves which are affected by a nuclear potential-well pull, defined by the classical potential term Q2. This depicts a quantum-relativistic model of transition dynamics for quantum waves, which is proposed to be described by the representation:hence requiring the non-self-adjoint third order differential operator:which is associated to the fluctuations in energy and the harmonic oscillator well-potential in the following relation:where the absolute value of the change in kinetic energy, EΔ is the energy difference between energy levels m and n and is the harmonic force-constant. The Hamiltonian for the attenuation of a wave under a well potential in a box is therefore proposed to:where the notation H stands for s ingle wave a ttenuation Hamiltonian.

The operator in (9) is remarkably simple and grants a simple Hamiltonian which describes attenuation in quantum physics for systems confined by a well-potential. Eqn. (9) can be expanded with additional terms such as electromagnetic field vector potential and be extended to describe wave behaviour and wave changes under various physical phenomena and hence an additional model of study in addition to statistical physics for wave properties. With this in mind, an algebraic outline of the operator properties is made, and the properties of the solutions are ultimately derived numerically in order to study this Hamiltonian as a potential candidate for pinning or quasi-pinning methods for quantum algorithms and quantum measurements of waves and signals.

Operator Properties

In order to conceive the properties of the modified Schrödinger equation, we wish to study the attributes and the allowed operations of the inherent Hamiltonian (9) in Hilbert space . This is done by the following set of propositions and proofs.

PROPOSITION I.

H is linear.

Proof. For H to be linear, the following relation must be satisfied:

Inserting for H in the left-hand side of the equation, and considering the operator under the Heisenberg picture, with and the potential coefficient set to unity:

By the sum rule:it follows that:giving

For the right-hand side of the equation, we have:which yields:and

(14) is equal to (17), hence, H is linear.

PROPOSITION II.

H is symmetric for a set of orthogonal functions in a subdomain in : .

Proof.

For H to be symmetric in , the given relation follows:With H given explicitly and using the Heisenberg picture, with :where the left side of the equal sign is equal to:We consider the two terms on the right hand side as and :

Solving for A and its two Hermitian counterparts with integration by parts applied with and ψ* = v:and the equivalent for the right hand side in (19):which yields (19) fully written out:which is valid only for orthogonal functions . We therefore prove:hence, A is symmetric.

Solving for B in (23):which is also valid only for orthogonal functions . We therefore provehence B is symmetric. With (A + B) ⊆ H, H is symmetric.

PROPOSITION III.

H is an unbounded operator in . H is continuous in the compact interval , where , and .

Proof. A proof for unbounded operators is generally related to the unmet criteria for boundedness, where, for H = (A + B) to be bounded, the following relation must be satisfied:

Applying (A) on the given function, ψ which is square integrable in L2[0,2π], infinitely differentiable and is thus element of :gives with . Using the Heisenberg picture, with :and,hence, does not satisfy condition (25) and (26) and is therefore unbounded.

Applying (28) on (29), the same accounts for B:hence B is unbounded. H = A + B, whence H is an unbounded operator.

PROPOSITION IV.

H is variant under time-reversal and spatial reflections.

Proof.

For H to variant, then the following condition applies:withandwe use A and B as above, where , , B = x2 and . The potential term is treated as a Hermitian function, and by following the Hermitian relationship: f(x)* = f(−x), the minus term appears in (36). We use as before, the Heisenberg picture with integration by parts as in (21):andand for B:With 〈ψ, Aψ〉 ≠ 〈A*ψ, ψ〉 and 〈ψ, Bψ〉 ≠ 〈B*ψ, ψ〉, we have that

H is variant under spatial and temporal transformations. As we consider a mono-dimensional attenuation model in (9), the variance is not expected to change the outcome as the absolute value of the change in energy is used (given in (8)).

PROPOSITION V.

H satisfies the Born-Sommerfeld condition in where where {ψ1, ψ2, ..., ψ} ∈ Y and .

Proof.

For an operator to satisfy the Born-Sommerfeld condition, it must be Hermitian and symmetric, which is proved in proposition II. H is therefore Hermitian and probability is conserved. By this, the domain of orthogonal functions is maximal, and hence no boundary conditions have been neglected for the physical problem described by (8).

PROPOSITION VI.

H is a positive operator.

Proof Any positive operator has a positive inner-product, 〈Hψ, ψ*〉. Using the infinetly differentiable and square integrable trial function , we get:hence H is a positive operator yielding positive expectation values, when acting on a Hermitian pair of infinitely differentiable functions, similar to the standard Schrödinger Hamiltonian.

Spectral properties of H

Spectrum of H

For H to have a real spectrum, it must be closed in its domain. We have seen from the previous section that H is an unbounded operator, and hence it’s domain cannot be defined everywhere in . However, H can be defined in a subdomain domain which is composed of two sets of elements, Y1 and Y2, where , hence and .

For H to be closed, the following must be satisfied:andwhich is met only for the case of Hermitian sets Y1 and Y2. H is therefore closed in a subset , and thus H has a real and non-negative spectrum (see proof VI above) within the sub-domain . This set of orthogonal functions forms a closed domain , where H is compact. This domain is is the cone in space of principal eigenvalues of the operator and form the spectrum of real averages derived from the solutions of equation (8), via the inner product defined in the closed domain of orthogonal functions (see proof II and V).

If H is not closed on , then its set of eigenvalues σ(H) is complex. Inversely, as H is closed on , its set of eigenvalues is real and the derivation of averages yields real expectation values. H has therefore a complex point spectrum σp(Hsa) outside the sub-domain and a real point spectrum within the sub-domain . As stated above, the sub-domain on is composed of orthogonal eigenfunctions. In other words, H is an operator that can be defined in various quantum mechanical integrals (i.e. inner product for expectation values or the quantum Lagrangian for exchange statistics) for orthogonal solutions.

His a spectral operator

H is, as most unbounded operators, a spectral operator, a definition which attributes to it special properties, with emphasis on the following:which implies that its continuity has a limit from the left. This gives the particular advantage in that it can be represented as a Riemann-Stieltjes integral, or it can be written by Stone’s theorem[16] in the form (following the Heisenberg picture):which makes the time-dimension, which is strongly continuous, adaptive to the operator H. This allows to solve the equation to obtain the state at any subsequent time. This forms a foundation for a time-dependent/-independent mode of analysis of eigenfunctions of H.

Finally, as Eq. (45) is assured, we can further outline by Stone’s theorem that every one-parameter group of unitary transformations admits the spectral representationwhere {E} is a spectral family such that . Hence the solutions of the Hamiltonian H invoked in the modified Schrödinger equation in Eq. (8), are composed of terms with complex exponentials with generalized trigonometric polynomials.

Relationship of H with other quantum operators

The operator H is considered for commutation with the Schrödinger operator in the harmonic potential picture. To begin with, we can inspect the Schrödinger Hamiltonian and the modified Hamiltonian for commutation:which gives

The same outcome is derived for the full Schrödinger operator in the harmonic potential picture with the operator in Eq. (9). By Eq. (49) it can be derived that the operator in Eq. (7) commutes with the kinetic energy and the momentum of a particle. Hence, the attenuation of a kinetic wave, can be defined simultaneously as the momentum and the kinetic energy of a quantum wave.

Results

The initial analysis of the modified Schrödinger equation in Eq. (8) is performed in the Heisenberg picture, with , and with the value of |EΔ| set for simplicity to a constant, 1, and the coefficient of the harmonic oscillator potential is equally considered at unity. The differential equation studied numerically is therefore,where ψ = ψ1(x) + iψ2(x) contains ψ1(x) and ψ2(x) respectively the imaginary and real part of ψ. By (50) we obtain the equationand thus the system of differential equationswhich are studied numerically using AsymptoticDSolveValue and NDSolve in Mathematica under specific initial conditions (I.C.). The first set of I.C. yields a wavefunction that starts a unity at the zero-point, which is an initial condition conventionally used for harmonic-oscillator solutions to the Schrödinger equation[17]. The derivative is set to zero, in order to indicate a zeroth-order momentum at the start, similar to s-orbitals. The last term, y″(0) is considered over several levels, given that the behaviour of the particle’s kinetic energy in a transition between two states m and n is not measurable by conventional quantum mechanical operators[7], and hence we represent it by step-wise incremental kinetic energy values. We consider therefore the simple initial conditions for evaluating the pattern of the change in kinetic energy of the wavefunction with the following I.C.s:

Approximated series solutions depending on increment in kinetic energy

The Mathematica code AsymptoticDSolveValue was used to derive approximated solutions in series, for the problem given in (52). The solutions represent the functions that describe the energy transitions between energy levels following initial conditions in (53) with a fixed momentum value of zero, and transitions from a kinetic energy factor of 1, 10, 15 and so forth. These functions are approximated polynomials and have immobile local maxima and minima independent on the initial conditions reflecting kinetic energy change (higher levels of n). The real part of the polynomial solution is given bywhere k is the level of kinetic energy given in (53), hence k can be 1, 10, 15, …, k. All levels have the same stationary points in the wavefunction, respectively −5.88 ∧ 2.04. The imaginary part is equally given as:with stationary points at 2.60 ∧ 4.67.

The stationary point in Fig. 1 indicate turning points in the plots, where the change in kinetic energy approaches zero. Hence, by following the series solutions given above, one undulation appears to occur in the process of attenuation of kinetic energy for a wave during a transition from level m to n. This suggests that attenuation in kinetic energy does not evolve in a constant fashion during a change from one level to the next, but in an undular manner.

Figure 1

The approximated series solution of the modified Schroedinger equation (52) for the first I.C. given in (53).

Numerical solutions depending on increment in kinetic energy

The numerical solutions derived using NDSolve depended equally to the series solutions, by the same set of initial conditions in (53) for the problem given in (52). Increments in kinetic energy gave no change in the pattern of the wavefunction (Fig. 2), similar to the series approximations with stationary equilibrium points. This indicates that the process of transition between levels occurs in an equal fashion for any defined interval of potential energy ΔE.

Figure 2

The numerical solution of the modified Schroedinger equation (52) for the initial conditions from Eq. (53). The same plot resulted for all I.C.

Numerical solutions depending on energy level transition

As a third part of the numerical analysis, it is equally critical to determine the effect of the width of the excitation gap between energy levels on the solutions of Eq. (52). The lowest energy transitions have larger gaps than the transitions which reach the classical limits, hence the behaviour of the solutions to (52) where the energy gap is defined by variable value ΔE:is investigated in Mathematica[18] and shown in Fig. 3, where ΔE is modulated from a low value (narrow excitation gap - classical limit) to a high value (wide excitation gap - quantum level).

Figure 3

The numerical solutions of the modified Schroedinger equation (52) depending on variations of ΔE from A (large value of ΔE - quantum level) towards the classical limit, D (small value of ΔE).

Approximated series solutions depending on increment in momentum

In order to study the pattern of energy transition for different start-momenta (i.e s-, p-, d- and f- and g- orbitals), the initial conditions:are used to form series solutions to the problem in Eq. (52) by differences in momentum of the waves. The real components of the series solutions following increment in momentum (subscript) are:

The imaginary components of the solutions are given as:

The real components have the following nodes with the x-axis, shown in Table 1 (where the attenuation of energy is zero - stable states).

Table 1

Stationary points of series solutions (zero attenuation).

Momentum	ψ(0)1	ψ(0)2	ψ(0)3
1	−0.92	2.45	5.06
2	−14.86	2.41	4.79
3	2.36	5.68	23.20
4	−33.19	2.37	4.61
5	−43.27	−3.49	2.19

In an equal fashion, the imaginary components have variable stationary points (data not shown). The variation of stationary points suggests that the attenuation proceeds by different undulation patterns for waves of different momentum. The stationary points are positions where the attenuation is zero, given that the wavefunction is rectilinear. A zero attenuation indicates an intermediate stationary state in change of energy. The plots of the 5 states of Eqs 58 and 59 are shown in Fig. 4.

Figure 4

The approximated series solutions of the modified Schroedinger equation (52) by increment in momentum (A-E refers to momenta in Table 1) in the initial conditions.

Numerical solutions depending on increment in momentum

The values of the momentum in the initial conditions are, equally to the previous section, studied numerically. The plots show a trend where an increment in momentum in the initial conditions yields a real part of the solutions that evolves from having an undular behaviour within a well-like parabolic curve towards forming a hyperbolic curve, similar to the hyperbolic trigonometric function sinh(x) (Fig. 5), while the inverse occurs for the imaginary part. Accounting for that the imaginary part is related to the differential operator, and the real part is closely associated to the effect from the potential term x2, it suggests that zero momentum (i.e. s-orbitals) allows the wave to attenuate energy in a dynamic undular fashion (with several stationary points in attenuation), when absorbing energy and undergoing a transition to the next level of momentum (i.e. s- to p-orbital). Higher transitions (i.e. from p- to d-, or d- to f-) yield instead a smoother attenuation of energy from one level to the next.

Figure 5

The numerical solutions of the modified Schroedinger equation (52) by increment in momentum (A-E refers to momenta in Table 1) in the initial conditions.

Discussion

The first point to consider for discussion is the resulting existence of stationary points in the attenuation pattern of a wave, when undergoing a transition between energy levels. The first series-wave functions in (54) which can be attributed to describe the attenuation of energy between any level of kinetic energy with fixed momentum (n = 1 → n = 2, n = 1 ← n = 2, n = 40 → n = 41 or n = 40 ← n = 41) show the particular property of having two preserved stationary points occurring for any transition. Noting that the transitions only differ by start-level of kinetic energy, the existence of two stationary points where the attenuation of energy is zero, defines that any transition proceeds in an undular manner, independent of the state. This suggests that the absorption or the release of a photon does not proceed linearly, but harmonically. The implications of this feature are that “quasi-states” can occur between energy levels, and that energy is exchanged bidirectionally, rather than uni-directionally. A second point of interest for discussion from the results are the implications of the momentum of the wave when discerning between transitions. The series solutions in Eq. (58) and also the numerical solutions (Figs 2 and 4), show that attenuation of energy for wave with differing momentum proceed in different patterns. The first point of interest lies in the variation of the distribution of the stationary points in (58). These stationary points, where the attenuation of the wave is zero, appear in narrower intervals for sharp orbitals, slightly wider interval for principal orbitals and wide intervals for diffuse and fundamental orbitals. In other words, attenuation of energy for s-orbitals (i.e. 1 s to 2 s) can be pinned to proceed with a high “frequency” in attenuation of energy. If the orbitals are more evolved in angular momentum (i.e. d orbitals), the attenuation proceeds more smoothly, with wider infinitesimal time-periods between stationary phases of attenuation. A third point of discussion is the property of the third order differential operator, which is a non-self-adjoint operator. Non-self-adjoint operators in quantum mechanics are less common, however some researchers have studied non-self-adjoint Hamiltonians in context with particular physical problems. Hokkyo[19] studied the role of adjoint wavefunctions for non-self-adjoint Hamiltonians in nuclear reactor models. Hokkyo used the non-self-adjoint representation to study the solutions for unstable states which can be exemplified by non-relativistic mechanisms such as the Gamow theory of α-decay and the unstable V particle in the NVfJ model[19]. Davies[20] indicate that the resolvent norms of non-self-adjoint quantum operators with complex potentials becomes very large even though the complex variable is far from the spectrum of the operator. This may imply the existence of non-trivial or simply incorrect eigenvalues the more . In this study we have considered the Eq. (8) under the Heisenberg picture, where for computational limitations, however the operator yields real eigenvalues (see section on Spectrum properties) and should therefore not suffer for complex eigenvalues. However, the operator is not self-adjoint, and has no upper or lower limit, contrary to the regular Schrödinger operator which has a lower limit. This may however not suffer for physical inconsistencies, and rather represent a challenge in determining its spectrum. At last, one wishes to ask the following question: How do these results serve the scientific community? There are two answers to this question. The first is the educational perspective and scientific value of conceiving that all processes of energy-exchange occur most-likely by wave-like patterns. The second answer is the value to technology development. Although infinitesimal time-periods reign the attenuation process of waves, the differentiation of attenuation between waves of unequal momentum may be a relevant property of quantum waves for pinning and quasi-pinning of states in the development.

Conclusions

A new wave model has been derived for waves affected by a well potential by the algebraic relations of the momentum and position operators. The attenuation Hamiltonian includes the imaginary third order differential operator, which - in similarity to the second order differential operator from the regular Schrödinger equation - has algebraic properties that allow for orthogonal functions to be derived analytically. The attenuation Hamiltonian yields a wave-attenuation model with non-quantized solutions which describe the change of the kinetic energy of waves subjected to a well potential, during transitions from one energy level to the next. This Hamiltonian can describe the release and absorption of photons and be an important operator for pinning and quasi-pinning algorithms.

Methods

The numerical analysis and derivation of approximated series solutions is conducted using NDSolve and AsymptoticDSolveValue with the Mathematica Development Platform[18]. The  Mathematica protocol is available at Fjordforsk A.S. website, http://www.fjordforsk.no/sci_rep_1860.zip.