Sandpile groups of generalized de Bruijn and Kautz graphs and circulant matrices over finite fields

A maximal minor $M$ of the Laplacian of an $n$-vertex Eulerian digraph $\Gamma$ gives rise to a finite group $\mathbb{Z}^{n-1}/\mathbb{Z}^{n-1}M$ known as the sandpile (or critical) group $S(\Gamma)$ of $\Gamma$. We determine $S(\Gamma)$ of the generalized de Bruijn graphs $\Gamma=\mathrm{DB}(n,d)$ with vertices $0,\dots,n-1$ and arcs $(i,di+k)$ for $0\leq i\leq n-1$ and $0\leq k\leq d-1$, and closely related generalized Kautz graphs, extending and completing earlier results for the classical de Bruijn and Kautz graphs. Moreover, for a prime $p$ and an $n$-cycle permutation matrix $X\in\mathrm{GL}_n(p)$ we show that $S(\mathrm{DB}(n,p))$ is isomorphic to the quotient by $\langle X\rangle$ of the centralizer of $X$ in $\mathrm{PGL}_n(p)$. This offers an explanation for the coincidence of numerical data in sequences A027362 and A003473 of the OEIS, and allows one to speculate upon a possibility to construct normal bases in the finite field $\mathbb{F}_{p^n}$ from spanning trees in $\mathrm{DB}(n,p)$.Comment: I+24 page

An Efficient Block Circulant Preconditioner For Simulating Fracture Using Large Fuse Networks

{\it Critical slowing down} associated with the iterative solvers close to the critical point often hinders large-scale numerical simulation of fracture using discrete lattice networks. This paper presents a block circlant preconditioner for iterative solvers for the simulation of progressive fracture in disordered, quasi-brittle materials using large discrete lattice networks. The average computational cost of the present alorithm per iteration is $O(rs log s) + delops$, where the stiffness matrix ${\bf A}$ is partioned into $r$-by-$r$ blocks such that each block is an $s$-by-$s$ matrix, and $delops$ represents the operational count associated with solving a block-diagonal matrix with $r$-by-$r$ dense matrix blocks. This algorithm using the block circulant preconditioner is faster than the Fourier accelerated preconditioned conjugate gradient (PCG) algorithm, and alleviates the {\it critical slowing down} that is especially severe close to the critical point. Numerical results using random resistor networks substantiate the efficiency of the present algorithm.Comment: 16 pages including 2 figure

The Dynamics of Charges Induced by a Charged Particle Traversing a Dielectric Slab

We studied the dynamics of surfacea and wake charges induced by a charged particle traversing a dielectric slab. It is shown that after the crossing of the slab first boundary, the induced on the slab surface charge (image charge) is transformed into the wake charge, which overflows to the second boundary when the particle crosses it. It is also shown, that the polarization of the slab is of an oscillatory nature, and the net induced charge in a slab remains zero at all stages of the motion.Comment: 12 pages, 1 figur

Phase transition in the Higgs model of scalar dyons

In the present paper we investigate the phase transition "Coulomb--confinement" in the Higgs model of abelian scalar dyons -- particles having both, electric $e$ and magnetic $g$, charges. It is shown that by dual symmetry this theory is equivalent to scalar fields with the effective squared electric charge e^{*2}=e^2+g^2. But the Dirac relation distinguishes the electric and magnetic charges of dyons. The following phase transition couplings are obtained in the one--loop approximation: \alpha_{crit}=e^2_{crit}/4\pi\approx 0.19, \tilde\alpha_{crit}=g^2_{crit}/4\pi\approx 1.29 and \alpha^*_{crit}\approx 1.48.Comment: 16 pages, 2 figure

Excitation of surface plasmon-polaritons in metal films with double periodic modulation: anomalous optical effects

We perform a thorough theoretical analysis of resonance effects when an arbitrarily polarized plane monochromatic wave is incident onto a double periodically modulated metal film sandwiched by two different transparent media. The proposed theory offers a generalization of the theory that had been build in our recent papers for the simplest case of one-dimensional structures to two-dimensional ones. A special emphasis is placed on the films with the modulation caused by cylindrical inclusions, hence, the results obtained are applicable to the films used in the experiments. We discuss a spectral composition of modulated films and highlight the principal role of ``resonance'' and ``coupling'' modulation harmonics. All the originating multiple resonances are examined in detail. The transformation coefficients corresponding to different diffraction orders are investigated in the vicinity of each resonance. We make a comparison between our theory and recent experiments concerning enhanced light transmittance and show the ways of increasing the efficiency of these phenomena. In the appendix we demonstrate a close analogy between ELT effect and peculiarities of a forced motion of two coupled classical oscillators.Comment: 24 pages, 17 figure

Radiating Shear-Free Gravitational Collapse with Charge

We present a new shear free model for the gravitational collapse of a spherically symmetric charged body. We propose a dissipative contraction with radiation emitted outwards. The Einstein field equations, using the junction conditions and an ansatz, are integrated numerically. A check of the energy conditions is also performed. We obtain that the charge delays the black hole formation and it can even halt the collapse.Comment: 22 pages, 9 figures. It has been corrected several typos and included several references. Accepted for publication in GR