Linear Systems Research Articles

Variational quantum algorithm for the Poisson equation based on the banded Toeplitz systems

Abstract For solving the Poisson equation it is usually possible to discretize it into solving the corresponding linear system $Ax=b$. Variational quantum algorithms (VQAs) for \textcolor{red}{the discreted} Poisson equation have been studied before. We give a VQA based on the banded Toeplitz systems for solving the Poisson equation with respect to the structural features of matrix $A$. In detail, we decompose the matrix $A$ and $A^2$ into a linear combination of the corresponding banded Toeplitz matrix and sparse matrices with only a few non-zero elements. For the one-dimensional Poisson equation with different boundary conditions and the $d$-dimensional Poisson equation with Dirichlet boundary conditions, \textcolor{red}{the number of decomposition terms is less than the work in [Phys. Rev. A108, 032418 (2023)]}. Based on the decomposition of the matrix, we design quantum circuits that \textcolor{red}{evaluate efficiently} the cost function. \textcolor{red}{Additionally, numerical simulation verifies the feasibility of the proposed algorithm. In the end,} the VQAs for linear systems of equations and matrix–vector multiplications with $K$-banded Teoplitz matrix $T_n^K$ \textcolor{red}{are} given, where $T_n^K\in R^{n\times n}$ and $K\in O(ploy\log n)$.

Free vibration and supersonic flutter analyses of a sandwich cylindrical shell with CNT-reinforced honeycomb core integrated with piezoelectric layers

In this research, the free vibration and supersonic flutter analyses of a cylindrical sandwich shell with a reinforced honeycomb core (RHC) integrated with piezoelectric face sheets are investigated. The core is a honeycomb structure enriched with fibers and carbon nanotubes (CNTs). The layers of the shell are modeled based on the Cooper-Naghdi shell theory and the continuity conditions between the layers are considered. The aerodynamic pressure of the supersonic flow is modeled using the linear piston theory. Hamilton’s principle is used to derive the governing equations. A semi-analytical solution is presented including an exact solution in the circumferential direction followed by a numerical solution in the longitudinal direction via the differential quadrature method (DQM). The effects of several parameters on the frequencies, the corresponding damping coefficients, and the critical aerodynamic pressure are investigated. It is concluded that improving the stiffness of the honeycomb core with the CNTs does not necessarily improve the aeroelastic stability of the shell. The presented results can be used in the design and analysis of aerospace structures.

Amplitude variation with azimuth direct inversion for effective pressure-sensitivity parameters in horizontal transversely isotropic media

Effective pressure prediction is of great significance for fracturability evaluation in unconventional reservoirs. Quantitative characterization of effective pressure remains a great challenge in the inversion prediction of subsurface media. Our goal is to develop a novel Bayesian amplitude varying with azimuth (AVAZ) inversion method to estimate fracture weaknesses and effective pressure for a horizontal transversely isotropic (HTI) medium. Using Schoenberg’s linear slip model and pore space compressibility theory, an anisotropic saturated stiffness matrix directly characterized by effective pressure is derived based on a weak-contrast interface separating two weakly anisotropic half-spaces of HTI medium. Based on the perturbation stiffness matrix and scattering theory, we derive a PP-wave reflection coefficient approximation directly characterized by fracture weaknesses and effective pressure-sensitive parameters, which establishes a quantitative relationship between seismic reflection coefficient and effective pressure. The analysis of accuracy and sensitivity is performed to verify the reliability and applicability of the constructed reflection coefficient. And then, the azimuth-amplitude-difference forward solver is constructed based on seismic reflection coefficients to eliminate the influence of isotropic background perturbations on seismic reflection coefficients and reduce the condition number of the forward solver. Finally, we present a novel Bayesian azimuth-amplitude-difference AVAZ inversion method for transversely isotropic media constrained by Cauchy and Gaussian regularization to directly estimate the fracture weaknesses and effective pressure-sensitive parameters. The application of synthetic data and field data verifies that the inversion results are in good agreement with the actual data, indicating that the method can effectively realize the stable prediction of anisotropic parameters and effective pressure-sensitive parameters of complex media.

