Showing 1–50 of 59 results for author: Pandey, A

Exploring Algorithmic Solutions for the Independent Roman Domination Problem in Graphs

Authors: Kaustav Paul, Ankit Sharma, Arti Pandey

Abstract: Given a graph $G=(V,E)$, a function $f:V\to \{0,1,2\}$ is said to be a \emph{Roman Dominating function} if for every $v\in V$ with $f(v)=0$, there exists a vertex $u\in N(v)$ such that $f(u)=2$. A Roman Dominating function $f$ is said to be an \emph{Independent Roman Dominating function} (or IRDF), if $V_1\cup V_2$ forms an independent set, where $V_i=\{v\in V~\vert~f(v)=i\}$, for… ▽ More Given a graph $G=(V,E)$, a function $f:V\to \{0,1,2\}$ is said to be a \emph{Roman Dominating function} if for every $v\in V$ with $f(v)=0$, there exists a vertex $u\in N(v)$ such that $f(u)=2$. A Roman Dominating function $f$ is said to be an \emph{Independent Roman Dominating function} (or IRDF), if $V_1\cup V_2$ forms an independent set, where $V_i=\{v\in V~\vert~f(v)=i\}$, for $i\in \{0,1,2\}$. The total weight of $f$ is equal to $\sum_{v\in V} f(v)$, and is denoted as $w(f)$. The \emph{Independent Roman Domination Number} of $G$, denoted by $i_R(G)$, is defined as min$\{w(f)~\vert~f$ is an IRDF of $G\}$. For a given graph $G$, the problem of computing $i_R(G)$ is defined as the \emph{Minimum Independent Roman Domination problem}. The problem is already known to be NP-hard for bipartite graphs. In this paper, we further study the algorithmic complexity of the problem. In this paper, we propose a polynomial-time algorithm to solve the Minimum Independent Roman Domination problem for distance-hereditary graphs, split graphs, and $P_4$-sparse graphs. △ Less

Submitted 12 July, 2024; v1 submitted 4 July, 2024; originally announced July 2024.

Musielak-Orlicz-Sobolev embeddings: Necessary and Sufficient Conditions

Authors: Ankur Pandey, Nijjwal Karak

Abstract: In this paper we study the necessary and sufficient conditions on domain for Musielak-Orlicz-Sobolev embedding of the space $W^{1,Φ(\cdot,\cdot)}(Ω)$ where $Φ(x,t):=t^{p(x)}{(\log(e+t))}^{q(x)}.$ In this paper we study the necessary and sufficient conditions on domain for Musielak-Orlicz-Sobolev embedding of the space $W^{1,Φ(\cdot,\cdot)}(Ω)$ where $Φ(x,t):=t^{p(x)}{(\log(e+t))}^{q(x)}.$ △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 15 pages

On Generalized Transmuted Lifetime Distribution

Authors: Alok Kumar Pandey, Alam Ali, Ashok Kumar Pathak

Abstract: This article presents a new class of generalized transmuted lifetime distributions which includes a large number of lifetime distributions as sub-family. Several important mathematical quantities such as density function, distribution function, quantile function, moments, moment generating function, stress-strength reliability function, order statistics, Rényi and q-entropy, residual and reversed… ▽ More This article presents a new class of generalized transmuted lifetime distributions which includes a large number of lifetime distributions as sub-family. Several important mathematical quantities such as density function, distribution function, quantile function, moments, moment generating function, stress-strength reliability function, order statistics, Rényi and q-entropy, residual and reversed residual life function, and cumulative information generating function are obtained. The methods of maximum likelihood, ordinary least square, weighted least square, Cramér-von Mises, Anderson Darling, and Right-tail Anderson Darling are considered to estimate the model parameters in a general way. Further, a well-organized Monte Carlo simulation experiments have been performed to observe the behavior of the estimators. Finally, two real data have also been analyzed to demonstrate the effectiveness of the proposed distribution in real-life modeling. △ Less

Submitted 19 May, 2024; originally announced May 2024.

Comments: 26 pages, 8 figures

MSC Class: 60E05; 62F10; 62E15; 65C05; 33B20

A note on the equivalence of Gromov boundary and metric boundary

Authors: Vasudevarao Allu, Abhishek Pandey

Abstract: In this paper, we introduce the concept of quasihyperbolically visible spaces. As a tool, we study the connection between the Gromov boundary and the metric boundary. In this paper, we introduce the concept of quasihyperbolically visible spaces. As a tool, we study the connection between the Gromov boundary and the metric boundary. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: 12 pages. Comments are welcome. arXiv admin note: text overlap with arXiv:2306.03815

MSC Class: 30L10; 30L99; 30C65; 51F30; 53C23

Record-based transmuted unit omega distribution: different methods of estimation and applications

Authors: Ashok Kumar Pathak, Mohd. Arshad, Alok Kumar Pandey, Alam Ali

Abstract: Dombi et al. (2019) introduced a three parameter omega distribution and showed that its asymptotic distribution is the Weibull model. We propose a new record-based transmuted generalization of the unit omega distribution by considering Balakrishnan and He (2021) approach. We call it the RTUOMG distribution. We derive expressions for some statistical quantities, like, probability density function,… ▽ More Dombi et al. (2019) introduced a three parameter omega distribution and showed that its asymptotic distribution is the Weibull model. We propose a new record-based transmuted generalization of the unit omega distribution by considering Balakrishnan and He (2021) approach. We call it the RTUOMG distribution. We derive expressions for some statistical quantities, like, probability density function, distribution, hazard function, quantile function, moments, incomplete moments, inverted moments, moment generating function, Lorenz curve, and Bonferroni curve of the proposed distribution. The numerical values of various measures of central tendency and coefficient of skewness and kurtosis are also presented. Concepts of stochastic ordering and some results related to ordered statistics of the RTUOMG distribution are discussed. The parameters of the RTUOMG distribution are estimated using five distinct estimators. Additionally, the Monte Carlo simulations are performed to assess the performance of these estimators. Finally, two real data sets are analyzed to demonstrate the utility of the RTUOMG distribution. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 26 pages, 14 figures

MSC Class: 60E05; 62F10; 62E15; 65C05; 33C05

Exact recovery in Gaussian weighted stochastic block model and planted dense subgraphs: Statistical and algorithmic thresholds

Authors: Aaradhya Pandey, Sanjeev Kulkarni

Abstract: In this paper, we study the exact recovery problem in the Gaussian weighted version of the Stochastic block model with two symmetric communities. We provide the information-theoretic threshold in terms of the signal-to-noise ratio (SNR) of the model and prove that when SNR $<1$, no statistical estimator can exactly recover the community structure with probability bounded away from zero. On the oth… ▽ More In this paper, we study the exact recovery problem in the Gaussian weighted version of the Stochastic block model with two symmetric communities. We provide the information-theoretic threshold in terms of the signal-to-noise ratio (SNR) of the model and prove that when SNR $<1$, no statistical estimator can exactly recover the community structure with probability bounded away from zero. On the other hand, we show that when SNR $>1$, the Maximum likelihood estimator itself succeeds in exactly recovering the community structure with probability approaching one. Then, we provide two algorithms for achieving exact recovery. The Semi-definite relaxation as well as the spectral relaxation of the Maximum likelihood estimator can recover the community structure down to the threshold value of 1, establishing the absence of an information-computation gap for this model. Next, we compare the problem of community detection with the problem of recovering a planted densely weighted community within a graph and prove that the exact recovery of two symmetric communities is a strictly easier problem than recovering a planted dense subgraph of size half the total number of nodes, by establishing that when the same SNR$< 3/2$, no statistical estimator can exactly recover the planted community with probability bounded away from zero. More precisely, when $1 <$ SNR $< 3/2$ exact recovery of community detection is possible, both statistically and algorithmically, but it is impossible to exactly recover the planted community, even statistically, in the Gaussian weighted model. Finally, we show that when SNR $>2$, the Maximum likelihood estimator itself succeeds in exactly recovering the planted community with probability approaching one. We also prove that the Semi-definite relaxation of the Maximum likelihood estimator can recover the planted community structure down to the threshold value of 2. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Community detection in the hypergraph stochastic block model and reconstruction on hypertrees

Authors: Yuzhou Gu, Aaradhya Pandey

Abstract: We study the weak recovery problem on the $r$-uniform hypergraph stochastic block model ($r$-HSBM) with two balanced communities. In this model, $n$ vertices are randomly divided into two communities, and size-$r$ hyperedges are added randomly depending on whether all vertices in the hyperedge are in the same community. The goal of weak recovery is to recover a non-trivial fraction of the communit… ▽ More We study the weak recovery problem on the $r$-uniform hypergraph stochastic block model ($r$-HSBM) with two balanced communities. In this model, $n$ vertices are randomly divided into two communities, and size-$r$ hyperedges are added randomly depending on whether all vertices in the hyperedge are in the same community. The goal of weak recovery is to recover a non-trivial fraction of the communities given the hypergraph. Pal and Zhu (2021); Stephan and Zhu (2022) established that weak recovery is always possible above a natural threshold called the Kesten-Stigum (KS) threshold. For assortative models (i.e., monochromatic hyperedges are preferred), Gu and Polyanskiy (2023) proved that the KS threshold is tight if $r\le 4$ or the expected degree $d$ is small. For other cases, the tightness of the KS threshold remained open. In this paper we determine the tightness of the KS threshold for a wide range of parameters. We prove that for $r\le 6$ and $d$ large enough, the KS threshold is tight. This shows that there is no information-computation gap in this regime and partially confirms a conjecture of Angelini et al. (2015). On the other hand, we show that for $r\ge 5$, there exist parameters for which the KS threshold is not tight. In particular, for $r\ge 7$, the KS threshold is not tight if the model is disassortative (i.e., polychromatic hyperedges are preferred) or $d$ is large enough. This provides more evidence supporting the existence of an information-computation gap in these cases. Furthermore, we establish asymptotic bounds on the weak recovery threshold for fixed $r$ and large $d$. We also obtain a number of results regarding the broadcasting on hypertrees (BOHT) model, including the asymptotics of the reconstruction threshold for $r\ge 7$ and impossibility of robust reconstruction at criticality. △ Less

