Logarithmic pinpricks in wavefunctions

Waves in the plane, punctured by excision of a small disk with radius much smaller than the wavelength, can be modified by being forced to vanish on the boundary of the disk. Such waves exhibit a logarithmically thin ‘pinprick’, and logarithmically weak oscillations persisting far away. As the radius vanishes, these modifications become asymptotically invisible. Examples are punctured plane waves, and a punctured unit disk; in the latter case, the pinprick causes a logarithmic shift in the eigenvalues. It is conjectured that the plane can be densely covered with asymptotically invisible pinpricks, and that there are analogous phenomena in higher dimensions. The curious phenomenon of pinpricks is not hard to understand, and would be worth presenting in graduate courses on waves.


Introduction
In the anatomy of waves, zeros are important.Examples are Dirichlet boundary conditions, nodal lines in 2D and surfaces in 3D for real waves, and their counterparts in complex scalar waves, namely points and lines of phase singularity.The aim here is to explore a simple type of zero.
Consider 2D for definiteness, and consider waves ( ) y r that vanish at points in the plane x y , .Point zeros of several types have already been studied.In complex scalar waves, they occur naturally, as phase singularities (also called wave dislocations or wave vortices [1][2][3][4][5][6][7][8]). In the Šeba quantum billiard [9], they can be created by a particular boundary condition corresponding to a Dirac delta potential, whose significance is that it can introduce wave chaos in systems where the classical dynamics is not chaotic-for example for free motion in a punctured rectangle.In the Aharonov-Bohm effect [10,11], a line of magnetic flux generates a different type of zero, around which there is a phase change [12].Mathematically, the Šeba and Aharonov-Bohm zeros can be regarded as self-adjoint extensions of the quantum Hamiltonian in the punctured plane [13].In addition, plane-wave superpositions forced to vanish at points on a boundary [14] constitute an effective numerical method for computing mode wavefunctions and eigenvalues.
The zero considered here is different, and much simpler.A disk of small radius d, centred on a chosen point, is excised, and y is required to vanish on its boundary, i.e. on a small circle of radius d, and the influence of this condition as d  0 is studied.The result is that the zero differs from the ambient ('unperturbed') wave inside a logarithmically thin tube, outside which there are logarithmically weak ripples.The limit d  0 is singular, and can be regarded as a special case of the Šeba billiard boundary condition.In the limit d = 0, the zero is infinitely thin; it is asymptotically invisible, and we call it a 'logarithmic pinprick'.The limit d  0 suggests application of perturbation theory.But this is not straightforward because the limit is a singular perturbation, so it is clearer to work with explicit solutions.
Section 2 gives the general theory of the pinprick.Section 3 gives two examples: section 3.1 describes a pinprick in a plane wave; section 3.2 describes a pinprick in a circular billiard (a disk with Dirichlet boundary conditions), where the eigenvalues (energy levels) are logarithmically shifted.The concluding section 4 discusses some extensions: peppering the plane densely with pinpricks, and analogous effects in higher dimensions; also, a note about limits.
The pinprick phenomenon is not immediately obvious, and the underlying analysis is not difficult.Therefore it would be an instructive exercise, worth presenting to students at the graduate level.

General pinprick theory
A sufficiently general setting is: free waves in the plane x y , , satisfying the Helmholtz equation with wavelength l: We excise a disk of radius d l  centred on a chosen point r , 0 and require y to vanish on its boundary (figure 1): If the ambient wave, i.e. the wave before imposing the zero condition, is ( ) y r , 0 the zero condition can be most simply accommodated, sufficiently accurately, by adding a small multiple of the Bessel function of the second kind, Y : This is an exact solution of the linear equation (2.1) because it is the sum of two exact solutions.But it is an approximation to the boundary condition (2.2).On the circle ( ) 0 0 which is of order d for d l  and therefore negligible.This is further justified in the appendix, by showing that the discrepancy can be eliminated by a series of corrections, each corresponding to a much thinner pinprick and weaker accompanying oscillations.Close to r , 0 i.e. for small d k , we can replace both Bessel functions by their small-argument limiting forms (see equation (10.8.2) of [15], including the logarithmic and constant terms), so where, in terms of the Euler constant g = 0.57722, ( )

2.5
This is the asymptotic form of y close to the pinprick, and our main result, illustrated in figure 1. Far from r , 0 i.e. away from the pinprick, the Bessel function 3) oscillates, indicating persistent undulations decorating the ambient wave.In the next section we will see that for small d this influence is weak; we also describe a slight generalisation.

