Next Article in Journal
Our Peculiar Motion Inferred from Number Counts of Mid Infra Red AGNs and the Discordance Seen with the Cosmological Principle
Next Article in Special Issue
No-Slip Boundary Conditions for Electron Hydrodynamics and the Thermal Casimir Pressure
Previous Article in Journal
Painlevé–Gullstrand form of the Lense–Thirring Spacetime
Previous Article in Special Issue
Measurement of the Casimir Force between 0.2 and 8 μm: Experimental Procedures and Comparison with Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Van der Waals Interactions of Moving Particles with Surfaces of Cylindrical Geometry

Nanoscale Physics Group, Kabardino-Balkarian State University, 360004 Nalchik, Russia
Universe 2021, 7(4), 106; https://doi.org/10.3390/universe7040106
Submission received: 17 March 2021 / Revised: 6 April 2021 / Accepted: 13 April 2021 / Published: 19 April 2021
(This article belongs to the Special Issue The Casimir Effect: From a Laboratory Table to the Universe)

Abstract

:
General nonrelativistic theory has been developed and the expressions obtained for the tangential (dissipative) and radial (conservative) image forces and van der Waals forces (vdW) acting on charged and neutral particles when they move parallel to the axis of a cylinder with circular cross-section, or in the space between coaxial cylinders. Numerical calculations of vdW forces have been performed for metal (Au) and dielectric (Si) materials of cylinders (filaments) and Cs atoms at velocities ~ 10 7 m / s . A remarkable result is that in the case of metal cylinders (atomic filaments and chains) the dynamic vdW potential can be repulsive for certain values of the velocity–distance parameter and the characteristic atomic frequency. In the case of a Si material, the dynamic vdW potential increases relative to the static one (by modulus) when the velocity–distance parameter V ω 0 / R changes from zero to ~ 1.3 and then tends to zero.

1. Introduction

Atom–surface interaction is a long-explored research problem of physics and adjacent areas. A particular case is the van der Waals interaction (vdW) caused by the zero-point and thermal fluctuations of the electromagnetic field and its material sources. The surface curvature is known to affect the spectrum of surface excitations (plasmons) and, accordingly, vdW forces of interaction of particles with the surfaces at distances of 1 to 10 nm [1,2,3,4]. Calculations of the vdW and Casimir–Polder potentials in cylindrical configurations (even in static case) invoke great interest [5,6] due to advances in nanotechnology, since it became possible to measure the Casimir and Casimir–Polder forces with increased precision (see Ref. [7] for a review). However, currently there are practically no relevant studies of nonequilibrium vdW interactions in cylindrical geometry. Dynamic corrections to conservative vdW forces in this case were first considered in [8]. More general calculations of the conservative and dissipative vdW forces acting on particles moving parallel to the generatrix of a cylindrical surface were carried out in [9,10]. In addition to the general theoretical importance associated with nonequilibrium vdW forces in the systems with cylindrical symmetry, studying the electromagnetic and fluctuation-electromagnetic interactions of moving charged and neutral atomic particles with cylindrical surfaces has been given a great impetus since the discovery of carbon nanotubes in 1991 [11], when a lot of theoretical teams started to work on the theoretical modelling and computer simulation of ion channeling through carbon nanotubes and capillary structures [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]. However, the role of the image forces for charged particles and vdW forces for neutral ones is underestimated due to an insufficient knowledge of the role of surface curvature at distances exceeding atomic dimensions. An interesting, but still not fully explored, feature of the image forces [16] and vdW forces [28,29,30,31] of moving particles is their velocity dependence. In particular, a particle moving near a surface, along with an attractive velocity-dependent force (an image force or vdW force in the case of charged and neutral particles, respectively), experiences a tangential dissipative (stopping) force (friction force). Both of these forces affect the particle dynamics. In contrast to the case of plane geometry [28,29,30,31], in systems with cylindrical symmetry, the velocity and curvature effects are much more complex [8,9,10,13]. Meanwhile, in the standard description of the dynamics (channeling) of charged and neutral particles in nanotubes and the capillary systems, only pairwise short-range interaction potentials between the particles and the atoms of walls are used, or continuous potentials obtained by averaging the pair potentials [13,14,15,16,17,18,19,20,21]. At channel diameters of the order of several nanometers and more, the contribution of short-range pair interactions in the main part of the cross section of the channels is close to zero, and the main role will be played by the image forces for charged particles and vdW forces (for neutral ones) with the walls of the channels. Due to this, the development of an adequate theoretical basis required in future profound studies of the channeling of charged (neutral) particles through nanotube and capillary structures is an urgent problem.
In this article, the particles are assumed to move nonrelativistically parallel to the generatrix of a concave/convex cylindrical surface, or in the space between coaxial cylindrical shells. The properties of surface materials are described by a local dielectric permittivity. Within the framework of electromagnetic fluctuation theory, we calculate the tangential (dissipative) and radial (conservative) image and vdW forces acting on particles, and the rate of their heating (for nanoparticles) when (in the general case) the particles and the cylinder walls have different local temperatures. In the limit of an infinitely large radius of cylinders, the results are reduced to those in the case of a flat surface or a plane-parallel gap. Numerical calculations of the radial and tangential vdW forces are carried out in the case of the metal (Au) and dielectric (Si) materials of the cylinders (atomic filaments and chains) and neutral Cs atoms at velocities of ~107 m/s. The role of pair interactions and nonlocal effects is briefly touched upon. An intriguing result following from numerical calculations in Section 6 is that dynamic vdW potentials can be repulsive at high velocities, for certain combinations of resonance wavelengths and distances.

2. General Theory

Following [9,31], we recall the main details of the calculation of the vdW force (both its conservative and dissipative parts) and the heating rate for a neutral particle moving parallel to the generatrix of a cylindrical surface (convex and concave).
The case of a charged particle is also treated as a constituent point of the theory. Though the velocity is nonrelativistic ( V c ) and the retardation effects can be neglected, it can reach the values of order ~ 10   V B (with V B being the Bohr velocity). Figure 1 shows the cylindrical coordinate systems used and different configurations to be analyzed: convex cylindrical surface (a), cylindrical channel (b), and the case of particle motion between two coaxial cylindrical shells (c). The wall surface material is characterized by a local dielectric constant ε(ω), and the particle is characterized by the dipole polarizability α(ω). We use a cylindrical coordinate system with coordinates ( r , φ , z ) , where the z axis is parallel to the symmetry axis of the cylinders (cylindrical shells) and their generators. The radial distance r is counted off the axes of the cylinders. We further assume that an external charged or neutral particle with density ρ ( r , φ , z , t ) moves along a classical trajectory with constant velocity V localized in the vacuum region outside (Figure 1a), inside (Figure 1b) or between the cylinders (Figure 1c). The instantaneous particle coordinates are ( R , 0 , z = V t ) . Within the framework of the nonrelativistic approximation V << c, the main task is to find the electric field E on a moving particle at the point of its instantaneous location. Since = ϕ , with ϕ being the electric potential, the initial problem reduces to solving the Poisson equation Δ Φ = 4 π ρ   .

2.1. Green’s Functions of Poisson’s Equation

Written in cylindrical coordinates , φ , z , the Poisson equation takes the form
( 2 r 2 + 1 r r + 1 r 2 2 φ 2 + 2 z 2 ) Φ ( r , φ , z , t ) = 4 π ρ ( r , φ , z , t )
According to the symmetry of the problem, we expand both parts of (1) into the Fourier integral over the frequency ω and the wave vector k (the direction of k is parallel to the z axis) and into a Fourier series over the angle φ :  
Φ ( r , φ , z , t ) = 1 2 π n = n = + + d ω 2 π + d k 2 π u n ( ω k ; r ) exp ( i ( k z ω t + n φ ) )
ρ ( r , φ , z , t ) = 1 2 π n = n = + + d ω 2 π + d k 2 π ρ n ( ω k ; r ) exp ( i ( k z ω t + n φ ) )
Substituting (2), (3) into (1) yields
( d 2 d r 2 + 1 r d d r n 2 r 2 k 2 ) u n ( ω k , r ) = 4 π ρ n ( ω k ; r )
Without the loss of generality, let us consider the case shown in Figure 1a. Then, the solution to (4) should meet the boundary conditions at the surface r = a of the cylinder:
u n ( ω k ; r = a + 0 ) = u n ( ω k ; r = a 0 ) ,   d u n d r r = a + 0 = ε ( ω ) d u n d r r = a 0 .
It is worth noting that Equation (4) is valid both inside and outside the cylinder, since the problem is electrostatic.
To proceed further, we first find the Green’s function G n ( r , r ) of Poisson’s Equation (4)
( d 2 d r 2 + 1 r d d r n 2 r 2 k 2 ) G n ( r , r ) = δ ( r r )
with additional conditions
G n ( r = r + 0 , r ) = G n ( r = r 0 , r ) ,   d G n d r r = r + 0 d G n d r r = r 0 = 1 ,
and Equation (5). Then, the general solution to inhomogeneous Equation (4) that meets boundary conditions (5) and (7) reads
u n ( ω k ; r ) = 4 π   d r G n ( r , r ) ρ n ( ω k ; r ) .
For r r , Equation (6) is a modified Bessel equation with two linearly independent solutions I n ( k r ) and K n ( k r ) . We find G n ( r , r ) in Appendix A. Equations (A2), (A4), (A6) and (A7) correspond to cases (a), (b), (c) shown in Figure 1.

