Relations between transfer matrices and numerical stability analysis to avoid the $\Omega d$ problem

Abstract

The transfer matrix method is usually employed to study problems described by $N$ equations of matrix Sturm-Liouville (MSL) kind. In some cases a numerical degradation (the so called $\Omega d$ problem) appears thus impairing the performance of the method. We present here a procedure that can overcome this problem in the case of multilayer systems having piecewise constant coefficients. This is performed by studying the relations between the associated transfer matrix and other transfer matrix variants. In this way it was possible to obtain the matrices which can overcome the $\Omega d$ problem in the general case and then in problems which are particular cases of the general one. In this framework different strategies are put forward to solve different boundary condition problems by means of these numerically stable matrices. Numerical and analytic examples are presented to show that these stable variants are more adequate than other matrix methods to overcome the $\Omega d$ problem. Due to the ubiquity of the MSL system, these results can be applied to the study of many elementary excitations in multilayer structures.

23 pages, 2 figures

References in corpus (2)

• Wave impedance matrices for cylindrically anisotropic radially inhomogeneous elastic solids
• Stable methods to solve the impedance matrix for radially inhomogeneous cylindrically anisotropic structures

Cited by in corpus (8)

• Anomalous energy transport in laminates with exceptional points
• Dynamic homogenization of composite and locally resonant flexural systems
• Hybrid matrix method for stable numerical analysis of the propagation of Dirac electrons in gapless bilayer graphene superlattices
• Quantum thermostatted disordered systems and sensitivity under compression
• Time-Energy Uncertainty Limit for spin-related wavepacket evolution
• Surveying the Multicomponent Scattering Matrix: Unitarity and Symmetries
• Modeling Nonlinear Optics with the Transfer Matrix Method
• Enhancement of the thermoelectric properties in bilayer graphene structures induced by Fano resonances