On the coefficients of asymptotic expansion for the harmonic number by Ramanujan

Let \(H_n=\sum _{k=1}^{n}\frac{1}{k}\) be the \(n\)th harmonic number. In this paper, we establish a new asymptotic expansion of \(H_n\). By using the result obtained, we derive a recurrence relation for determining the coefficients of Ramanujan’s asymptotic expansion for the harmonic number. Also, we establish asymptotic expansion of \(H_n\) in terms of \(n(n+1)+1/3\).     

On some constacyclic codes over $$\mathbb {Z}_{4}\left[ u\right] /\left\langle u^{2}-1\right\rangle $$, their $$\mathbb {Z}_4$$Z4 images,and new codes

In this paper, we study \(\lambda \)-constacyclic codes over the ring \(R=\mathbb {Z}_4+u\mathbb {Z}_4\) where \(u^{2}=1\), for \(\lambda =3+2u\) and \(2+3u\). Two new Gray maps from R to \(\mathbb {Z}_4^{3}\) are defined with the goal of obtaining new linear codes over \(\mathbb {Z}_4\). The Gray images of \(\lambda \)-constacyclic codes over R are determined. We then conducted a computer search and obtained many \(\lambda \)-constacyclic codes over R whose \(\mathbb {Z}_4\)-images have better parameters than currently best-known linear codes over \(\mathbb {Z}_4\).     

Standard Jordan partitions with three parameters

We extend previous work on standard two-parameter Jordan partitions by Barry (Commun Algebra 43:4231–4246, 2015) to three parameters. Let \(J_r\) denote an \(r \times r\) matrix with minimal polynomial \((t-1)^r\) over a field F of characteristic p. For positive integers \(n_1\), \(n_2\), and \(n_3\) satisfying \(n_1 \le n_2 \le n_3\), the Jordan canonical form of the \(n_1 n_2 n_3 \times n_1 n_2 n_3\) matrix \(J_{n_1} \otimes J_{n_2} \otimes J_{n_3}\) has the form \(J_{\lambda _1} \oplus J_{\lambda _2} \oplus \cdots \oplus J_{\lambda _m}\) where \(\lambda _1 \ge \lambda _2 \ge \cdots \ge \lambda _m>0\) and \(\sum _{i=1}^m \lambda _i=n_1 n_2 n_3\). The partition \(\lambda (n_1,n_2,n_3:p)=(\lambda _1, \lambda _2,\ldots , \lambda _m)\) of \(n_1 n_2 n_3\), which depends on \(n_1\), \(n_2\), \(n_3\), and p, will be called a Jordan partition. We will define what we mean by a standard Jordan partition and give necessary and sufficient conditions for its existence.     

Extremal eigenvalues and eigenvectors of deformed Wigner matrices

We consider random matrices of the form \(H = W + \lambda V, \lambda \in {\mathbb {R}}^+\), where \(W\) is a real symmetric or complex Hermitian Wigner matrix of size \(N\) and \(V\) is a real bounded diagonal random matrix of size \(N\) with i.i.d. entries that are independent of \(W\). We assume subexponential decay of the distribution of the matrix entries of \(W\) and we choose \(\lambda \sim 1\), so that the eigenvalues of \(W\) and \(\lambda V\) are typically of the same order. Further, we assume that the density of the entries of \(V\) is supported on a single interval and is convex near the edges of its support. In this paper we prove that there is \(\lambda _+\in {\mathbb {R}}^+\) such that the largest eigenvalues of \(H\) are in the limit of large \(N\) determined by the order statistics of \(V\) for \(\lambda >\lambda _+\). In particular, the largest eigenvalue of \(H\) has a Weibull distribution in the limit \(N\rightarrow \infty \) if \(\lambda >\lambda _+\). Moreover, for \(N\) sufficiently large, we show that the eigenvectors associated to the largest eigenvalues are partially localized for \(\lambda >\lambda _+\), while they are completely delocalized for \(\lambda <\lambda _+\). Similar results hold for the lowest eigenvalues.     

Computing the Permanent of (Some) Complex Matrices

We present a deterministic algorithm, which, for any given \(0< \epsilon < 1\) and an \(n \times n\) real or complex matrix \(A=\left( a_{ij}\right) \) such that \(\left| a_{ij}-1 \right| \le 0.19\) for all \(i, j\) computes the permanent of \(A\) within relative error \(\epsilon \) in \(n^{O\left( \ln n -\ln \epsilon \right) }\) time. The method can be extended to computing hafnians and multidimensional permanents.     

