Search

Two approaches are developed for the study of the bound states of a one-dimensional Dirac equation with the potential consisting of N δ-function centers. One of these uses Green’s function method. This method is applicable to a finite number N of δ-point centers, reducing the bound state problem to finding the energy eigenvalues from the determinant of a 2 N × 2 N matrix. The second approach starts with the matrix for a single delta-center that connects the two-sided boundary conditions for this center. This connection matrix is obtained from the squeezing limit of a piecewise constant approximation of the delta-function. Having then the connection matrices for each center, the transmission matrix for the whole system is obtained by multiplying the one-center connection matrices and the free transfer matrices between neighbor centers. An equation for bound state energies is derived in terms of the elements of the total transfer matrix. Within both approaches, the transcendental equations for bound state energies are derived, the solutions to which depend on the strength of delta-centers and the distance between them, and this dependence is illustrated by numerical calculations. The bound state energies for the potentials composed of one, two, and three delta-centers ( N = 1, 2, 3) are computed explicitly. The principle of strength additivity is analyzed in the limits as the delta-centers merge at a single point or diverge to infinity.

In the framework of Higgs inflation, we consider the electromagnetic field nonminimally coupled to gravity via the parity-preserving $\propto RF^{2}$ and parity-violating $\propto RF\tilde{F}$ terms. Using the perturbation theory to the leading order in these couplings, we study the generation of the electromagnetic field during the inflation stage. We derive the effective kinetic and axial coupling functions arising in the Einstein frame in the usual metric as well as Palatini formulations of gravity. For both formulations, we determine the power spectrum, energy density, and helicality of the generated electromagnetic fields for different values of the nonminimal coupling constants. Finally, we estimate the maximal present-day magnitude of the magnetic field as $10^{-14}-10^{-15}\,$G with the correlation length of order 10 pc., 23 pages, 10 figures

We study the explosive production of gauge fields during axion inflation in a novel gradient expansion formalism that describes the time evolution of a set of bilinear electromagnetic functions in position space. Based on this formalism, we are able to simultaneously account for two important effects that have thus far been mostly treated in isolation: (i) the backreaction of the produced gauge fields on the evolution of the inflaton field and (ii) the Schwinger pair production of charged particles in the strong gauge-field background. This allows us to show that the suppression of the gauge-field production due to the Schwinger effect can prevent the backreaction in scenarios in which it would otherwise be relevant. Moreover, we point out that the induced current, $\boldsymbol{J} = \sigma \boldsymbol{E}$, also dampens the Bunch-Davies vacuum fluctuations deep inside the Hubble horizon. We describe this suppression by a new parameter $\Delta$ that is related to the time integral over the conductivity $\sigma$ which hence renders the description of the entire system inherently nonlocal in time. Finally, we demonstrate how our formalism can be used to construct highly accurate solutions for the mode functions of the gauge field in Fourier space, which serves as a starting point for a wealth of further phenomenological applications, including the phenomenology of primordial perturbations and baryogenesis., Comment: 24 pages, 7 figures

In order to describe magnetogenesis during inflation in the kinetic coupling model, we utilize a gradient expansion which is based on the fact that only long-wavelength (superhorizon) modes undergo amplification. For this purpose, we introduce a set of functions (bilinear combinations of electromagnetic fields with an arbitrary number of spatial curls) satisfying an infinite chain of equations. Apart from the usual mode enhancement due to interaction with the inflaton, these equations also take into account the fact that the number of relevant modes constantly grows during inflation. Truncating this chain, we show that even with a relatively small number of equations, it is possible to describe the electric and magnetic energy densities with a few percent accuracy during the whole inflation stage. We arrive at this conclusion for different types of coupling functions (increasing, decreasing, and nonmonotonic) in the regime with strong backreaction and its absence., 22 pages, 14 figures

We investigate collisional decay of the axial charge in an electron-photon plasma at temperatures 10 MeV - 100 GeV. We demonstrate that the decay rate of the axial charge is first order in the fine-structure constant $\Gamma_{\rm flip}\propto \alpha m_{e}^{2}/T$ and thus orders of magnitude greater than the naive estimate which has been in use for decades. This counterintuitive result arises through infrared divergences regularized at high temperature by environmental effects. The decay of axial charge plays an important role in the problems of leptogenesis and cosmic magnetogenesis., Comment: 7 pages, 3 figures + Supplementary Material (3 pages)