A quantile-based block Kaczmarz algorithm for solving large consistent linear systems

A quantile-based block Kaczmarz algorithm for solving large consistent linear systems

Static output feedback control for positive 2D discrete systems in Roesser model

Abstract This paper investigates the output stabilization problem for multi-input multi-output positive 2D systems described by a linear discrete-time Roesser system. This kind of systems have the property that the horizontal and vertical states take non-negative values whenever the initial boundaries are non-negative. Conditions for the existence of desired static output feedback controllers guaranteeing the resultant closed-loop system to be positive and asymptotically stable are obtained. Also, the synthesis of static output feedback controllers, including the requirement of positivity of the controllers, is solved. These results are then extended to the case of uncertain plants. All of the proposed conditions are expressed in terms of Linear Programming, which can be easily verified by using standard numerical software. Finally, several numerical examples are included to illustrate the validity and the effectiveness of the developed results.

Stability Analysis of Matrices and Single-Parameter Matrix Families on Special Regions

This paper proposes an efficient method for determining the D-stability of matrices using Kelley’s cutting-plane method. The study examines the D-stability of matrices in symmetric regions of the complex plane, defined by quadratic matrix inequalities (QMI) and polynomial functions. Using Kelley’s cutting-plane method, we present an iterative algorithm for determining the D-stability of a given matrix. Further, the robustness of single-parameter matrix families is analyzed, and a method is proposed for considering their D-stability. Through examples, we indicate the possible applicability of the proposed approach in addressing stability problems in linear systems.

3D Controlled‐source electromagnetic modelling in anisotropic media using secondary potentials and a cascadic multigrid solver

AbstractQuantitative interpretation of the data from controlled‐source electromagnetic methods, whether via forward modelling or inversion, requires solving a considerable number of forward problems, and multigrid methods are often employed to accelerate the solving process. In this study, a new extrapolation cascadic multigrid method is employed to solve the large sparse complex linear system arising from the finite element approximation of Maxwell's equations using secondary potentials. The equations using secondary potentials are discretized by the classic nodal finite element method on nonuniform rectilinear grids. The resulting linear systems are solved by the extrapolation cascadic multigrid method with a new prolongation operator and preconditioned Stabilized bi‐conjugate gradient method smoother. High‐order interpolation and global extrapolation formulas are utilized to construct the multigrid prolongation operator. The extrapolation cascadic multigrid method with the new prolongation operator is easier to implement and more flexible in application than the original one. Finally, several synthetic examples including layered models, models with anisotropic anomalous bodies or layers, are used to validate the accuracy and efficiency of the proposed method. Numerical results show that the extrapolation cascadic multigrid method improves the efficiency of 3D controlled‐source electromagnetic forward modelling a lot, compared with traditional iterative solvers and some state‐of‐the‐art methods or software (e.g., preconditioned flexible generalized minimal residual method, emg3d) in the considered models and grid settings. The efficiency benefit is more evident as the number of unknowns increases, and the proposed method is more efficient at low frequencies. The extrapolation cascadic multigrid method can also be used to solve systems of equations arising from related applications, such as induction logging, airborne electromagnetic, etc.

Impulse Controllability for Singular Hybrid Coupled Systems

This study examines the concept of impulse controllability within singular hybrid coupled systems through the utilisation of decentralised proportional plus derivative (P-D) output feedback. By employing the Differential Mean Value Theorem, the nonlinear model can be converted into a linear parameter-varying large-scale system. Our analysis leads to the establishment of algebraic conditions that are both necessary and sufficient for the existence of a decentralised P-D output feedback controller that can guarantee impulse controllability in these complex systems. Moreover, we address the issue of admissibility within these systems by employing matrix trace inequalities. We present a novel sufficient condition for impulse controllability, which offers a new perspective on addressing this challenging problem. To validate our findings, we present numerical examples that demonstrate the effectiveness of the proposed methodologies in practice.

Generalized linear response theory for pumped systems and its application to transient optical properties

Generalized linear response theory for pumped systems and its application to transient optical properties