Submitted 11 June, 2024; v1 submitted 9 February, 2024; originally announced February 2024.

Making Distribution State Estimation Practical: Challenges and Opportunities

Authors: Frederik Geth, Marta Vanin, Werner Van Westering, Terese Milford, Amritanshu Pandey

Abstract: In increasingly digitalized and metered distribution networks, state estimation is generally recognized as a key enabler of advanced network management functionalities. However, despite decades of research, the real-life adoption of state estimation in distribution systems remains sporadic. This systematization of knowledge paper discusses the cause for this while comparing industrial and academic… ▽ More In increasingly digitalized and metered distribution networks, state estimation is generally recognized as a key enabler of advanced network management functionalities. However, despite decades of research, the real-life adoption of state estimation in distribution systems remains sporadic. This systematization of knowledge paper discusses the cause for this while comparing industrial and academic experiences and reviewing well- and less-established research directions. We argue that to make distribution system state estimation more practical and applicable in the field, new perspectives are needed. In particular, research should move away from conventional approaches and embrace generalized problem specifications and more comprehensive workflows. These, in turn, require algorithm advancements and more general mathematical formulations. We discuss lines of work to enable the delivery of tangible research. △ Less

Submitted 12 November, 2023; originally announced November 2023.

Comments: 10 pages

Predictive Optimization of Hybrid Energy Systems with Temperature Dependency

Authors: Tanmay Mishra, Amritanshu Pandey, Mads R. Almassalkhi

Abstract: Hybrid Energy Systems (HES), amalgamating renewable sources, energy storage, and conventional generation, have emerged as a responsive resource for providing valuable grid services. Subsequently, modeling and analysis of HES has become critical, and the quality of grid services hedges on it. Currently, most HES models are temperature-agnostic. However, the temperature-dependent factors can signifi… ▽ More Hybrid Energy Systems (HES), amalgamating renewable sources, energy storage, and conventional generation, have emerged as a responsive resource for providing valuable grid services. Subsequently, modeling and analysis of HES has become critical, and the quality of grid services hedges on it. Currently, most HES models are temperature-agnostic. However, the temperature-dependent factors can significantly impact HES performance, necessitating advanced modeling and optimization techniques. With the inclusion of temperature-dependent models, the challenges and complexity of solving optimization problem increases. In this paper, the electro-thermal modeling of HES is discussed. Based on this model, a nonlinear predictive optimization framework is formulated. A simplified model is developed to address the challenges associated with solving nonlinear problems. Further, projection and homotopy approaches are proposed. In the homotopy method, the NLP is solved by incrementally changing the C-rating of the battery. Simulation-based analysis of the algorithms highlights the effects of different battery ratings, ambient temperatures, and energy price variations. Finally, comparative assessments with a temperature-agnostic approach illustrates the effectiveness of electro-thermal methods in optimizing HES. △ Less

Submitted 15 March, 2024; v1 submitted 1 November, 2023; originally announced November 2023.

Computation of Grundy dominating sequences in (co-)bipartite graphs

Authors: Boštjan Brešar, Arti Pandey, Gopika Sharma

Abstract: A sequence $S$ of vertices of a graph $G$ is called a dominating sequence of $G$ if $(i)$ each vertex $v$ of $S$ dominates a vertex of $G$ that was not dominated by any of the vertices preceding vertex $v$ in $S$, and $(ii)$ every vertex of $G$ is dominated by at least one vertex of $S$. The Grundy Domination problem is to find a longest dominating sequence for a given graph $G$. It has been known… ▽ More A sequence $S$ of vertices of a graph $G$ is called a dominating sequence of $G$ if $(i)$ each vertex $v$ of $S$ dominates a vertex of $G$ that was not dominated by any of the vertices preceding vertex $v$ in $S$, and $(ii)$ every vertex of $G$ is dominated by at least one vertex of $S$. The Grundy Domination problem is to find a longest dominating sequence for a given graph $G$. It has been known that the decision version of the Grundy Domination problem is NP-complete even when restricted to chordal graphs. In this paper, we prove that the decision version of the Grundy Domination problem is NP-complete for bipartite graphs and co-bipartite graphs. On the positive side, we present a linear-time algorithm that solves the Grundy Domination problem for chain graphs, which form a subclass of bipartite graphs. △ Less

Submitted 16 October, 2023; originally announced October 2023.

Comments: 18 pages, 3 figures

MSC Class: 05C69; 05C65; 05C85

Conformal Warped Product Submersion

Authors: Harmandeep Kaur, Abhishek Pandey, Gauree Shanker

Abstract: In this paper, the concept of Riemannian warped product submersion is generalized to the conformal case. We introduce the notion of conformal warped product submersion. It is a submersion between warped product manifolds that preserves angles between the horizontal vectors. The fundamental tensors of submersion are derived for conformal warped product submersion. In this paper, the concept of Riemannian warped product submersion is generalized to the conformal case. We introduce the notion of conformal warped product submersion. It is a submersion between warped product manifolds that preserves angles between the horizontal vectors. The fundamental tensors of submersion are derived for conformal warped product submersion. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 9 pages

MSC Class: 53C15; 53C18; 53C20; 53B25

A characterization of $b$-generalized skew derivations on lie ideal of a prime ring

Authors: Ashutosh Pandey, Mani Shankar Pandey

Abstract: Let $R$ be a prime ring of characteristic different from $2$, $U$ be its Utumi quotient ring, $C$ be its extended centroid and $L$ be a non-central Lie ideal of $R$. Suppose $F$ and $G$ are two non-zero $b$-generalized skew derivations of $R$ associated with the same automorphism $α$ such that $$puF(u) + F(u)uq = G(u^2),\ \text{with} \ p + q \notin C,~~ \text{for all $u \in L$. }$… ▽ More Let $R$ be a prime ring of characteristic different from $2$, $U$ be its Utumi quotient ring, $C$ be its extended centroid and $L$ be a non-central Lie ideal of $R$. Suppose $F$ and $G$ are two non-zero $b$-generalized skew derivations of $R$ associated with the same automorphism $α$ such that $$puF(u) + F(u)uq = G(u^2),\ \text{with} \ p + q \notin C,~~ \text{for all $u \in L$. }$$ This paper aims to study the complete structure of the $\mathrm{b}$-generalized skew derivations $F$ and $G$. △ Less

Submitted 23 June, 2023; originally announced June 2023.

MSC Class: 16N60; 16W25

Visible quasihyperbolic geodesics

Authors: Vasudevarao Allu, Abhishek Pandey

Abstract: In this article, we initiate the study of visibility property in the context of quasihyperbolic metric. Let $Ω$ be a bounded domain in $\mathbb{R}^n$ and $\partial_{Euc}Ω$ be its Euclidean boundary. We say that the pair $\{p,q\}\subset\partial_{Euc}Ω$ has visible quasihyperbolic geodesic if there exist open neighborhoods $U_p$ of $p$ and $U_q$ of $q$ and a compact set $K_{\{p,q\}}\subset Ω$ such t… ▽ More In this article, we initiate the study of visibility property in the context of quasihyperbolic metric. Let $Ω$ be a bounded domain in $\mathbb{R}^n$ and $\partial_{Euc}Ω$ be its Euclidean boundary. We say that the pair $\{p,q\}\subset\partial_{Euc}Ω$ has visible quasihyperbolic geodesic if there exist open neighborhoods $U_p$ of $p$ and $U_q$ of $q$ and a compact set $K_{\{p,q\}}\subset Ω$ such that for any quasihyperbolic geodesic which joins points in $U_p$ and $U_q$ intersects with $K$. If a domain has visible quasihyperbolic geodesic for every pair of points, we say that it is a visibility domain. The major part of this paper is devoted to provide a rich collection of visibility domains. In this direction we provide a general visibility criteria for a domain to be a visibility domain, and using this we obtain the visibility of Uniform domains, John domains and domains satisfying quasihyperbolic boundary conditions. We also study the visibility of bounded hyperbolic domains with respect to hyperbolic and quasihyperbolic metric. Further, we explore the relation between Gromov hyperbolicity and visibility. As an application of visibility we prove a sufficient condition for the continuous extension of quasihyperbolic isometries and quasi-isometries. Further, motivated by the work of Bharali and Zimmer [{\it Trans. Amer. Math. Soc. (2022)}] and using the end compactification we define the visibility of unbounded domains and study its relationship with Gromov hyperbolicity. Finally, we provide some useful examples of visibility domains which are neither John domain nor domains satisfying quasihyperbolic boundary conditions. △ Less