By using the first-principles approach, we derive a system of three quantum kinetic equations governing the production and evolution of charged scalar particles by an electric field in an expanding universe. Analyzing the ultraviolet asymptotic behavior of the kinetic functions, we found the divergent parts of the electric current and the energy-momentum tensor of the produced particles and determined the corresponding counterterms. The renormalized system of equations is used to study the generation of electromagnetic fields during and after inflation in the kinetic coupling model $\mathcal{L}_{\rm EM}=-(1/4)f^{2}(\phi)F_{\mu\nu}F^{\mu\nu}$ with the Ratra coupling function $f=\exp(\beta\phi/M_{p})$. It is found that the electric current of created particles is retarded with respect to the electric field. This leads to an oscillatory behavior of both quantities in agreement with the results obtained previously in phenomenological kinetic and hydrodynamical approaches., Comment: 22 pages, 4 figures

Using the kinetic approach, we study the impact of the charged particles dynamics due to the Schwinger effect on the electric field evolution during inflation. As a simple model of the electric field generation, we consider the kinetic coupling of the electromagnetic field to the inflaton via the term $f^{2}(\phi)F_{\mu\nu}F^{\mu\nu}$ with the Ratra coupling function $f=\exp(\beta\phi/M_{p})$. The production of charged particles is taken into account in the Boltzmann kinetic equation through the Schwinger source term. Produced particles are thermalized due to collisions which we model by using the collision integral in the self-consistent relaxation time approximation. We found that the current of created particles exhibits a non-Markovian character and cannot be described by a simple Ohm's law relation $j\propto E$. On the contrary, the electric current as well as the electric field are oscillatory functions of time with decreasing amplitudes and a phase difference due to the ballistic motion of charged carriers. Our qualitative results are checked by using a hydrodynamic approach. Deriving a closed system of equations for the number, current, and energy densities of charged particles and determining its solution, we find a good agreement with the results obtained in the kinetic approach., Comment: 17 pages, 4 figures

We study the electrical conductivity of hot Abelian plasma containing scalar charge carriers in the leading logarithmic order in coupling constant $\alpha$ using the Boltzmann kinetic equation. The leading contribution to the collision integral is due to the M{\o}ller and Bhabha scattering of scalar particles with a singular cross section in the region of small momentum transfer. Regularizing this singularity by taking into account the hard thermal loop corrections to the propagators of intermediate particles, we derive the second order differential equation which determines the kinetic function. We solve this equation numerically and also use a variational approach in order to find a simple analytical formula for the conductivity. It has the standard parametric dependence on the coupling constant $\sigma\approx 2.38\, T/(\alpha \log\alpha^{-1})$ with the prefactor taking a somewhat lower value compared to the fermionic case. Finally, we consider the general case of hot Abelian plasma with an arbitrary number of scalar and fermionic particle species and derive the simple analytical formula for its conductivity., Comment: 36 pages, 2 figures, 4 tables

We study the generation of electromagnetic fields during inflation when the conformal invariance of Maxwell's action is broken by the kinetic coupling $f^{2}(\phi)F_{\mu\nu}F^{\mu\nu}$ of the electromagnetic field to the inflaton field $\phi$. We consider the case where the coupling function $f(\phi)$ decreases in time during inflation and, as a result, the electric component of the energy density dominates over the magnetic one. The system of equations which governs the joint evolution of the scale factor, inflaton field, and electric energy density is derived. The backreaction occurs when the electric energy density becomes as large as the product of the slow-roll parameter $\epsilon$ and inflaton energy density, $\rho_{E}\sim \epsilon \rho_{\rm inf}$. It affects the inflaton field evolution and leads to the scale-invariant electric power spectrum and the magnetic one which is blue with the spectral index $n_{B}=2$ for any decreasing coupling function. This gives an upper limit on the present-day value of observed magnetic fields below $10^{-22}\,{\rm G}$. It is worth emphasizing that since the effective electric charge of particles $e_{\rm eff}=e/f$ is suppressed by the coupling function, the Schwinger effect becomes important only at the late stages of inflation when the inflaton field is close to the minimum of its potential. The Schwinger effect abruptly decreases the value of the electric field, helping to finish the inflation stage and enter the stage of preheating. It effectively produces the charged particles, implementing the Schwinger reheating scenario even before the fast oscillations of the inflaton. The numerical analysis is carried out in the Starobinsky model of inflation for the powerlike $f\propto a^{\alpha}$ and Ratra-type $f=\exp(\beta\phi/M_{p})$ coupling functions., Comment: 21 pages, 8 figures

