Induced vacuum magnetic flux in quantum spinor matter in the background of a topological defect in two-dimensional space

A topological defect in the form of the Abrikosov-Nielsen-Olesen vortex is considered as a gauge-flux-carrying tube that is impenetrable for quantum matter. The relativistic spinor matter field is quantized in the vortex background in $2+1$-dimensional conical space-time which is a section orthogonal to the vortex axis; the most general set of boundary conditions ensuring the impenetrability of the vortex core is employed. We find the induced vacuum current circulating around the vortex and the induced vacuum magnetic field strength pointing along the vortex axis. The requirement of finiteness and physical plausibility for the total induced vacuum magnetic flux allows us to restrict the variety of admissible boundary conditions. The dependence of the results on the transverse size of the vortex, as well as on the vortex flux and the parameter of conicity, is elucidated.


Introduction
Spontaneous breakdown of continuous symmetries can give rise to topological defects with rather interesting properties. A topological defect in three-dimensional space, which is characterized by the nontrivial first homotopy group, is known as the Abrikosov-Nielsen-Olesen (ANO) vortex [1,2]. The vortex is described classically in terms of a spin-zero (Higgs) field that condenses and a spin-one field corresponding to the spontaneously broken gauge group; the former field is coupled to the latter one in the minimal way with constantẽ cond . Single valuedness of the condensate field and finiteness of the vortex energy implement that the vortex flux is related toẽ cond , Φ = dxV(x) = 2π/ẽ cond , (1.1) where V(x) is the vector potential of the spin-one gauge field, and the integral is over a path enclosing the vortex tube once (natural units = c = 1 are used). As some amount of energy (mass) is stored in the core of a topological defect, this core is a source of gravity. Such a source in the case of the linear ANO vortex makes the spatial region outside the vortex core to be conical, i.e., with the deficit angle equal to 8πĢ M: the squared length element in the outer region is Ģ is the gravitational constant, and M is the linear density of mass stored in the core. The transverse size of the vortex core is determined by the correlation length, and quantity M is of order of the inverse correlation length squared. Since constant Ģ is of the order of the Planck length squared, the effects of conicity, which are characterized by the value of the deficit angle, are negligible for vortices in ordinary superconductors. However topological defects of the type of ANO vortices may arise in a field that is seemingly rather different from condensed matter physics -in cosmology and high-energy physics. This was realized by Kibble [3,4] and Vilenkin [5,6] (see also [7]), and, from the beginning of the 1980s, such topological defects are known under the name of cosmic strings. Cosmic strings with the thickness of the order of the Planck length are definitely ruled out by astrophysical observations, but there remains a room for cosmic strings with the thickness that is more than 3 orders larger than the Planck length (see, e.g., [8]), although the direct evidence for their existence is lacking. A recent development in material science also provides an unexpected link between condensed matter and high-energy physics, which is caused to a large extent by the experimental discovery of graphene -a two-dimensional crystalline allotrope formed by a monolayer of carbon atoms [9]. Low-energy electronic excitations in graphene are characterized by dispersion law which is the same as that for Dirac fermions in relativistic field theory, with the only distinction that the velocity of light is changed to the Fermi velocity, see [10,11]. It is well established by now that a sheet of graphene is always corrugated and covered by ripples which can be either intrinsic or induced by roughness of a substrate. A single topological defect (disclination) warps a sheet of graphene, rolling it into a nanocone which is similar to the transverse section of a spatial region out of a cosmic string; carbon nanocones with deficit angles equal to N d π/3 (N d = 1, 2, 3, 4, 5, i.e. ν = 6 5 , 3 2 , 2, 3, 6) were observed experimentally, see [12,13] and references therein. Moreover, theory also predicts saddle-like nanocones with the deficit angle taking negative values unbounded from below, N d = −1, −2, −3, ..., −∞, i.e. ν = 6 7 , 3 4 , 2 3 , ..., 0, which can be regarded as corresponding to cosmic strings with negative mass density. Note that nanoconical structures may arise as well in a diverse set of condensed matter systems known as the two-dimensional Dirac materials, ranging from honeycomb crystalline allotropes (silicene and germanene [14], phosphorene [15]) to high-temperature cuprate superconductors [16] and topological insulators [17].
While considering the effect of the ANO vortex on the vacuum of quantum matter, the following two circumstances should be kept in mind. First, the phase with broken symmetry exists outside the vortex core and the vacuum is to be defined there. Hence, the quantum matter field does not penetrate inside the core, obeying a boundary condition at its side edge. Second, the impact of the ANO vortex on quantum matter is through a vector potential of the vortex-forming spin-one field, and the quantum matter field is assumed to couple to this vector potential in the minimal way with coupling constantẽ. Hence, the ANO vortex flux has no effect on the surrounding matter in the framework of classical theory, and such an effect, if it exists, is of a purely quantum nature. This phenomenon should be understood as a quantum-field-theoretical manifestation of the famous Aharonov-Bohm effect [18] and is characterized by the periodic dependence on the value of the vortex flux, Φ (1.1), with the period equal to London flux quantum 2π/ẽ.
A crucial task in the study of the effect of the ANO vortex on the vacuum of quantum matter is to elucidate the dependence on a boundary condition at the edge of the vortex core. It seems reasonable to start from the most general set of mathematically admissible boundary conditions and then, after obtaining the outcoming effect, to restrict this set by physically motivated arguments. Another task is to elucidate the dependence on the transverse size of the vortex core. These two tasks will be thoroughly scrutinized and solved in the course of the present study by considering a somewhat simplified case of two-dimensional space 1 being the transverse section of a three-dimensional spatial region out of the ANO vortex.
It should be noted that the current, the condensate, and the energy-momentum tensor that are induced in the vacuum of quantum relativistic spinor matter were considered in the above-described context in [20,21,22]. However, a particular boundary condition was employed, and the issue of a dependence of the results on the choice of admissible boundary conditions remained undisclosed.
The current and the magnetic field strength, as well as the energy density and the Casimir force, that are induced in the vacuum of quantum relativistic scalar matter at ν = 1 in a space of arbitrary dimension were considered in the above-described context in [23,24,25,26]. In these studies the Dirichlet boundary condition was employed; a physical motivation herein is in the assumption of a perfect reflection of quantum matter from the edge of the vortex core.
In the case of quantum relativistic spinor matter, neither the Dirichlet nor the Neumann boundary condition is admissible. A physically motivated demand is the absence of the matter flux across the boundary. In 2 + 1-dimensional spacetime with a connected boundary, this demand yields a one-parameter family of boundary conditions; see Section 4 below. Employing such boundary conditions, we shall find the induced vacuum current and the induced vacuum magnetic field strength; further physical arguments will be shown to remove an ambiguity in the choice of boundary conditions.
In the next section we define the current and the magnetic field that are induced in the vacuum of quantum relativistic spinor matter in the background of the ANO vortex of nonzero transverse size. In Section 3 we present the complete set of solutions to the Dirac equation that is relevant to the problem considered. In Section 4 we choose boundary conditions ensuring the absence of the matter flux across the edge of the vortex core. The induced vacuum current is obtained in Section 5. In Section 6 we consider the induced vacuum magnetic field and its total flux with the use of both analytical and numerical methods. Finally, the results are summarized and discussed in Section 7. Some details in the derivation of the expression for the induced vacuum current are given in Appendices A and B. The case of the infinitely thin vortex is reviewed in Appendix C. The results for massless quantum spinor matter are presented in Appendix D.
2 Preliminaries and definition of physical characteristics of the vacuum The operator of the second-quantized spinor field is presented as where a † E and a E (b † E and b E ) are the spinor particle (antiparticle) creation and destruction operators obeying the anticommutation relations, ψ E (x) is the solution to the stationary Dirac equation, and symbol denotes summation over the discrete part and integration (with a certain measure) over the continuous part of the energy spectrum; ground state |vac is conventionally defined by relation In the case of the ANO vortex background, the Dirac hamiltonian takes form where ω is the spin connection corresponding to conical space (1.2). The current that is induced in the vacuum is given by expression (sgn(u) is the sign function, sgn(u) = ±1 at u ≷ 0). The magnetic field strength, B I (x), is also induced in the vacuum, as a consequence of the Maxwell equation, where the electromagnetic coupling constant, e, differs in general fromẽ. The total flux of the induced vacuum magnetic field is Since the vacuum of quantum matter exists outside the ANO vortex core, as was already emphasized, an issue of the choice of boundary conditions at the edge of the core is of primary concern. Turning to this issue, let us note first, that (2.4) is not enough to define the hamiltonian operator rigorously in a mathematical sense. To define an operator in an unambiguous way, one has to specify its domain of definition. Let the set of functions ψ be the domain of definition of operator H and the set of functionsψ be the domain of definition of its adjoint, operator H † . Then the operator is Hermitian (or symmetric in mathematical parlance), is valid; here functions ψ(x) andψ(x) are defined in space X with boundary ∂X. It is evident that condition (2.9) can be satisfied by imposing different boundary conditions for ψ andψ. But, a nontrivial task is to find a possibility that a boundary condition forψ is the same as that for ψ; then the domain of definition of H † coincides with that of H, and operator H is self-adjoint (for a review of the Weyl-von Neumann theory of self-adjoint operators see [27,28] with j(x) given by (2.5). If boundary ∂X is connected, then (2.11) is reduced to n · j| x∈∂X = 0, (2.12) where n is the unit normal that may be chosen as pointing outward to X. The last relation signifies the impenetrability of ∂X; i.e., the matter field is confined to X.
In the present paper we consider the vacuum polarization effects in 2 + 1-dimensional spacetime, which is a section orthogonal to the ANO vortex axis, i.e., at z = . The irreducible representation of the Clifford algebra is chosen in such a way that the Dirac matrices in flat 2 + 1-dimensional spacetime take form where σ 1 , σ 2 , and σ 3 are the Pauli matrices (a transition to another inequivalent irreducible representation can be made by changing the sign of β). In the background of the ANO vortex, the only one component of the vector potential and the spin connection is nonvanishing: 14) and the Dirac hamiltonian takes form (Z is the set of integer numbers), we present the Dirac equation as a system of two first-order differential equations for radial functions: [|u|] is the integer part of quantity u (i.e., the integer that is less than or equal to u), and { |u| } = u − [|u|] is the fractional part of quantity u, 0 ≤ { |u| } < 1. Using (2.16) and (2.17), one gets j r = 0, and the only component of the induced vacuum current, is independent of the angular variable. The induced vacuum magnetic field strength is directed along the vortex axis, with total flux where it is assumed without a loss of generality that the vortex core has the form of a tube of radius r 0 .