Submitted 18 August, 2023; v1 submitted 6 June, 2023; originally announced June 2023.

Comments: 51 pages, 2 figures. A revised version with significant additions. For instance, visibility of unbounded domains

MSC Class: Primary 30F45; 30L10; 30L99; 30C65. Secondary 51F99; 53C22

Transverse spectral instabilities in Rotation-Modified Kadomtsev-Petviashvili equation and related models

Authors: Bhavna, Ashish Kumar Pandey, Anastassiya Semenova

Abstract: The rotation modified Kadomtsev Petviashvili equation which is also known as the Kadomtsev Petviashvili Ostrovsky equation, describes the gradual wave field diffusion in the transverse direction to the direction of the propagation of the wave in a rotating frame of reference. This equation is a generalization of the Ostrovsky equation additionally having weak transverse effects. We investigate tra… ▽ More The rotation modified Kadomtsev Petviashvili equation which is also known as the Kadomtsev Petviashvili Ostrovsky equation, describes the gradual wave field diffusion in the transverse direction to the direction of the propagation of the wave in a rotating frame of reference. This equation is a generalization of the Ostrovsky equation additionally having weak transverse effects. We investigate transverse instability and stability of small periodic traveling waves of the Ostrovsky equation with respect to either periodic or square-integrable perturbations in the direction of wave propagation and periodic perturbations in the transverse direction of motion in the rotation-modified Kadomtsev Petviashvili equation. We also study transverse stability or instability in generalized rotation modified KP equation by taking dispersion term as general and quadratic and cubic nonlinearity. As a consequence, we obtain transverse stability or instability in two dimensional generalization of RMBO eqution, Ostrovsky Gardner equation, Ostrovsky fKdV equation, Ostrovsky mKdV equation, Ostrovsky ILW equation, Ostrovsky Whitham etc. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Modulational Instability in the Ostrovsky Equation and Related Models

Authors: Bhavna, Mathew A. Johnson, Ashish Kumar Pandey

Abstract: We study the modulational instability of small-amplitude periodic traveling wave solutions in a generalized Ostrovsky equation. Specifically, we investigate the invertibility of the associated linearized operator in the vicinity of the origin and derive a modulational instability index that depends on the dispersion and nonlinearity. We then show that the small-amplitude periodic traveling waves o… ▽ More We study the modulational instability of small-amplitude periodic traveling wave solutions in a generalized Ostrovsky equation. Specifically, we investigate the invertibility of the associated linearized operator in the vicinity of the origin and derive a modulational instability index that depends on the dispersion and nonlinearity. We then show that the small-amplitude periodic traveling waves of the Ostrovsky equation exhibit modulational instability if the wavenumber is greater than a critical value which agrees with previously found numerical results in \cite{Whitfield2014Rotation-inducedWaves} both qualitatively and quantitatively. We also study the effects of surface tension on modulational instability using the index. △ Less

Submitted 14 May, 2023; originally announced May 2023.

Bayesian Analysis of Generalized Hierarchical Indian Buffet Processes for Within and Across Group Sharing of Latent Features

Authors: Lancelot Fitzgerald James, Juho Lee, Abhinav Pandey

Abstract: Bayesian nonparametric hierarchical priors provide flexible models for sharing of information within and across groups. We focus on latent feature allocation models, where the data structures correspond to multisets or unbounded sparse matrices. The fundamental development in this regard is the Hierarchical Indian Buffet process (HIBP), devised by Thibaux and Jordan (2007). However, little is know… ▽ More Bayesian nonparametric hierarchical priors provide flexible models for sharing of information within and across groups. We focus on latent feature allocation models, where the data structures correspond to multisets or unbounded sparse matrices. The fundamental development in this regard is the Hierarchical Indian Buffet process (HIBP), devised by Thibaux and Jordan (2007). However, little is known in terms of explicit tractable descriptions of the joint, marginal, posterior and predictive distributions of the HIBP. We provide explicit novel descriptions of these quantities, in the Bernoulli HIBP and general spike and slab HIBP settings, which allows for exact sampling and simpler practical implementation. We then extend these results to the more complex setting of hierarchies of general HIBP (HHIBP). The generality of our framework allows one to recognize important structure that may otherwise be masked in the Bernoulli setting, and involves characterizations via dynamic mixed Poisson random count matrices. Our analysis shows that the standard choice of hierarchical Beta processes for modeling across group sharing is not ideal in the classic Bernoulli HIBP setting proposed by Thibaux and Jordan (2007), or other spike and slab HIBP settings, and we thus indicate tractable alternative priors. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Comments: This is an extensive re-write and extension of arXiv:2103.11407 where variations of the results for the HIBP (but not HHIBP) were established

MSC Class: 60C05; 60G09; 60G57; 60E99

Complexity of total dominator coloring in graphs

Authors: Michael A. Henning, Kusum, Arti Pandey, Kaustav Paul

Abstract: Let $G=(V,E)$ be a graph with no isolated vertices. A vertex $v$ totally dominate a vertex $w$ ($w \ne v$), if $v$ is adjacent to $w$. A set $D \subseteq V$ called a total dominating set of $G$ if every vertex $v\in V$ is totally dominated by some vertex in $D$. The minimum cardinality of a total dominating set is the total domination number of $G$ and is denoted by $γ_t(G)$. A total dominator col… ▽ More Let $G=(V,E)$ be a graph with no isolated vertices. A vertex $v$ totally dominate a vertex $w$ ($w \ne v$), if $v$ is adjacent to $w$. A set $D \subseteq V$ called a total dominating set of $G$ if every vertex $v\in V$ is totally dominated by some vertex in $D$. The minimum cardinality of a total dominating set is the total domination number of $G$ and is denoted by $γ_t(G)$. A total dominator coloring of graph $G$ is a proper coloring of vertices of $G$, so that each vertex totally dominates some color class. The total dominator chromatic number $χ_{td}(G)$ of $G$ is the least number of colors required for a total dominator coloring of $G$. The Total Dominator Coloring problem is to find a total dominator coloring of $G$ using the minimum number of colors. It is known that the decision version of this problem is NP-complete for general graphs. We show that it remains NP-complete even when restricted to bipartite, planar and split graphs. We further study the Total Dominator Coloring problem for various graph classes, including trees, cographs and chain graphs. First, we characterize the trees having $χ_{td}(T)=γ_t(T)+1$, which completes the characterization of trees achieving all possible values of $χ_{td}(T)$. Also, we show that for a cograph $G$, $χ_{td}(G)$ can be computed in linear-time. Moreover, we show that $2 \le χ_{td}(G) \le 4$ for a chain graph $G$ and give characterization of chain graphs for every possible value of $χ_{td}(G)$ in linear-time. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Comments: V1, 18 pages, 1 figure

Cosecure Domination: Hardness Results and Algorithm

Authors: Kusum, Arti Pandey

Abstract: For a simple graph $G=(V,E)$ without any isolated vertex, a cosecure dominating set $D$ of $G$ satisfies the following two properties (i) $S$ is a dominating set of $G$, (ii) for every vertex $v \in S$ there exists a vertex $u \in V \setminus S$ such that $uv \in E$ and $(S \setminus \{v\}) \cup \{u\}$ is a dominating set of $G$. The minimum cardinality of a cosecure dominating set of $G$ is calle… ▽ More For a simple graph $G=(V,E)$ without any isolated vertex, a cosecure dominating set $D$ of $G$ satisfies the following two properties (i) $S$ is a dominating set of $G$, (ii) for every vertex $v \in S$ there exists a vertex $u \in V \setminus S$ such that $uv \in E$ and $(S \setminus \{v\}) \cup \{u\}$ is a dominating set of $G$. The minimum cardinality of a cosecure dominating set of $G$ is called cosecure domination number of $G$ and is denoted by $γ_{cs}(G)$. The Minimum Cosecure Domination problem is to find a cosecure dominating set of a graph $G$ of cardinality $γ_{cs}(G)$. The decision version of the problem is known to be NP-complete for bipartite, planar, and split graphs. Also, it is known that the Minimum Cosecure Domination problem is efficiently solvable for proper interval graphs and cographs. In this paper, we work on various important graph classes in an effort to reduce the complexity gap of the Minimum Cosecure Domination problem. We show that the decision version of the problem remains NP-complete for circle graphs, doubly chordal graphs, chordal bipartite graphs, star-convex bipartite graphs and comb-convex bipartite graphs. On the positive side, we give an efficient algorithm to compute the cosecure domination number of chain graphs, which is an important subclass of bipartite graphs. In addition, we show that the problem is linear-time solvable for bounded tree-width graphs. Further, we prove that the computational complexity of this problem varies from the domination problem. △ Less

Submitted 25 February, 2023; originally announced February 2023.

Comments: V1, 19 pages, 2 figures

Generalized skew derivations on ideal with engel conditions

Authors: Ashutosh Pandey, Balchand Prajapati

Abstract: Let R be a prime ring of characteristic different from 2, U be the Utumi quotient ring of R and C be the extended centroid of R. Let F be a generalized skew derivation on R, I be a non-zero ideal of R. Then we give the complete structure of F satisfying certain conditions. Let R be a prime ring of characteristic different from 2, U be the Utumi quotient ring of R and C be the extended centroid of R. Let F be a generalized skew derivation on R, I be a non-zero ideal of R. Then we give the complete structure of F satisfying certain conditions. △ Less

