Diagonal Function Research Articles

Thermoacoustics of multi-burner combustors with plenum and chamber cross-talk

This paper presents a thermoacoustic model in the frequency domain for multi-burner combustors composed of multiple plenum-nozzle-chamber modules with cross-talk at both plenums and chambers. The thermoacoustic model is given as a multi-input-multi-output transfer function matrix from heat fluctuation input vector to velocity fluctuation output vector. The discrete rotational symmetry of a multi-burner combustor is exploited to obtain a simple, diagonal transfer function matrix. The resonance equation, whose roots are resonances, is given as a set of scalar equations. The explicit form of the resonance equation reveals how combustor and operating parameters influence the resonances, and which parameters are significant for the longitudinal-azimuthal coupling and the plenum-chamber coupling. The diagonalization also allows us to show constructively that every mode shape is a Bloch wave and those mode shapes can characterize all resonances in terms of longitudinal and azimuthal mode numbers. At a low frequency, the resonance equation can be approximated as biquadratic and quadratic equations for a typical combustor geometry and operating conditions. A bunch of low-frequency resonances, peculiar to multi-burner combustors with cross-talk, are interpreted as azimuthal Helmholtz modes in which plenums/chambers and cross-talk ducts serve as the cavities and necks, respectively. Our combustor model is combined with a diagonal flame transfer function matrix for combustion instability analysis and we develop a set of scalar stability equations whose roots determine combustion instability. This stability equation reveals how combustion instability depends on the combustor geometry and operating parameters. Several numerical examples are presented to validate our thermoacoustic model and theoretical findings.

QCD sum rules for parameters of the B -meson distribution amplitudes

We obtain new estimates for the parameters $\lambda_{E}^2$, $\lambda_H^2$ and their ratio $\mathcal{R} = \lambda_{E}^2/\lambda_H^2$, which appear in the second moments of the $B$-meson light-cone distribution amplitudes defined in the heavy-quark effective field theory. The computation is based on two-point QCD sum rules for the diagonal correlation function and includes all contributions up to mass dimension seven in the operator-product expansion. For the ratio we get $\mathcal{R} = (0.1 \pm 0.1)$ with $\lambda_H^2 = (0.15 \pm 0.05) \, \text{GeV}^2$ and $\lambda_E^2 = (0.01 \pm 0.01) \, \text{GeV}^2$.

The Study of Monotonic Core Functions and Their Use to Build RNS Number Comparators

A non-positional residue number system (RNS) enjoys particularly efficient implementation of addition and multiplication, but non-modular arithmetic operations in RNS-like number comparison are known to be difficult. In this paper, a new technique for designing comparators of RNS numbers represented in an arbitrary moduli set is presented. It is based on using the core function for which it was shown that it must be monotonic to allow for RNS number comparison. The conditions of the monotonicity of the core function were formulated, which also ensured the minimal range of the core function (essential to obtain the best characteristics of the comparator). The best choice is a core function in which only one coefficient corresponding to the largest modulus is set to 1 whereas all other coefficients are set to 0. It is also shown that the already known diagonal function is nothing else but the special case of the core function with all coefficients set to 1. Performance evaluation suggests that the new comparator uses less hardware and in some cases also introduces smaller delay than its counterparts based on diagonal function. The potential applications of the new comparator include some recently developed homomorphic encryption algorithms implemented using RNS.

RNS Number Comparator Based on a Modified Diagonal Function

Number comparison has long been recognized as one of the most fundamental non-modular arithmetic operations to be executed in a non-positional Residue Number System (RNS). In this paper, a new technique for designing comparators of RNS numbers represented in an arbitrary moduli set is presented. It is based on a newly introduced modified diagonal function, whose strictly monotonic properties make it possible to replace the cumbersome operations of finding the remainder of the division by a large and awkward number with significantly simpler computations involving only a power of 2 modulus. Comparators of numbers represented in sample RNSs composed of varying numbers of moduli and offering different dynamic ranges, designed using various methods, were synthesized for the 65 nm technology. The experimental results suggest that the new circuits enjoy a delay reduction ranging from over 11% to over 75% compared to the fastest circuits designed using existing methods. Moreover, it is achieved using less hardware, the reduction of which reaches over 41%, and is accompanied by significantly reduced power-consumption, which in several cases exceeds 100%. Therefore, it seems that the presented method leads to the design of the most efficient current hardware comparators of numbers represented using a general RNS moduli set.

Three-dimensional slowly rotating black hole in Einstein-power-Maxwell theory

A three-dimensional slowly rotating black hole solution in the presence of negative cosmological constant in the Einstein-power-Maxwell theory is studied. It is shown that in the small rotation limit, the electric field, diagonal metric function and thermodynamic properties are the same as for static case, whereas the small rotation gives in addition a nondiagonal metric function and magnetic field which are also small. For these functions cased by rotation of black hole, exact integral solution and analytic asymptotic solution were obtained.

Stability and stabilisation of heterogeneous switched systems with arbitrary switching constraints

Heterogeneous switched systems are a class of special hybrid systems containing different dimensional subsystems. To research the stability and stabilisation of heterogeneous switched systems, a novel approach is presented to design static output feedback controllers for asymptotically stabilising heterogeneous continuous-time and discrete-time switched systems under arbitrary switching rules. This approach, which has less computation complexity, includes a main diagonal square block common Lyapunov function for stability analysis of heterogeneous switched systems and a main diagonal square block common homogeneous polynomial Lyapunov function for stabilisation of heterogeneous switched systems. In the study of stabilisation of heterogeneous switched systems, a two-step iterative method including two convex optimisation functions is established by sums of squares to express existing properties condition and solvable properties condition for determining block controller. At last, a practical example and two numerical examples are presented to show the feasibility of the proposed results.

The soft quark Sudakov

There has been recent interest in understanding the all loop structure of the subleading power soft and collinear limits, with the goal of achieving a systematic resummation of subleading power infrared logarithms. Most of this work has focused on subleading power corrections to soft gluon emission, whose form is strongly constrained by symmetries. In this paper we initiate a study of the all loop structure of soft fermion emission. In mathcal{N} = 1 QCD we perform an operator based factorization and resummation of the associated infrared logarithms using the formalism introduced in [1], and prove that they exponentiate into a Sudakov due to their relation to soft gluon emission. We verify this result through explicit calculation to mathcal{O}left({alpha}_s^3right) . We show that in QCD, this simple Sudakov exponentiation is violated by endpoint contributions proportional to (CA−CF)n which contribute at leading logarithmic order. Combining our mathcal{N} = 1 result and our calculation of the endpoint contributions to mathcal{O}left({alpha}_s^3right) , we conjecture a result for the soft quark Sudakov in QCD, a new all orders function first appearing at subleading power, and give evidence for its universality. Our result, which is expressed in terms of combinations of cusp anomalous dimensions in different color representations, takes an intriguingly simple form and also exhibits interesting similarities to results for large-x logarithms in the off diagonal splitting functions.

Classification of Moduli Sets for Residue Number System With Special Diagonal Functions

The paper presents algorithms for the generation of Residue Number System (RNS) triples with $SQ=2^{k}-1$ and quadruples with $SQ=2^{k}$ for some k. Triples and quadruples allow us to design efficient hardware implementations of non-modular operations in RNS such as division, sign detection, comparison of numbers, reverse conversion with using of a diagonal function from requiring division with the remainder by the diagonal module SQ. Division with a remainder in the general case is the most complex arithmetic operation in computer technology. However, the consideration of special cases can significantly simplify this operation and increase the efficiency of hardware implementation. We show that there are 5573 good RNS triples (2301 even and 2372 odd) with elements less than 10 000, as the values of SQ vary from $2^{5}-1$ to $2^{27}-1$ . In contrast, RNS quadruples with $SQ=2^{k}$ seem to be quite rare. Restricting our search to sums of the elements in a quadruple less than 4000 we find that exactly 31 such quadruples exist. Their values of SQ vary between 220 and 230 with always even exponent. We suggest the measure of RNS balance and find perfectly balanced RNS among triples according to this measure. We demonstrate the advantages of more balanced quadruples by means of hardware implementation.

Generalized Spectrum of Second Order Differential Operators

We analyze the spectrum of the operator $\Delta^{-1} [\nabla \cdot (K\nabla u)]$, where $\Delta$ denotes the Laplacian and $K=K(x,y)$ is a symmetric tensor. Our main result shows that this spectrum can be derived from the spectral decomposition $K=Q \Lambda Q^T$, where $Q=Q(x,y)$ is an orthogonal matrix and $\Lambda=\Lambda(x,y)$ is a diagonal matrix. More precisely, provided that $K$ is continuous, the spectrum equals the convex hull of the ranges of the diagonal function entries of $\Lambda$. The involved domain is assumed to be bounded and Lipschitz, and both homogeneous Dirichlet and homogeneous Neumann boundary conditions are considered. We study operators defined on infinite dimensional Sobolev spaces. Our theoretical investigations are illuminated by numerical experiments, using discretized problems. The results presented in this paper extend previous analyses which have addressed elliptic differential operators with scalar coefficient functions. Our investigation is motivated by both preconditioning issues (efficient numerical computations) and the need to further develop the spectral theory of second order PDEs (core analysis).

Optimal decisions and expected values in two player zero sum games with diagonal game matrixes-explicit functions, general proofs and effects of parameter estimation errors

In this paper the two player zero sum games with diagonal game matrixes TPZSGD are analyzed Many important applications of this particular class of games are found in military decision problems in customs and immigration strategies and police work Explicit functions are derived that give the optimal frequences of different decisions and the expected results of relevance to the different decision makers Arbitrary numbers of decision alternatives are covered It is proved that the derived optimal decision frequency formulas correspond to the unique optimization results of the two players It is proved that the optimal solutions for both players always lead to a unique completely mixed strategy Nash equilibrium For each player the optimal frequency of a particular decision is strictly greater than and strictly less than With comparative statics analyses the directions of the changes of optimal decision frequences and expected game values as functions of changes in different parameter values are determined The signs of the optimal changes of the decision frequences of the different players are also determined as functions of risk in different parameter values Furthermore the directions of changes of the expected optimal value of the game are determined as functions of risk in the different parameter values Finally some of the derived formulas are used to confirm earlier game theory results presented in the literature It is demonstrated that the new functions can be applied to solve common military problems

Dual switched positive systems – a less conservative condition for diagonal quadratic stability

ABSTRACTThe key result of the paper exploits the duality of arbitrary switching positive linear systems, in order to derive a sufficient condition for the existence and construction of diagonal quadratic copositive Lyapunov functions. This condition is less conservative than a result also relying on duality, recently reported in the literature; our condition operates with milder inequalities than those required by the existing result. We prove that the algebraic relaxation of these inequalities is related to a relaxation of the primal and dual systems’ properties, referring to the trajectory growth rates, or, equivalently, to the spectra of the associated column representatives. The existing result is incorporated into the new one as a particular case, and the advantages offered by the latter are illustrated by a relevant numerical example. Our developments focus on the continuous-time case, so as to permit direct comparisons with the existing result; some brief discussion suggests how the proposed approach can be applied, mutatis mutandis, to discrete-time dynamics.

Connection formulas for the generalized Ising correlation functions

We derive and prove the connection formulas for the generalized diagonal Ising model correlation functions.

Construction of Residue Number System Using Hardware Efficient Diagonal Function

The residue number system (RNS) is a non-positional number system that allows one to perform addition and multiplication operations fast and in parallel. However, because the RNS is a non-positional number system, magnitude comparison of numbers in RNS form is impossible, so a division operation and an operation of reverse conversion into a positional form containing magnitude comparison operations are impossible too. Therefore, RNS has disadvantages in that some operations in RNS, such as reverse conversion into positional form, magnitude comparison, and division of numbers are problematic. One of the approaches to solve this problem is using the diagonal function (DF). In this paper, we propose a method of RNS construction with a convenient form of DF, which leads to the calculations modulo 2 n , 2 n − 1 or 2 n + 1 and allows us to design efficient hardware implementations. We constructed a hardware simulation of magnitude comparison and reverse conversion into a positional form using RNS with different moduli sets constructed by our proposed method, and used different approaches to perform magnitude comparison and reverse conversion: DF, Chinese remainder theorem (CRT) and CRT with fractional values (CRTf). Hardware modeling was performed on Xilinx Artix 7 xc7a200tfbg484-2 in Vivado 2016.3 and the strategy of synthesis was highly area optimized. The hardware simulation of magnitude comparison shows that, for three moduli, the proposed method allows us to reduce hardware resources by 5.98–49.72% in comparison with known methods. For the four moduli, the proposed method reduces delay by 4.92–21.95% and hardware costs by twice as much by comparison to known methods. A comparison of simulation results from the proposed moduli sets and balanced moduli sets shows that the use of these proposed moduli sets allows up to twice the reduction in circuit delay, although, in several cases, it requires more hardware resources than balanced moduli sets.

Ion Impact Induced Ultrafast Electron Dynamics in Finite Graphene‐Type Hubbard Clusters

Strongly correlated systems of fermions have an interesting phase diagram arising from the Hubbard gap. Excitation across the gap leads to the formation of doubly occupied lattice sites (doublons) which offers interesting electronic and optical properties. Moreover, when the system is driven out of equilibrium interesting collective dynamics may arise that are related to the spatial propagation of doublons. Here, a novel mechanism that was recently proposed by the authors [Balzer et al., Phys. Rev. Lett. 121, 267602 (2018)] is verified by exact diagonalization and nonequilibrium Green functions (NEGF) simulations – fermionic doublon creation by the impact of energetic ions. The formation of a nonequilibrium steady state with homogeneous doublon distribution is reported. The effect should be particularly important for correlated finite systems, such as graphene nanoribbons, and directly observable with fermionic atoms in optical lattices. It is demonstrated that doublon formation and propagation in correlated lattice systems can be accurately simulated with NEGF. In addition to two‐time results, single‐time results within the generalized Kadanoff–Baym ansatz (GKBA) with Hartree–Fock propagators (HF‐GKBA) is presented. Finally systematic improvements of the GKBA that use correlated propagators (correlated GKBA) and a correlated initial state are discussed.

Reflected maxmin copulas and modeling quadrant subindependence

Copula models have become popular in different applications, including modeling shocks, in view of their ability to describe better the dependence concepts in stochastic systems. The class of maxmin copulas was recently introduced by Omladič and Ružić [21]. It extends the well known classes of Marshall–Olkin and Marshall copulas by allowing the external shocks to have different effects on the two components of the system. By a reflection (flip) in one of the variables we introduce a new class of bivariate copulas called reflected maxmin (RMM) copulas. We explore their properties and show that symmetric RMM copulas relate to general RMM copulas similarly as semilinear copulas relate to Marshall copulas. We transfer that relation also to maxmin copulas. We also characterize possible diagonal functions of symmetric RMM copulas.

The reproducing kernel of $\mathcal H^2$ and radial eigenfunctions of the hyperbolic Laplacian

In the paper we characterize the reproducing kernel $\mathcal {K}_{n,h}$ for the Hardy space $\mathcal {H}^2$ of hyperbolic harmonic functions on the unit ball $\mathbb {B}$ in $\mathbb {R}^n$. Specifically we prove that \[ \mathcal {K}_{n,h}(x,y) = \sum _{\alpha =0}^\infty S_{n,\alpha }(\lvert x\rvert )S_{n,\alpha }(\lvert y\rvert ) Z_\alpha (x,y), \] where the series converges absolutely and uniformly on $K\times \mathbb {B}$ for every compact subset $K$ of $\mathbb {B}$. In the above, $S_{n,\alpha }$ is a hypergeometric function and $Z_\alpha $ is the reproducing kernel of the space of spherical harmonics of degree α. In the paper we prove that \[ 0\le \mathcal K_{n,h}(x,y) \le \frac {C_n}{(1-2\langle x,y\rangle + \lvert x \rvert^2 \lvert y \rvert^2)^{n-1}}, \] where $C_n$ is a constant depending only on $n$. It is known that the diagonal function $\mathcal K_{n,h}(x,x)$ is a radial eigenfunction of the hyperbolic Laplacian $\varDelta_h $ on $\mathbb{B} $ with eigenvalue $\lambda _2 = 8(n-1)^2$. The result for $n=4$ provides motivation that leads to an explicit characterization of all radial eigenfunctions of $\varDelta_h $ on $\mathbb{B} $. Specifically, if $g$ is a radial eigenfunction of $\varDelta_h $ with eigenvalue $\lambda _\alpha = 4(n-1)^2\alpha (\alpha -1)$, then \[ g(r) = g(0) \frac {p_{n,\alpha }(r^2)}{(1-r^2)^{(\alpha -1)(n-1)}}, \] where $p_{n,\alpha }$ is again a hypergeometric function. If α is an integer, then $p_{n,\alpha }(r^2)$ is a polynomial of degree $2(\alpha -1)(n-1)$.

Algebra of Countably Functions and Theorems of Completeness

An algebraic approach to the theory of countable functions is given. Compositions (superpositions) of functions are used instead of recursions. Arithmetic and analytic algorithms are defined. All closed sets are founded. Mathematically precise definitions of logic algorithms with quantifiers of existence and universality are given. Logic algorithm for fast-growing function is built as example. Classification of functions is given. There are non-computable functions. These functions are fictitious (useless) and their set is continual. The set of computable functions is countable. Incompleteness of disjunction and negation, conjunction and negation, of Pierce, Sheffer and diagonal Webb functions is proved. The completeness of the set of one-place functions and any all-valued essential function (Slupecki theorem) is proved for computable functions. Existence of generators of all computable functions is proved too.

Doublon Formation by Ions Impacting a Strongly Correlated Finite Lattice System.

Strongly correlated systems of fermions have a number of exciting collective properties. Among them, the creation of a lattice that is occupied by doublons, i.e., two quantum particles with opposite spins, offers interesting electronic properties. In the past a variety of methods have been proposed to control doublon formation, both, spatially and temporally. Here, a novel mechanism is proposed and verified by exact diagonalization and nonequilibrium Green functions simulations-fermionic doublon creation by the impact of energetic ions. We report the formation of a nonequilibrium steady state with homogeneous doublon distribution. The effect should be particularly important for strongly correlated finite systems, such as graphene nanoribbons, and directly observable with fermionic atoms in optical lattices.

Construction of a Correct Algorithm and Spatial Neural Network for Recognition Problems with Binary Data

Conditions under which it is possible to design a correct algorithm and a six-level spatial neural network reproducing the computations performed by this algorithm for recognition problems with binary data (Ω-regular problems) are found. A distinctive feature of this network is the use of diagonal activation functions in its internal layers, which significantly simplify intermediate computations in the inner and outer loops. Given an Ω-regular problem, the network sequentially computes the rows of the classification matrix for the test sample objects. The computational process (i.e., the inner loop) for each test object consists inside the elementary 3-level network (i.e., μ-block) of a single iteration determined by a single object of the training set. The proposed approach to the neural network construction does not rely on the conventional approach based on the minimization of a functional; rather, it is based on the operator theory developed by Zhuravlev for solving recognition and classification problems.

Coherent deeply virtual Compton scattering off He4

Coherent Deeply virtual Compton scattering off the $^4$He nucleus is studied in impulse approximation. A convolution formula for the nuclear Generalized Parton Distribution (GPD) is derived in terms of the $^4$He non-diagonal spectral function and of the GPD of the struck nucleon. A model of the nuclear non diagonal spectral function, based on the momentum distribution corresponding to the AV18 nucleon-nucleon interaction, is used in the actual calculation. Typical impulse approximation results are reproduced, in proper limits, for the nuclear form factor and for nuclear parton distributions. The nuclear generalized parton distribution and the Compton form factor are evaluated using, as nucleonic ingredient, a well known generalized parton distribution model. An overall very good agreement is found with the data recently published by the EG6 experiment at the Jefferson Laboratory (JLab). More refined nuclear calculations are addressed and will be necessary for the expected improved accuracy of the next generation of experiments at JLab, with the 12 GeV electron beam and high luminosity.