Solution to the Dirac equation
The solution to the system of equations, (2.18), is given in terms of cylindrical functions.
Let us define where l = n − n c , and where l ′ = −n + n c ; here J ρ (u) and Y ρ (u) are the Bessel and Neumann functions of order ρ, k = √ E 2 − m 2 . In the case of ν > 1 and 0 < F < 1 2 1 − 1 ν 1 2 (1 − ν) < G < 0 , the complete set of solutions to (2.18) is given by In the case of ν > 1 and 1 2 (1 + 1 ν ) < F < 1 1 < G < 1 2 (1 + ν) , the complete set of solutions to (2.18) is given by One can note that both upper and lower components of each mode consist of two terms: one (given by the Bessel function) is vanishing and another one (given by the Neumann function) is diverging in the limit of r → 0.
In the case of ν ≥ 1 and 1 , there is a peculiar mode corresponding to n = n c . This mode can be composed either from the pair of columns or from the pair of columns both terms in the first variant have divergent upper components, whereas both terms in the second variant have divergent lower components. Instead of these variants we choose the following form: at r → 0; they are irregular but square integrable. The latter circumstance requires a quest for a self-adjoint extension, and the Weyl-von Neumann theory of deficiency indices (see [27,28]) has to be employed. In the case of ν ≥ 1 and 1 2 (1 − 1 ν ) < F < 1 2 (1 + 1 ν ), as well as in the case of 1 2 ≤ ν < 1 and 1 2 1 ν − 1 < F < 1 2 3 − 1 ν , when there is one irregular mode, the deficiency index is (1,1), and the one-parametric family of self-adjoint extensions can be introduced with the use of condition where Θ is the self-adjoint extension parameter [29,30]. In view of relations the modes are orthonormalized as the modes corresponding to the continuous spectrum: In addition, there is a bound state at cos Θ < 0 with energy E BS in the gap between the continuums, −m < E BS < m. Its mode is , (3.15) and the value of its energy is determined from relation Γ(u) is the Euler gamma function, and K ρ (u) is the Macdonald function of order ρ.