Submitted 31 January, 2023; originally announced January 2023.

Neargeodesics in Gromov hyperbolic John domains in Banach spaces

Authors: Vasudevarao Allu, Abhishek Pandey

Abstract: In this paper, we prove that neargeodesics in Gromov hyperbolic John domains in Banach space are cone arcs. This result gives an improvement of a result of Li [Theorem 1, Int. J. Math. 25 (2014)]. In this paper, we prove that neargeodesics in Gromov hyperbolic John domains in Banach space are cone arcs. This result gives an improvement of a result of Li [Theorem 1, Int. J. Math. 25 (2014)]. △ Less

Submitted 14 June, 2024; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: 14 pages. Revised version accepted for publication in Annales Fennici Mathematici

MSC Class: 30C65; 30L10; 30F45; 30C20

Actionable Three-Phase Infeasibility Optimization with Varying Slack Sources

Authors: Elizabeth Foster, Timothy McNamara, Amritanshu Pandey, Larry Pileggi

Abstract: Modern distribution grids that include numerous distributed energy resources (DERs) and battery electric vehicles (BEVs) will require simulation and optimization methods that can capture behavior under infeasible operating scenarios to assess reliability. A three-phase infeasibility analysis (TPIA) localizes and identifies power deficient areas in distribution feeders via a non-convex optimization… ▽ More Modern distribution grids that include numerous distributed energy resources (DERs) and battery electric vehicles (BEVs) will require simulation and optimization methods that can capture behavior under infeasible operating scenarios to assess reliability. A three-phase infeasibility analysis (TPIA) localizes and identifies power deficient areas in distribution feeders via a non-convex optimization that injects and subsequently minimizes slack sources, subject to AC network constraints. In this paper, we extend the TPIA framework by introducing operational bounds to ensure realistic, actionable solutions. We incorporate current, reactive power, and susceptance slack sources to model real-world assets, and discuss their potential use cases. We show that the voltage-bounded TPIA formulations provide actionable solutions for realistic networks of up to 5360 nodes where power flow simulations either fail or return low-voltage solutions. We demonstrate reactive power compensation using the slack susceptance formulation on an infeasible test case. △ Less

Submitted 16 February, 2023; v1 submitted 6 December, 2022; originally announced December 2022.

b-generalized skew derivations acting as 2-Jordan multiplier on multilinear polynomials in prime rings

Authors: Mani Shankar Pandey, Ashutosh Pandey

Abstract: Let R be a prime ring of characteristic not equal to 2, U be its Utumi quotient ring and C be the extended centroid of R. Let φbe a multilinear polynomial over C, which is not central valued on R and F, G be two b-generalized skew derivations on R. The purpose of this article is to describe all possible forms of the b-generalized skew derivations F and G satisfying the identity… ▽ More Let R be a prime ring of characteristic not equal to 2, U be its Utumi quotient ring and C be the extended centroid of R. Let φbe a multilinear polynomial over C, which is not central valued on R and F, G be two b-generalized skew derivations on R. The purpose of this article is to describe all possible forms of the b-generalized skew derivations F and G satisfying the identity $F (u)u + uG(u) = G(u ^2)$, for all u \in {φ(ζ) | ζ = (ζ1 . . . , ζn) \inR^n}. Consequently, we discuss the cases when this identity acts as Jordan derivation and 2- Jordan multiplier on prime rings △ Less

Submitted 31 January, 2023; v1 submitted 6 December, 2022; originally announced December 2022.

MSC Class: 16N60; 16W25

Proof of The Generalized Zalcman Conjecture for Initial Coefficients of Univalent Functions

Authors: Vasudevarao Allu, Abhishek Pandey

Abstract: Let $\mathcal{S}$ denote the class of analytic and univalent ({\it i.e.}, one-to-one) functions $f(z)= z+\sum_{n=2}^{\infty}a_n z^n$ in the unit disk $\mathbb{D}=\{z\in \mathbb{C}:|z|<1\}$. For $f\in \mathcal{S}$, Ma proposed the generalized Zalcman conjecture that $$|a_{n}a_{m}-a_{n+m-1}|\le (n-1)(m-1),\,\,\,\mbox{ for } n\ge2,\, m\ge 2,$$ with equality only for the Koebe function… ▽ More Let $\mathcal{S}$ denote the class of analytic and univalent ({\it i.e.}, one-to-one) functions $f(z)= z+\sum_{n=2}^{\infty}a_n z^n$ in the unit disk $\mathbb{D}=\{z\in \mathbb{C}:|z|<1\}$. For $f\in \mathcal{S}$, Ma proposed the generalized Zalcman conjecture that $$|a_{n}a_{m}-a_{n+m-1}|\le (n-1)(m-1),\,\,\,\mbox{ for } n\ge2,\, m\ge 2,$$ with equality only for the Koebe function $k(z) = z/(1 - z)^2$ and its rotations. In this paper using the properties of holomorphic motion and Krushkal's Surgery Lemma \cite{Krushkal-1995}, we prove the generalized Zalcman conjecture when $n=2$, $m=3$ and $n=2$, $m=4$. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: 14 pages. arXiv admin note: text overlap with arXiv:2209.10595

MSC Class: 30C45; 30C50

Support points of some classes of analytic and univalent functions

Authors: Vasudevarao Allu, Abhishek Pandey

