Polynomial Expansion Research Articles

Chebyshev polynomial approach to Loschmidt echo: Application to quench dynamics in two-dimensional quasicrystals.

The understanding of quantum phase transitions in disordered or quasicrystal media is a central issue in condensed matter physics. In this paper we investigate localization properties of the two-dimensional Aubry-André model. We find that the system exhibits self-duality for the transformation between position and momentum spaces at a critical quasiperiodic potential, leading to an energy-independent Anderson transition. Most importantly, we present the implementation of an efficient and accurate algorithm based on the Chebyshev polynomial expansion of the Loschmidt echo, which characterizes the nonequilibrium dynamics of quantum quenched quasiperiodic systems. We analytically prove that the system under quench dynamics displays dynamical quantum phase transitions and further provide numerical verification by computing the polynomial expansion of the Loschmidt echo. Our results may provide insight into the realization of electronic transport in experiments.

Longitudinal failure criterion considering local fiber misplacement for UD composites

A new longitudinal failure criterion for unidirectional fibre reinforced (UD) composites is proposed, incorporating the contribution of local fibre misalignment to fibre failure. The failure function was established based on the quadratic polynomial expansion of stress components in the local coordinate system parallel to the misaligned fibres, with coefficients rationally determined based on the basic strength parameters of UD composites. The local stress components were used to determine the failure status associated with fibre compression, while the tensile load along fibres was assumed to have negligible influence on the development of the fibre misalignment. The predicted failure envelopes agreed well with the experimental data for different stress combinations, particularly in the regimes where the fibre failure was influenced by transverse stresses. This validates the predictive ability of the proposed criteria and highlights the importance of the contribution of matrix-failure-related stresses to the fibre-dominated failure behaviour of fibre reinforced composites.

POLYNOMIAL SOLUTIONS FOR THE KOLMOGOROV–WIENER PREDICTION OF MODELED SMOOTHED HEAVY-TAIL PROCESS

Nowadays telecommunication traffic in systems with data packet transfer is considered as a heavy-tail random process. In a couple of rather simple models traffic is considered to be stationary one. In our recent papers we generated modeled heavy-tail data, which is based on the smoothing of the fractional Gaussian noise. In particular, the applicability if the continuous Kolmogorov–Wiener filter to the prediction of such data was investigated, the corresponding Wiener–Hopf integral equation was solved on the basis of the truncated Walsh function expansion. However, a question occurs – may another truncated orthogonal function expansion be applied to the problem under consideration? So, the corresponding investigation may be an actual question. In our recent papers we investigated theoretical fundamentals of the Kolmogorov–Wiener filter construction for different models, in particular, on the basis of the truncated polynomial expansion method and on the basis of the truncated trigonometric Fourier series expansion method. In this paper we restrict ourselves to the investigation of the applicability of the truncated polynomial expansion method to the problem under consideration, the corresponding method is based on the Сhebyshev polynomials of the first kind. The applicability of another polynomial or trigonometric expansions to the problem under consideration may be discussed in other papers. The aim of the work is to investigate the applicability of the Galerkin method based on the Chebyshev polynomials of the first kind to the Kolmogorov–Wiener prediction of smoothed heavy-tail data. The methodology consists in the solving of the Wiener–Hopf integral equation on the basis of the truncated polynomial expansion method which is based on the Chebyshev polynomials of the first kind. The scientific novelty consists in the proof of the fact that the Galerkin method based on the Chebyshev polynomials of the first kind may be applied to the Kolmogorov–Wiener prediction of smoothed heavy-tail data. The conclusions are as follows. The truncated polynomial expansion method based on the Chebyshev polynomials of the first kind may give reliable results in the framework of the Kolmogorov–Wiener prediction of smoothed heavy-tail data.

On the direct quantization of Maxwell field