Effect of Natural Ionizing Radiation on Health Indicators in Region with Monazite Sand in Brazil

Abstract Objective This article compares the occurrence of cancer and germinal damage in the city of Guarapari-Espírito Santo (ES), an area with high natural ionizing radiation, with other coastal cities. Material and Methods The evaluated cities were: Guarapari (ES), Campos (Rio de Janeiro), Rio Grande (Rio Grande do Sul), and Ilhéus (Bahia); the following factors were considered: mortality rate % from tumors (2007–2017), mortality rate % from tumors up to 15 years of age (2007–2017), hospital morbidity rate % from tumors up to 15 years of age (2008–2018), and hospital morbidity rate % from congenital malformation and chromosomal anomalies (2008–2018). Radiometric surveys were conducted by the Applied Physics Group at Universidade Federal do Espírito Santo. Results The natural radiation levels recorded in Guarapari resulted in accumulated dose between 3.65 and 10.95 mSv/year, and 1 mSv in the other cities. The highest average cancer mortality rates were: Rio Grande = 22.4%; Guarapari = 17.6%; Campos = 16.7%; and Ilhéus = 11.8%. The proportional hospital mortality and morbidity rates for cancer up to 15 years of age are as follows: Ilhéus = 3.36 and 5.87%; Rio Grande = 0.79 and 7.38%; Guarapari = 0.64 and 7.25%; and Campos = 0.39 and 9.13%. The hospital morbidity rate due to congenital malformations and chromosomal anomalies was 0.72% for Campos, 0.63% for Rio Grande, 0.62% for Guarapari, and 0.43% for Ilhéus. Conclusion There was no increase in cases of cancer or germ damage in Guarapari. These results indicated a dose threshold in the induction of these damages, contradicting the current linear no-threshold theory.

Constructing a model of the axis form in a S-shaped riser of a cultivator paw

When cultivating the soil, a force of resistance to the movement of cultivator paw acts on it. It has a variable value and causes a moment of force applied to the riser of the paw. Under the action of the moment, the elastic axis of the riser changes its shape. This affects the position of the paw in the soil. The form of an S-shaped riser whose elastic axis consists of two circle arcs has been considered in this study. During cultivator operation, one part of the riser bends, increasing the curvature of the elastic axis, and the other, on the contrary, unbends, that is, its curvature decreases. The modeling of the shape of the elastic axis of the paw riser is based on the theory of resistance of materials, according to which the curvature of the elastic axis of the cantilevered band is directly proportional to the applied moment and inversely proportional to its stiffness. If the shape of the cross-section of the riser along its entire length is unchanged and the properties of the metal are also the same, then the stiffness is constant. In the case of small deflections of the band, the linear theory of bending is used, but the deflections in the riser are significant, so the nonlinear theory has been used for this case. At the same time, it is taken into account that the elastic axis of the riser already has an initial curvature, the sign of which changes after passing through the point of connection of the component arcs. To model the shape of the elastic axis of the S-shaped paw riser, the deformation of the arcs of the circles that form this paw was calculated separately. Numerical integration methods were used to find the shape of the deformed elastic axes of both parts of the riser. They were connected into a whole and a deformed elastic axis of the S-shaped riser was obtained. Two variants of the riser with different lengths of their elastic axis, but the same height and the same angle of entry into the soil, were considered. The combination of the component arcs of the riser shows that, under the action of the same force, the deviation from the specified movement depth for one riser is 2 cm, and for the other – 4 cm.

Existence of a mild solution and approximate controllability for fractional random integro-differential inclusions with non-instantaneous impulses

This paper investigates the existence and approximate controllability (ACA) of fractional neutral-type stochastic differential inclusions (NTSDIs) characterized by non-instantaneous impulses within a separable Hilbert space (HS) framework. Employing the Atangana–Baleanu–Caputo (ABC) derivative, we transform the system into an equivalent fixed-point (FP) problem through an integral operator. Subsequently, the Bohnenblust–Karlin FP theorem is leveraged to establish existence results. By assuming ACA of the corresponding linear system, we derive sufficient conditions for the ACA of the nonlinear stochastic impulsive control system. Our analysis relies on concepts from stochastic analysis, fractional calculus, FP theory, semigroup theory, and the theory of multivalued maps (MVMs). The theoretical findings are illustrated through a concrete example.

