Harish Murali
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Harish Murali has a total of 7 co-authors. Here are the top co-authors by number of shared publications:
- Pedro Vieira
Perimeter Institute for Theoretical Physics
3 shared publications - Andrea L. Guerrieri
City St. George’s, Univ. of London
2 shared publications - Martin Kruczenski
Department of Physics and Astronomy and PQSEI † Purdue University 47907 West Lafayette IN USA
2 shared publications - Vladimir Kazakov
Laboratoire de Physique de l’École Normale Supérieure, CNRS, Université PSL, Sorbonne Universités
2 shared publications - Adolfo Holguin
Trinity College Dublin
1 shared publication - João Penedones
Institute of Physics, Ecole Polytechnique Fédérale de Lausanne (EPFL)
1 shared publication
Latest Publications
- Open Access
- 9 April 2026
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2026 (4) , 69
(Un)solvable matrix models for bps correlators
A bstract: We propose and study a family of complex matrix models computing the protected two- and three-point correlation functions in $$ \mathcal{N}=4 $$ SYM. Our description allows us to directly relate the eigenvalue density of the matrix model for “Huge” operators with ∆ ~ N 2 to the shape of droplets in the dual Lin-Lunin-Maldacena (LLM) geometry. We demonstrate how to determine the eigenvalue distribution for various choices of operators such as those of exponential, character, or coherent state type, which then allows us to efficiently compute one-point functions of light chiral primaries in generic LLM backgrounds. In particular, we successfully match the results for light probes with the supergravity calculations of Skenderis and Taylor. We provide a large N formalism for one-point functions of “Giant” probes, such as operators dual to giant graviton branes in LLM backgrounds, and explicitly apply it for particular backgrounds. We also explicitly compute the correlator of three huge half-BPS operators of exponential type and stacks of determinant operators by reducing them to the known matrix model problems such as the Potts or O ( n ) model on random planar graphs. Finally, we point out a curious relation between the correlators of $$ \frac{1}{4} $$ -BPS and $$ \frac{1}{8} $$ -BPS coherent state operators and the Eguchi-Kawai reduction of the Principal Chiral Model in 2 D and 3 D correspondingly.- 9 April 2026
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2026 (4) , 69
Open Access(Un)solvable matrix models for bps correlators
A bstract: We propose and study a family of complex matrix models computing the protected two- and three-point correlation functions in $$ \mathcal{N}=4 $$ SYM. Our description allows us to directly relate the eigenvalue density of the matrix model for “Huge” operators with ∆ ~ N 2 to the shape of droplets in the dual Lin-Lunin-Maldacena (LLM) geometry. We demonstrate how to determine the eigenvalue distribution for various choices of operators such as those of exponential, character, or coherent state type, which then allows us to efficiently compute one-point functions of light chiral primaries in generic LLM backgrounds. In particular, we successfully match the results for light probes with the supergravity calculations of Skenderis and Taylor. We provide a large N formalism for one-point functions of “Giant” probes, such as operators dual to giant graviton branes in LLM backgrounds, and explicitly apply it for particular backgrounds. We also explicitly compute the correlator of three huge half-BPS operators of exponential type and stacks of determinant operators by reducing them to the known matrix model problems such as the Potts or O ( n ) model on random planar graphs. Finally, we point out a curious relation between the correlators of $$ \frac{1}{4} $$ -BPS and $$ \frac{1}{8} $$ -BPS coherent state operators and the Eguchi-Kawai reduction of the Principal Chiral Model in 2 D and 3 D correspondingly. - Open Access
- 5 February 2026
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2026 (2) , 74
Universality of heavy operators in matrix models
A bstract: In large N theories with a gravity dual, generic heavy operators should be dual to black holes in the bulk. The microscopic details of such operators should then be irrelevant in the low energy theory. We look for such universality in the strong coupling limit of a very simple two matrix model — the Hoppe model. Using analytics as well as Monte Carlo simulations, we show that there exists a universal black hole regime where the eigenvalue densities are given by parabolas and the correlation functions of probes in these backgrounds are completely determined by a few parameters. An important feature of strong coupling in this model is that the matrices commute and one can define joint eigenvalue distributions which also exhibit universality. These two results extend the beautiful findings of Berenstein, Hanada and Hartnoll [1]. Not all heavy operators are universal and at strong coupling there is a sharp phase boundary between the universal and non universal regimes (Of course this should not be confused with the universality of eigenvalue spacing in matrix models). Moreover, in the non universal phase, we also find an interesting phenomenon we call Abelianization where some eigenvalues run off to infinity, reminiscent of heavy dual giant gravitons in $$ \mathcal{N}=4 $$ SYM.- 5 February 2026
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2026 (2) , 74
Open AccessUniversality of heavy operators in matrix models
A bstract: In large N theories with a gravity dual, generic heavy operators should be dual to black holes in the bulk. The microscopic details of such operators should then be irrelevant in the low energy theory. We look for such universality in the strong coupling limit of a very simple two matrix model — the Hoppe model. Using analytics as well as Monte Carlo simulations, we show that there exists a universal black hole regime where the eigenvalue densities are given by parabolas and the correlation functions of probes in these backgrounds are completely determined by a few parameters. An important feature of strong coupling in this model is that the matrices commute and one can define joint eigenvalue distributions which also exhibit universality. These two results extend the beautiful findings of Berenstein, Hanada and Hartnoll [1]. Not all heavy operators are universal and at strong coupling there is a sharp phase boundary between the universal and non universal regimes (Of course this should not be confused with the universality of eigenvalue spacing in matrix models). Moreover, in the non universal phase, we also find an interesting phenomenon we call Abelianization where some eigenvalues run off to infinity, reminiscent of heavy dual giant gravitons in $$ \mathcal{N}=4 $$ SYM. - Open Access
- 17 September 2025
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2025 (9) , 1-63
Huge BPS operators and fluid dynamics in $ \mathcal{N}=4 $ SYMIn the bulk dual of holography, huge operators correspond to sources so heavy that they fully backreact on the space-time geometry. Here we study the correlation function of three such huge operators when they are given by 1/2 BPS operators in $$ \mathcal{N}=4 $$ SYM theory, dual to IIB Strings in AdS 5 × S 5. We unveil simple matrix model representations for these correlators which we can sometimes solve analytically. For general huge operators, we translate these matrix model expressions into a 1 + 1 dimensional hydrodynamical fluid problem. This fluid is integrable thus unveiling a novel integrable sector of the AdS/CFT duality in a full fledged gravitational regime, very far from the usual free string planar regime where integrability reigns supreme. We explain how an adiabatic deformation method can be developed to yield the solution to an integrable discrete formulation of these fluids — the rational Calogero-Moser Model — so we can access the general three point correlation functions of generic huge 1/2-BPS operators. Everything will be done on the gauge theory side of the duality. It would be fascinating to find the holographic dual of these matrix models and fluids.- 17 September 2025
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2025 (9) , 1-63
Open AccessHuge BPS operators and fluid dynamics in $ \mathcal{N}=4 $ SYMIn the bulk dual of holography, huge operators correspond to sources so heavy that they fully backreact on the space-time geometry. Here we study the correlation function of three such huge operators when they are given by 1/2 BPS operators in $$ \mathcal{N}=4 $$ SYM theory, dual to IIB Strings in AdS 5 × S 5. We unveil simple matrix model representations for these correlators which we can sometimes solve analytically. For general huge operators, we translate these matrix model expressions into a 1 + 1 dimensional hydrodynamical fluid problem. This fluid is integrable thus unveiling a novel integrable sector of the AdS/CFT duality in a full fledged gravitational regime, very far from the usual free string planar regime where integrability reigns supreme. We explain how an adiabatic deformation method can be developed to yield the solution to an integrable discrete formulation of these fluids — the rational Calogero-Moser Model — so we can access the general three point correlation functions of generic huge 1/2-BPS operators. Everything will be done on the gauge theory side of the duality. It would be fascinating to find the holographic dual of these matrix models and fluids. - Open Access
- 13 June 2023
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2023 (6) , 1-32
Where is M-theory in the space of scattering amplitudes?
We use the S-matrix bootstrap to carve out the space of unitary, analytic, crossing symmetric and supersymmetric graviton scattering amplitudes in nine, ten and eleven dimensions. We extend and improve the numerical methods of our previous work in ten dimensions. A key new tool employed here is unitarity in the celestial sphere. In all dimensions, we find that the minimal allowed value of the Wilson coefficient α, controlling the leading correction to maximal supergravity, is very close but not equal to the minimal value realized in Superstring theory or M-theory. This small difference may be related to inelastic effects that are not well described by our numerical extremal amplitudes. Although α has a unique value in M-theory, we found no evidence of an upper bound on α in 11D.
- 13 June 2023
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2023 (6) , 1-32
Open AccessWhere is M-theory in the space of scattering amplitudes?
We use the S-matrix bootstrap to carve out the space of unitary, analytic, crossing symmetric and supersymmetric graviton scattering amplitudes in nine, ten and eleven dimensions. We extend and improve the numerical methods of our previous work in ten dimensions. A key new tool employed here is unitarity in the celestial sphere. In all dimensions, we find that the minimal allowed value of the Wilson coefficient α, controlling the leading correction to maximal supergravity, is very close but not equal to the minimal value realized in Superstring theory or M-theory. This small difference may be related to inelastic effects that are not well described by our numerical extremal amplitudes. Although α has a unique value in M-theory, we found no evidence of an upper bound on α in 11D.
- Open Access
- 12 April 2021
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2021 (4) , 1-27
The R-matrix bootstrap for the 2d O(N) bosonic model with a boundary
The S-matrix bootstrap is extended to a 1+1d theory with O(N) symmetry and a boundary in what we call the R-matrix bootstrap since the quantity of interest is the reflection matrix (R-matrix). Given a bulk S-matrix, the space of allowed R-matrices is an infinite dimensional convex space from which we plot a two dimensional section given by a convex domain on a 2d plane. In certain cases, at the boundary of the domain, we find vertices corresponding to integrable R-matrices with no free parameters. In other cases, when there is a one-parameter family of integrable R-matrices, the whole boundary represents integrable theories. We also consider R-matrices which are analytic in an extended region beyond the physical cuts, thus forbidding poles (resonances) in that region. In certain models, this drastically reduces the allowed space of R-matrices leading to new vertices that again correspond to integrable theories. We also work out the dual problem, in particular in the case of extended analyticity, the dual function has cuts on the physical line whenever unitarity is saturated. For the periodic Yang-Baxter solution that has zero transmission, we computed the R-matrix initially using the bootstrap and then derived its previously unknown analytic form.
- 12 April 2021
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2021 (4) , 1-27
Open AccessThe R-matrix bootstrap for the 2d O(N) bosonic model with a boundary
The S-matrix bootstrap is extended to a 1+1d theory with O(N) symmetry and a boundary in what we call the R-matrix bootstrap since the quantity of interest is the reflection matrix (R-matrix). Given a bulk S-matrix, the space of allowed R-matrices is an infinite dimensional convex space from which we plot a two dimensional section given by a convex domain on a 2d plane. In certain cases, at the boundary of the domain, we find vertices corresponding to integrable R-matrices with no free parameters. In other cases, when there is a one-parameter family of integrable R-matrices, the whole boundary represents integrable theories. We also consider R-matrices which are analytic in an extended region beyond the physical cuts, thus forbidding poles (resonances) in that region. In certain models, this drastically reduces the allowed space of R-matrices leading to new vertices that again correspond to integrable theories. We also work out the dual problem, in particular in the case of extended analyticity, the dual function has cuts on the physical line whenever unitarity is saturated. For the periodic Yang-Baxter solution that has zero transmission, we computed the R-matrix initially using the bootstrap and then derived its previously unknown analytic form.
- Open Access
The R-matrix bootstrap for the 2d O(N) bosonic model with a boundary
The S-matrix bootstrap is extended to a 1+1d theory with $O(N)$ symmetry and a boundary in what we call the R-matrix bootstrap since the quantity of interest is the reflection matrix (R-matrix). Given a bulk S-matrix, the space of allowed R-matrices is an infinite dimensional convex space from which we plot a two dimensional section given by a convex domain on a 2d plane. In certain cases, at the boundary of the domain, we find vertices corresponding to integrable R-matrices with no free parameters. In other cases, when there is a one-parameter family of integrable R-matrices, the whole boundary represents integrable theories. We also consider R-matrices which are analytic in an extended region beyond the physical cuts, thus forbidding poles (resonances) in that region. In certain models, this drastically reduces the allowed space of R-matrices leading to new vertices that again correspond to integrable theories. We also work out the dual problem, in particular in the case of extended analyticity, the dual function has cuts on the physical line whenever unitarity is saturated. For the periodic Yang-Baxter solution that has zero transmission, we computed the R-matrix initially using the bootstrap and then derived its previously unknown analytic form.Open AccessThe R-matrix bootstrap for the 2d O(N) bosonic model with a boundary
The S-matrix bootstrap is extended to a 1+1d theory with $O(N)$ symmetry and a boundary in what we call the R-matrix bootstrap since the quantity of interest is the reflection matrix (R-matrix). Given a bulk S-matrix, the space of allowed R-matrices is an infinite dimensional convex space from which we plot a two dimensional section given by a convex domain on a 2d plane. In certain cases, at the boundary of the domain, we find vertices corresponding to integrable R-matrices with no free parameters. In other cases, when there is a one-parameter family of integrable R-matrices, the whole boundary represents integrable theories. We also consider R-matrices which are analytic in an extended region beyond the physical cuts, thus forbidding poles (resonances) in that region. In certain models, this drastically reduces the allowed space of R-matrices leading to new vertices that again correspond to integrable theories. We also work out the dual problem, in particular in the case of extended analyticity, the dual function has cuts on the physical line whenever unitarity is saturated. For the periodic Yang-Baxter solution that has zero transmission, we computed the R-matrix initially using the bootstrap and then derived its previously unknown analytic form.