Self-adjointness and choice of boundary conditions
The Dirac hamiltonian operator in the background of the ANO vortex of nonzero radius is valid, see (2.8) and (2.9), and sets of functions ψ andψ coincide. Ergo, the quest is for a boundary condition in the form where K is a matrix (element of the Clifford algebra) which without a loss of generality can be chosen to be Hermitian and has to obey conditions [K, α r ] + = 0 (4.3) and One can simply go through four linearly independent elements of the Clifford algebra in 2 + 1-dimensional spacetime which is a section orthogonal to the ANO vortex axis, and find with real coefficients satisfying c 2 1 + c 2 2 = 1. Using obvious parametrization we finally obtain K = iβα r e −iθα r . Thus, boundary condition (4.2) with K given by (4.7) is the most general boundary condition ensuring the self-adjointness of the Dirac hamiltonian operator in the background of the ANO vortex of nonzero radius r 0 in transverse section z = const, and parameter θ can be interpreted as the self-adjoint extension parameter. Value θ = 0 corresponds to the MIT bag boundary condition, which was proposed long ago as the condition ensuring the confinement of the matter field [40]. However, it should be comprehended that a condition with an arbitrary value of θ ensures the confinement equally as well as that with θ = 0. Imposing boundary condition (4.2) with matrix K (4.7) on the solution to the Dirac equation, ψ E (x) (2.17), we obtain the condition for the modes: which allows us to determine their coefficients: tan(µ (∨) ρ ) = cos θ 2 + π 4 (m + E)Y ρ (kr 0 ) − sin θ 2 + π 4 kY ρ−1 (kr 0 ) − cos θ 2 + π 4 (m + E)J ρ (kr 0 ) + sin θ 2 + π 4 kJ ρ−1 (kr 0 ) , (4.10) tan(µ ρ ) = cos θ 2 + π 4 kJ 1−ρ (kr 0 ) + sin θ 2 + π 4 (m − E)J −ρ (kr 0 ) − cos θ 2 + π 4 kJ ρ−1 (kr 0 ) + sin θ 2 + π 4 (m − E)J ρ (kr 0 ) .
Comparing the case of a vortex of nonzero transverse size with that of an infinitely thin one, we conclude that in the first case the total hamiltonian is extended with the use of the only one self-adjoint extension parameter, whereas in the second case several partial hamiltonians are extended, and the number of self-adjoint extension parameters can be zero (no need for extension, the operator is essentially self-adjoint), one, four, etc. The values of the self-adjoint extension parameters in the second case can be fixed from the first case by limiting procedure r 0 → 0 [41]. The nonpeculiar modes (ρ > 1) in this limit become regular and independent of θ, since, as was already noted, The peculiar modes (0 < ρ < 1) in this limit become irregular and square integrable, and , when there is one peculiar mode. Inserting the contribution of the appropriate modes (see (3.1), (3.2), (3.5), (3.6), and (4.12)) to (2.20), we obtain .