2.2. Image Forces on a Moving Charged Particle

The charge density ρ of a moving particle with a charge Z e has the form [32]
ρ ( r , φ , z , t ) = Z e r δ ( r R ) δ ( φ ) δ ( z V t ) .
Using (9), the Fourier-transformed density ρ n ( ω k ; r ) is given by
ρ n ( ω k ; r ) = 2 π Z e r δ ( r R ) δ ( ω k V ) .
Next, substituting (10) along with (A2), (A4) and (A6) into (8), and leaving only the terms corresponding to the induced potentials of walls, proportional to Δ n and Δ n , yields
u n i n ( ω k ; r ) = 8 π 2 Z e δ ( ω k V ) · { Δ n K n ( k r ) K n ( k R ) ,   Δ n I n ( k r ) I n ( k R ) , 1 D [ Δ n , 1 K n ( k r ) K n ( k R ) + Δ n , 2 I n ( k r ) I n ( k R ) ]
Here, and in the similar expressions below (namely, (14), (16), (30)) the upper, middle and lower expressions after the curly brace correspond to cases (a), (b), (c) in Figure 1, namely external configuration ( r > a ), internal configuration ( 0 < r < a ), and configuration between coaxial cylinders, ( a 1 < r < a 2 ). The functions Δ n , Δ n and Δ n , 1 ,   Δ n , 2 are defined in (A3), (A5) and after Equation (A7), D = 1 Δ n , 1 Δ n , 2 .
The induced Fourier-components E r , n i n ( ω k ; r ) of the electric field E ( r , φ , z , t ) are obtained using Equation (2) and identity E = ϕ
E r , n i n ( ω k ; r ) = d d r u n ( ω k ; r ) ,   E φ , n i n ( ω k ; r ) = i n r u n ( ω k ; r ) E z , n i n ( ω k ; r ) = i k u n ( ω k ; r )   .
The force on a charged particle is then F = Z e E i n with the electric field E i n obtained by the direct Fourier-transform using (12) and (11). Substituting Fourier-components (12) into the integral expression for E i n taken with instantaneous coordinates ( R , 0 , V t ) of the particle, integrating by frequency ω and using the analytical properties of Δ n , Δ n (even real components and odd imaginary ones) one obtains
F R = 2 ( Z e ) 2 π n = 0 d k k f n ( R ) ( k R , k V ) ,
f n ( R ) ( x , y ) = { K n ( x ) K n ( x ) Re   Δ n ( y ) I n ( x ) I n ( x ) Re   Δ n ( y ) K n ( x ) K n ( x ) Re ( Δ n , 1 ( y ) D ) + I n ( x ) I n ( x ) Re ( Δ n , 2 ( y ) D )
F z = 2 ( Z e ) 2 π n = 0 d k k f n ( z ) ( k R , k V ) ,
f n ( z ) ( x , y ) = { K n 2 ( x ) Im   Δ n ( y ) I n 2 ( x ) Im   Δ n ( y ) K n 2 ( x ) Im ( Δ n , 1 ( y ) D ) + I n 2 ( x ) Im ( Δ n , 2 ( y ) D ) .
In (14), the Bessel functions with a prime denote their derivatives with respect to the total argument x = k R . The force component F φ turns out to be zero, as expected from the symmetry. Equations (13)–(16) completely coincide with [16], where the authors considered the case of a charged particle moving inside a cylindrical channel.

2.3. VdW Forces on a Moving Neutral Particle

In the case of a neutral particle (an atom) with a fluctuating dipole moment d ( t ) , the charge density is ρ = div P with P ( r , φ , z , t ) = 1 r δ ( r R ) δ ( φ ) δ ( z V t ) d ( t ) being the polarization vector. Then, the Fourier-transform ρ n ( ω k ; r ) has the form
ρ n ( ω k ; r ) = [ δ ( r R ) r d r ( ω ) + i n r 2 δ ( r R ) d φ ( ω ) + i k r δ ( r R ) d z ( ω ) ] ,
where ω = ω k V . The induced Fourier-component u n i n ( ω k ; r ) in this case is calculated by subsituting (27) and the induced components of G n ( r , r ) from (A2), (A4), (A6), (A7) into (8). The resulting expressions are as follows.
(i)
external configuration, r > a (Figure 1a)
u n i n ( ω k ; r ) = 4 π Δ n K n ( k r ) [ K n ( k R ) | k | d r ( ω ) + i n R K n ( k R ) d φ ( ω ) + i k K n ( k R ) d z ( ω ) ]
(ii)
internal configuration, 0 < r < a (Figure 1b)
u n i n ( ω k ; r ) = 4 π Δ n I n ( k r ) [ I n ( k R ) | k | d r ( ω ) + i n R I n ( k R ) d φ ( ω ) + i k I n ( k R ) d z ( ω ) ]
(iii)
configuration between coaxial cylinders, a 1 < r < a 2 (Figure 1c)
u n i n ( ω k ; r ) = 4 π Δ n , 1 D K n ( k r ) [ K n ( k R ) | k | d r ( ω ) + i n R K n ( k R ) d φ ( ω ) + i k K n ( k R ) d z ( ω ) ] + 4 π Δ n , 2 D I n ( k r ) [ I n ( k R ) | k | d r ( ω ) + i n R I n ( k R ) d φ ( ω ) + i k I n ( k R ) d z ( ω ) ]
One can see that Equation (20) is a sum of (18) and (19) taken with the factor D 1 = ( 1 Δ n , 1 Δ n , 2 ) 1 . To find vdW forces acting on a particle, we use the expressions [31,33]
F r = r ( d s p E i n ) + r ( d i n E s p ) ,
F z = z ( d s p E i n ) + z ( d i n E s p ) ,
with standard definitions of the induced and spontaneous components of the dipole moment and the electric field. Instead of using (21), when calculating the radial (conservative) force component F r , it is simpler to find first the conservative potential (vdW energy)
U = 1 2 ( d s p E i n ) 1 2 ( d i n E s p )
Having determined U , the radial force F r is then calculated by F r = U / r .
Consider the first terms in (21) and (23), performing the Fourier-expansions for d ( t ) and E i n ( r , φ , z , t ) :
d s p ( t ) =   d ω 2 π d s p ( ω ) exp ( i ω t )
E i n ( r , φ , z , t ) = 1 2 π n = d ω 2 π d k 2 π E n i n ( ω k ; r ) exp ( i ( k z ω t + n φ ) )
Substituting (24), (25) into (22), (23) yields
F z ( 1 ) = 1 ( 2 π ) 4     d ω d ω d k ( i k ) d s p ( ω ) E n i n ( ω k ; R ) exp ( i ( ω + ω k V ) t )
U ( 1 ) = 1 2 ( 2 π ) 4     d ω d ω d k d s p ( ω ) E n i n ( ω k ; R ) exp ( i ( ω + ω k V ) t )
If not specified, the summation by n in (25)–(27) and in what follows is carried out from to + .
The induced components of the electric field are calculated by inserting (18)–(20) into (19). The appeared correlators of the dipole moment are found using the fluctuation-dissipation relation [34]
d i s p ( ω ) d k s p ( ω k V ) = 2 π δ i k δ ( ω + ω k V ) α ( ω k V ) coth ( ω k V ) 2 T 1 .
For all configurations shown in Figure 1, the correlators in (26), (27) can be written in the form
d s p ( ω ) E n i n ( ω k ; R ) = 8 π 2 δ ( ω + ω k V ) α ( ω k V ) coth ( ω k V ) 2 T 1 k 2 Z n ( k R , ω ) ,
where Z n ( k R , ω ) is given by
Z n ( k R , ω ) = { S n ( k R ) Δ n ( ω ) S n ( k R ) Δ n ( ω ) S n ( k R ) Δ n , 1 ( ω ) 1 Δ n , 1 ( ω ) Δ n , 2 ( ω ) + S n ( k R ) Δ n , 2 ( ω ) 1 Δ n , 1 ( ω ) Δ n , 2 ( ω )
S n ( x ) = K n 2 ( x ) ( n 2 x 2 + 1 ) + ( K n ( x ) 2 ,
S n ( x ) = I n 2 ( x ) ( n 2 x 2 + 1 ) + ( I n ( x ) ) 2 .
Note that in (31), (32) we assume that x = k R . Substituting (29) into (26), (27), integrating by frequency ω and making use of the transformation of integrals with allowance for the analytical properties of the functions under the integral yields
F z ( 1 ) = π 2 0 d ω + d k k 3 α ( ω + ) coth ( ω + 2 T 1 ) n = Im Z n ( k R , ω )
U ( 1 ) = 2 π 2 0 d ω + d k k 2 α ( ω + ) coth ( ω 2 T 1 ) n = Re Z n ( k R , ω )
where ω + = ω + k V ; α ( ω + ) and α ( ω + ) in what follows are the imaginary and real parts of polarizability α ( ω + ) .
To find the contribution from spontaneous field E s p to (21)–(23), we expand E s p into the Fourier integral in frequency and wave vector at the particle location point r 0 ( t ) = ( R , 0 , V t ) :
E s p ( r 0 ( t ) ) =   d ω 2 π d k 2 π E ω k s p ( R ) exp ( i ( ω t k V ) )
and express the induced dipole moment d i n ( t ) through E s p ( r 0 ( t ) ) :
d i n ( t ) =   d ω 2 π d k 2 π α ( ω k V ) E ω k s p ( R ) exp ( i ( ω k V ) t )
Substituting (35), (37) into (22), (23) yields
F z ( 2 ) = 1 ( 2 π ) 4   d ω d ω d k d k ( i k ) α ( ω k V ) E ω k s p ( R ) E ω k s p ( R ) exp ( i ( ω + ω k V ) t )
U ( 2 ) = 1 2 ( 2 π ) 4   d ω d ω d k d k α ( ω k V ) E ω k s p ( R ) E ω k s p ( R ) exp ( i ( ω + ω k V ) t )
Based on the theory of stationary electromagnetic fluctuations [34,35,36], the correlators in (37), (38) are worked out using the fluctuation-dissipation relation (see Appendix B):
E ω k ( R ) E ω k ( R ) = 2 ( 2 π ) 2 δ ( ω + ω ) δ ( k + k ) cot h ω 2 T n = Im   Z n ( k R , ω )
Substituting (39) into (37), (38), integrating by ω , k and simplifying yields
F z ( 2 ) = π 2 0 d ω + d k k 3 α ( ω + ) coth ( ω 2 T 2 ) n = Im   Z n ( k R , ω )
U ( 2 ) = 2 π 2 0 d ω + d k k 2 α ( ω + ) coth ( ω 2 T 2 ) n = Im   Z n ( k R , ω )
Summing (33) with (40) and (34) with (41) results in
F z = π 2 0 d ω + d k k 3 α ( ω + ) [ coth ( ω 2 T 2 ) coth ( ω + 2 T 1 ) ] n = Im   Z n ( k R , ω )
U = 2 π 2 0 d ω + d k k 2 [ α ( ω + ) coth ( ω 2 T 2 ) n = Im Z n ( k R , ω ) + α ( ω + ) coth ( ω + 2 T 1 ) n = Re   Z n ( k R , ω ) ]
Using (43), the radial force F R = U / R is then given by
F R = 2 π 2 0 d ω + d k k 3 [ α ( ω + ) coth ( ω 2 T 2 ) n = Im   Z n ( k R , ω ) + α ( ω + ) coth ( ω + 2 T 1 ) n = Re   Z n ( k R , ω ) ]
where Z n ( k R , ω ) denotes the derivative of Z n ( k R , ω ) with respect to its first argument.
Expressions (42)–(44) describe all configurations in Figure 1, when relating (a), (b), (c) to the first, second and third lines in the definition of Z n ( k R , ω ) in (30). Moreover, using (42)–(44), one can also retrieve all known results [31] for friction forces, dynamic vdW forces and interaction potentials in planar configurations. This is shown in Section 3.

2.4. Heating Rate of a Neutral Particle

The heating rate d Q / d t of a nanoparticle is calculated quite similar to vdW forces. The starting equation has the form [31,33]
d Q / d t = ( d ˙ sp E in ) + ( d ˙ in E sp ) ,
where the overdots for the dipole moments stand for time derivatives. Further calculations are very similar to those performed in Section 2.3. The resulting expression is simply obtained when substituting ( ω + ) for k under the integrand in (42):
d Q / d t = π 2 0 d ω + d k k 2 ω + α ( ω + ) [ coth ( ω 2 T 2 ) coth ( ω + 2 T 1 ) ] n = Im Z n ( k R , ω )

3. Transition from Cylindrical to Plane-Geometry Surfaces

Using (42)–(44), one can retrieve all known results for planar configurations. We consider, as an example, the friction force F z . In the coordinates chosen in [31], the relevant plane-geometry result is
F x = π 2 0 d ω + d k x k x + d k y exp ( 2 k z ) α ( ω + ) [ coth ( ω 2 T 2 ) coth ( ω 2 T 1 ) ] Δ ( ω )
In (47), in contrast to (42), the particle moves along the x -direction, ω + = ω + k x V , k = k x 2 + k y 2 , z is the distance to surface, and Δ ( ω ) = ε ( ω ) 1 ε ( ω ) + 1 is the p-wave reflection factor in the electrostatic limit. Performing the integration over k y in Equation (47) yields
F z = π 2 0 d ω + d k k 3 α ( ω + ) [ 2 K 0 ( 2 k z ) + K 1 ( 2 k z ) k z ] [ coth ( ω 2 T 2 ) coth ( ω + 2 T 1 ) ] Δ ( ω )
Equation (48) now has a form very similar to (42). Using (30), (31) and (A3), it is easy to check that
n = Im Z n ( k R , ω ) = n = [ K n 2 ( x ) ( n 2 x 2 + 1 ) + ( d K n ( x ) d x ) 2 ] I n ( y ) K n ( y ) Im ( ε 1 ε [ I n ( y ) K n ( y ) / I n ( y ) K n ( y ) ] )
where x = k R ,   y = k a . Taking the limit a ,   R ,   R a = z = const , and making use of the asymptotic relations K n ( x ) ~ π / 2 x exp ( x ) ,   K n ( x ) K n ( x ) , I n ( x )   I n ( x ) ~ 1 / 2 π x exp ( x )   , we see that the factor Im   ( ) in (49) reduces to Im   Δ ( ω ) = Δ ( ω ) . Finally, making use of the relations [4] (for x ,   y ,   x y = z )
n = K n 2 ( x ) I n ( y ) K n ( y ) = K 0 ( 2 ( x y ) ) ,   n = n 2 x 2 K n 2 ( x ) I n ( y ) K n ( y ) = K 1 ( 2 ( x y ) ) 2 ( x y )
  n = ( K n ( x ) ) 2 I n ( y ) K n ( y ) = K 0 ( 2 ( x y ) ) + K 1 ( 2 ( x y ) ) 2 ( x y ) ,  
we find
n = Im Z n ( k R , ω ) = [ 2 K 0 ( 2 k z ) + K 1 ( 2 k z ) k z ] Im ( ε 1 ε + 1 )
Therefore, Equation (42) reduces to plane-geometry result (47). The identities (50) also hold valid for concave cylindrical surfaces when replacing K n ( x ) I n ( x ) and x y y x = z . In this case, Z n ( k R , ω ) is defined by the second-line in (30). Using (30), we can also check the result [31] for a particle moving in a plane-parallel gap. In the same way, one can verify the transformation properties of Equations (43) and (44).

4. VdW Energy and Stopping (Friction) Force

For a more detailed study of dynamic effects, we represent (43) in the form
U = U ( 0 ) + U ( 1 ) ,
U ( 0 ) = π 2 0 d ω + d k k 2 cot h ( ω 2 T 2 ) Im [ α ( ω + ) n = 0 ϵ n Z n ( k R , ω ) ]
U ( 1 ) = π 2 0 d ω + d k k 2 α ( ω + ) [ coth ( ω 2 T 1 ) coth ( ω 2 T 2 ) ] n = 0 ϵ n Re   Z n ( k R , ω )
where ϵ n = 1 for n 0 , and ϵ n = 1/2 for n = 0 , and the property Z n ( x , y ) = Z n ( x , y ) is taken into account. Here, U ( 0 ) is the part of the vdW energy that remains nonzero at V = 0 regardless of temperature T 2 . The behavior of U ( 1 ) is more complicated. At T 1 = T 2 = 0 , taking into account the analytic properties of functions in (53) and making use of the identity lim T 0 coth ( x / T ) = sign ( x ) , one obtains
U ( 0 ) = π 2 0 d ξ + d k k 2 Im [ i α ( i ξ + k V ) n = 0 ϵ n Z n ( k R , i ξ ) ]
U ( 1 ) = 2 π 2 0 d k k 2 0 k V d ω α ( ω ) n = 0 ϵ n Re   Z n ( k R , ω )
where ω = ω k V . In the case V = 0 , Equation (55) coincides with the well-known static results [1,2,3,4,5], whereas one obtains U ( 1 ) 0 only for V 0 . To proceed further, we consider the oscillator model of atomic polarizability
α ( ω ) = α ( 0 ) ω 0 2 ω 0 2 ω 2 + i 0 sign ( ω ) = P α ( 0 ) ω 0 2 ω 0 2 ω 2 + i π α ( 0 ) ω 0 2 [ δ ( ω ω 0 ) δ ( ω + ω 0 ) ]
with α ( 0 ) and ω 0 being the static polarizability and resonance frequency. Substituting (58) into (56), (57) and simplifying, yields
U ( 0 ) = 2 α ( 0 ) ω 0 2 π 2 0 d ω 0 + d k k 2 ( ω 0 2 + ω 2 k 2 V 2 ) [ ( ω 0 2 + ω 2 k 2 V 2 ) 2 + 4 ω 2 k 2 V 2 ] n = 0 ϵ n Z n ( k R , i ω )
U ( 1 ) = ω 0 α ( 0 ) π ω 0 / V d k k 2 n = 0 ϵ n Re   Z n ( k R , k V ω 0 )
In the case of a flat surface at V 0 , Equations (56) and (59) reduce to [28,29,30]. The plane-geometry equivalent of (60) was first obtained in Ref. [30]. The analysis performed in [30] showed that U ( 0 ) has a V 2 poportional (attractive) correction to the energy at small V . At larger velocities, U ( 0 ) reaches a maximum and then tends to zero. An intriguing fact (in the case of a flat surface) is that U ( 1 ) can be positive at a certain choice of the atomic frequency ω 0 and the plasma frequency of the metal, but the resulting potential remains attractive. However, for metallized cylindrical filaments, channels and atomic chains, as we will see in Section 5, the total potential can be repulsive. In contrast to this, for dielectric (nonconducting) surfaces the resulting vdW potential proves to be attractive (as usual).
For V = 0 and T 2 0 , using the standard transformation of the integration contour to the complex plane for U ( 0 ) , Equation (53) takes the known form (Equation (22) in [3])
U ( 0 ) = 4 T 2 π m = 0 ϵ m α ( i ξ m ) 0 + d k k 2 n = 0 ϵ n Z n ( k R , i ξ m )
with ξ m = 2 π m T 2 .
As for the stopping force (42) at T 1 = T 2 = 0 , we obtain
F z = 2 π 2 0 d k k 3 0 k V d ω α ( ω k V ) n = 0 ϵ n Im   Z n ( k R , ω )
Inserting (58) into (62) yields
F z = ω 0 α ( 0 ) π ω 0 / V d k k 3 n = 0 ϵ n Im   Z n ( k R , k V ω 0 )
It is worthwhile noting that function Z n ( k R , ω ) has an exponentially decaying asymptotics at large k R . Then, according to (63), the force F z becomes negligibly small at V ω 0 R . If ω 0 ~ 10 15   s 1 (for atoms) and R ~ 10   n m , this implies V 10 7   m / s . The corresponding quantum friction limit [37] for atoms is retrieved when using Equation (61) with a more precise definition of the dressed atomic polarizability [38]. We neglect here corrections with respect to the local equilibrium case and spatial dispersion [39,40,41,42].
The case of the friction force F z acting on nanoparticles with a linear velocity dependence (at nonzero temperature and V ω 0 R ) follows from (42) when expanding α ( ω + ) and cot h ( ω + 2 T 1 ) in powers of ( k V ) , retaining the first-order expansion terms:
F z = 4 V π 2 0 d ω { α ( ω ) ω [ coth ( ω 2 T 2 ) coth ( ω 2 T 1 ) ] + α ( ω ) ( ω ) coth ( ω 2 T 1 ) } 0 d k k 4 n = 0 ϵ n Im Z n ( k R , ω )
For T 1 = T 2 = T , Equation (64) reduces to
F z = 4 V π 2 0 d ω α ( ω ) sin h 2 ( ω / 2 T ) 0 d k k 4 n = 0 ϵ n Im Z n ( k R , ω )
The validity of (64), (65) is restricted by the velocities V 3.9 · 10 6 m / s for R = 10 nm , T = 300 K , since the characteristic frequency here is ω 0 ~ k B T / . At higher velocities and (or) R, T values, the nonretarded relativistic consideration is required.

5. Interactions with Cylindrical Filaments and Atomic Chains

In the case of configuration shown in Figure 1a, as follows from the general formulas the main contribution to the k   –integrals comes from the wave vectors k 1 / R . Therefore, for an atom at large distance R compared to the radius a of the cylinder, we have a R and k a 1 .   Therefore , we can expand Δ n in powers of k a . Then, taking into account the asymptotics of cylindrical functions for small values of the argument ( z 1 ) [43]
I 0 ( z ) 1 + z 2 4   ,   I n ( z ) 1 Γ ( n + 1 ) ( z 2 ) n ,   n = 1 , 2 ,
K 0 ( z ) ln ( γ z 2 ) , K n ( z ) Γ ( n ) 2 ( 2 z ) n ,   n = 1 , 2
where γ = e 0.577 (with 0.577 being the Euler–Mascheroni constant), it is sufficient to retain only the first two terms in the sums by n , including Δ 0 ( ω ) and Δ 1 ( ω ) . Using (11) and (66), (67) yields
Δ 0 ( ω ) = ( k a ) 2 4 ε 1 1 + ε ( k a ) 2 2 ln ( 2 / γ k a )
Δ 1 ( ω ) = ( k a ) 2 2 ε 1 ε + 1
For n 2   , other functions Δ n ( ω ) and Z n ( k R , ω ) in (49) can be neglected since they contain higher powers of k a .
In the case of a dielectric filament (except strongly polar ones with large ε ), one can neglect the term with ε in (68). Then, the sums in Equations (54)–(65) take the form
ϵ n Z n = 0.5 Z 0 + Z 1 = ( k a ) 2 2 ( ε 1 ) [ 1 2 R 1 ( x ) + 1 ( ε + 1 ) R 2 ( x ) ] ,
R 1 ( x ) = K 0 2 ( x ) + K 1 2 ( x ) ,
R 2 ( x ) = K 1 2 ( x ) + K 1 2 ( x ) x 2 + ( d K 1 ( x ) / d x ) 2 = R 1 ( x ) + 2 ( K 1 2 ( x ) x 2 + K 0 ( x ) K 1 ( x ) x ) .
Note that in (70)–(72) we use x = k R . For V = 0 , making use of the integral table [43]
0 d x x a 1 K m ( x ) K n ( x ) = 2 a 3 Γ ( a ) Γ ( a + m + n 2 ) Γ ( a + m n 2 ) Γ ( a m + n 2 ) Γ ( a m n 2 )
we obtain from (56) or (61) (in the limit T 2 0 ) the known result [3,44]
U ( R ) = 9 128 R 3 ( a R ) 2 0 d ω α ( i ω ) ( ε ( i ω ) 1 ) ( ε ( i ω ) + 7 ) ( ε ( i ω ) + 1 )
The dynamic potentials (59), (60) can be worked out using a single-oscillator model
ε ( ω ) = a + b ω 1 2 ω 1 2 ω 2 i γ ω
with parameters a , b , ω 1 , γ chosen to fit the UV range of the spectrum when calculating potentials (59), (60), and the IR range when calculating friction forces (64) and (65).
When computing the interaction of an atom with a cylindrical metal filament, the contributions from (68) and (69) are essentially different since we should retain the whole denominator in (68).
In the static case ( V = 0 ) , using the plasma limit of the Drude permittivity ε ( i ω ) = 1 + ω p 2 / ω 2 with allowance for the condition ( a / R ) 2 ln ( 1.123 R / a ) 1 (here 1.123 = 2 / e 0.577 ) taking into account only the term n = 0 in (59), and changing the order of integration (with the ω integration being the first), one obtains
U ( R ) = α ( 0 ) a 2 2 π R 4   0 d x x 3 R 1 ( x ) ln 1 2 ( 1.123 x R a ) ω p α ( 0 ) a 2 2 π R 4 ln 1 / 2 ( 1.123 R / a )
This is exactly the well-known result (Equation (43) in [44]) obtained in the electrostatic approximation, using another method. The distance dependence in (76) is a surprising feature of the vdW interaction caused by the nonadditivity and plasmon spectra of a conducting filament. In the limit of high temperatures, retaining in (61) the term with m = n = 0 , we retrieve another result of [44]
U ( R ) = π α ( 0 ) T 8 R 3 ln ( 1.123 R / a )
For a metal filament, dynamic potentials (59), (60) can be simplified further. The resulting expression U = U ( 0 ) + U ( 1 ) we cast in the form
U ( R , V ) = ω 0 α ( 0 ) 8 π R 3 ( a R ) 2 ( ω p ω 0 ) 2 ( U 1 + U 2 + 2 U 3 + 2 U 4 )
The contributions U 1 4 are defined by Equations (A16)–(A19) in Appendix C. The terms U 1 , U 2 correspond to terms n = 0 , 1 in (59), and U 3 , U 4 –to the similar ones in (60). The static result (76) follows from (A16).
When calculating stopping force F z by means of (63), we use the Drude model ε ( ω ) = 1 ω p 2 / ( ω 2 + i γ ω ) . As shown in in Appendix C, Equation (63) takes the form
F z = γ α ( 0 ) 4 π R 4 ( a R ) 2 ( ω p ω 0 ) 2 ( F 1 + 2 F 2 )
with F 1 and F 2 defined in (A21) and (A22).
We now pass to the interactions of an atom with atomic chains. For dielectric chains with atomic spacing d ,   the transition is very simple, using the limit of rarified medium ε ( ω ) 1 = 4 π n α 2 ( ω ) , where n 0 , α 2 ( ω ) is the polarizability of an atom of the string and n = 1 / π a 2 d . Substituting these relations into (68), (69) yields
  ϵ n Z n = 2 k 2 d α 2 ( ω ) [ K 0 2 ( x ) + K 1 2 ( x ) + K 1 2 ( x ) x 2 + K 0 ( x ) K 1 ( x ) x ]
Using (79) and (57) with parameters α 2 ( 0 ) and ω 02 of the atomic chain, and ω 01 and α 1 ( 0 ) of the projectile atom, Equations (60), (61) and (63) take the form
U ( 0 ) = ω 01 α 1 ( 0 ) α 2 ( 0 ) 2 π d R 5 ( ω 02 ω 01 ) 2 0 d x x 4 [ R 1 ( x ) + R 2 ( x ) ] f 1 ( p 1 , b 1 x )
U ( 1 ) = ω 01 α 1 ( 0 ) α 2 ( 0 ) π d R 5 ( ω 02 ω 01 ) 2 P 1 / b 1 d x x 4 [ ( b 1 x 1 ) 2 p 1 2 ] [ R 1 ( x ) + R 2 ( x ) ]
F z = ω 01 α 1 ( 0 ) α 2 ( 0 ) 2 d R 5 ( ω 02 R V ) φ 1 ( ( ω 01 + ω 02 ) R / V )
where f 1 ( x , y ) is defined by (A15), φ 1 ( x ) = x 5 ( R 1 ( x ) + R 2 ( x ) ) , b 1 = V / ω 01 R , p 1 = ω 02 / ω 01 , and one should take the principal value of the integral in (82).
For a metal string, the rarefied-gas approximation is not valid due to remaining cooperative effects, giving rise to certain effective parameters of a one-dimensional plasma. Therefore, to describe the interactions between a projectile atom and a metallized atomic chain, a more adequate way is to use (78) and (79) with a d / 2 and the appropriate plasma parameters of the atomic chain. In the case of undamped plasmons ( γ 0 + ) , the expressions (A21) and (A22) can be simplified further and (79) reduced to (using a = d / 2 ) :
F z ω 0 α ( 0 ) 4 2 R 4 ln 1 2 ( 2.246 R a ) ( a R ) ( ω p ω 0 ) φ 1 ( b ) R 1 ( b ) ω 0 α ( 0 ) 4 R 4 ( a R ) 2 ( ω p ω 0 ) 2 φ 1 ( ( 1 + p ) / b ) R 2 ( ( 1 + p ) / b )
where b 1 = b p φ 1 2 ( 2 R / d ) , b = V / ω 0 R (see Appendix C). The second term in (84) is nearly d / R 2 times (or more) less than the first one, since (usually) ( 1 + p ) / b > b (note that φ 1 ( x ) decays exponentially fast at x 1 ). The more gentle dependence on R   of the first term in (84), is a result of one-dimensional plasmon dispersion, similar to (76) [44].

6. Numerical Analysis

Due to the plethora of new analytical results obtained above, we will restrict ourselves by demonstrating only a few numerical examples (Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6) of the dynamic vdW potentials and stopping forces, namely the interactions of neutral atoms (in particular, Cs) with metal and dielectric filaments and atomic chains of Au, Si, and interactions with the walls of a capillary metallized with gold.
Figure 2 shows the ratio between the dynamic vdW potential U d and the modulus of the static potential U s shown as a function of the reduced velocity–distance parameter V / ω o R for an atom in parallel motion to a metal filament, ω o is the characteristic atomic frequency, ω p is the metal plasma frequency. The calculations were carried out using Equation (78) with a / R = 0.1 and ω p / ω 0 = 0.3 ,   1 ,   3 (red, blue and green lines). The static potentials U s   in Figure 2 and in Figure 4 and Figure 5 were computed by the corresponding formulas for U d assuming V = 0 . It should be noted that in this velocity range, the calculations of U d are nearly independent on the choice of γ . The nonzero value of γ simplifies calculating the integral U 4 in (78) by using (A19). For γ = 0 , the result is the same, but one should take the principal value of this integral.
The remarkable finding is that the dynamic vdW potential can be positive (repulsive) (panel (b)), though the coupling constant is an order (or more) less than that of the static potential U s . At moderate velocity–distance parameters V / ω o R ~ 0.5 , the dynamic potential is attractive, and the coupling constant is 1.5 times higher than the static one.
Figure 3 plots the dynamic vdW potential U d for a Cs atom as a function of distance R to a metal chain of Au atoms (Equation (78) with parameters a = d / 2 = 0.2   nm ,   ω p = 9   eV ,   γ = 0.03   eV , corresponding to gold. For a comparison, shown is the string continuous potential in the Moliere approximation for the Cs-Au interatomic potential U ( r ) = Z 1 Z 2 e 2 r 1 Φ ( r / a s ) (dashed green), with Z 1 and Z 2 being the charges of the nuclei (55 and 79 for Cs and Au). Here, Φ ( r / a s ) is the screening function and a s = 0.04683 ( Z 1 + Z 2 ) 2 / 3 nm–the screening length [45]. The red and blue lines correspond to velocities of 3 · 10 7 and 1.5 · 10 7 m / s . For the Cs atom, we used in (57) the values α ( 0 ) = 0.057   nm 3 and ω 0 = 3.89   eV . Parameters ω p and γ are the same as in Figure 2. The peak structure of the red and blue curves at R < 0.5   nm is likely due to the violation of the used approximation a / R 1 . Figure 3 also demonstrates the repulsive dynamic vdW potential. In the range 0 < V / ω o R < 0.6 , the vdW interaction increases in absolute value relative to the static one (Figure 2a). This is similar to the case of a flat surface [30], but the repulsive effect was not found there. The observed feature is the pure result of the cylindrical symmetry and conductive properties of the metal filaments and atomic chains. Mathematically, the repulsion appears due to the negative contribution to the integral (A14) from the term y 2 in the numerator, which is proportional to V 2 . This leads to a change in sign in (59) and is enhanced by an increase (in absolute value) of the negative contribution to the interaction resulting from (6).
It is worthwhile noting that repulsive Casimir and Casimir–Polder forces known in so far were obtained only in static configurations of two parallel semi-infinite dielectric plates if they are separated by a dielectric fluid that is intermediate between those of the dielectric plates [46].
On the contrary, the interactions of the Cs atom with the dielectric filament and the chain (the Si chain was chosen as an example) turn out to be attractive (as in the static case). This can be seen from Figure 4 showing the ratio between the dynamic potential U d and the absolute values of the static potential | U s | for the Cs atom and the atomic chain of Si. The red and blue lines and their sum (green) were calculated using Equations (81) and (82). Both the U ( 0 ) and U ( 1 ) contributions have the negative sign, but the coupling constant of the main distance dependence ~ R 5 increases in the range 0 < V / ω 01 R < 3 ( ω 01 = 3.89   eV here means the atomic frequency of Cs).
The repulsive vdW interaction also manifests itself in the case of the motion of Cs atoms inside a capillary with Au metallized walls. The results are plotted in Figure 5 with the same color convention as in Figure 4: the red and blue lines correspond to the contributions U ( 0 ) and U ( 1 ) , calculated by Equations (58) and (59). The green line and the inset show the resulting potential. All parameters for the Cs atom and gold are the same as those used when plotting Figure 3. The change in the sign of the vdW potential in this case occurs close to the value V / ω o R = 1.5 . The toothlike structure results from the contributions of multiple harmonics in the sums of (58) and (59).
Finally, we calculated the stopping force (Equation (63)) for the Cs atom in a capillary with Au metallized walls (Figure 6).
The peak structure in Figure 6 is caused by the excitation of plasmons. This is similar to the case of the excitation of plasmons by charged particles in cylindrical channels [16].
The numerical examples presented here are illustrative, and additional analysis is required to clarify many of the subtle points of the vdW dynamic interactions that are beyond the scope of this article. Among them are the presence of many multiple resonances of the dielectric response, the transition region from the vdW to the short-range interactions with atoms of walls, the nonlocality of the dielectric permittivity, which increases at small separations, etc.

7. Conclusions

In this paper, we refined some points and developed further the nonrelativistic theory of vdW interactions; namely, we calculated the attractive potential and stopping (friction) forces between neutral ground-state atoms and cylindrical surfaces of different configurations. As a constituent part of the theory, the electromagnetic image forces for the charged particles in the systems under consideration were also calculated. The transition from the case of a cylindrical surface to a flat surface is demonstrated and the obtained expressions for image forces and vdW forces and potentials are shown to completely coincide with all known results in plane and cylindrical geometries. The formulas describing the configurations with convex, concave cylindrical surfaces and two coaxial cylinders with a particle between them are represented in a universal form. A new remarkable result is that the dynamic vdW potential can be positive (repulsive) in the case of metallized atomic chains, filaments and capillaries, providing a steering effect for the channeled neutral atoms.

Funding

This research received no external funding.

Acknowledgments

This article is dedicated to the blessed memory of Arthur A. Kyasov, who made a fundamental contribution at the initial stage of formulating ideas that were further developed in this research, to whom I owe happy decades of dense friendship and fruitful cooperation.

Conflicts of Interest

The author declare no conflict of interest with any parties.

Abbreviations

The following abbreviation is used in the manuscript.
vdW van der Waals.

Appendix A

In the case of a particle outside cylinder, we seek G n ( r , r ) in the form
G n ( r , r ) = { D n I n ( k r ) , 0 r a A n K n ( k r ) + B n I n ( k r ) , a < r r C n K n ( k r ) , r > r
with A n ,   B n ,   C n ,   D n being the constants. Having determined them using (A1), (5) and (7), one obtains
G n ( r , r ) = { r [ Δ n K n ( k r ) I n ( k r ) K n ( k a ) I n ( k a ) K n ( k r ) I n ( k r ) ] , 0 < r a r [ Δ n K n ( k r ) K n ( k r ) K n ( k r ) I n ( k r ) ] , a < r r r [ Δ n K n ( k r ) K n ( k r ) I n ( k r ) K n ( k r ) ] , r > r
Δ n ( ω , k ) = ( ε 1 ) I n ( k a ) I n ( k a ) ε K n ( k a ) I n ( k a ) I n ( k a ) K n ( k a )
Here, K n ( x ) = d K n ( x ) / d x , I n ( x ) = d I n ( x ) / d x , and arguments ω in ε and ω ,   k in Δ n , Δ n (see below and in the main text) are omitted for brevity, or ω is retained when writing Δ n , Δ n . In the case when the particle is inside the cylindrical channel, the Green’s function is calculated similar to (A1). The resulting expressions are
G n ( r , r ) = { r [ Δ n I n ( k r ) K n ( k r ) K n ( k a ) / I n ( k a ) I n ( k r ) K n ( k r ) ] , r > a r [ Δ n I n ( k r ) I n ( k r ) K n ( k r ) I n ( k r ) ] , 0 < r r r [ Δ n I n ( k r ) I n ( k r ) I n ( k r ) K n ( k r ) ] , r < r a
Δ n ( ω , k ) = ( ε 1 ) K n ( k a ) K n ( k a ) ε I n ( k a ) K n ( k a ) K n ( k a ) I n ( k a )
In (A2), (A4), the terms including Δ n and Δ n correspond to the induced potential of the cylinder walls, whereas the other terms–to the “bare” electric potential of a point charge 1 / 4 π .
In the case when the particle is between two coaxial cylinders with inner and outer radii a 1 and a 2 , the Green’s function is given by
(i)
a 1 < r r :
G n ( r , r ) = r I n ( k r ) K n ( k r ) + r D [ Δ n , 1 K n ( k r ) K n ( k r ) + Δ n , 2 I n ( k r ) I n ( k r ) ] r D Δ n , 1 Δ n , 2 [ I n ( k r ) K n ( k r ) + K n ( k r ) I n ( k r ) ]
(ii)
r < r a 2
G n ( r , r ) = r K n ( k r ) I n ( k r ) + r D [ Δ n , 1 K n ( k r ) K n ( k r ) + Δ n , 2 I n ( k r ) I n ( k r ) ] r D [ Δ n , 1 K n ( k r ) K n ( k r ) + Δ n , 2 I n ( k r ) I n ( k r ) ]
Here, D = 1 Δ n , 1 Δ n , 2 , and Δ n , 1 , Δ n , 2 are defined in (A3) and (A5) taking ε = ε 1 ,   a = a 1 and ε = ε 2 ,   a = a 2   , respectively (see Figure 1c). Using (8), (A2), (A4), (A6) and (A7) it is easy to consider particular cases.

Appendix B

Within the theory of stationary electromagnetic fluctuations [34,35,36], as applied to our case, the correlator of spontaneous electric fields generated by a dielectric medium, is
E ω k ( r ) E ω k E ( r ) = ( 2 π ) 2 δ ( ω + ω ) δ ( k + k ) ω 2 c 2 cot h ω 2 T Im D i i ( ω k , r , r )
where D i k ( ω k , r , r ) is the retarded Green’s function of a photon (assuming that the medium is homogeneous, isotropic and nonmagnetic), and the summation over i is implied in (A8) ( i = r , φ , z ) . At fixed values of m and r we have [35]
D i k ( ω k , r , r ) = E l ( ω k , r ) = l ϕ ω k i n ( r )
where ϕ ω k i n ( r ) is the Fourier-transform of the induced potential produced by the point electric dipole d l ( ω ) = ( c 2 / ω 2 ) δ l m located at point r . To find this potential, one should replace the quantities d l ( ω ) ( l = r , φ , z ) in (18)–(20) by the above components. Using (20) to obtain ϕ ω k i n ( r ) and substituting the result into (A9) yields
D r r ( ω k , r , R ) = ( 2 c 2 ω 2 ) n = e i n φ k D 2 [ Δ n , 1 K n ( k r ) K n ( k R ) + Δ n , 2 I n ( k r ) I n ( k R ) ]
D φ φ ( ω k , r , R ) = ( 2 c 2 ω 2 ) n = e i n φ n 2 D R 2 [ Δ n , 1 K n ( k r ) K n ( k R ) + Δ n , 2 I n ( k r ) I n ( k R ) ]
D z z ( ω k , r , R ) = ( 2 c 2 ω 2 ) n = e i n φ k D 2 [ Δ n , 1 K n ( k r ) K n ( k R ) + Δ n , 2 I n ( k r ) I n ( k R ) ]
where D = 1 Δ n , 1 Δ n , 2 . For other configurations shown in Figure 1, one should take either Δ n , 1 = 0 ., D = 1 (Figure 1a), or Δ n , 2 = 0 , D = 1 (Figure 1b). Substituting (A10)–(A12) taken at φ = 0 ,   r = r = R into (A8), we find a unified fluctuation-dissipation relation, which is valid for all considered cases
E ω k ( R ) E ω k ( R ) = 2 ( 2 π ) 2 δ ( ω + ω ) δ ( k + k ) cot h ω 2 T n = Im Z n ( k R , ω )
where Z n ( k R , ω ) is defined in (30).

Appendix C

We use the plasma model ε ( i ω ) = 1 + ω p 2 / ω 2 to calculate (59) and Drude model ε ( ω ) = 1 ω p 2 / ( ω 2 + i γ ω ) to calculate (60), introducing the parametrization y = ω / ω 0 , p = ω p / ω 0 ,   τ = γ / ω 0 , b = V / ω 0 R , x = k R . Changing the order of integration in (58) with allowance for (68) and (69), the integration by ω in (59) is carried out, yielding the integral of the form
f 1 ( x , y ) 4 π 0 d z ( 1 + z 2 y 2 ) ( z 2 + x 2 ) [ z 2 + ( y 1 ) 2 ] [ z 2 + ( y + 1 ) 2 ]
Working out integral (A14) by making use of a partial fraction expansion explicitly yields
f 1 ( x , y ) = 1 x 2 ( y + 1 ) 2 sign ( y 1 ) x 2 ( y 1 ) 2 + 2 ( 1 x 2 y 2 ) x [ x 2 ( y 1 ) 2 ] [ x 2 ( y + 1 ) 2 ]
Using (A15), (71) and (72) and the previously used parametrization, the expressions for U 1 and U 2 in (78) take the form
U 1 = 0 d x x 4 R 1 ( x ) f 1 ( p φ 1 / 2 ( a x / R ) , b x )
U 2 = 0 d x x 4 R 2 ( x ) f 1 ( p / 2 , b x )
where φ ( x ) = 0.5 x 2 ln ( 1.123 / x ) . The numerical coefficient 1.123 stems from 2 / γ in (67). With the same notation, expressions for U 3 and U 4 in (78) reduce to
U 3 = 1 / b d x x 4 R 1 ( x ) f 2 ( b x 1 , p φ 1 / 2 ( a x / R ) , τ )
U 4 = 1 / b d x x 4 R 2 ( x ) f 2 ( b x 1 , p / 2 , τ )
where R 1 ( x ) ,   R 2 ( x ) are defined in Equations (71) and (72) and f 2 ( x , y , z ) reads
f 2 ( x , y , z ) = x 2 y 2 ( x 2 y 2 ) 2 + x 2 z 2
Using the same parametrization in (63), expressions for F 1 and F 2 in (79) take the form
F 1 = 1 / b d x x 5 R 1 ( x ) ( b x 1 ) [ ( b x 1 ) 2 p 2 φ ( a x / R ) ] 2 + τ 2 ( b x 1 ) 2
F 2 = 1 / b d x x 5 R 2 ( x ) ( b x 1 ) [ ( b x 1 ) 2 p 2 ] 2 + τ 2 ( b x 1 ) 2

References

  1. Schmeits, M.; Lucas, M. Physical adsorption and surface plasmons. Surf. Sci. 1977, 64, 176–196. [Google Scholar] [CrossRef]
  2. Schmeits, M.; Lucas, M. Physical adsorption and surface plasmons. Prog. Surf. Sci. 1983, 14, 1–51. [Google Scholar] [CrossRef]
  3. Nabutovskii, V.M.; Belosludov, V.R.; Korotkikh, A.M. Interaction potential between small neutral particle and spherical or cylindrical surfaces. Zh. Eksp. Teor. Fiz. 1979, 77, 700–706. [Google Scholar]
  4. Marvin, A.M.; Toigo, F. VdW interaction between a point particle and a metallic surface. 1. Theory. Phys. Rev. A 1982, 25, 782–802. [Google Scholar] [CrossRef]
  5. Eberlein, C.; Zietal, R. Force on a neutral atom near conducting microstructures. Phys. Rev. A 2007, 75, 32516. [Google Scholar] [CrossRef] [Green Version]
  6. Bezerra, V.B.; de Mello, B.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Saharian, A.A. Exact Casimir-Polder potential between a particle and an ideal metal cylindrical shell and the proximity force approximation. Eur. Phys. J. C 2011, 71, 1614. [Google Scholar] [CrossRef] [Green Version]
  7. Bordag, M.; Klimchitskaya, G.L.; Mohideen, U.; Mostepanenko, V.M. Advances in the Casimir Effect; Oxford University Press: Oxford, UK, 2009. [Google Scholar]
  8. Kyasov, A.A.; Dedkov, G.V. Electromagnetic interactions of moving particles with cylindrical surfaces. Surf. Sci. 2001, 491, 124. [Google Scholar] [CrossRef]
  9. Dedkov, G.V.; Kyasov, A.A. Nonrelativistic theory of electromagnetic forces on particles and nanoprobes moving near a surface. Phys. Low-Dim. Struct. 2003, 1/2, 1–86. [Google Scholar]
  10. Kyasov, A.A. Fluctuation-Electromagnetic Interaction of Moving Particles with Surfaces. Ph.D. Thesis, Kabardino-Balkarian State University, Nalchik, Russia, 2004. (In Russian). [Google Scholar]
  11. Iijima, S. Helical microtubules of graphitic carbon. Nature 1991, 354, 56–58. [Google Scholar] [CrossRef]
  12. Klimov, V.V.; Letokhov, V.S. Hard X-ray radiation emitted by a charged particle moving in a carbon nanotube. Phys. Lett. A 1996, 222, 424–428. [Google Scholar] [CrossRef]
  13. Dedkov, G.V. Fullerene nanotubes can be used when transporting gamma-quanta, neutrons, ion beams and radiation from relativistic particles. Nucl. Instr. Meth. Phys. Res. 1998, 143, 584–590. [Google Scholar] [CrossRef]
  14. Gevorgian, L.A.; Ispirian, K.A.; Ispirian, R.K. High-energy particle channeling in nanotubes. Nucl. Instr. Meth. Phys. Res. 1998, 145, 155–159. [Google Scholar] [CrossRef]
  15. Zhevago, N.K.; Glebov, V.I. Channeling of fast charged and neutral particles in nanotubes. Phys. Lett. A 1998, 250, 360–368. [Google Scholar] [CrossRef]
  16. Arista, N.R.; Fuentes, M.A. Interaction of charged particles with surface plasmons in cylindrical channels in solids. Phys. Rev. B 2001, 63, 165401. [Google Scholar] [CrossRef] [Green Version]
  17. Biryukov, V.M.; Bellucci, S. Nanotube diameter optimal for channeling of of high-energy particle beam. Phys. Lett. B 2002, 542, 111–115. [Google Scholar] [CrossRef] [Green Version]
  18. Greenenko, A.A.; Shul’ga, N.F. Fast ion passing through straight and bent nanotubes. Nucl. Instr. Meth. Phys. Res. B 2003, 205, 767–772. [Google Scholar] [CrossRef]
  19. Krasheninnikov, A.A.; Nordlund, K. Multiwalled carbon nanotubes as apertures and conduits for energetic ions. Phys. Rev. B 2005, 71, 245408. [Google Scholar] [CrossRef] [Green Version]
  20. Mowbray, D.J.; Chung, S.; Miškovič, Z.L. Dynamic interactions of fast ions with carbon nanotubes. Nucl. Instr. Meth. Phys. Res. B 2005, 230, 142. [Google Scholar] [CrossRef]
  21. Artru, X.; Fomin, S.P.; Shul’ga, N.F.; Ispirian, K.A.; Zhevago, N.K. Carbon nanotubes and fullerites in high-energy and X-ray physics. Phys. Rep. 2005, 412, 89. [Google Scholar] [CrossRef]
  22. Borka, D.; Petrovič, S.; Neškovič, N. Influence of the dynamical image potential on the rainbows in ion channeling through short carbon nanotubes. Nucl. Instr. Meth. Phys. Res. B 2006, 73, 62902. [Google Scholar] [CrossRef]
  23. Mišcovič, Z.L. Ion Channeling through carbon nanotubes. Radiat. Eff. Defects Solids 2007, 162, 185–205. [Google Scholar] [CrossRef]
  24. Borka, D.; Mowbray, D.J.; Miškcovič, Z.L.; Petrovič, S.; Nescovič, N. Donut and dynamic polarization effects in proton channeling through nanotubes. New. J. Phys. 2010, 12, 43021. [Google Scholar] [CrossRef]
  25. Stolterfoht, N. Simulation and analysis of ion guiding through a nanocapillary in insulating polymers. Phys. Rev. A 2013, 87, 12902. [Google Scholar] [CrossRef]
  26. Stepanov, A.V.; Filippov, G.M. Channeling of low energy atomic particles in carbon nanotubes with heterojunctions. Nucl. Instr. Meth. Phys. Res. 2017, 402, 263. [Google Scholar] [CrossRef]
  27. Borka, D.; Borka, V.I. Channeling of protons through radial deformed double wall carbon nanotubes. Atoms 2019, 7, 88. [Google Scholar] [CrossRef] [Green Version]
  28. Ferrell, T.L.; Ritchie, R.H. Dynamical and geometrical effects on the physisorption of atoms. Phys. Rev. A 1980, 21, 1305. [Google Scholar] [CrossRef]
  29. Annett, J.F.; Echenique, P.M. VdW interaction between an atom and a surface at finite separations. Phys. Rev. B 1986, 34, 6853–6859. [Google Scholar] [CrossRef]
  30. Dedkov, G.V.; Kyasov, A.A. Dynamical atom-surface interaction. Surf. Sci. 2011, 605, 1077–1081. [Google Scholar] [CrossRef]
  31. Dedkov, G.V.; Kyasov, A.A. Fluctuation-electromagnetic interaction under dynamic and thermal nonequilibrium conditions. Physics-Uspekhi 2017, 60, 1–27. [Google Scholar] [CrossRef]
  32. Jackson, J.D. Classical Electrodynamics; Wiley: New York, NY, USA, 1975. [Google Scholar]
  33. Buhmann, S.Y. Dispersion Forces II. Many-Body Effects, Excited Atoms, Finite Temperature and QUANTUM Friction; Springer Tracts in Modern Physics Volume 248; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  34. Landau, L.D.; Lifshitz, E.M. Statistical Physics Part I; Pergamon: Oxford, UK, 1980. [Google Scholar]
  35. Lifshitz, E.M.; Pitaevskii, L.P. Statistical Physics Part II; Pergamon: Oxford, UK, 1980. [Google Scholar]
  36. Barash, Y.S. Van der Waals Forces; Nauka: Moscow, Russia, 1988. (In Russian) [Google Scholar]
  37. Pendry, J.B. Shearing the vacuum-quantum friction. J. Phys. C 1997, 9, 10301–10320. [Google Scholar] [CrossRef]
  38. Høye, J.S.; Brevik, I.; Milton, K.A. Casimir friction between polarizable particle and half-space with radiation damping at zero temperature. J. Phys. A 2015, 48, 365004. [Google Scholar] [CrossRef] [Green Version]
  39. Intravaia, F.; Behunin, R.O.; Henkel, C.; Busch, K.; Dalvit, D.A.R. Failure of local thermal equilibrium in quantum friction. Phys. Rev. Lett. 2016, 117, 100402. [Google Scholar] [CrossRef] [Green Version]
  40. Intravaia, F.; Oelschläger, M.; Reiche, D.; Dalvit, D.A.R.; Busch, K. Quantum rolling friction. Phys. Rev. Lett. 2019, 123, 120401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Reiche, D.; Busch, K.; Intravaia, F. Nonadditive enhancement of nonequilibrium atom-surface interactions. Phys. Rev. Lett. 2020, 124, 193603. [Google Scholar] [CrossRef] [PubMed]
  42. Reiche, D.; Dalvit, D.A.R.; Busch, K.; Intravaia, F. Spatial dispersion in atom-surface quantum friction. Phys. Rev. B 2017, 95, 155448. [Google Scholar] [CrossRef] [Green Version]
  43. Abramovitz, M.; Stegun, I.A. Handbook of Mathematical Functions; Dover: New York, NY, USA, 1972. [Google Scholar]
  44. Barash, Y.S.; Kyasov, A.A. Interaction potential for two filaments and for an atom interacting with a filament. Sov. Phys. JETP 1989, 68, 39–45. [Google Scholar]
  45. Gemmell, D.S. Channeling and related effects in the motion of charged particlesthrough crystals. Rev. Mod. Phys. 1974, 46, 129–236. [Google Scholar] [CrossRef]
  46. Milton, K.A.; Abalo, E.K.; Parashar, P.; Pourtolami, N.; Brevik, I.; Ellingsen, S.A. Repulsive Casimir and Casimir-Polder forces. J. Phys. A 2012, 45, 374006. [Google Scholar] [CrossRef]
Figure 1. Cylindrical configurations and coordinate systems: a particle outside a cylinder (a); inside cylindrical channel (b); between two coaxial cylinders (c). R is the position of the particle in the cross-section of the cylinder and V its velocity, always parallel to the cylinder axis.
Figure 1. Cylindrical configurations and coordinate systems: a particle outside a cylinder (a); inside cylindrical channel (b); between two coaxial cylinders (c). R is the position of the particle in the cross-section of the cylinder and V its velocity, always parallel to the cylinder axis.
Universe 07 00106 g001
Figure 2. Comparison of the dynamic U d and static U s vdW potentials (Equation (78) with (A16)–(A19)) for an atom moving parallel to a metal filament at a / R = 0.1 and γ / ω 0 = 0.005 .     The   red ,   blue   and   green   lines   correspond   to   ω p / ω 0 = 0.3 ,   1 ,   3 . The temperature T = 0 is assumed. Figure (b) is an extension of (a) to greater values of V / ω o R .
Figure 2. Comparison of the dynamic U d and static U s vdW potentials (Equation (78) with (A16)–(A19)) for an atom moving parallel to a metal filament at a / R = 0.1 and γ / ω 0 = 0.005 .     The   red ,   blue   and   green   lines   correspond   to   ω p / ω 0 = 0.3 ,   1 ,   3 . The temperature T = 0 is assumed. Figure (b) is an extension of (a) to greater values of V / ω o R .
Universe 07 00106 g002
Figure 3. Dynamic vdW potential (Equation (78) with (A16)–(A19)) for the Cs atom moving parallel to a metal chain of Au atoms as a function of separation R . Figure (b) is an extension of (a). The red and blue lines correspond to velocities of 3 · 10 7 and 1.5 · 10 7 m / s . The dotted lines show the continuous velocity-independent string potential in the Moliere approximation for the Cs—Au interaction potential [45]. The temperature T = 0 is assumed.
Figure 3. Dynamic vdW potential (Equation (78) with (A16)–(A19)) for the Cs atom moving parallel to a metal chain of Au atoms as a function of separation R . Figure (b) is an extension of (a). The red and blue lines correspond to velocities of 3 · 10 7 and 1.5 · 10 7 m / s . The dotted lines show the continuous velocity-independent string potential in the Moliere approximation for the Cs—Au interaction potential [45]. The temperature T = 0 is assumed.
Universe 07 00106 g003
Figure 4. Dynamic vdW potential (Equations (81) and (82) of a Cs atom with a chain of Si atoms, normalized to its absolute value for an atom at rest, plotted as a function of parameter V ω 0 / R at T = 0 . We used the values ω 01   =   3.89   eV and ω 02   =   8.15   eV for Cs and Si atoms, respectively, and interatomic spacing d   = 0.44   nm for the chain of Si atoms.
Figure 4. Dynamic vdW potential (Equations (81) and (82) of a Cs atom with a chain of Si atoms, normalized to its absolute value for an atom at rest, plotted as a function of parameter V ω 0 / R at T = 0 . We used the values ω 01   =   3.89   eV and ω 02   =   8.15   eV for Cs and Si atoms, respectively, and interatomic spacing d   = 0.44   nm for the chain of Si atoms.
Universe 07 00106 g004
Figure 5. Comparison of the dynamic and static vdW potentials for the Cs atom moving inside a capillary with Au metallized walls ( T = 0 ) . The red and blue lines correspond to Equations (58) and (59), green line and inset show the resulting potential. A fixed value a / R = 1.2 was used.
Figure 5. Comparison of the dynamic and static vdW potentials for the Cs atom moving inside a capillary with Au metallized walls ( T = 0 ) . The red and blue lines correspond to Equations (58) and (59), green line and inset show the resulting potential. A fixed value a / R = 1.2 was used.
Universe 07 00106 g005
Figure 6. The stopping force of a Cs atom moving inside a capillary of 5 nm radius with Au metallized walls, plotted as a function of the distance from the capillary axis. The red, blue and green lines correspond to velocities of 3 · 10 7 ,   1.5 · 10 7 and 6 · 10 6   m / s . The temperature T = 0 is assumed.
Figure 6. The stopping force of a Cs atom moving inside a capillary of 5 nm radius with Au metallized walls, plotted as a function of the distance from the capillary axis. The red, blue and green lines correspond to velocities of 3 · 10 7 ,   1.5 · 10 7 and 6 · 10 6   m / s . The temperature T = 0 is assumed.
Universe 07 00106 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dedkov, G.V. Van der Waals Interactions of Moving Particles with Surfaces of Cylindrical Geometry. Universe 2021, 7, 106. https://doi.org/10.3390/universe7040106

AMA Style

Dedkov GV. Van der Waals Interactions of Moving Particles with Surfaces of Cylindrical Geometry. Universe. 2021; 7(4):106. https://doi.org/10.3390/universe7040106

Chicago/Turabian Style

Dedkov, George V. 2021. "Van der Waals Interactions of Moving Particles with Surfaces of Cylindrical Geometry" Universe 7, no. 4: 106. https://doi.org/10.3390/universe7040106

APA Style

Dedkov, G. V. (2021). Van der Waals Interactions of Moving Particles with Surfaces of Cylindrical Geometry. Universe, 7(4), 106. https://doi.org/10.3390/universe7040106

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop