Search

We numerically demonstrate the excitation of leaky spin waves (SWs) guided along a ferromagnetic stripe by an obliquely incident SW beam on the thin film edge placed below the stripe. During propagation, leaky waves emit energy back to the layer in the form of plane waves and several laterally shifted parallel SW beams. This resonance excitation, combined with interference effects of the reflected and re-emitted waves, results in the magnonic Woods anomaly and significant increase of the Goos-Hanchen shift magnitude. Hence, we provide a unique platform to control SW reflection and to transfer SWs from a 2D platform into the 1D guiding mode that can be used to form a transdimensional magnonic router.

Efficient numerical methods are required for the design of optimised devices. In magnonics, the primary computational tool is micromagnetic simulations, which solve the Landau-Lifshitz equation discretised in time and space. However, their computational cost is high, and the complexity of their output hinders insight into the physics of the simulated system, especially in the case of multimode propagating wave-based devices. We propose a finite-element modal method allowing an efficient solution of the scattering problem for dipole-exchange spin waves propagating perpendicularly to the magnetisation direction. The method gives direct access to the scattering matrix of the whole system and its components. We extend the formula for the power carried by a magnetostatic mode in the Damon-Eshbach configuration to the case with exchange, allowing the scattering coefficients to be normalised to represent the fraction of the input power transferred to each output channel. We apply the method to the analysis of spin-wave scattering on a basic functional block of magnonic circuits, consisting of a resonator dynamically coupled to a thin film. The results and the method are validated by comparison with micromagnetic simulations., Comment: Revised version

Using frequency-domain finite element calculations cross-checked with micromagnetic simulations, we demonstrate that the phase of spin waves reflected from an interface between a permalloy film and a bilayer can be controlled by changing dimensions of the bilayer. Treating the bilayer formed by the permalloy film and a ferromagnetic stripe as a segment of a multi-mode waveguide, we show that spin-wave Fabry-Perot resonances of one of its modes are responsible for the high sensitivity of the phase of reflected waves to stripe width and the stripe-film separation. Thus, the system is a unique realization of a fully magnonic Gires-Tournois interferometer based on a two-modes resonator, which can be treated as a magnonic counterpart of a metasurface, since it enables manipulation of the phase of spin waves at subwavelength distances. Knowledge gained from these calculations might be used to design magnonic devices such as flat lenses or magnetic particle detectors.

Operator products occur naturally in a range of regularized boundary integral equation formulations. However, while a Galerkin discretisation only depends on the domain space and the test (or dual) space of the operator, products require a notion of the range. In the boundary element software package Bempp we have implemented a complete operator algebra that depends on knowledge of the domain, range and test space. The aim was to develop a way of working with Galerkin operators in boundary element software that is as close to working with the strong form on paper as possible while hiding the complexities of Galerkin discretisations. In this paper, we demonstrate the implementation of this operator algebra and show, using various Laplace and Helmholtz example problems, how it significantly simplifies the definition and solution of a wide range of typical boundary integral equation problems.

In recent years there have been tremendous advances in the theoretical understanding of boundary integral equations for Maxwell problems. In particular, stable dual pairing of discretisation spaces have been developed that allow robust formulations of the preconditioned electric field, magnetic field and combined field integral equations. Within the BEM++ boundary element library we have developed implementations of these frameworks that allow an intuitive formulation of the typical Maxwell boundary integral formulations within a few lines of code. The basis of these developments is an efficient and robust implementations of Calder\'on identities together with a product algebra that hides and automates most technicalities involved in assembling Galerkin boundary integral equations. In this paper we demonstrate this framework and use it to derive very simple and robust software formulations of the standard preconditioned electric field, magnetic field and regularised combined field integral equations for Maxwell.