(5.4)
In Appendix A, the summation in (5.2) is performed, yielding where I ρ (u) is the modified Bessel function of order ρ (p and N are the positive integers, δ ω, ω ′ is the Kronecker symbol, δ ω, ω ′ = 0 at ω ′ = ω and δ ω, ω = 1), while (5.3) is transformed to the following expression: and In Appendix B, j ϕ (5.5) is transformed to the following expression: Summing (5.5), (5.6) and (5.9), we obtain the final form for the induced vacuum current and express it in terms of F instead of G (see (2.19)), The analysis in Appendix A is sufficient to consider cases when there are no peculiar modes. In the case of ν > 1 and 0 < F < 1 In the case of ν > 1 and 1 Note that both (5.13) and (5.14) consist of two parts: one (with a factor of m/(2π) 2 ) is independent of r 0 , and another one (with a factor of r/π 2 ) is vanishing in the limit of It is instructive to present result (5.11) in the same way; evidently, j (a) ϕ (r) then coincides with the current that is obtained by imposing the condition of minimal irregularity in the case of an infinitely thin vortex [31,32], see (3.11) and (4.16). We obtain for the decomposition of (5.11) according to (5.15): and and the use is made of relations It should be noted that the r 0 -independent part of the current in cases 1 whereas it depends on θ if θ = ±π/2, see (5.17). The latter is distinct from cases of the absence of peculiar modes, when the r 0 -independent part of the current is always independent of θ, see the first four lines in (5.13) and (5.14). Note also that limits F → 1/2 and r 0 → 0 in general do not commute. Indeed, we obtain a discontinuity at F = 1/2, if limit r 0 → 0 is taken first, where the use is made of relation (5.27) When the order of limits is reversed, then j which follows from particular cases of (5.24) and (5.25), Adding to (5.28) and taking limit r 0 → 0, we get The limits do commute in special cases only: ) the discontinuity at F = 1/2 is absent in these cases.