Simultaneous fault-tolerant consensus control and disturbance suppression for multi-agent systems with polynomial form fault under switching topology

For continuous-time linear leader-follower multi-agent systems with polynomial faults, the problem of fault-tolerant state consensus under switching topologies is studied. firstly, the average dwell time method is introduced, and observers are designed for different topologies to observe faults and disturbances. Secondly, according to the observed estimated value, the corresponding fault-tolerant control consensus protocol is proposed considering the switching situation of the communication topology. Finally, both theoretical analysis and simulation examples verify that the proposed design scheme can compensate faults and disturbance rejection.

Fractional Einstein–Gauss–Bonnet Scalar Field Cosmology

Our paper introduces a new theoretical framework called the Fractional Einstein–Gauss–Bonnet scalar field cosmology, which has important physical implications. Using fractional calculus to modify the gravitational action integral, we derived a modified Friedmann equation and a modified Klein–Gordon equation. Our research reveals non-trivial solutions associated with exponential potential, exponential couplings to the Gauss–Bonnet term, and a logarithmic scalar field, which are dependent on two cosmological parameters, m and α0=t0H0 and the fractional derivative order μ. By employing linear stability theory, we reveal the phase space structure and analyze the dynamic effects of the Gauss–Bonnet couplings. The scaling behavior at some equilibrium points reveals that the geometric corrections in the coupling to the Gauss–Bonnet scalar can mimic the behavior of the dark sector in modified gravity. Using data from cosmic chronometers, type Ia supernovae, supermassive Black Hole Shadows, and strong gravitational lensing, we estimated the values of m and α0, indicating that the solution is consistent with an accelerated expansion at late times with the values α0=1.38±0.05, m=1.44±0.05, and μ=1.48±0.17 (consistent with Ωm,0=0.311±0.016 and h=0.712±0.007), resulting in an age of the Universe t0=19.0±0.7 [Gyr] at 1σ CL. Ultimately, we obtained late-time accelerating power-law solutions supported by the most recent cosmological data, and we proposed an alternative explanation for the origin of cosmic acceleration other than ΛCDM. Our results generalize and significantly improve previous achievements in the literature, highlighting the practical implications of fractional calculus in cosmology.

Cavity-magnon-polariton spectroscopy of strongly hybridized electro-nuclear spin excitations in LiHoF4

We first present a formalism that incorporates the input-output formalism and the linear response theory to employ cavity-magnon-polariton coupling as a spectroscopic tool for investigating strongly hybridized electro-nuclear spin excitations. A microscopic relation between the generalized susceptibility and the scattering parameter |S11|\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$|S_{11}|$$\\end{document} in strongly hybridized cavity-magnon-polariton systems has been derived without resorting to semi-classical approximations. The formalism is then applied to both analyze and simulate a specific systems comprising a model quantum Ising magnet (LiHoF4\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\hbox {LiHoF}_4$$\\end{document}) and a high-finesse 3D re-entrant cavity resonator. Quantitative information on the electro-nuclear spin states in LiHoF4\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\hbox {LiHoF}_4$$\\end{document} is extracted, and the experimental observations across a broad parameter range were numerically reproduced, including an external magnetic field traversing a quantum critical point. The method potentially opens a new avenue not only for further studies on the quantum phase transition in LiHoF4\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\hbox {LiHoF}_4$$\\end{document} but also for a wide range of complex magnetic systems.

A least-squares Fourier frame method for nonlocal diffusion models on arbitrary domains

