BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Govind Menon

Govind Menon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
nlin.AOgr-qcmath.NAmath.PRastro-ph.HEmath.APnlin.SIastro-ph.COcond-mat.stat-mechMachine Learningmath-phmath.DGmath.DSmath.MP

Catalog footprint

What is connected

20works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2026arXiv

An entropy formula for the Deep Linear Network

We study the Riemannian geometry of the Deep Linear Network (DLN) as a foundation for a thermodynamic description of the learning process. The main tools are the use of group actions to analyze overparametrization and the use of Riemannian submersion from the space of parameters to the space of observables. The foliation of the balanced manifold in the parameter space by group orbits is used to define and compute a Boltzmann entropy. We also show that the Riemannian geometry on the space of observables defined in [2] is obtained by Riemannian submersion of the balanced manifold. The main technical step is an explicit construction of an orthonormal basis for the tangent space of the balanced manifold using the theory of Jacobi matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Null and Force-Free Electromagnetic Field

This paper describes the electrodynamics of a null and force-free field in completely geometric terms. As was previously established in \cite{Menon_FF20}, solutions to force-free electrodynamics are governed by the existence of certain special types of foliations of spacetime. Here the nature of the foliations in a coordinate-free formalism in the null case is prescribed. All of the general results are illustrated by constructing a null, force-free electrodynamic field in a Friedmann-Lemaitre-Robertson- Walker (FLRW) spacetime.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Smoothed Analysis for the Conjugate Gradient Algorithm

The purpose of this paper is to establish bounds on the rate of convergence of the conjugate gradient algorithm when the underlying matrix is a random positive definite perturbation of a deterministic positive definite matrix. We estimate all finite moments of a natural halting time when the random perturbation is drawn from the Laguerre unitary ensemble in a critical scaling regime explored in Deift et al. (2016). These estimates are used to analyze the expected iteration count in the framework of smoothed analysis, introduced by Spielman and Teng (2001). The rigorous results are compared with numerical calculations in several cases of interest.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Energy landscapes for the self-assembly of supramolecular polyhedra

We develop a mathematical model for the energy landscape of polyhedral supramolecular cages recently synthesized by self-assembly [Sun et al., Science 2010]. Our model includes two essential features of the experiment: (i) geometry of the organic ligands and metallic ions; and (ii) combinatorics. The molecular geometry is used to introduce an energy that favors square-planar vertices (modeling $\mathrm{Pd}^{2+}$ ions) and bent edges with one of two preferred opening angles (modeling boomerang-shaped ligands of two types). The combinatorics of the model involve $2$-colorings of edges of polyhedra with $4$-valent vertices. The set of such $2$-colorings, quotiented by the octahedral symmetry group, has a natural graph structure, and is called the combinatorial configuration space. The energy landscape of our model is the energy of each state in the combinatorial configuration space. The challenge in the computation of the energy landscape is a combinatorial explosion in the number of $2$-colorings of edges. We describe sampling methods based on the symmetries of the configurations and connectivity of the configuration graph. When the two preferred opening angles encompass the geometrically ideal angle, the energy landscape exhibits a very low-energy minimum for the most symmetric configuration at equal mixing of the two angles, even when the average opening angle does not match the ideal angle.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Force-free Currents and the Newman-Penrose Tetrad of a Kerr Black Hole: Exact Local Solutions

In a previous article we derived a class of solutions to the force-free magnetosphere in a Kerr background. Here, the streaming surface, defined by constant values of the toriodal component of the electromagnetic vector potential $A$, were generated by constant values of $θ$. The electromagnetic current vector flowed along the in-falling principle null geodesic vector of the geometry. Subsequently, we generalized this to obtain an out-going principle null geodesic vector as well. In this article, we derive solutions that are complimentary to the above mentioned criteria. Namely, here the solution has a streaming surface generated by spheres of constant radial coordinate $r$, and the current vector is generated by linear combinations of $m$ and $m^\star$, the remaining bases vectors in the Newman-Penrose null tetrad.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Local, Non-Geodesic, Timelike Currents in the Force-Free Magnetosphere of a Kerr Black Hole

In this paper, we use previously developed exact solutions to present some of the curious features of a force-free magnetosphere in a Kerr background. More precisely, we obtain a hitherto unseen timelike current in the force-free magnetosphere that does not flow along a geodesic. The electromagnetic field in this case happens to be magnetically dominated. This too is a feature that has entered the literature for the first time. Changing the sign of a single parameter in our solutions generates a spacelike current that creates an electromagnetic field that is electrically dominated.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the condition number of the critically-scaled Laguerre Unitary Ensemble