Induced vacuum magnetic field and its flux
In the case of ν ≥ 1 and 1 , we obtain the following expression for B I (2.21): where Σ(w, v) is given by (5.12), while C (v) is given by (5.23). In the case of ν > 1 and 0 < F < 1 2 1 − 1 ν , the induced vacuum magnetic field is given by (6.2) with B (a) I given by the right-hand side of (6.3) and B (b) I given by the right-hand side of (6.6). In the case of ν > 1 and 1 2 1 + 1 ν < F < 1, the induced vacuum magnetic field is given by (6.2) with B (a) I given by the right-hand side of (6.4) and B (b) I given by the right-hand side of (6.7).
Turning to the total flux of the induced vacuum magnetic field, see (2.22), we present it as Φ I = Φ We obtain in the case of ν ≥ 1 and 1 2 1 − 1 ν < F < 1 2 1 + 1 ν , or 1 2 ≤ ν < 1 and and In the case of ν > 1 and 0 < F < 1 is given by the right-hand side of (6.12). In the case of ν > 1 and 1 is given by the right-hand side of (6.13). As follows from (6.12)-(6.14), Φ   18) in Appendix A, contour C 0 in the complex z plane is depicted in Fig. 8. Consequently, we get B (a) and Only a simple pole at z = 0 of the integrands contributes to the integrals in (6.20) and (6.21). Calculation of its residue yields Considering the case of 1 2 < ν ≤ 1 and 1 2 1 ν − 1 < F < 1 2 3 − 1 ν at F = 1/2, we obtain that are given by the right-hand sides of (6.22) and (6.23), respectively. Note that lim r 0 →0 Φ (a) I θ =±π/2 is discontinuous at F = 1/2 and its limiting values are independent of ν: ; (6.24) this is clearly a consequence of (5.26).
As for the remaining part of the total flux, Φ I , it can be presented as provided that the following condition holds: The case of F = 1/2 needs a special comment, since, due to relation (5.31), the current in this case is odd function of θ. Whereas the current and, consequently, the induced magnetic field with its flux vanish at θ = 0, they can be nonvanishing if discontinuous at θ = π. Indeed, periodicity in θ with period 2π, together with oddness in θ, Namely this is obtained from the appropriate formulas at F = 1/2 and θ = π ± : As a consequence, we obtain and in particular, cf. (6.15), lim In the case of F = 1/2 continuity in θ is maintained, and the induced vacuum current in this case takes form Consequently, we obtain the following expressions for the induced vacuum magnetic field, (qr ′ ) , F = 1/2, (6.36) and its flux, (qr) , F = 1/2. (6.37) With the use of relations (see [42]) and the Schläfli contour integral representation for the Macdonald function, 2 K ρ (w) = 1 4i sin(ρπ) C 0 dz e w cosh z+ρz , the integration over r can be performed. As a result, we obtain the following representation for the induced vacuum magnetic flux (6.39) in particular, , is positive at F > 1/2 and negative at F < 1/2; its absolute value increases with the increase of ν.
is positive and e −1 m Φ I | F = 1 2 , θ=π − is negative, being of the same absolute value that is independent of ν, see (6.32). Continuity of the results in θ is broken at θ = π and F = 1/2 only.
A more detailed analysis of the behavior of the induced vacuum magnetic flux can be obtained with the use of numerical computations. Taking F = 0.7 and ν = 1, we calculate the dimensionless flux, e −1 m Φ I | θ= π 2 ∓ π 2 , as a function of mr 0 , see Fig. 2. In the case of θ = 0, this function decreases with the increase of mr 0 , becoming vanishingly small ( 10 −7 ) at mr 0 ≥ 1. On the contrary, in the case of θ = π, this function increases at no allowance with the increase of mr 0 .
The dimensionless flux in the case of θ = 0 at several values of ν, as well as of mr 0 , is presented as a function of F in Fig. 3. As mr 0 increases, the absolute value of this function decreases as compared to the value at mr 0 = 0, becoming negligible in the vicinity of F = 1/2. However, the vicinity is shrinked as ν increases (this is also demonstrated by Fig. 4), and the flux at mr 0 ≥ 1 can equal its value at mr 0 = 0 for sufficiently large values of ν, unless F = 1/2.
The dimensionless flux in the case of θ = π at several values of ν, as well as of mr 0 , is presented as a function of F in Fig. 5. In the case of 1 2 < ν ≤ 2 the absolute value of this function increases with the increase of mr 0 , see Figs. 5a and 5b. However, in the case of ν > 2, the increase takes place in the vicinity of F = 1/2 and, otherwise, there is a decrease, see Figs. 5c and 5d (this is also demonstrated by Fig. 6a). Note that the flux at large values of mr 0 fails to depend on ν (lines corresponding to different values of ν merge together) at least in the case of 1 2 < ν ≤ 4, see Fig. 6b.