Abstract: Let $\mathcal{A}$ denote the class of analytic functions in the unit disk $\mathbb{D}:=\{z\in\mathbb{C}:|z|<1\}$ satisfying $f(0)=0$ and $f'(0)=1$. Let $\mathcal{U}$ be the class of functions $f\in\mathcal{A}$ satisfying $$\left|f'(z)\left(\frac{z}{f(z)}\right)^2-1 \right|< 1 \quad\mbox{ for } z\in\mathbb{D},$$ and $\mathscr{G}$ denote the class of functions $f\in \mathcal{A}$ satisfying… ▽ More Let $\mathcal{A}$ denote the class of analytic functions in the unit disk $\mathbb{D}:=\{z\in\mathbb{C}:|z|<1\}$ satisfying $f(0)=0$ and $f'(0)=1$. Let $\mathcal{U}$ be the class of functions $f\in\mathcal{A}$ satisfying $$\left|f'(z)\left(\frac{z}{f(z)}\right)^2-1 \right|< 1 \quad\mbox{ for } z\in\mathbb{D},$$ and $\mathscr{G}$ denote the class of functions $f\in \mathcal{A}$ satisfying $${\rm Re\,}\left(1+\frac{zf''(z)}{f'(z)}\right)>-\frac{1}{2} \quad\mbox{ for } z\in\mathbb{D}.$$ △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: 10 pages

MSC Class: 30C45; 30C50

Journal ref: Comptes Rendus. Mathématique, Volume 359 (2021) no. 4, pp. 465-473

On the generalized Zalcman conjecture

Authors: Vasudevarao Allu, Abhishek Pandey

Abstract: Let $\mathcal{S}$ denote the class of analytic and univalent ({\it i.e.}, one-to-one) functions $ f(z)= z+\sum_{n=2}^{\infty}a_n z^n$ in the unit disk $\mathbb{D}=\{z\in \mathbb{C}:|z|<1\}$. For $f\in \mathcal{S}$, In 1999, Ma proposed the generalized Zalcman conjecture that $$|a_{n}a_{m}-a_{n+m-1}|\le (n-1)(m-1),\,\,\,\mbox{ for } n\ge2,\, m\ge 2,$$ with equality only for the Koebe function… ▽ More Let $\mathcal{S}$ denote the class of analytic and univalent ({\it i.e.}, one-to-one) functions $ f(z)= z+\sum_{n=2}^{\infty}a_n z^n$ in the unit disk $\mathbb{D}=\{z\in \mathbb{C}:|z|<1\}$. For $f\in \mathcal{S}$, In 1999, Ma proposed the generalized Zalcman conjecture that $$|a_{n}a_{m}-a_{n+m-1}|\le (n-1)(m-1),\,\,\,\mbox{ for } n\ge2,\, m\ge 2,$$ with equality only for the Koebe function $k(z) = z/(1 - z)^2$ and its rotations. In the same paper, Ma \cite{Ma-1999} asked for what positive real values of $λ$ does the following inequality hold? \begin{equation}\label{conjecture} |λa_na_m-a_{n+m-1}|\le λnm -n-m+1 \,\,\,\,\, (n\ge 2, \,m\ge3). \end{equation} Clearly equality holds for the Koebe function $k(z) = z/(1 - z)^2$ and its rotations. In this paper, we prove the inequality (\ref{conjecture}) for $λ=3, n=2, m=3$. Further, we provide a geometric condition on extremal function maximizing (\ref{conjecture}) for $λ=2,n=2, m=3$. △ Less

Submitted 13 April, 2024; v1 submitted 21 September, 2022; originally announced September 2022.

Comments: Revised version, accepted for publication in Annali di Matematica Pura ed Applicata (1923 -), 11 Pages

MSC Class: 30C45; 30C50

Transverse Spectral Instabilities in Konopelchenko-Dubrovsky Equation

Authors: Bhavna, Ashish Kumar Pandey, Sudhir Singh

Abstract: We study the transverse spectral stability of the one-dimensional small-amplitude periodic traveling wave solutions of the (2+1)-dimensional Konopelchenko-Dubrovsky (KD) equation. We show that these waves are transversely unstable with respect to two-dimensional perturbations that are periodic in both directions with long wavelength in the transverse direction. We also show that these waves are tr… ▽ More We study the transverse spectral stability of the one-dimensional small-amplitude periodic traveling wave solutions of the (2+1)-dimensional Konopelchenko-Dubrovsky (KD) equation. We show that these waves are transversely unstable with respect to two-dimensional perturbations that are periodic in both directions with long wavelength in the transverse direction. We also show that these waves are transversely stable with respect to perturbations which are either mean-zero periodic or square-integrable in the direction of the propagation of the wave and periodic in the transverse direction with finite or short wavelength. We discuss the implications of these results for special cases of the KD equation - namely, KP-II and mKP-II equations. △ Less

Submitted 1 April, 2022; v1 submitted 23 March, 2022; originally announced March 2022.

Complexity of Paired Domination in AT-free and Planar Graphs

Authors: Vikash Tripathi, Ton Kloks, Arti Pandey, Kaustav Paul, Hung-Lung Wang

Abstract: For a graph $G=(V,E)$, a subset $D$ of vertex set $V$, is a dominating set of $G$ if every vertex not in $D$ is adjacent to atleast one vertex of $D$. A dominating set $D$ of a graph $G$ with no isolated vertices is called a paired dominating set (PD-set), if $G[D]$, the subgraph induced by $D$ in $G$ has a perfect matching. The Min-PD problem requires to compute a PD-set of minimum cardinality. T… ▽ More For a graph $G=(V,E)$, a subset $D$ of vertex set $V$, is a dominating set of $G$ if every vertex not in $D$ is adjacent to atleast one vertex of $D$. A dominating set $D$ of a graph $G$ with no isolated vertices is called a paired dominating set (PD-set), if $G[D]$, the subgraph induced by $D$ in $G$ has a perfect matching. The Min-PD problem requires to compute a PD-set of minimum cardinality. The decision version of the Min-PD problem remains NP-complete even when $G$ belongs to restricted graph classes such as bipartite graphs, chordal graphs etc. On the positive side, the problem is efficiently solvable for many graph classes including intervals graphs, strongly chordal graphs, permutation graphs etc. In this paper, we study the complexity of the problem in AT-free graphs and planar graph. The class of AT-free graphs contains cocomparability graphs, permutation graphs, trapezoid graphs, and interval graphs as subclasses. We propose a polynomial-time algorithm to compute a minimum PD-set in AT-free graphs. In addition, we also present a linear-time $2$-approximation algorithm for the problem in AT-free graphs. Further, we prove that the decision version of the problem is NP-complete for planar graphs, which answers an open question asked by Lin et al. (in Theor. Comput. Sci., $591 (2015): 99-105$ and Algorithmica, $ 82 (2020) :2809-2840$). △ Less

Submitted 10 December, 2021; originally announced December 2021.

Three-Phase Infeasibility Analysis for Distribution Grid Studies

Authors: Elizabeth Foster, Amritanshu Pandey, Larry Pileggi

Abstract: With the increase of distributed energy resources in the distribution grid, planning to ensure sufficient infrastructure and resources becomes critical. Planning at the distribution level is limited by the complexities of optimizing unbalanced systems. In this paper we develop a three-phase infeasibility analysis that identifies weak locations in a distribution network. This optimization is formul… ▽ More With the increase of distributed energy resources in the distribution grid, planning to ensure sufficient infrastructure and resources becomes critical. Planning at the distribution level is limited by the complexities of optimizing unbalanced systems. In this paper we develop a three-phase infeasibility analysis that identifies weak locations in a distribution network. This optimization is formulated by adding slack current sources at nodes in the system and minimizing their norm subject to distribution power flow constraints. Through this analysis we solve instances of power flow that would otherwise be infeasible and diverge. Under conditions when power flow is feasible, our approach is equivalent to standard three-phase power flow; however, for cases where power flow fails, the nonzero slack injection currents compensate for missing power to make the grid feasible. Since an uncountable number of injected currents can provide feasibility, we further explore the optimization formulation that best fits the solution objective through use of both a least squares and an L1 norm objective. Our L1 norm formulation localizes power deficient locations through its inherent sparsity. We show the efficacy of this approach on realistic unbalanced testcases up to 8500 nodes and for a scenario with a high penetration of electric vehicles. △ Less

Submitted 5 May, 2022; v1 submitted 21 October, 2021; originally announced October 2021.

Two-Stage Homotopy Method to Incorporate Discrete Control Variables into AC-OPF

Authors: Timothy McNamara, Amritanshu Pandey, Aayushya Agarwal, Larry Pileggi

Abstract: Alternating-Current Optimal Power Flow (AC-OPF) is an optimization problem critical for planning and operating the power grid. The problem is traditionally formulated using only continuous variables. Typically, control devices with discrete-valued settings, which provide valuable flexibility to the network and improve resilience, are omitted from AC-OPF formulations due to the difficulty of integr… ▽ More Alternating-Current Optimal Power Flow (AC-OPF) is an optimization problem critical for planning and operating the power grid. The problem is traditionally formulated using only continuous variables. Typically, control devices with discrete-valued settings, which provide valuable flexibility to the network and improve resilience, are omitted from AC-OPF formulations due to the difficulty of integrality constraints. We propose a two-stage homotopy algorithm to solve the AC-OPF problem with discrete-valued control settings. This method does not rely on prior knowledge of control settings or other initial conditions. The first stage relaxes the discrete settings to continuous variables and solves the optimization using a robust homotopy technique. Once the solution has been obtained using relaxed models, second homotopy problem gradually transforms the relaxed settings to their nearest feasible discrete values. We test the proposed algorithm on several large networks with switched shunts and adjustable transformers and show it can outperform a similar state-of-the-art solver. △ Less

Submitted 14 October, 2021; originally announced October 2021.

Comments: Under review: submitted for consideration for 22nd Power Systems Computation Conference

On the density of visible lattice points along polynomials

Authors: Sneha Chaubey, Ashish Kumar Pandey

Abstract: Recently, the notion of visibility from the origin has been generalized by viewing lattice points through curved lines of sights, where the family of curves considered are $y=mx^k$, $k\in\mathbb{N}$. In this note, we generalize the notion of visible lattice points for a given polynomial family of curves passing through the origin, and study the density of visible lattice points for this family. Th… ▽ More Recently, the notion of visibility from the origin has been generalized by viewing lattice points through curved lines of sights, where the family of curves considered are $y=mx^k$, $k\in\mathbb{N}$. In this note, we generalize the notion of visible lattice points for a given polynomial family of curves passing through the origin, and study the density of visible lattice points for this family. The density of visible lattice points for family of curves $y=mx^k$, $k\in\mathbb{N}$ is well understood as one has nice arithmetic interpretations in terms of a generalized gcd function, which seems to be absent for general polynomial families. We pose "Visibility density conjecture" regarding the density of visible lattice points for polynomial families passing through the origin, and show some numerical results supporting the conjecture. We obtain a lower bound on the density for a class of quadratic families. In addition, we discuss some ideas of the proof of the conjecture for the simplest quadratic family. △ Less

Submitted 17 September, 2021; originally announced September 2021.

A linear-time algorithm for semitotal domination in strongly chordal graphs

Authors: Vikash Tripathi, Arti Pandey, Anil Maheshwari

Abstract: In a graph $G=(V,E)$ with no isolated vertex, a dominating set $D \subseteq V$, is called a semitotal dominating set if for every vertex $u \in D$ there is another vertex $v \in D$, such that distance between $u$ and $v$ is at most two in $G$. Given a graph $G=(V,E)$ without isolated vertices, the Minimum Semitotal Domination problem is to find a minimum cardinality semitotal dominating set of… ▽ More In a graph $G=(V,E)$ with no isolated vertex, a dominating set $D \subseteq V$, is called a semitotal dominating set if for every vertex $u \in D$ there is another vertex $v \in D$, such that distance between $u$ and $v$ is at most two in $G$. Given a graph $G=(V,E)$ without isolated vertices, the Minimum Semitotal Domination problem is to find a minimum cardinality semitotal dominating set of $G$. The semitotal domination number, denoted by $γ_{t2}(G)$, is the minimum cardinality of a semitotal dominating set of $G$. The decision version of the problem remains NP-complete even when restricted to chordal graphs, chordal bipartite graphs, and planar graphs. Galby et al. in [6] proved that the problem can be solved in polynomial time for bounded MIM-width graphs which includes many well known graph classes, but left the complexity of the problem in strongly chordal graphs unresolved. Henning and Pandey in [20] also asked to resolve the complexity status of the problem in strongly chordal graphs. In this paper, we resolve the complexity of the problem in strongly chordal graphs by designing a linear-time algorithm for the problem. △ Less

Submitted 5 September, 2021; originally announced September 2021.

Transverse spectral instability in generalized Kadomtsev-Petviashvili equation

Authors: Bhavna, Atul Kumar, Ashish Kumar Pandey

Abstract: We study transverse stability and instability of one-dimensional small-amplitude periodic traveling waves of a generalized Kadomtsev-Petviashvili equation with respect to two-dimensional perturbations, which are either periodic or square-integrable in the direction of the propagation of the underlying one-dimensional wave and periodic in the transverse direction. We obtain transverse instability r… ▽ More We study transverse stability and instability of one-dimensional small-amplitude periodic traveling waves of a generalized Kadomtsev-Petviashvili equation with respect to two-dimensional perturbations, which are either periodic or square-integrable in the direction of the propagation of the underlying one-dimensional wave and periodic in the transverse direction. We obtain transverse instability results in KP-fKdV, KP-ILW, and KP-Whitham equations. Moreover, assuming the spectral stability of one-dimensional wave with respect to one-dimensional square-integrable periodic perturbations, we obtain transverse stability results in aforementioned equations. △ Less

Submitted 31 March, 2022; v1 submitted 1 September, 2021; originally announced September 2021.

A Novel Bivariate Generalized Weibull Distribution with Properties and Applications

Authors: Ashok Kumar Pathak, Mohd. Arshad, Qazi J. Azhad, Mukti Khetan, Arvind Pandey

Abstract: Univariate Weibull distribution is a well-known lifetime distribution and has been widely used in reliability and survival analysis. In this paper, we introduce a new family of bivariate generalized Weibull (BGW) distributions, whose univariate marginals are exponentiated Weibull distribution. Different statistical quantiles like marginals, conditional distribution, conditional expectation, produc… ▽ More Univariate Weibull distribution is a well-known lifetime distribution and has been widely used in reliability and survival analysis. In this paper, we introduce a new family of bivariate generalized Weibull (BGW) distributions, whose univariate marginals are exponentiated Weibull distribution. Different statistical quantiles like marginals, conditional distribution, conditional expectation, product moments, correlation and a measure component reliability are derived. Various measures of dependence and statistical properties along with ageing properties are examined. Further, the copula associated with BGW distribution and its various important properties are also considered. The methods of maximum likelihood and Bayesian estimation are employed to estimate unknown parameters of the model. A Monte Carlo simulation and real data study are carried out to demonstrate the performance of the estimators and results have proven the effectiveness of the distribution in real-life situations △ Less

Submitted 26 July, 2021; originally announced July 2021.

Point Singularities in Incompatible Elasticity

Authors: Animesh Pandey, Anurag Gupta

Abstract: The equations of stress equilibrium and strain compatibility/incompatibility are discussed for fields with point singularities in a planar domain. The sufficiency (or insufficiency) of the smooth maps, obtained by restricting the singular fields to the domain away from the singularities, in completely characterizing the equations of equilibrium and compatibility/incompatibility over the entire dom… ▽ More The equations of stress equilibrium and strain compatibility/incompatibility are discussed for fields with point singularities in a planar domain. The sufficiency (or insufficiency) of the smooth maps, obtained by restricting the singular fields to the domain away from the singularities, in completely characterizing the equations of equilibrium and compatibility/incompatibility over the entire domain, is established and illustrated with examples. The uniqueness of the solution to the stress problem of incompatible linear elasticity, allowing for singular fields, is proved. The uniqueness fails when the problem is considered solely in terms of the restricted maps. As applications of our framework, a general stress solution, in response to point supported body force and defect fields, is derived and a generalized notion of the force acting on a defect is developed. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Fast accurate approximation of convolutions with weakly singular kernel and its applications

Authors: Awanish Kumar Tiwari, Ambuj Pandey, Jagabandhu Paul, Akash Anand

Abstract: In this article, we present an $O(N \log N)$ rapidly convergent algorithm for the numerical approximation of the convolution integral with radially symmetric weakly singular kernels and compactly supported densities. To achieve the reduced computational complexity, we utilize the Fast Fourier Transform (FFT) on a uniform grid of size $N$ for approximating the convolution. To facilitate this and ma… ▽ More In this article, we present an $O(N \log N)$ rapidly convergent algorithm for the numerical approximation of the convolution integral with radially symmetric weakly singular kernels and compactly supported densities. To achieve the reduced computational complexity, we utilize the Fast Fourier Transform (FFT) on a uniform grid of size $N$ for approximating the convolution. To facilitate this and maintain the accuracy, we primarily rely on a periodic Fourier extension of the density with a suitably large period depending on the support of the density. The rate of convergence of the method increases with increasing smoothness of the periodic extension and, in fact, approximations exhibit super-algebraic convergence when the extension is infinitely differentiable. Furthermore, when the density has jump discontinuities, we utilize a certain Fourier smoothing technique to accelerate the convergence to achieve the quadratic rate in the overall approximation. Finally, we apply the integration scheme for numerical solution of certain partial differential equations. Moreover, we apply the quadrature to obtain a fast and high-order Nystöm solver for the solution of the Lippmann-Schwinger integral equation. We validate the performance of the proposed scheme in terms of accuracy as well as computational efficiency through a variety of numerical experiments. △ Less

Submitted 8 July, 2021; originally announced July 2021.

Comments: 22 pages, 18 figures

MSC Class: 65R10; 65R20; 44A35; 35J05; 78A45; 31A10

High-frequency instabilities of the Ostrovsky equation

Authors: Bhavna, Atul Kumar, Ashish Kumar Pandey

Abstract: We study spectral stability of small amplitude periodic traveling waves of the Ostrovsky equation. We prove that these waves exhibit spectral instabilities arising from a collision of pair of non-zero eigenvalues on the imaginary axis when subjected to square integrable perturbations on the whole real line. We also list all such collisions between pair of eigenvalues on the imaginary axis and do a… ▽ More We study spectral stability of small amplitude periodic traveling waves of the Ostrovsky equation. We prove that these waves exhibit spectral instabilities arising from a collision of pair of non-zero eigenvalues on the imaginary axis when subjected to square integrable perturbations on the whole real line. We also list all such collisions between pair of eigenvalues on the imaginary axis and do a Krein signature analysis. △ Less

Submitted 5 July, 2021; originally announced July 2021.

A Neumann type problem on an unbounded domain in the Heisenberg group

Authors: Ashutosh Pandey, Mukund Madhav Mishra, Shivani Dubey

Abstract: We discuss the wellposedness of the Neumann problem on a half-space for the Kohn-Laplacian in the Heisenberg group. We then construct the Neumann function and explicitly represent the solution of the associated inhomogeneous problem. We discuss the wellposedness of the Neumann problem on a half-space for the Kohn-Laplacian in the Heisenberg group. We then construct the Neumann function and explicitly represent the solution of the associated inhomogeneous problem. △ Less

Submitted 27 May, 2021; originally announced May 2021.

MSC Class: 31B20; 35H20; 35N15; 45B05; 65N80

Posterior distributions for Hierarchical Spike and Slab Indian Buffet processes

Authors: Lancelot F. James, Juho Lee, Abhinav Pandey

Abstract: Bayesian nonparametric hierarchical priors are highly effective in providing flexible models for latent data structures exhibiting sharing of information between and across groups. Most prominent is the Hierarchical Dirichlet Process (HDP), and its subsequent variants, which model latent clustering between and across groups. The HDP, may be viewed as a more flexible extension of Latent Dirichlet A… ▽ More Bayesian nonparametric hierarchical priors are highly effective in providing flexible models for latent data structures exhibiting sharing of information between and across groups. Most prominent is the Hierarchical Dirichlet Process (HDP), and its subsequent variants, which model latent clustering between and across groups. The HDP, may be viewed as a more flexible extension of Latent Dirichlet Allocation models (LDA), and has been applied to, for example, topic modelling, natural language processing, and datasets arising in health-care. We focus on analogous latent feature allocation models, where the data structures correspond to multisets or unbounded sparse matrices. The fundamental development in this regard is the Hierarchical Indian Buffet process (HIBP), which utilizes a hierarchy of Beta processes over J groups, where each group generates binary random matrices, reflecting within group sharing of features, according to beta-Bernoulli IBP priors. To encompass HIBP versions of non-Bernoulli extensions of the IBP, we introduce hierarchical versions of general spike and slab IBP. We provide explicit novel descriptions of the marginal, posterior and predictive distributions of the HIBP and its generalizations which allow for exact sampling and simpler practical implementation. We highlight common structural properties of these processes and establish relationships to existing IBP type and related models arising in the literature. Examples of potential applications may involve topic models, Poisson factorization models, random count matrix priors and neural network models △ Less

Submitted 21 March, 2021; originally announced March 2021.

Comments: 4 figures

MSC Class: 60C05; 60G09 (Primary); 60G57; 60E99 (Secondary)

Incremental Model Building Homotopy Approach for Solving Exact AC-Constrained Optimal Power Flow

Authors: Amritanshu Pandey, Aayushya Agarwal, Larry Pileggi

Abstract: Alternating-Current Optimal Power Flow (AC-OPF) is framed as a NP-hard non-convex optimization problem that solves for the most economical dispatch of grid generation given the AC-network and device constraints. Although there are no standard methodologies for obtaining the global optimum for the problem, there is considerable interest from planning and operational engineers in finding a local opt… ▽ More Alternating-Current Optimal Power Flow (AC-OPF) is framed as a NP-hard non-convex optimization problem that solves for the most economical dispatch of grid generation given the AC-network and device constraints. Although there are no standard methodologies for obtaining the global optimum for the problem, there is considerable interest from planning and operational engineers in finding a local optimum. Nonetheless, solving for the local optima of a large AC-OPF problem is challenging and time-intensive, as none of the leading non-linear optimization toolboxes can provide any timely guarantees of convergence. To provide robust local convergence for large complex systems, we introduce a homotopy-based approach that solves a sequence of primal-dual interior point problems. We utilize the physics of the grid to develop the proposed homotopy method and demonstrate the efficacy of this approach on U.S. Eastern Interconnection sized test networks. △ Less

Submitted 1 November, 2020; originally announced November 2020.

Comments: To be published in Proceedings of Hawaii International Conference on System Sciences-54, Hawaii, 2021

Fourier smoothed pre-corrected trapezoidal rule for solution of Lippmann-Schwinger integral equation

Authors: Ambuj Pandey, Akash Anand

Abstract: For the numerical solution of the Lippmann-Schwinger equation, while the pre-corrected trapezoidal rule converges with high-order for smooth compactly supported densities, it exhibits only the linear convergence in the case of discontinuity in material properties across the interface. In this short article, we propose a Nyström solver based on "Fourier smoothed pre-corrected trapezoidal rule" th… ▽ More For the numerical solution of the Lippmann-Schwinger equation, while the pre-corrected trapezoidal rule converges with high-order for smooth compactly supported densities, it exhibits only the linear convergence in the case of discontinuity in material properties across the interface. In this short article, we propose a Nyström solver based on "Fourier smoothed pre-corrected trapezoidal rule" that converges with second order for such scattering problems while maintaining the computational complexity of $O(N \log N)$. Moreover, the method is not only very simple to implement, it is also applicable to problems with geometrically complex inhomogeneities including those with corners and cusps. We present a variety of numerical experiments including comparative studies with competing approaches reported in [J. Comput. Phys., 200(2) (2004), 670--694] by Bruno and Hyde, and in [J. Fourier Anal. Appl., 11(4) (2005), 471-487 ] by Andersson and Holst to exemplify its performance in terms of speed and accuracy. This Fourier smoothed numerical integration scheme can also be adapted to other problems of interest where the convolution integral with discontinuous density is required to be computed. △ Less

Submitted 13 July, 2020; originally announced July 2020.

Zalcman Conjecture for certain analytic and univalent functions

Authors: Vasudevarao Allu, Abhishek Pandey

Abstract: Let $\mathcal{A}$ denote the class of analytic functions in the unit disk $\mathbb{D}$ of the form $f(z)= z+\sum_{n=2}^{\infty}a_n z^n$ and $\mathcal{S}$ denote the class of functions $f\in\mathcal{A}$ which are univalent ({\it i.e.}, one-to-one). In 1960s, L. Zalcman conjectured that $|a_n^2-a_{2n-1}|\le (n-1)^2$ for $n\ge 2$, which implies the famous Bieberbach conjecture $|a_n|\le n$ for… ▽ More Let $\mathcal{A}$ denote the class of analytic functions in the unit disk $\mathbb{D}$ of the form $f(z)= z+\sum_{n=2}^{\infty}a_n z^n$ and $\mathcal{S}$ denote the class of functions $f\in\mathcal{A}$ which are univalent ({\it i.e.}, one-to-one). In 1960s, L. Zalcman conjectured that $|a_n^2-a_{2n-1}|\le (n-1)^2$ for $n\ge 2$, which implies the famous Bieberbach conjecture $|a_n|\le n$ for $n\ge 2$. For $f\in \mathcal{S}$, Ma \cite{Ma-1999} proposed a generalized Zalcman conjecture $$|a_{n}a_{m}-a_{n+m-1}|\le (n-1)(m-1) $$ for $n\ge 2, m\ge 2$. Let $\mathcal{U}$ be the class of functions $f\in\mathcal{A}$ satisfying $$ \left|f'(z)\left(\frac{z}{f(z)}\right)^2-1 \right|< 1 \quad\mbox{ for } z\in\mathbb{D}. $$ and $\mathcal{F}$ denote the class of functions $f\in \mathcal{A}$ satisfying ${\rm Re\,}(1-z)^{2}f'(z)>0$ in $\mathbb{D}$. In the present paper, we prove the Zalcman conjecture and generalized Zalcman conjecture for the class $\mathcal{U}$ using extreme point theory. We also prove the Zalcman conjecture and generalized Zalcman conjecture for the class $\mathcal{F}$ for the initial coefficients. △ Less

Submitted 13 June, 2020; originally announced June 2020.

Comments: 15 pages

MSC Class: 30C45; 30C50

How many zeros of a random sparse polynomial are real?

Authors: Gorav Jindal, Anurag Pandey, Himanshu Shukla, Charilaos Zisopoulos

Abstract: We investigate the number of real zeros of a univariate $k$-sparse polynomial $f$ over the reals, when the coefficients of $f$ come from independent standard normal distributions. Recently Bürgisser, Ergür and Tonelli-Cueto showed that the expected number of real zeros of $f$ in such cases is bounded by $O(\sqrt{k} \log k)$. In this work, we improve the bound to $O(\sqrt{k})$ and also show that th… ▽ More We investigate the number of real zeros of a univariate $k$-sparse polynomial $f$ over the reals, when the coefficients of $f$ come from independent standard normal distributions. Recently Bürgisser, Ergür and Tonelli-Cueto showed that the expected number of real zeros of $f$ in such cases is bounded by $O(\sqrt{k} \log k)$. In this work, we improve the bound to $O(\sqrt{k})$ and also show that this bound is tight by constructing a family of sparse support whose expected number of real zeros is lower bounded by $Ω(\sqrt{k})$. Our main technique is an alternative formulation of the Kac integral by Edelman-Kostlan which allows us to bound the expected number of zeros of $f$ in terms of the expected number of zeros of polynomials of lower sparsity. Using our technique, we also recover the $O(\log n)$ bound on the expected number of real zeros of a dense polynomial of degree $n$ with coefficients coming from independent standard normal distributions. △ Less

Submitted 6 November, 2019; originally announced November 2019.

Variety Membership Testing, Algebraic Natural Proofs, and Geometric Complexity Theory

Authors: Markus Bläser, Christian Ikenmeyer, Vladimir Lysikov, Anurag Pandey, Frank-Olaf Schreyer

Abstract: We study the variety membership testing problem in the case when the variety is given as an orbit closure and the ambient space is the set of all 3-tensors. The first variety that we consider is the slice rank variety, which consists of all 3-tensors of slice rank at most $r$. We show that the membership testing problem for the slice rank variety is $\NP$-hard. While the slice rank variety is a un… ▽ More We study the variety membership testing problem in the case when the variety is given as an orbit closure and the ambient space is the set of all 3-tensors. The first variety that we consider is the slice rank variety, which consists of all 3-tensors of slice rank at most $r$. We show that the membership testing problem for the slice rank variety is $\NP$-hard. While the slice rank variety is a union of orbit closures, we define another variety, the minrank variety, expressible as a single orbit closure. Our next result is the $\NP$-hardness of membership testing in the minrank variety, hence we establish the $\NP$-hardness of the orbit closure containment problem for 3-tensors. Algebraic natural proofs were recently introduced by Forbes, Shpilka and Volk and independently by Grochow, Kumar, Saks and Saraf. Bläser et al. gave a version of an algebraic natural proof barrier for the matrix completion problem which relies on $\coNP \subseteq \exists \BPP$. It implied that constructing equations for the corresponding variety should be hard. We generalize their approach to work with any family of varieties for which the membership problem is $\NP$-hard and for which we can efficiently generate a dense subset. Therefore, a similar barrier holds for the slice rank and the minrank varieties, too. This allows us to set up the slice rank and the minrank varieties as a test-bed for geometric complexity theory (GCT). We determine the stabilizers of the tensors that generate the orbit closures of the two varieties and prove that these tensors are almost characterized by their symmetries. We prove several nontrivial equations for both the varieties using different GCT methods. Many equations also work in the regime where membership testing in the slice rank or minrank varieties is $\NP$-hard. We view this as a promising sign that the GCT approach might indeed be successful. △ Less

Submitted 6 November, 2019; originally announced November 2019.

The extended xgamma distribution

Authors: Mahendra Saha, Abhimanyu Singh Yadav, Arvind Pandey, Shivanshi Shukla, Sudhansu S Maiti

Abstract: This article aims to introduced a new distribution named as extended xgamma (EXg) distribution. This generalization is derived from xgamma distribution (Xg), a special finite mixture of exponential and gamma distributions [see, Sen et al. ($2016$)]. Some important statistical properties, viz., survival characteristics, moments, mean deviation and random number generation have been derived. Further… ▽ More This article aims to introduced a new distribution named as extended xgamma (EXg) distribution. This generalization is derived from xgamma distribution (Xg), a special finite mixture of exponential and gamma distributions [see, Sen et al. ($2016$)]. Some important statistical properties, viz., survival characteristics, moments, mean deviation and random number generation have been derived. Further, maximum likelihood estimation for the estimation of the unknown parameters have also been discussed for the complete sample. The application of the proposed model has been illustrated through a real data set and observed that the proposed model might be taken as an better alternative to some well known lifetime distributions. △ Less

Submitted 30 August, 2019; originally announced September 2019.

Comments: 8 pages, 2 figures

Direct/iterative hybrid solver for scattering by inhomogeneous media

Authors: Oscar P. Bruno, Ambuj Pandey

Abstract: This paper presents a fast high-order method for the solution of two-dimensional problems of scattering by penetrable inhomogeneous media, with application to high-frequency configurations containing (possibly) discontinuous refractivities. The method relies on a hybrid direct/iterative combination of 1)~A differential volumetric formulation (which is based on the use of appropriate Chebyshev diff… ▽ More This paper presents a fast high-order method for the solution of two-dimensional problems of scattering by penetrable inhomogeneous media, with application to high-frequency configurations containing (possibly) discontinuous refractivities. The method relies on a hybrid direct/iterative combination of 1)~A differential volumetric formulation (which is based on the use of appropriate Chebyshev differentiation matrices enacting the Laplace operator) and, 2)~A second-kind boundary integral formulation. The approach enjoys low dispersion and high-order accuracy for smooth refractivities, as well as second-order accuracy (while maintaining low dispersion) in the discontinuous refractivity case. The solution approach proceeds by application of Impedance-to-Impedance (ItI) maps to couple the volumetric and boundary discretizations. The volumetric linear algebra solutions are obtained by means of a multifrontal solver, and the coupling with the boundary integral formulation is achieved via an application of the iterative linear-algebra solver GMRES. In particular, the existence and uniqueness theory presented in the present paper provides an affirmative answer to an open question concerning the existence of a uniquely solvable second-kind ItI-based formulation for the overall scattering problem under consideration. Relying on a modestly-demanding scatterer-dependent precomputation stage (requiring in practice a computing cost of the order of $O(N^α)$ operations, with $α\approx 1.07$, for an $N$-point discretization), together with fast ($O(N)$-cost) single-core runs for each incident field considered, the proposed algorithm can effectively solve scattering problems for large and complex objects possibly containing strong refractivity contrasts and discontinuities. △ Less