We review the theoretical and experimental results connected with the electron states in two-dimensional Dirac systems paying a special attention to the atomic collapse in graphene. Two-electron bound states of a Coulomb impurity are considered too. A rather subtle role of a magnetic field in the supercritical charge problem in graphene is discussed. The electron states in the field of two equally charged impurities are studied and the conditions for supercritical instability to occur are determined. It is shown that the supercriticality of novel type is realized in gapped graphene with two unlikely charged impurities. For sufficiently large charges of impurities, it is found that the wave function of the occupied electron bound state of the highest energy changes its localization from the negatively charged impurity to the positively charged one as the distance between the impurities increases. The specifics of the atomic collapse in bilayer graphene is considered and it is shown that the atomic collapse in this material is not related to the phenomenon of the fall-to-center.

By assuming the kinetic coupling $f^{2}(\phi)FF$ of the effective inflaton field $\phi$ with the electromagnetic field, we explore magnetogenesis during the inflation and preheating stages in the Starobinsky model. We consider the case of the exponential coupling function $f(\phi)=\exp(\alpha\phi/M_{p})$ and show that for $\alpha \sim 12-15$ it is possible to generate the large scale magnetic fields with strength $\gtrsim 10^{-15}$ Gauss at the present epoch. The spectrum of generated magnetic fields is blue with the spectral index $n=1+s$, $s>0$. We have found that for the relevant values of the coupling parameter, $\alpha=12-15$, model avoids the back-reaction problem for all relevant modes., Comment: 10 pages, 5 figures

The electron states in the field of a charged impurity in graphene in a magnetic field are studied numerically. It is shown that a charged impurity removes the degeneracy of Landau levels converting them into bandlike structures. As the charge of impurity grows, the repulsion of sublevels of different Landau levels with the same value of orbital momentum takes place leading to the redistribution of the wave function profiles of these sublevels near the impurity. By studying the polarization effects, it is shown in agreement with the recent experiments that the effective charge of impurity can be very effectively tuned by chemical potential. If the chemical potential is situated inside a Landau level, then the charge of impurity is strongly diminished. In addition, the polarization function in this case has a peak at zero momentum, which leads to the sign-changing oscillations of the screened potential as a function of distance. If the chemical potential lies between the Landau levels, then the screened potential does not change sign, the screening is minimal, and the charged impurity can strongly affect the electron spectrum., Comment: 11 pages, 7 figures; final version published in PRB

We reveal a new type of supercritical behavior in gapped graphene with two oppositely charged impurities by studying the two-dimensional Dirac equation for quasiparticles with the Coulomb potential regularized at small distances accounting the lattice effects. By utilizing the variational Galerkin--Kantorovich method, we show that for supercritical electric dipole the wave function of the electron bound state changes its localization from the negatively charged impurity to the positively charged one as the distance between the impurities changes. Such a migration of the wave function corresponds to the electron and hole spontaneously created from the vacuum in bound states screening the positively and negatively charged impurities of the supercritical electric dipole, respectively. We generalize our results to a particle-hole asymmetric case, where the charges of impurities differ in signs and absolute values and demonstrate that the necessary energetic condition for the supercriticality of novel type to occur is that the energy levels of single positively and negatively charged impurities traverse together the energy distance separating the upper and lower continua. The robustness of the supercriticality of novel type is confirmed by the study of an exactly solvable 1D problem of the Dirac equation with the square well and barrier potential modeling an electric dipole potential., Comment: 15 pages, 13 Figures, and 3 video files in Supplemental Material. arXiv admin note: text overlap with arXiv:1506.08379