We consider the Laguerre Unitary Ensemble (aka, Wishart Ensemble) of sample covariance matrices $A = XX^*$, where $X$ is an $N \times n$ matrix with iid standard complex normal entries. Under the scaling $n = N + \lfloor \sqrt{ 4 c N} \rfloor$, $c > 0$ and $N \rightarrow \infty$, we show that the rescaled fluctuations of the smallest eigenvalue, largest eigenvalue and condition number of the matrices $A$ are all given by the Tracy--Widom distribution ($β= 2$). This scaling is motivated by the study of the solution of the equation $Ax=b$ using the conjugate gradient algorithm, in the case that $A$ and $b$ are random: For such a scaling the fluctuations of the halting time for the algorithm are empirically seen to be universal.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Coarse graining, dynamic renormalization and the kinetic theory of shock clustering

We demonstrate the utility of the equation free methodology developed by one of the authors (I.G.K) for the study of scalar conservation laws with disordered initial conditions. The numerical scheme is benchmarked on exact solutions in Burgers turbulence corresponding to Levy process initial data. For these initial data, the kinetics of shock clustering is described by Smoluchowski's coagulation equation with additive kernel. The equation free methodology is used to develop a particle scheme that computes self-similar solutions to the coagulation equation, including those with fat tails.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Universality in Numerical Computations with Random Data. Case Studies

The authors present evidence for universality in numerical computations with random data. Given a (possibly stochastic) numerical algorithm with random input data, the time (or number of iterations) to convergence (within a given tolerance) is a random variable, called the halting time. Two-component universality is observed for the fluctuations of the halting time, i.e., the histogram for the halting times, centered by the sample average and scaled by the sample variance, collapses to a universal curve, independent of the input data distribution, as the dimension increases. Thus, up to two components, the sample average and the sample variance, the statistics for the halting time are universally prescribed. The case studies include six standard numerical algorithms, as well as a model of neural computation and decision making. A link to relevant software is provided in for the reader who would like to do computations of his'r own.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How long does it take to compute the eigenvalues of a random symmetric matrix?

We present the results of an empirical study of the performance of the QR algorithm (with and without shifts) and the Toda algorithm on random symmetric matrices. The random matrices are chosen from six ensembles, four of which lie in the Wigner class. For all three algorithms, we observe a form of universality for the deflation time statistics for random matrices within the Wigner class. For these ensembles, the empirical distribution of a normalized deflation time is found to collapse onto a curve that depends only on the algorithm, but not on the matrix size or deflation tolerance provided the matrix size is large enough (see Figure 4, Figure 7 and Figure 10). For the QR algorithm with the Wilkinson shift, the observed universality is even stronger and includes certain non-Wigner ensembles. Our experiments also provide a quantitative statistical picture of the accelerated convergence with shifts.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical solution of Dyson Brownian motion and a sampling scheme for invariant matrix ensembles

The Dyson Brownian Motion (DBM) describes the stochastic evolution of $N$ points on the line driven by an applied potential, a Coulombic repulsion and identical, independent Brownian forcing at each point. We use an explicit tamed Euler scheme to numerically solve the Dyson Brownian motion and sample the equilibrium measure for non-quadratic potentials. The Coulomb repulsion is too singular for the SDE to satisfy the hypotheses of rigorous convergence proofs for tamed Euler schemes. Nevertheless, in practice the scheme is observed to be stable for time steps of O(1/N^2) and to relax exponentially fast to the equilibrium measure with a rate constant of O(1) independent of $N$. Further, this convergence rate appears to improve with $N$ in accordance with O(1/N) relaxation of local statistics of the Dyson Brownian motion. This allows us to use the Dyson Brownian motion to sample $N\times N$ Hermitian matrices from the invariant ensembles. The computational cost of generating $M$ independent samples is O(MN^4) with a naive scheme, and O(MN^3\log N) when a fast multipole method is used to evaluate the Coulomb interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Complete integrability of shock clustering and Burgers turbulence

We consider scalar conservation laws with convex flux and random initial data. The Hopf-Lax formula induces a deterministic evolution of the law of the initial data. In a recent article, we derived a kinetic theory and Lax equations to describe the evolution of the law under the assumption that the initial data is a spectrally negative Markov process. Here we show that: (i) the Lax equations are Hamiltonian and describe a principle of least action on the Markov group that is in analogy with geodesic flow on $SO(N)$; (ii) the Lax equations are completely integrable and linearized via a loop-group factorization of operators; (iii) the associated zero-curvature equations can be solved via inverse scattering. Our results are rigorous for $N$-dimensional approximations of the Lax equations, and yield formulas for the limit $N \to \infty$. The main observation is that the Lax equations are a $N \to \infty$ limit of a Markovian variant of the $N$-wave model. This allows us to introduce a variety of methods from the theory of integrable systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Jet Formation in the magnetospheres of supermassive black holes: analytic solutions describing energy loss through Blandford-Znajek processes