Summary and discussion
In the present paper, we have studied the impact of boundary conditions at the edge of the ANO vortex core on the vacuum polarization effects in quantum relativistic spinor matter in two-dimensional space. The most general boundary condition that is compatible with the self-adjointness of the energy operator in first-quantized theory, Dirac hamiltonian (2.15), is (see (4.2) and (4.7)) where θ is the self-adjoint extension parameter. This condition is also the most general one ensuring the impenetrability of the vortex core edge, i.e., the confinement of the matter field to the region out of the vortex core. We find that a current circulating in the vacuum around the vortex is given by expression (5.11) in the case of ν ≥ 1 and 1 2 1 − 1 ν < F < 1 2 1 + 1 ν , or 1 2 ≤ ν < 1 and 1 2 1 ν − 1 < F < 1 2 3 − 1 ν ; it is given by expression (5.13) in the case of ν > 1 and 0 < F < 1 2 1 − 1 ν and by expression (5.14) in the case of ν > 1 and 1 2 1 + 1 ν < F < 1. At large distances from the vortex, r → ∞, the b) a) n=3 n=4 n = 1 Figure 6: The dimensionless induced flux at θ = π as a function of F : a) mr 0 = 0 (solid lines) and mr 0 = 10 −1 (dashed lines), b) mr 0 = 0 (solid lines) and mr 0 = 5 (dashed line).

current decreases exponentially as
whereas it decreases as 1/r in the case of massless spinor matter; see Appendix D.
As a manifestation of the Aharonov-Bohm effect, the current is periodic in the value of the vortex flux, Φ; i.e., it depends on F and not on n c , see (2.19); moreover, it changes sign under simultaneous changes F → 1 − F and θ → −θ, see (5.31). One can introduce the charge conjugation transformation with the vortex flux changing its sign, The boundary condition for a conjugated wave function differs from (7.1): By requiring invariance of the boundary condition under such a charge conjugation, one restricts the values of the self-adjoint extension parameter to θ = π 2 ∓ π 2 . In the latter case the current changes sign under change F → 1 − F ; i.e., it is odd, as well as periodic in the value of the vortex flux. Consequently, it vanishes at θ = 0 and F = 1/2, while it is discontinuous in θ at θ = π and F = 1/2, see (6.30).
It is appropriate here to discuss the dependence on the transverse size of the vortex core and the limiting procedure as this size tends to zero, r 0 → 0. For this task it is instructive to decompose the current into r 0 -independent, j (a) ϕ , and r 0 -dependent, j (b) ϕ , pieces, see (5.15) -(5.22); the r 0 -dependent piece vanishes at r 0 → 0. It should be noted that, in the case of the infinitely thin (r 0 = 0) vortex, the Dirac hamiltonian is essentially self-adjoint for ν > 1 and either 0 < F < 1 2 1 − 1 ν or 1 2 1 + 1 ν < F < 1; otherwise, there emerges the self-adjoint extension problem with one, or four, or more parameters. One self-adjoint extension parameter, Θ, appears for ν ≥ 1 and 1 2 1 − 1 ν < F < 1 2 1 + 1 ν , or for 1 2 ≤ ν < 1 and 1 The results for ν ≥ 1 and 0 < F < 1, as well as for 1 2 ≤ ν < 1 and 1 2 1 ν − 1 < F < 1 2 3 − 1 ν , are comprehensively presented in Appendix C. The value of Θ can be fixed by limiting procedure r 0 → 0 starting from the r 0 > 0 case. In this way, the condition of minimal irregularity in the r 0 = 0 case is obtained in the form of (4.16). If this condition is supplemented with the requirements of invariance under charge conjugation (7.2) and continuity in θ, then it takes the form of (4.16) with θ = 0 at F = 1/2; namely in the latter form it was first proposed in [31,32].
As a consequence of the Maxwell equation, the magnetic field strength is also induced in the vacuum, pointing along the vortex axis; the relevant expressions in the case of the most general boundary condition, (7.1), are given by (6.1)-(6.9). This allows us to consider the total magnetic flux which is induced in the vacuum. As follows from our numerical analysis, the latter is finite at θ = π 2 ∓ π 2 only, but otherwise, it is divergent. Thus, the physical condition that the induced vacuum magnetic flux be finite corresponds to the requirement of invariance under charge conjugation (7.2). The flux for boundary conditions maintaining the charge conjugation invariance is given at F = 1/2 by expression (6.38), it vanishes at θ = 0 and F = 1/2, and it is nonvanishing and discontinuous in θ at θ = π and F = 1/2, see (6.32). The flux in the case of r 0 = 0 is discontinuous at F = 1/2; moreover, its absolute value at θ = π and F = 1/2, see (6.33), is twice as large as its absolute value at θ = π in the limit F → 1/2, see (6.24). The case of an infinitely thin vortex is an idealization that is inappropriate to physics reality, since, as has already been noted in Introduction, the transverse size of the vortex, r 0 , is of the order of the correlation length. In the case of r 0 > 0, the differences in behavior of the flux at θ = 0 and at θ = π are comprehensively illustrated by Figs. 2-6. Whereas the flux at θ = 0 decreases in its absolute value with the increase of r 0 , the flux at θ = π in general increases at no allowance in its absolute value with the increase of r 0 (although there is a moderate decrease in vicinities of F = 0 and F = 1 at ν > 2). Such a behavior of the flux, as that at θ = π, can hardly be regarded as physical. Quantity r −1 0 is identified with the energy scale of spontaneous symmetry breaking, i.e., the mass of the corresponding Higgs particle. It looks rather unlikely that a topological defect influences the surrounding quantum matter with the matter particle mass, m, exceeding the Higgs particle mass, m cond ; the more unlikely is the unrestricted growth of this influence with the increase of quotient m/m cond . The influence of a topological defect on the surrounding quantum matter at m cond > m looks much more physically plausible. Namely this situation is realized in the case of quantum scalar matter obeying the Dirichlet boundary condition at the vortex edge, see [23,24,25,26].
We conclude that, although we have solved the problem analytically with the use of the most general set of boundary conditions ensuring the impenetrability of the vortex core, the analysis of the behavior of the induced vacuum magnetic flux allows us to restrict the variety of admissible boundary conditions. The requirement of the flux to be finite, which is equivalent to the requirement of invariance under charge conjugation (7.2), restricts the values of the self-adjoint extension parameter to θ = π 2 ∓ π 2 . The further requirement of physical plausibility of the finite flux behavior, which is equivalent to the formal requirement of continuity in the dependence on θ, yields θ = 0. As long as the transverse size of the vortex is taken into account, the induced vacuum effects at θ = 0 are continuous in F and vanishing at F = 1/2. At this point we would like to emphasize the crucial distinction between the cases of massive and massless quantum spinor matter. The latter case requires an introduction of the maximal size of the system, r max . We discover that, for r max ≫ r 0 (in conformance to the reality), the induced vacuum effects for both θ = 0 and θ = π are physically plausible; moreover, they coincide, being independent of the transverse size of the vortex, see (D.15) and (D. 16) in Appendix D. Because of this distinction, the results in the massless case are discontinuous at F = 1/2 with a jump which is independent of ν (see (D.19)). Note in this respect that the current and the magnetic field that are induced in the vacuum of quantum scalar matter under the Dirichlet boundary condition are continuous in F and vanishing at F = 1/2 even in the r 0 = 0 case, see [43,44,45]. In contrast to this, the emergence of a peculiar mode in the solution to the Dirac equation (see Section 3) leads to a discontinuity of the results at F = 1/2. For massive quantum spinor matter, the discontinuity is present in the r 0 = 0 case, but is eliminated by a choice of the physically motivated boundary condition in the r 0 > 0 case. The discontinuity nonetheless stays on for massless quantum spinor matter. It would be interesting to perform a similar study for other characteristics of the vacuum, for instance, for the induced vacuum energy-momentum tensor. and and it is implied that µ (∧) ρ and µ (∨) ρ , determined by relations (4.9) and (4.10), depend on √ k 2 instead of k. The integral over real k can be transformed into the integral over a contour circumventing anticlockwise the positive imaginary semiaxis in the complex k plane. It is evident in the case of j (1) ϕ that the latter contour is reduced to a contour circumventing a part of the positive imaginary semiaxis, see Fig. 7. As a result, we get which can be decomposed as and Using the Schläfli contour integral representation, where contours C + , C − , and C 0 in the complex z plane are shown in Fig. 8. The vertical segments of contours C + , C − and C 0 are infinitesimally close to the imaginary axis, not coinciding with it in order to avoid simple poles of the integrand at z = 0 and z = ±2ipπ/ν (p is the positive integer, 1 ≤ p ≤ [|ν/2|]). Contour C 0 circumvents poles out of the origin (existing at ν ≥ 2), whereas the contribution of the pole at the origin (existing at ν = 0) is explicitly separated in (A.11). Substituting (A.11) into (A.9), we change integration variable y → v = y(qr) −2 and take a derivative, . . . Integrating over q and v, we get 1 + 1 2mr sin(pπ/ν) e −2mr sin(pπ/ν) cos(2Gpπ/ν) where the finite sum over integer p and the last term with the Kronecker δ symbol (N is the positive integer) are due to a contribution of simple poles on the imaginary axis out of the origin.
In a similar way we calculate the sum in (A.10): . (A.14) Substituting (A.14) into (A.10) and integrating over q and y, we get As a result, we obtain the following expression for j (1) ϕ (A.7): we can present (A.16) as given by (A.13), (A.15), and (A.18), respectively. Using the latter relation, we finally obtain (5.5).
Turning now to j  As a consequence of (A.5) and (A.26), we obtain (5.6).

Appendix B: contribution of peculiar modes to the current
Similarly to that in the beginning of Appendix A, one can obtain With the help of (A.1) and (B.1), j ϕ (5.4) is presented as and it is implied that µ 1−G , as determined by (4.11), depends on √ k 2 instead of k. The integral over real k can be transformed into the integral over a contour circumventing anticlockwise the positive imaginary semiaxis in the complex k plane. The latter contour is reduced to a contour consisting of two parts: one encircling a simple pole emerging at cos θ < 0 and another one circumventing a part of the positive imaginary semiaxis, see Fig. 9. The contribution of the pole cancels out the last term in (B.2), and for the remaining part we get .

Appendix D: case of massless quantum spinor matter
We present here the results for the case of massless quantum spinor matter in the background of the ANO vortex of nonzero transverse size.