Submitted 28 July, 2023; v1 submitted 12 July, 2019; originally announced July 2019.

Lagrangian coherent sets in turbulent Rayleigh-Bénard convection

Authors: Christiane Schneide, Martin Stahn, Ambrish Pandey, Oliver Junge, Péter Koltai, Kathrin Padberg-Gehle, Jörg Schumacher

Abstract: Coherent circulation rolls and their relevance for the turbulent heat transfer in a two-dimensional Rayleigh--Bénard convection model are analyzed. The flow is in a closed cell of aspect ratio four at a Rayleigh number ${\rm Ra}=10^6$ and at a Prandtl number ${\rm Pr}=10$. Three different Lagrangian analysis techniques based on graph Laplacians -- distance spectral trajectory clustering, time-aver… ▽ More Coherent circulation rolls and their relevance for the turbulent heat transfer in a two-dimensional Rayleigh--Bénard convection model are analyzed. The flow is in a closed cell of aspect ratio four at a Rayleigh number ${\rm Ra}=10^6$ and at a Prandtl number ${\rm Pr}=10$. Three different Lagrangian analysis techniques based on graph Laplacians -- distance spectral trajectory clustering, time-averaged diffusion maps and finite-element based dynamic Laplacian discretization -- are used to monitor the turbulent fields along trajectories of massless Lagrangian particles in the evolving turbulent convection flow. The three methods are compared to each other and the obtained coherent sets are related to results from an analysis in the Eulerian frame of reference. We show that the results of these methods agree with each other and that Lagrangian and Eulerian coherent sets form basically a disjoint union of the flow domain. Additionally, a windowed time-averaging of variable interval length is performed to study the degree of coherence as a function of this additional coarse graining which removes small-scale fluctuations that cause trajectories to disperse quickly. Finally, the coherent set framework is extended to study heat transport. △ Less