The solution gap of the Brezis–Nirenberg problem on the hyperbolic space

We consider the positive solutions of the nonlinear eigenvalue problem \(-\Delta _{\mathbb {H}^n} u = \lambda u + u^p, \) with \(p=\frac{n+2}{n-2}\) and \(u \in H_0^1(\Omega ),\) where \(\Omega \) is a geodesic ball of radius \(\theta _1\) on \(\mathbb {H}^n.\) For radial solutions, this equation can be written as an ordinary differential equation having n as a parameter. In this setting, the problem can be extended to consider real values of n. We show that if \(2<n<4\) this problem has a unique positive solution if and only if \(\lambda \in \left( n(n-2)/4 +L^*\,,\, \lambda _1\right) .\) Here \(L^*\) is the first positive value of \(L = -\ell (\ell +1)\) for which a suitably defined associated Legendre function \(P_{\ell }^{-\alpha }(\cosh \theta ) >0\) if \(0 < \theta <\theta _1\) and \(P_{\ell }^{-\alpha }(\cosh \theta _1)=0,\) with \(\alpha = (2-n)/2\).     

Eigenvalue Counting Function for Robin Laplacians on Conical Domains

We study the discrete spectrum of the Robin Laplacian \(Q^{\Omega }_\alpha \) in \(L^2(\Omega )\), \(u\mapsto -\Delta u, \quad D_n u=\alpha u \text { on }\partial \Omega \), where \(D_n\) is the outer unit normal derivative and \(\Omega \subset {\mathbb {R}}^{3}\) is a conical domain with a regular cross-section \(\Theta \subset {\mathbb {S}}^2\), n is the outer unit normal, and \(\alpha >0\) is a fixed constant. It is known from previous papers that the bottom of the essential spectrum of \(Q^{\Omega }_\alpha \) is \(-\alpha ^2\) and that the finiteness of the discrete spectrum depends on the geometry of the cross-section. We show that the accumulation of the discrete spectrum of \(Q^\Omega _\alpha \) is determined by the discrete spectrum of an effective Hamiltonian defined on the boundary and far from the origin. By studying this model operator, we prove that the number of eigenvalues of \(Q^{\Omega }_\alpha \) in \((-\infty ,-\alpha ^2-\lambda )\), with \(\lambda >0\), behaves for \(\lambda \rightarrow 0\) as

$$\begin{aligned} \dfrac{\alpha ^2}{8\pi \lambda } \int _{\partial \Theta } \kappa _+(s)^2\mathrm {d}s +o\left( \frac{1}{\lambda }\right) , \end{aligned}$$

where \(\kappa _+\) is the positive part of the geodesic curvature of the cross-section boundary.

     

Iteration of certain exponential-like meromorphic functions

The dynamics of functions \(f_\lambda (z)= \lambda \frac{\mathrm{e}^{z}}{z+1}\ \text{ for }\ z\in \mathbb {C}, \lambda >0\) is studied showing that there exists \(\lambda ^* > 0\) such that the Julia set of \(f_\lambda \) is disconnected for \(0< \lambda < \lambda ^*\) whereas it is the whole Riemann sphere for \(\lambda > \lambda ^*\). Further, for \(0< \lambda < \lambda ^*\), the Julia set is a disjoint union of two topologically and dynamically distinct completely invariant subsets, one of which is totally disconnected. The union of the escaping set and the backward orbit of \(\infty \) is shown to be disconnected for \(0<\lambda < \lambda ^*\) whereas it is connected for \(\lambda > \lambda ^*\). For complex \(\lambda \), it is proved that either all multiply connected Fatou components ultimately land on an attracting or parabolic domain containing the omitted value of the function or the Julia set is connected. In the latter case, the Fatou set can be empty or consists of Siegel disks. All these possibilities are shown to occur for suitable parameters. Meromorphic functions \(E_n(z) =\mathrm{e}^{z}(1+z+\frac{z^2}{2!}+\cdots +\frac{z^n}{n!})^{-1}\), which we call exponential-like, are studied as a generalization of \(f(z)=\frac{\mathrm{e}^{z}}{z+1}\) which is nothing but \(E_1(z)\). This name is justified by showing that \(E_n\) has an omitted value 0 and there are no other finite singular value. In fact, it is shown that there is only one singularity over 0 as well as over \(\infty \) and both are direct. Non-existence of Herman rings are proved for \(\lambda E_n \).     

Derivations Vanishing Identities Involving Generalized Derivations and Multilinear Polynomial in Prime Rings

Let R be a prime ring of characteristic different from 2 with Utumi quotient ring U and extended centroid C, \(f(x_1,\ldots ,x_n)\) be a multilinear polynomial over C, which is not central valued on R. Suppose that d is a non-zero derivation of R, F and G are two generalized derivations of R such that \(d\{F(u)u-uG^2(u)\}=0\) for all \(u\in f(R)\). Then one of the following holds:

1. (i)
   there exist \(a, b, p\in U\), \(\lambda \in C\) such that \(F(x)=\lambda x+bx+xa^2\), \(G(x)=ax\), \(d(x)=[p, x]\) for all \(x\in R\) with \([p, b]=0\) and \(f(x_1,\ldots , x_n)^2\) is central valued on R;
    
2. (ii)
   there exist \(a, b, p\in U\) such that \(F(x)=ax\), \(G(x)=xb\), \(d(x)=[p,x]\) for all \(x\in R\) and \(f(x_1,\ldots , x_n)^2\) is central valued on R with \([p, a-b^2]=0\);
    
3. (iii)
   there exist \(a\in U\) such that \(F(x)=xa^2\) and \(G(x)=ax\) for all \(x\in R\);
    
4. (iv)
   there exists \(a\in U\) such that \(F(x)=a^2x\) and \(G(x)=xa\) for all \(x\in R\) with \(a^2\in C\);
    
5. (v)
   there exist \(a, p\in U\), \(\lambda , \alpha , \mu \in C\) such that \(F(x)=\lambda x-a^2x\), \(G(x)=xa\) and \(d(x)=[p,x]\) for all \(x\in R\) with \(a^2=\mu -\alpha p\) and \(\alpha p^2+(\lambda -2\mu ) p\in C\);
    
6. (vi)
   there exist \(a\in U\), \(\lambda \in C\) such that R satisfies \(s_4\) and either \(F(x)=\lambda x+xa^2\), \(G(x)=ax\) or \(F(x)=\lambda x-a^2x\), \(G(x)=xa\) for all \(x\in R\).
    

     

Relative Difference Sets and Hadamard Matrices from Perfect Quaternionic Arrays

Let \(G=\mathbf{C}_{n_1}\times \cdots \times \mathbf{C}_{n_m}\) be an abelian group of order \(n=n_1\dots n_m\), where each \(\mathbf{C}_{n_t}\) is cyclic of order \(n_t\). We present a correspondence between the (4n, 2, 4n, 2n)-relative difference sets in \(G\times Q_8\) relative to the centre \(Z(Q_8)\) and the perfect arrays of size \(n_1\times \dots \times n_m\) over the quaternionic alphabet \(Q_8\cup qQ_8\), where \(q=(1+i+j+k)/2\). In view of this connection, for \(m=2\) we introduce new families of relative difference sets in \(G\times Q_8\), as well as new families of Williamson and Ito Hadamard matrices with G-invariant components.     

Some results on generalized strong external difference families

A generalized strong external difference family (briefly \((v, m; k_1,\dots ,k_m; \lambda _1,\dots ,\lambda _m)\)-GSEDF) was introduced by Paterson and Stinson in 2016. In this paper, we give some nonexistence results for GSEDFs. In particular, we prove that a \((v, 3;k_1,k_2,k_3; \lambda _1,\lambda _2,\lambda _3)\)-GSEDF does not exist when \(k_1+k_2+k_3< v\). We also give a first recursive construction for GSEDFs and prove that if there is a \((v,2;2\lambda ,\frac{v-1}{2};\lambda ,\lambda )\)-GSEDF, then there is a \((vt,2;4\lambda ,\frac{vt-1}{2};2\lambda ,2\lambda )\)-GSEDF with \(v>1\), \(t>1\) and \(v\equiv t\equiv 1\pmod 2\). Then we use it to obtain some new GSEDFs for \(m=2\). In particular, for any prime power q with \(q\equiv 1\pmod 4\), we show that there exists a \((qt, 2;(q-1)2^{n-1},\frac{qt-1}{2};(q-1)2^{n-2},(q-1)2^{n-2})\)-GSEDF, where \(t=p_1p_2\dots p_n\), \(p_i>1\), \(1\le i\le n\), \(p_1, p_2,\dots ,p_n\) are odd integers.     

Diagonalization of the Finite Hilbert Transform on Two Adjacent Intervals