In the continuum low energy model, we calculate the one-loop dynamical polarization functions in ABC-stacked (rhombohedral) $n$-layer graphene in a magnetic field. Neglecting the trigonal warping effects, they are derived as functions of wavevector and frequency at finite chemical potential, temperature, band gap, and the width of Landau levels. The analytic results are given in terms of digamma functions and generalized Laguerre polynomials and have the form of double sums over Landau levels. Various particular limits for polarization functions (static, clean, etc.) are discussed. The intralayer and interlayer screened Coulomb potentials are numerically calculated as functions of momentum and frequency., 17 pages, 7 figures, small changes, to match the version published in PRB

The supercritical instability in a system of two identical charged impurities in gapped graphene described in the continuous limit by the two-dimensional Dirac equation has been studied. The case where the charge of each impurity is subcritical, but their sum exceeds the critical value calculated in the version with a single Coulomb center, is considered. Using the developed variational method, the dependence of the critical distance $R_{\rm cr}$ between the impurities on their total charge is calculated. The $R_{\rm cr}$-value is found to grow as the total impurity charge increases and the quasiparticle band gap decreases. The results of calculations are compared with those obtained in earlier researches., 10 pages, 2 figures

We study the supercritical instability in gapped graphene with two charged impurities separated by distance R using the two-dimensional Dirac equation for electron quasiparticles. Attention is paid to a situation when charges of impurities are subcritical, whereas their total charge exceeds a critical one. The critical distance R_{cr} in the system of two charged centers is defined as that at which the electron bound state with the lowest energy reaches the boundary of the lower continuum. A variational calculation of the critical distance R_{cr} separating the supercritical (RR_{cr}) regimes is carried out. It is shown that the critical distance R_{cr} increases as the quasiparticle gap decreases. The energy and width of a quasistationary state as functions of the distance between two impurities are derived in the quasiclassical approximation., Comment: 12 pages, 1 figure, comments and references added, to match the published version

We study the Dirac equation for quasiparticles in gapped graphene with two oppositely charged impurities by using the technique of linear combination of atomic orbitals and variational Galerkin--Kantorovich method. We show that for sufficiently large charges of impurities the wave function of the occupied electron bound state of the highest energy changes its localization from the negatively charged impurity to the positively charged one as the distance between the impurities increases. This migration of the electron wave function of supercritical electric dipole is a generalization of the familiar phenomenon of the atomic collapse of single charged impurity to the case where electron-hole pairs are spontaneously created from vacuum in bound states with charge impurities thus partially screening them., 8 pages, 6 Figures, and video; submitted to Europhysics Letters

In the low-energy two-band as well as four-band continuum models, we study the supercritical instability in gapped bilayer graphene in the field of a charged impurity. It is found that the screening effects are crucially important in bilayer graphene. If they are neglected, then the critical value for the impurity charge as the lowest-energy bound state dives into the lower continuum tends to zero as the gap $\Delta$ vanishes. If the screened Coulomb interaction is considered, then the critical charge tends to a finite value for $\Delta \to 0$. The different scalings of the kinetic energy of quasiparticles and the Coulomb interaction with respect to the distance to the charged impurity ensure that the wave function of the electron bound state does not shrink toward the impurity as its charge increases. This results in the absence of the fall-to-center phenomenon in bilayer graphene although the supercritical instability is realized., Comment: 15 pages, 6 figures

We study the process of inflationary magnetogenesis in the natural single-field inflation model, whose parameters are chosen in accordance with the recent observations by the Planck collaboration [1]. The conformal invariance of the Maxwell action is broken by a kinetic coupling with the inflaton field by means of the coupling function as a power of the scale factor, I(phi) proportional to a(alpha), and alpha < 0 is used in order to avoid the strong coupling problem. For such.., the electric component of the energy density dominates over the magnetic one and, for alpha less than or similar to -2.2, it causes a strong back-reaction, which can spoil inflation and terminate the enhancement of the magnetic field. It is found that the magnetic fields generated without back-reaction problem cannot exceed similar to 10(-20) G at the present epoch, and their spectrum has a blue tilt.