It is widely accepted that, on ensemble average, the transmission T of guided modes decays exponentially with the waveguide length L due to small imperfections, leading to the important figure of merit defined as the attenuation-rate coefficient alpha = -/L. In this letter, we evidence that the exponential-damping law is not valid in general for periodic monomode waveguides, especially as the group velocity decreases. This result that contradicts common beliefs and experimental practices aiming at measuring alpha is supported by a theoretical study of light transport in the limit of very small imperfections, and by numerical results obtained for two waveguide geometries that offer contrasted damping behaviours.

This article concerns the design of antireflection structures which, placed on a photonic crystal surface, significantly diminish the fraction of energy lost to reflected waves. After a review of the classes of these structures proposed to date, a new method is presented in detail for the design of antireflection gratings operating in a wide range of angles of incidence. The proposed algorithm is illustrated by means of several examples, showing the advantages and limitations., Comment: Submitted to Phys. Rev. B

We propose a new class of compact integrated optical circulators providing a large isolation level while maintaining a straightforward technological feasibility. Their layout is based on a nonreciprocal radial Bragg cavity composed of concentric magneto-optical rings. The circulator ports are standard rib waveguides, butt-coupled to the cavity by cutting through its outer rings. The device is specifically designed for operation in a uniform external magnetic field. Using a coupled-mode description of the complete cavity/waveguide-port system, we explore the rich behaviour of cavity circulators in presence of varying levels of direct port-to-port coupling. We demonstrate numerically a strongly miniaturized two-dimensional cavity circulator, with a total footprint of less than $(10\lambda)^2$, achieving a 20-dB isolation level at telecom frequencies over a bandwidth of 130 GHz. The device is found to be very tolerant with respect to fabrication imperfections. We finish with an outlook on three-dimensional versions of the proposed nonreciprocal cavities., Comment: Revised

We present an extension of the multiple-scattering calculational method to the case of gyrotropic media.

By exploiting the concepts of magnetic group theory we show how unidirectional behavior can be obtained in two-dimensional magneto-photonic crystals (MOPhC) with uniform magnetization. This group theory approach generalizes all previous investigations of one-way MOPhCs including those based on the use of antiparallel magnetic domains in the elementary crystal cell. Here, the theoretical approach is illustrated for one MOPhC example where unidirectional behavior is obtained by appropriately lowering the geometrical symmetry of the elementary motifs. One-way transmission is numerically demonstrated for a photonic crystal slice., Comment: 4 pages, 3 figures, submitted to Phys. Rev. Lett

We propose a technique of compensating the spurious reflections implied by the multiple-scattering (MS) method, commonly used for analyzing finite photonic crystal (PC) systems, to obtain exact values of characteristic parameters, such as reflection and transmission coefficients, of PC functional elements. Rather than a modification of the MS computational algorithm, our approach involves postprocessing of results obtained by the MS method. We derive analytical formulas for the fields excited in a finite system, taking explicitly into account the spurious reflections occurring at the artificial system boundaries. The intrinsic parameters of the investigated functional element are found by fitting the results of MS simulations to those obtained from the formulas derived. Devices linked with one and two semi-infinite waveguides are analyzed explicitly; possible extensions of the formalism to more complex circuits are discussed as well. The accuracy of the proposed method is tested in a number of systems; the results of our calculations prove to be in good agreement with those obtained independently by other authors., Comment: 8 pages, 5 figures

We propose a quantitative model explaining the mechanism of light collimation by leaky surface modes that propagate on a corrugated surface around the output of a photonic crystal waveguide. The dispersion relation of these modes is determined for a number of surface terminations. Analytical results obtained on the basis of the model are compared to those of rigorous numerical simulations. Maximum collimation is shown to occur at frequency values corresponding to excitation of surface modes whose wave number retains a nonzero real part., Comment: 6 pages, 7 figures. Version 2: corrected sign of k_x' (sections 4-6, fig. 2), minor clarifications in section 2. Version 3: significant changes, including reformulation of the model using the theory of aperture antennas, as well as extended discussion of the accuracy of the model