In this paper, we apply the generalized integration constants (GCI) method (Belhadi 2023 https://arxiv.org/abs/2303.08236), in field theory to quantize Maxwell and the Klein–Gordon free fields. The study is performed in both position and momentum spaces, to obtain the equal-time Dirac brackets among the fields and their conjugate momenta. The idea is to compute the brackets near the initial instant using the Taylor polynomial expansion, and then deduce directly their expressions at any later time. In the case of the Maxwell field, the interdependence of the field components (constraints) requires the use of the Helmholtz theorem to separate the transversal and longitudinal parts. Our work finishes with the study of the O(3) nonlinear sigma model using the GCI approach.

Investigation of Viscoelastic Guided Wave Properties in Anisotropic Laminated Composites Using a Legendre Orthogonal Polynomials Expansion-Assisted Viscoelastodynamic Model.

This study investigates viscoelastic guided wave properties (e.g., complex-wavenumber-, phase-velocity-, and attenuation-frequency relations) for multiple modes, including different orders of antisymmetric, symmetric, and shear horizontal modes in viscoelastic anisotropic laminated composites. To obtain those frequency-dependent relations, a guided wave characteristic equation is formulated based on a Legendre orthogonal polynomials expansion (LOPE)-assisted viscoelastodynamic model, which fuses the hysteretic viscoelastic model-based wave dynamics and the LOPE-based mode shape approximation. Then, the complex-wavenumber-frequency solutions are obtained by solving the characteristic equation using an improved root-finding algorithm, which leverages coefficient matrix determinant ratios and our proposed local tracking windows. To trace the solutions on the dispersion curves of different wave modes and avoid curve-tracing misalignment in regions with phase-velocity curve crossing, we presented a curve-tracing strategy considering wave attenuation. With the LOPE-assisted viscoelastodynamic model, the effects of material viscosity and fiber orientation on different guided wave modes are investigated for unidirectional carbon-fiber-reinforced composites. The results show that the viscosity in the hysteresis model mainly affects the frequency-dependent attenuation of viscoelastic guided waves, while the fiber orientation influences both the phase-velocity and attenuation curves. We expect the theoretical work in this study to facilitate the development of guided wave-based techniques for the NDT and SHM of viscoelastic anisotropic laminated composites.

Probing New physics with high-redshift quasars: axions and non-standard cosmology

The Hubble diagram of quasars, as candidates to “standardizable” candles, has been used to measure the expansion history of the Universe at late times, up to very high redshifts (z ∼ 7). It has been shown that this history, as inferred from the quasar dataset, deviates at ≳ 3σ level from the concordance (ΛCDM) cosmology model preferred by the cosmic microwave background (CMB) and other datasets. In this article, we investigate whether new physics beyond ΛCDM (BΛCDM) or beyond the Standard Model (BSM) could make the quasar data consistent with the concordance model. We first show that an effective redshift-dependent relation between the quasar UV and X-ray luminosities, complementing previous phenomenological work in the literature, can potentially remedy the discrepancy. Such a redshift dependence can be realized in a BSM model with axion-photon conversion in the intergalactic medium (IGM), although the preferred parameter space is in tension with various other astrophysical constraints on axions, at a level depending on the specific assumptions made regarding the IGM magnetic field. We briefly discuss a variation of the axion model that could evade these astrophysical constraints. On the other hand, we show that models beyond ΛCDM such as one with a varying dark energy equation of state (wCDM) or the phenomenological cosmographic model with a polynomial expansion of the luminosity distance, cannot alleviate the tension. The code for our analysis, based on emcee [1] and corner.py [2], is publicly available at github.com/ChenSun-Phys/high_z_candles.

Energy stability of magnetohydrodynamic flow in channels and ducts

We study the energy stability of pressure-driven laminar magnetohydrodynamic flow in a rectangular duct with a transverse homogeneous magnetic field and electrically insulating walls. For sufficiently strong fields, the laminar velocity distribution has a uniform core and convex Hartmann and Shercliff boundary layers on the walls perpendicular and parallel to the magnetic field. The problem is discretized by a double expansion in Chebyshev polynomials in the cross-stream coordinates. The linear eigenvalue problem for the critical Reynolds number depends on the streamwise wavenumber, Hartmann number and the aspect ratio. We consider the limits of small and large aspect ratios in order to compare with stability models based on one-dimensional base flows. For large aspect ratios, we find good numerical agreement with results based on the quasi-two-dimensional approximation. The lift-up mechanism dominates in the limit of a zero streamwise wavenumber and provides a linear dependence between the critical Reynolds and Hartmann numbers in the duct. As the aspect ratio is reduced away from unity, the duct results converge to Orr's original energy stability result for spanwise uniform perturbations imposed on the plane Poiseuille base flow. We also examine different possible symmetries of eigenmodes as well as the purely hydrodynamic case in the duct geometry.

Cluster cosmology redux: a compact representation for the halo mass function

ABSTRACT Groups and clusters of galaxies imprint coherent, arcminute-scale features across the spectrophotometric sky, especially optical-IR clusters of galaxies, spectral distortions in the cosmic microwave background, and extended sources of X-ray emission. The space–time density of the host dark matter halo population – the halo mass function (HMF) – is a common theoretical basis for modelling such observable features. We explore a compact representation – a dual-quadratic (DQ-HMF) form – that features readily interpretable parameters representing polynomial expansions of the space–time number density surface, first in terms of log-mass, then in redshift. The DQ-HMF form fits Mira-Titan N-body emulator expectations for halo masses $10^{13.7-14.5} \, h^{-1}\, {\rm M}_\odot$ over redshifts, 0.1 < z < 1.5 to within $\sim \! 5~{{\ \rm per\ cent}}$. We provide best-fitting parameters for a Planck 2018 cosmology and demonstrate model self-similarity in the Ωm−$\, \sigma _8$ plane. Convolving with a minimal mass–observable relation (MOR) yields closed-form expressions for counts, mean mass, and mass variance of cluster samples characterized by an observable property. Performing information-matrix forecasts of potential parameter constraints from existing and future surveys, we demonstrate the potential for percent-level constraints on model parameters by an LSST-like optical cluster survey of 300 000 clusters and a richness–mass variance of 0.32. Even better constraints could potentially be achieved by a survey with one-tenth the sample size but with a reduced selection property variance of 0.12. Potential benefits and extensions to the basic MOR parametrization are discussed.

On the orthogonality of Atkin-like polynomials and orthogonal polynomial expansion of generalized Faber polynomials

On the orthogonality of Atkin-like polynomials and orthogonal polynomial expansion of generalized Faber polynomials

Uncertainty Benchmarks for Time-Dependent Transport Problems

Verification solutions for uncertainty quantification (UQ) are presented for time-dependent transport problems where , the scattering ratio, is uncertain. The method of polynomial chaos expansions is employed for quick and accurate calculation of the quantities of interest (QoIs), and uncollided solutions are used to treat part of the uncertainty calculation analytically. We find that approximately six moments in the polynomial expansion are required to represent the solutions to these problems accurately. Additionally, the results show that if the uncertainty interval spans c = 1, which means it is uncertain whether the system is multiplying or not, the confidence interval will grow in time. Finally, since the QoI is a strictly increasing function, the percentile values are known and can be used to verify the accuracy of the expansion. These results can be used to test UQ methods for time-dependent transport problems.

Vibroacoustic analysis of submerged fluid-filled cylindrical shell

This study presents a comprehensive Chebyshev-Fourier spectral method for analyzing the free and forced vibroacoustic characteristics of submerged fluid-filled cylindrical shells in heavy fluids. The vibroacoustic interaction system of submerged cylindrical shells consists of three discrete regions: The Sanders shell theory is applied for the shell structure, whereas the wave equation is used to determine the fluid inside. The external fluid pressure on the shell is described using the Helmholtz boundary integral equation. The governing equations for the shell's motion, the formulation of the external fluid pressure on the fluid-structure-interaction (FSI) surface using first-class Chebyshev polynomial and Fourier series expansions, and the internal sound pressure employing two-dimensional Chebyshev polynomial and Fourier series expansions, are seamlessly integrated with Gauss-Chebyshev-Lobatto sampling techniques. Artificial springs are implemented to replicate the boundary conditions of the cylindrical shell. The precision of the proposed methodology is corroborated by comparisons with existing literature and finite element method (FEM) results. Furthermore, the research examines how the properties of the acoustical medium on the vibroacoustic responses of the cylindrical shell. Results indicate that the presence of air inside the shell slightly influence on its acoustic radiation. However, replacing the internal medium with water significantly alters the shell's original radiation characteristics.

Connecting descent and peak polynomials

A permutation $\sigma=\sigma_1 \sigma_2 \cdots \sigma_n$ has a descent at $i$ if $\sigma_i>\sigma_{i+1}$. A descent $i$ is called a peak if $i>1$ and $i-1$ is not a descent. The size of the set of all permutations of $n$ with a given descent set is a polynomials in $n$, called the descent polynomial. Similarly, the size of the set of all permutations of $n$ with a given peak set, adjusted by a power of $2$ gives a polynomial in $n$, called the peak polynomial. In this work we give a unitary expansion of descent polynomials in terms of peak polynomials. Then we use this expansion to give an interpretation of the coefficients of the peak polynomial in a binomial basis, thus giving a constructive proof of the peak polynomial positivity conjecture.

Optimizations of the Accelerated Axial Polynomial Method of Characteristics of the APOLLO3® Deterministic Neutron Transport Code

For some years now, the TDT (two- and three-dimensional transport) solver of the APOLLO3® deterministic neutron transport code has been able to perform lattice calculations on three-dimensional extruded and unstructured geometries. A polynomial expansion of the angular flux has been implemented to better describe the flux gradient axially to reduce the number of computational meshes required to reach a given accuracy. Then the polynomial approximation was extended to macroscopic cross sections to perform evolution calculations. Besides these transport schemes, synthetic acceleration has also been implemented, relying on double PN approximations of the angular flux on the boundaries of the spatial regions. The solver has already introduced several techniques to reduce the transport and memory footprint; for example, for the storage of the surfaces crossed by a trajectory or the classification of chords. In this paper, new optimizations are presented. One deals with how monomials of the polynomial basis are integrated along trajectories. Another one concerns the computation of the source term of the transmission equation in the case of polynomial cross sections. The last optimization exploits the fact that, along horizontal trajectories, the flux and the cross sections are constant to speed up the sweep algorithm. Calculations on 5 × 5 and 7 × 7 pressurized water reactor assemblies were performed to assess the gains of these recently developed strategies. The results show good improvements both in computing time and in memory footprint reductions.

Polynomial expansion method for full wave three-dimensional analysis of dielectric waveguides and periodic structures

In summary, the utilization of Legendre polynomial expansion in the modal analysis of stratified dielectric layers with doubly periodic permittivity profiles offers a departure from conventional methods. This novel approach, grounded in the analytical projection of Maxwell’s equations onto the Hilbert space defined by Legendre polynomials, results in well-behaved algebraic equations. These equations, in turn, facilitate the derivation of propagation constants and electromagnetic field profiles, circumventing issues related to numerical instability and oscillatory behavior. Moreover, the method's adaptability to extend its application to nonperiodic dielectric waveguides through the periodic repetition concept further underscores its versatility and potential impact in electromagnetic field analysis. Finally, to validate the proposed method, we conducted a comparative analysis of three standard test cases against previously reported results in the literature. The comparison showcased a high level of agreement, affirming the accuracy and efficacy of the presented approach.

Model order reduction for discrete time-delay systems based on Laguerre function expansion

We investigate model order reduction for discrete time-delay systems via orthogonal polynomial expansion in the frequency domain. The transfer function of systems is expanded in the framework of Laguerre function basis. We show that Laguerre coefficients of the states satisfy a linear system and reduced models generated by projection methods can preserve some Laguerre coefficients of the original systems. We prove that the subspace spanned by Laguerre coefficients is exactly a high order Krylov subspace, thereby leading to an efficient computation of projection matrices and a more accurate coefficient-matching property in the two-sided framework. Further the Laguerre expansion of systems is employed to enable an approximate but fast execution of balanced truncation for discrete time-delay systems. Specifically, Gramians are approximated based on the derived Laguerre coefficients, which circumvents the main bottleneck of balanced truncation methods. Finally, two numerical examples are simulated to demonstrate the feasibility and effectiveness of the proposed methods.

Riemann-Hilbert hierarchies for hard edge planar orthogonal polynomials

abstract: We obtain a full asymptotic expansion for orthogonal polynomials with respect to weighted area measure on a Jordan domain $\mathscr{D}$ with real-analytic boundary. The weight is fixed and assumed to be real-analytically smooth and strictly positive, and for any given precision $\varkappa$, the expansion holds with an $\mathrm{O}(N^{-\varkappa-1})$ error in $N$-dependent neighborhoods of the exterior region as the degree $N$ tends to infinity. The main ingredient is the derivation and analysis of Riemann-Hilbert hierarchies---sequences of scalar Riemann-Hilbert problems---which allows us to express all higher order correction terms in closed form. Indeed, the expansion may be understood as a Neumann series involving an explicit operator. The expansion theorem leads to a semiclassical asymptotic expansion of the corresponding hard edge probability wave function in terms of distributions supported on $\partial\mathscr{D}$.

Random Green’s Function Method for Large-Scale Electronic Structure Calculation

We report a linear-scaling random Green’s function (rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H2O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H2O clusters at T = 0 K can reach an error of ∼ 1 meV per H2O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.

Endpoint estimates for harmonic analysis operators associated with Laguerre polynomial expansions

AbstractIn this paper, we give a criterion to prove boundedness results for several operators from the Hardy‐type space to and also from to the space of functions of bounded mean oscillation , with respect to the probability measure on when is a multi‐index in . We shall apply it to establish endpoint estimates for Riesz transforms, maximal operators, Littlewood–Paley functions, multipliers of the Laplace transform type, fractional integrals, and variation operators in the Laguerre setting.

Orthogonal polynomial expansions for the valuation of options under the stochastic volatility models with stochastic correlation

This work provides a new method for pricing options under the generalized stochastic volatility models with Jacobi stochastic correlation. Our method is based on the observation that the generalized models belong to the class of polynomial diffusions and therefore the option prices can be efficiently computed via orthogonal polynomial expansions. We take the Heston and Schöbel-Zhu models with stochastic correlation as two specific examples and are able to derive the analytical formulas for the option prices. We also illustrate the accuracy of the proposed method through a number of numerical experiments.

On the Thermally Induced Interfacial Behavior of Thin Two-Dimensional Hexagonal Quasicrystal Films with an Adhesive Layer

The mechanical response of a quasicrystal thin film is strongly affected by an adhesive layer along the interface. In this paper, a theoretical model is proposed to study a thin two-dimensional hexagonal quasicrystal film attached to a half-plane substrate with an adhesive layer, which undergoes a thermally induced deformation. A perfect non-slipping contact condition is assumed at the interface by adopting the membrane assumption. An analytical solution to the problem is obtained by constructing governing integral–differential equations for both single and multiple films in terms of interfacial shear stresses that are reduced to a linear algebraic system via the series expansion of Chebyshev polynomials. The solution is compared to that without adhesive layers, and the effects of the aspect ratio of films, material mismatch, and the adhesive layer, as well as the interaction between films, are discussed in detail. It is found that the adhesive layer can soften the localized stress concentration. This study is instructive to the accurate safety assessment and functional design of a quasicrystal film system.