We continue the study of stability of solving the interior problem of tomography. The starting point is the Gelfand–Graev formula, which converts the tomographic data into the finite Hilbert transform (FHT) of an unknown function f along a collection of lines. Pick one such line, call it the x-axis, and assume that the function to be reconstructed depends on a one-dimensional argument by restricting f to the x-axis. Let \(I_1\) be the interval where f is supported, and \(I_2\) be the interval where the Hilbert transform of f can be computed using the Gelfand–Graev formula. The equation to be solved is \(\left. {\mathcal {H}}_1 f=g\right| _{I_2}\), where \({\mathcal {H}}_1\) is the FHT that integrates over \(I_1\) and gives the result on \(I_2\), i.e. \({\mathcal {H}}_1: L^2(I_1)\rightarrow L^2(I_2)\). In the case of complete data, \(I_1\subset I_2\), and the classical FHT inversion formula reconstructs f in a stable fashion. In the case of interior problem (i.e., when the tomographic data are truncated), \(I_1\) is no longer a subset of \(I_2\), and the inversion problems becomes severely unstable. By using a differential operator L that commutes with \({\mathcal {H}}_1\), one can obtain the singular value decomposition of \({\mathcal {H}}_1\). Then the rate of decay of singular values of \({\mathcal {H}}_1\) is the measure of instability of finding f. Depending on the available tomographic data, different relative positions of the intervals \(I_{1,2}\) are possible. The cases when \(I_1\) and \(I_2\) are at a positive distance from each other or when they overlap have been investigated already. It was shown that in both cases the spectrum of the operator \({\mathcal {H}}_1^*{\mathcal {H}}_1\) is discrete, and the asymptotics of its eigenvalues \(\sigma _n\) as \(n\rightarrow \infty \) has been obtained. In this paper we consider the case when the intervals \(I_1=(a_1,0)\) and \(I_2=(0,a_2)\) are adjacent. Here \(a_1 < 0 < a_2\). Using recent developments in the Titchmarsh–Weyl theory, we show that the operator L corresponding to two touching intervals has only continuous spectrum and obtain two isometric transformations \(U_1\), \(U_2\), such that \(U_2{\mathcal {H}}_1 U_1^*\) is the multiplication operator with the function \(\sigma (\lambda )\), \(\lambda \ge (a_1^2+a_2^2)/8\). Here \(\lambda \) is the spectral parameter. Then we show that \(\sigma (\lambda )\rightarrow 0\) as \(\lambda \rightarrow \infty \) exponentially fast. This implies that the problem of finding f is severely ill-posed. We also obtain the leading asymptotic behavior of the kernels involved in the integral operators \(U_1\), \(U_2\) as \(\lambda \rightarrow \infty \). When the intervals are symmetric, i.e. \(-a_1=a_2\), the operators \(U_1\), \(U_2\) are obtained explicitly in terms of hypergeometric functions.     

Group divisible designs with large block sizes