Plane wave
For an ambient plane wave in the x direction, with wavelength l = 1, punctured at the origin, i.e. Figure 2 shows the intensity   and slower and stronger in y than x.By contrast, (e), (f) show that the effect on the phase is comparatively weaker, illustrating the fact that pinpricks are not phase singularities (around which the phase would change by a nonzero multiple of p 2 ).Close to the zero, the limiting form of the intensity is (see 2.4) From this it follows that the intensity halfwidth of the pinprick has radius Far from the pinprick, large-argument Bessel asymptotics (section 10.17(i) of [15]) gives the asymptotic intensity oscillations along the x axis:  These agree to visual accuracy with the exact oscillations in figures 2(a) and (b) (comparison not shown).
For simplicity, an important aspect of the oscillations was not mentioned until now.In the pinprick formula (2.3), the replacement leaves the small d behaviour close to the pinprick unchanged, but alters the oscillations.For the particular choice = a 1, the factor is involving the Bessel function of the third kind.This choice corresponds to the oscillations representing outgoing radiation.In the appendix we will return to this generalisation.

Circular billiard
Consider azimuthally symmetric waves ( ) y r in the unit disk with Dirichlet boundary conditions, d-punctured at the centre = r 0. The ambient modes, nonvanishing at = r 0 and vanishing at = r 1, are the zero-order Bessel functions with zeros j ,    The shift in the eigenvalues, n can be easily calculated numerically, or perturbatively from (3.12): Figure 5 shows the fractional eigenvalue shifts as functions of d.
The formula (3.12) can be simplified using the approximation This is the asymptotics for  n 1, but it can be used even for = n 1 where the ratio is 0.9822; similarly, we can approximate (accurate to 2% even for = n 1).Thus the fractional shift of the d -punctured (pinpricked) mode n is The factor /n 1 shows that the eigenvalue shift for fixed small d gets smaller as n increases, as illustrated in figure 5 for = n 1 and = n 3.

Concluding remarks
The foregoing has considered a single pinprick, but it is clear that the plane can be d punctured at more than one place.Indeed, a dense peppering of the plane with zeros can be envisaged, which would be invisible in the limit d  0. How dense?I conjecture that in the limit as d  0 the pinpricks can form a set of fractal dimension < d 1; dimension = d 1 is excluded because this would include Dirichlet conditions along a line, which would certainly not be invisible.
For waves in three-dimensional space, asymptotically invisible pinpricks would form lines, arbitrarily curved, possibly closed, linked or knotted, and arbitrarily numerous.In addition, there can be point zeros, based on excising small spheres of radius d.O as in (3.5).Extension to higher space dimensions D is clear: asymptotically invisible pinpricks can be created, with codimensions =  C D 2, , .
Complementary to pinpricks would be positive spikes ('prickles') instead of zeros.To create these is simple: instead of the negative modification (2.3), add 1, to the ambient wave.As e decreases to zero, the prickles get asymptotically invisible.
A note about limits.The emphasis here is the pinprick limit d  0. Other limits need not commute with it.For example, if the d-excised disk is in a region of the plane with a rectangular perimeter, or if the plane is compactified into a 2-torus, the system is the Sinai billiard [16], whose geometrical ('classical') trajectories exhibit chaos: exponential separation of infinitesimally close initial trajectories as  ¥ t time .This obviously does not commute with d  0, because as d gets smaller the chaos gets longer to develop, and when d = 0 it never does.And in the wave version, the eigenvalues for finite d reflect the classical chaos by exhibiting random-matrix statistics [17][18][19].The d  0 limit also does not commute with the classical limit   0, which in this case corresponds to the short-wave limit l  0 of high excited states: as d gets smaller, it is necessary to search higher in the spectrum in order to see the random-matrix statistics (this is discussed in appendix C of [17]).
The boundary condition at d = r can be satisfied exactly for each n separately, by adding a contribution to ( ) J r . n The modification depends on other boundary conditions that the exact ( ) y r must satisfy.For a real wave, as in section 3.2, the modification is For the infinite plane, the modification can be any combination of incoming and outgoing waves, depending on the radiation condition.For the outgoing choice, it involves Bessel functions of the third kind, as described for = n 0 at the end of section 3.For d  1, the difference between the coefficients in (A.2) and (A.3) disappears; the small argument limiting forms of the Bessel functions [15] gives

O
Therefore the ¹ n 0 corrections correspond to thinner pinpricks, reflecting the familiar s wave scattering of scalar waves by infinitesimal obstacles.
For = n 0, the real and outgoing modifications are different for the  r 1 oscillations.In the example in section 3.1 we have chosen the real modification in equation (3.2), but the order of magnitude of the oscillations is the same for the outgoing case: both vanish as (( ( )) ) d -- log .

Figure 1 .
Figure 1.Logarithmic form of the pinprick close to the zero at | | -d = r r 0 .

3
values of d.The thinness of the pinprick is evident in (a), (b).
y axis:

(
in this case (2.3) is exact.)The eigenvalues k n 2 of the operator in (2.1), with the stated Dirichlet boundary condition at = r 1, are determined by

Figures 3 and 4 show the punctured modes for = n 1 and = n 3 . 6 6 Figure 4 .
Figures3 and 4show the punctured modes for = n 1 and = n 3. The pinpricks are clearly visible in figures 3(d) and 4(d), even for d = - 10 .6Away from the pinpricks, the modes are visibly altered for / d = 1 20, and hardly so for d = - 10 .6