In this paper, we provide exact solutions for the extraction of energy from a rotating black hole via both the electromagnetic Poynting flux and matter currents. By appropriate choice of a radially independent poloidal function $Λ(θ)$, we find solutions where the dominant outward energy flux is along the polar axis, consistent with a jet-like collimated outflow, but also with a weaker flux of energy along the equatorial plane. Unlike all the previously obtained solutions (Blandford & Znajek (1977), Menon & Dermer (2005)), the magnetosphere is free of magnetic monopoles everywhere.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

High Energy Radiation from Black Holes: A Summary

Bright gamma-ray flares observed from sources far beyond our Galaxy are best explained if enormous amounts of energy are liberated by black holes. The highest-energy particles in nature--the ultra-high energy cosmic rays--cannot be confined by the Milky Way's magnetic field, and must originate from sources outside our Galaxy. Here we summarize the themes of our book, "High Energy Radiation from Black Holes: Gamma Rays, Cosmic Rays, and Neutrinos", just published by Princeton University Press. In this book, we develop a mathematical framework that can be used to help establish the nature of gamma-ray sources, to evaluate evidence for cosmic-ray acceleration in blazars, GRBs and microquasars, to decide whether black holes accelerate the ultra-high energy cosmic rays, and to determine whether the Blandford-Znajek mechanism for energy extraction from rotating black holes can explain the differences between gamma-ray blazars and radio-quiet AGNs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Kinetic theory and Lax equations for shock clustering and Burgers turbulence

We study shock statistics in the scalar conservation law $\partial_t u + \partial_x f(u)=0$, $x \in \R$, $t>0$, with a convex flux $f$ and spatially random initial data. We show that the Markov property (in $x$) is preserved for a large class of random initial data (Markov processes with downward jumps and derivatives of Lévy processes with downward jumps). The kinetics of shock clustering is then described completely by an evolution equation for the generator of the Markov process $u(x,t)$, $x\in \R$. We present four distinct derivations for this evolution equation, and show that it takes the form of a Lax pair. The Lax equation admits a spectral parameter as in Manakov (1976), and has remarkable exact solutions for Burgers equation ($f(u)=u^2/2$). This suggests the kinetic equations of shock clustering are completely integrable.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Dynamics and self-similarity in min-driven clustering

We study a mean-field model for a clustering process that may be described informally as follows. At each step a random integer $k$ is chosen with probability $p_k$, and the smallest cluster merges with $k$ randomly chosen clusters. We prove that the model determines a continuous dynamical system on the space of probability measures supported in $(0,\infty)$, and we establish necessary and sufficient conditions for approach to self-similar form. We also characterize eternal solutions for this model via a Levy-Khintchine formula. The analysis is based on an explicit solution formula discovered by Gallay and Mielke, extended using a careful choice of time scale.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Timelike Geodesic Currents in the Stationary, Axisymmetric, Force-free Magnetosphere of a Kerr Black Hole

The structural properties of geodesic currents in an ambient Kerr background is studied from an analytical point of view. The geodesics in the congruence correspond to charged particles that carry energy and angular momentum from the black hole through the Blandford-Znajek mechanism. It is shown that the resulting magnetosphere naturally satisfies the Znajek regularity condition. Particular attention is paid here to the energy extracted by matter currents rather than by electromagnetic Poynting fluxes.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The scaling attractor and ultimate dynamics for Smoluchowski's coagulation equations

We describe a basic framework for studying dynamic scaling that has roots in dynamical systems and probability theory. Within this framework, we study Smoluchowski's coagulation equation for the three simplest rate kernels $K(x,y)=2$, $x+y$ and $xy$. In another work, we classified all self-similar solutions and all universality classes (domains of attraction) for scaling limits under weak convergence (Comm. Pure Appl. Math 57 (2004)1197-1232). Here we add to this a complete description of the set of all limit points of solutions modulo scaling (the scaling attractor) and the dynamics on this limit set (the ultimate dynamics). The main tool is Bertoin's Lévy-Khintchine representation formula for eternal solutions of Smoluchowski's equation (Adv. Appl. Prob. 12 (2002) 547--64). This representation linearizes the dynamics on the scaling attractor, revealing these dynamics to be conjugate to a continuous dilation, and chaotic in a classical sense. Furthermore, our study of scaling limits explains how Smoluchowski dynamics ``compactifies'' in a natural way that accounts for clusters of zero and infinite size (dust and gel).

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Universality classes in Burgers turbulence

We establish necessary and sufficient conditions for the shock statistics to approach self-similar form in Burgers turbulence with Lévy process initial data. The proof relies upon an elegant closure theorem of Bertoin and Carraro and Duchon that reduces the study of shock statistics to Smoluchowski's coagulation equation with additive kernel, and upon our previous characterization of the domains of attraction of self-similar solutions for this equation.

0 citations0 reviewsNo code linkedOpen access