Submitted 7 November, 2019; v1 submitted 8 July, 2019; originally announced July 2019.

Journal ref: Phys. Rev. E 100, 053103 (2019)

Multilevel Monte Carlo Acceleration of Seismic Wave Propagation under Uncertainty

Authors: Marco Ballesio, Joakim Beck, Anamika Pandey, Laura Parisi, Erik von Schwerin, Raul Tempone

Abstract: We interpret uncertainty in a model for seismic wave propagation by treating the model parameters as random variables, and apply the Multilevel Monte Carlo (MLMC) method to reduce the cost of approximating expected values of selected, physically relevant, quantities of interest (QoI) with respect to the random variables. Targeting source inversion problems, where the source of an earthquake is inf… ▽ More We interpret uncertainty in a model for seismic wave propagation by treating the model parameters as random variables, and apply the Multilevel Monte Carlo (MLMC) method to reduce the cost of approximating expected values of selected, physically relevant, quantities of interest (QoI) with respect to the random variables. Targeting source inversion problems, where the source of an earthquake is inferred from ground motion recordings on the Earth's surface, we consider two QoI that measure the discrepancies between computed seismic signals and given reference signals: one QoI, $\hbox{QoI}_E$, is defined in terms of the $L^2$-misfit, which is directly related to maximum likelihood estimates of the source parameters; the other, $\hbox{QoI}_W$, is based on the quadratic Wasserstein distance between probability distributions, and represents one possible choice in a class of such misfit functions that have become increasingly popular to solve seismic inversion in recent years. We simulate seismic wave propagation, including seismic attenuation, using a publicly available code in widespread use, based on the spectral element method. Using random coefficients and deterministic initial and boundary data, we present benchmark numerical experiments with synthetic data in a two-dimensional physical domain and a one-dimensional velocity model where the assumed parameter uncertainty is motivated by realistic Earth models. Here, the computational cost of the standard Monte Carlo method was reduced by up to 97% for $\hbox{QoI}_E$, and up to 78% for $\hbox{QoI}_W$, using a relevant range of tolerances. Shifting to three-dimensional domains is straight-forward and will further increase the relative computational work reduction. △ Less