We introduce a least-squares Fourier frame method for solving nonlocal diffusion models with Dirichlet volume constraint on arbitrary domains. The mathematical structure of a frame rather than a basis allows using a discrete least-squares approximation on irregular domains and imposing non-periodic boundary conditions. The method has inherited the one-dimensional integral expression of Fourier symbols of the nonlocal diffusion operator from Fourier spectral methods for any d spatial dimensions. High precision of its solution can be achieved via a direct solver such as pivoted QR decomposition even though the corresponding system is extremely ill-conditioned, due to the redundancy in the frame. The extension of AZ algorithm improves the complexity of solving the rectangular linear system to O(Nlog2⁡N) for 1d problems and O(N2log2⁡N) for 2d problems, compared with O(N3) of the direct solvers, where N is the number of degrees of freedom. We present ample numerical experiments to show the flexibility, fast convergence and asymptotical compatibility of the proposed method.

Limit Cycles in a Class of Continuous–Discontinuous Piecewise Linear Differential Systems with Three Zones

In this paper, we investigate the existence of limit cycles in piecewise linear systems separated by two parallel straight lines, where in one of these straight lines, the piecewise differential system is continuous, and in the other, it is discontinuous. Based on the methods of Poincaré map and first integral, we prove that these systems can have at most two limit cycles. Additionally, taking into account the type of singular points in the zones, we present examples with one or two limit cycles for all possible combinations.

Three double multi-step randomized extended Kaczmarz methods for solving large sparse inconsistent linear systems

Three double multi-step randomized extended Kaczmarz methods for solving large sparse inconsistent linear systems

Contributions to Generalized Oscillation Theory of Linear Hamiltonian Systems

In this paper we present several new contributions to the oscillation theory of linear differential equations, in particular of linear Hamiltonian systems, when the traditional Legendre condition is absent. Following our recent work (Discrete Contin. Dyn. Syst. 43(12):4139–4173, 2023), we introduce the multiplicity of a generalized right focal point and derive the corresponding local Sturmian separation theorem. We also examine the relation between the existence of finitely many generalized right focal points, or in the special case the nonexistence of generalized right focal points, with the Legendre condition. As the main tools we use new notions of the minimal multiplicities at a given point and the dual comparative index — an object from matrix analysis or differential geometry (Maslov index theory). Furthermore, we study local limit properties of the dual comparative index and the comparative index and apply them for deriving new oscillation results phrased in terms of the generalized right and left focal points. The investigation of the interplay between generalized right and left focal points leads to conditions characterizing the situation, when in the local Sturmian separation theorem the corresponding multiplicities attain the minimal possible value. This also provides a generalization of the concepts of the right and left proper focal point defined by Kratz (Analysis 23(2):163–183, 2003) and Wahrheit (Int. J. Differ. Equ. 2(2):221–244, 2007) to the setting, which does not impose the Legendre condition. The results are new even for completely controllable linear Hamiltonian systems, including the Sturm–Liouville differential equations of arbitrary even order.

A study of optical properties and electron energy loss spectra of ZnS by linear response theory

The linear response theory has been applied to study the optical properties and electron energy loss spectra of ZnS. The ground state is built using the Full-Potential Linearized Augmented Plane Wave method. Thereafter, the electronic properties such as band structure and density of states are calculated. Calculations of electron energy loss spectra and optical properties are performed on the top of these ground state electronic properties. The calculations are performed using three exchange–correlation kernels namely, adiabatic local density approximation, long range contribution, and the bootstrap. The random phase approximation is also considered. An approach proposed by us to find the material dependent parameter α is used in long range contribution kernel to capture more accurate optical properties. Three values of α viz 0.78, 1.09, and 1.22 are considered. It is observed that the α = 0.78 gives very good results. This proposed scheme is found to work very well for ZnS also. The effect of local field and electron–hole interaction on the spectra is examined. For example after incorporating local field effect the high frequency dielectric constant ε∞ [i.e. ε1E→0eV ], reduces by 9.19, 12.26, 13.51 and 11.02% for random phase approximation, adiabatic local density approximation, long range contribution and bootstrap kernel respectively. In EELS plasmon peak positions are found using the three kernels and random phase approximation. A good accord with experimental results is noted after incorporating local field effects. It is observed that local field effect causes redshift of prominent peak (i.e. plasmon peak) in electron energy loss spectra by 0.61 eV.