For positive integers n, k with \(3\le k\le n\), let \(X=\mathbb {F}_{2^n}\setminus \{0,1\}\), \({\mathcal {G}}=\{\{x,x+1\}:x\in X\}\), and \({\mathcal {B}}_k=\left\{ \{x_1,x_2,\ldots ,x_k\}\!\subset \!X:\sum \limits _{i=1}^kx_i=1,\ \sum \limits _{i\in I}x_i\!\ne \!1\ \mathrm{for\ any}\ \emptyset \!\ne \!I\!\subsetneqq \!\{1,2,\ldots ,k\}\right\} \). Lee et al. used the inclusion–exclusion principle to show that the triple \((X,{\mathcal {G}},{\mathcal {B}}_k)\) is a \((k,\lambda _k)\)-GDD of type \(2^{2^{n-1}-1}\) for \(k\in \{3,4,5,6,7\}\) where \(\lambda _k=\frac{\prod _{i=3}^{k-1}(2^n-2^i)}{(k-2)!}\) (Lee et al. in Des Codes Cryptogr,  https://doi.org/10.1007/s10623-017-0395-8, 2017). They conjectured that \((X,{\mathcal {G}},{\mathcal {B}}_k)\) is also a \((k,\lambda _k)\)-GDD of type \(2^{2^{n-1}-1}\) for any integer \(k\ge 8\). In this paper, we use a similar construction and counting principles to show that there is a \((k,\lambda _k)\)-GDD of type \((q^2-q)^{(q^{n-1}-1)/(q-1)}\) for any prime power q and any integers k, n with \(3\le k\le n\) where \(\lambda _k=\frac{\prod _{i=3}^{k-1}(q^n-q^i)}{(k-2)!}\). Consequently, their conjecture holds. Such a method is also generalized to yield a \((k,\lambda _k)\)-GDD of type \((q^{\ell +1}-q^{\ell })^{(q^{n-\ell }-1)/(q-1)}\) where \(\lambda _k=\frac{\prod _{i=3}^{k-1}(q^n-q^{\ell +i-1})}{(k-2)!}\) and \(k+\ell \le n+1\).     

Multiple positive solutions of the stationary Keller–Segel system

We consider the stationary Keller–Segel equation

$$\begin{aligned} {\left\{ \begin{array}{ll} -\Delta v+v=\lambda e^v, \quad v>0 \quad &{} \text {in }\Omega ,\\ \partial _\nu v=0 &{}\text {on } \partial \Omega , \end{array}\right. } \end{aligned}$$

where \(\Omega \) is a ball. In the regime \(\lambda \rightarrow 0\), we study the radial bifurcations and we construct radial solutions by a gluing variational method. For any given \(n\in \mathbb {N}_0\), we build a solution having multiple layers at \(r_1,\ldots ,r_n\) by which we mean that the solutions concentrate on the spheres of radii \(r_i\) as \(\lambda \rightarrow 0\) (for all \(i=1,\ldots ,n\)). A remarkable fact is that, in opposition to previous known results, the layers of the solutions do not accumulate to the boundary of \(\Omega \) as \(\lambda \rightarrow 0\). Instead they satisfy an optimal partition problem in the limit.

     

Existence of positive solution for nonlocal singular fourth order Kirchhoff equation with Hardy potential

This paper is concerned with the existence of positive solution to a class of singular fourth order elliptic equation of Kirchhoff type

$$\begin{aligned} \triangle ^2 u-\lambda M(\Vert \nabla u\Vert ^2)\triangle u-\frac{\mu }{\vert x\vert ^4}u=\frac{h(x)}{u^\gamma }+k(x)u^\alpha , \end{aligned}$$

under Navier boundary conditions, \(u=\triangle u=0\). Here \(\varOmega \subset {\mathbf {R}}^N\), \(N\ge 1\) is a bounded \(C^4\)-domain, \(0\in \varOmega \), h(x) and k(x) are positive continuous functions, \(\gamma \in (0,1)\), \(\alpha \in (0,1)\) and \(M:{\mathbf {R}}^+\rightarrow {\mathbf {R}}^+\) is a continuous function. By using Galerkin method and sharp angle lemma, we will show that this problem has a positive solution for \(\lambda > \frac{\mu }{\mu ^*m_0}\) and \(0<\mu <\mu ^*\). Here \(\mu ^*=\Big (\frac{N(N-4)}{4}\Big )^2\) is the best constant in the Hardy inequality. Besides, if \(\mu =0\), \(\lambda >0\) and h, k are Lipschitz functions, we show that this problem has a positive smooth solution. If \(h,k\in C^{2,\,\theta _0}(\overline{\varOmega })\) for some \(\theta _0\in (0,1)\), then this problem has a positive classical solution.

     

The Hardy–Schrödinger operator with interior singularity: the remaining cases

We consider the remaining unsettled cases in the problem of existence of energy minimizing solutions for the Dirichlet value problem \(L_\gamma u-\lambda u=\frac{u^{2^*(s)-1}}{|x|^s}\) on a smooth bounded domain \(\Omega \) in \({\mathbb {R}}^n\) (\(n\ge 3\)) having the singularity 0 in its interior. Here \(\gamma <\frac{(n-2)^2}{4}\), \(0\le s <2\), \(2^*(s):=\frac{2(n-s)}{n-2}\) and \(0\le \lambda <\lambda _1(L_\gamma )\), the latter being the first eigenvalue of the Hardy–Schrödinger operator \(L_\gamma :=-\Delta -\frac{\gamma }{|x|^2}\). There is a threshold \(\lambda ^*(\gamma , \Omega ) \ge 0\) beyond which the minimal energy is achieved, but below which, it is not. It is well known that \(\lambda ^*(\Omega )=0\) in higher dimensions, for example if \(0\le \gamma \le \frac{(n-2)^2}{4}-1\). Our main objective in this paper is to show that this threshold is strictly positive in “lower dimensions” such as when \( \frac{(n-2)^2}{4}-1<\gamma <\frac{(n-2)^2}{4}\), to identify the critical dimensions (i.e., when the situation changes), and to characterize it in terms of \(\Omega \) and \(\gamma \). If either \(s>0\) or if \(\gamma > 0\), i.e., in the truly singular case, we show that in low dimensions, a solution is guaranteed by the positivity of the “Hardy-singular internal mass” of \(\Omega \), a notion that we introduce herein. On the other hand, and just like the case when \(\gamma =s=0\) studied by Brezis and Nirenberg (Commun Pure Appl Math 36:437–477, 1983) and completed by Druet (Ann Inst H Poincaré Anal Non Linéaire 19(2):125–142, 2002), \(n=3\) is the critical dimension, and the classical positive mass theorem is sufficient for the merely singular case, that is when \(s=0\), \(\gamma \le 0\).     

Reducibility for a Class of Analytic Multipliers on the Dirichlet Space

It is proved that if \(\phi \) is a finite Blaschke product with four zeros, then \(M_\phi \) is reducible on the Dirichlet space with norm \(\Vert \ \Vert \) if and only if \(\phi =\phi _1\circ \phi _2\), where \(\phi _1, \phi _2\) are Blaschke products and \(\phi _2\) is equivalent to \(z^2\). Also, the same reducibility of \(M_\phi \) with finite Blaschke product \(\phi \) on the Dirichlet space under the equivalent norms \(\Vert \ \Vert _1\) and \(\Vert \ \Vert _0\) is given.     

Diophantine approximation with four squares and one <Emphasis Type="Italic">k</Emphasis>th power of primes

Let k be an integer with \(k\ge 3\) and \(\eta \) be any real number. Suppose that \(\lambda _1, \lambda _2, \lambda _3, \lambda _4, \mu \) are non-zero real numbers, not all of the same sign and \(\lambda _1/\lambda _2\) is irrational. It is proved that the inequality \(|\lambda _1p_1^2+\lambda _2p_2^2+\lambda _3p_3^2+\lambda _4p_4^2+\mu p_5^k+\eta |<(\max \ p_j)^{-\sigma }\) has infinitely many solutions in prime variables \(p_1, p_2, \ldots , p_5\), where \(0<\sigma <\frac{1}{16}\) for \(k=3,\ 0<\sigma <\frac{5}{3k2^k}\) for \(4\le k\le 5\) and \(0<\sigma <\frac{40}{21k2^k}\) for \(k\ge 6\). This gives an improvement of an earlier result.     

$$I_0$$ and combinatorics at $$\lambda ^+$$λ+

We investigate the compatibility of \(I_0\) with various combinatorial principles at \(\lambda ^+\), which include the existence of \(\lambda ^+\)-Aronszajn trees, square principles at \(\lambda \), the existence of good scales at \(\lambda \), stationary reflections for subsets of \(\lambda ^{+}\), diamond principles at \(\lambda \) and the singular cardinal hypothesis at \(\lambda \). We also discuss whether these principles can hold in \(L(V_{\lambda +1})\).     

Simultaneous power factorization in modules over Banach algebras

Let A be a Banach algebra with a bounded left approximate identity \(\{e_\lambda \}_{\lambda \in \Lambda }\), let \(\pi \) be a continuous representation of A on a Banach space X, and let S be a non-empty subset of X such that \(\lim _{\lambda }\pi (e_\lambda )s=s\) uniformly on S. If S is bounded, or if \(\{e_\lambda \}_{\lambda \in \Lambda }\) is commutative, then we show that there exist \(a\in A\) and maps \(x_n: S\rightarrow X\) for \(n\ge 1\) such that \(s=\pi (a^n)x_n(s)\) for all \(n\ge 1\) and \(s\in S\). The properties of \(a\in A\) and the maps \(x_n\), as produced by the constructive proof, are studied in some detail. The results generalize previous simultaneous factorization theorems as well as Allan and Sinclair’s power factorization theorem. In an ordered context, we also consider the existence of a positive factorization for a subset of the positive cone of an ordered Banach space that is a positive module over an ordered Banach algebra with a positive bounded left approximate identity. Such factorizations are not always possible. In certain cases, including those for positive modules over ordered Banach algebras of bounded functions, such positive factorizations exist, but the general picture is still unclear. Furthermore, simultaneous pointwise power factorizations for sets of bounded maps with values in a Banach module (such as sets of bounded convergent nets) are obtained. A worked example for the left regular representation of \(\mathrm {C}_0({\mathbb R})\) and unbounded S is included.