Submitted 5 September, 2019; v1 submitted 3 October, 2018; originally announced October 2018.

Improved convergence of fast integral equation solvers for acoustic scattering by inhomogeneous penetrable media with discontinuous material interface

Authors: Ambuj Pandey, Akash Anand

Abstract: In recent years, several fast solvers for the solution of the Lippmann-Schwinger integral equation that mathematically models the scattering of time-harmonic acoustic waves by penetrable inhomogeneous obstacles, have been proposed. While many of these fast methodologies exhibit rapid convergence for smoothly varying scattering configurations, the rate for most of them reduce to either linear or qu… ▽ More In recent years, several fast solvers for the solution of the Lippmann-Schwinger integral equation that mathematically models the scattering of time-harmonic acoustic waves by penetrable inhomogeneous obstacles, have been proposed. While many of these fast methodologies exhibit rapid convergence for smoothly varying scattering configurations, the rate for most of them reduce to either linear or quadratic when material properties are allowed to jump across the interface. A notable exception to this is a recently introduced Nyström scheme [J. Comput. Phys., 311 (2016), 258--274] that utilizes a specialized quadrature in the boundary region for a high-order treatment of the material interface. In this text, we present a solution framework that relies on the specialized boundary integrator to enhance the convergence rate of other fast, low order methodologies without adding to their computational complexity of $O(N \log N)$ for an $N$-point discretization. In particular, to demonstrate the efficacy of the proposed framework, we explain its implementation to enhance the order to convergence of two schemes, one introduced by Duan and Rokhlin [J. Comput. Phys., 228(6) (2009), 2152--2174] that is based on a pre-corrected trapezoidal rule while the other by Bruno and Hyde [J. Comput. Phys., 200(2) (2004), 670--694] which relies on a suitable decomposition of the Green's function via Addition theorem. In addition to a detailed description of these methodologies, we also present a comparative performance study of the improved versions of these two and the Nyström solver in [J. Comput. Phys., 311 (2016), 258--274] through a wide range of numerical experiments. △ Less

Submitted 1 June, 2018; originally announced June 2018.

Some Properties of Kenmotsu Manifolds Admitting a Semi-symmetric Non-metric Connection

Authors: S. K. Chaubey, A. C. Pandey, N. V. C. Shukla

Abstract: The aim of this paper is to study generalized recurrent, generalized Ricci-recurrent, weakly symmetric and weakly Ricci-symmetric Kenmotsu manifolds with respect to the semi-symmetric non-metric connection. The aim of this paper is to study generalized recurrent, generalized Ricci-recurrent, weakly symmetric and weakly Ricci-symmetric Kenmotsu manifolds with respect to the semi-symmetric non-metric connection. △ Less

Submitted 9 January, 2018; originally announced January 2018.

Comments: 13

MSC Class: $53C15$; $53B05$; $53C25$

Comparison of modulational instabilities in full-dispersion shallow water models

Authors: Ashish Kumar Pandey

Abstract: We study the modulational instability of a shallow water model, with and without surface tension, which generalizes the Whitham equation to include bi-directional propagation. Without surface tension, the small amplitude periodic traveling waves are modulationally unstable if their wave number is greater than a critical wave number predicting a Benjamin-Feir type instability and the result qualita… ▽ More We study the modulational instability of a shallow water model, with and without surface tension, which generalizes the Whitham equation to include bi-directional propagation. Without surface tension, the small amplitude periodic traveling waves are modulationally unstable if their wave number is greater than a critical wave number predicting a Benjamin-Feir type instability and the result qualitatively agrees with the shallow water model in [HP16]. With surface tension, the result qualitatively agrees with the physical problem except for the large surface tension limit which is accurately predicted by the shallow water model in [HP16]. We also compare the results with the Whitham and full-dispersion Camassa-Holm equations. △ Less

Submitted 1 August, 2017; originally announced August 2017.