About

Amitava Bhattacharjee is a Professor of Astrophysical Sciences at Princeton University. His office is located at 112 Nassau St., Princeton, NJ. His research focuses on astrophysical sciences, and he is a member of the Department of Astrophysical Sciences at Princeton. Further details about his specific research contributions are not provided in the page text.

Research topics

• Physics
• Mechanics
• Classical mechanics
• Computational physics
• Atomic physics

Selected publications

• Non-thermal particle acceleration in multi-species kinetic plasmas: universal power-law distribution functions and temperature inversion in the solar corona

  arXiv (Cornell University) · 2026-01-06

  preprintOpen accessSenior author

  Non-thermal power-law distribution functions are ubiquitous in astrophysical, space, and laboratory kinetic plasmas, but their origin remains unclear. A related puzzle is the temperature inversion of the solar corona. We show that these phenomena are deeply connected by developing a self-consistent quasilinear theory for electromagnetically driven, unmagnetized kinetic plasmas. The theory yields a multi-species Fokker-Planck equation with drive-induced diffusion from direct acceleration by broad-band turbulent or narrow-band wave-like fields, indirect acceleration by excited waves, and Balescu-Lenard diffusion/drag from Debye-scale fluctuations and Coulomb collisions. For a super-Debye turbulent electric-field spectrum, $|{\bf E}_{\bf k}|^2\propto k^{-α}$, electrons and ions relax toward a universal $f(v)\propto v^{-5}$, or $N(E)\propto E^{-2}$, attractor, equivalent to the high-energy tail of a $κ=1.5$ distribution, when $α\ge5$. This universality follows from Debye screening: large-scale fields accelerate unscreened fast particles but not screened slow ones. For shallower spectra, $α<5$, the tail scales as $v^{-α}$; incomplete relaxation and anisotropy also break universality. Anisotropic wave drives yield branch- and spectrum-dependent exponents. Because collisions cannot decelerate suprathermal particles, the tails resist Maxwellianization. In the solar atmosphere, such tails may be generated by chromospheric convection or nanoflares despite collisional and radiative losses. Direct wave heating energizes electrons through Landau-resonant interactions with whistler and electron-cyclotron waves, while ions may be accelerated by turbulent ambipolar fields. Resulting $κ\simeq1.5$--$3$ distributions produce an abrupt upper-chromosphere/lower-corona transition and velocity-filtration-driven inverted profiles, yielding coronal temperatures $\sim10^6\,{\rm K}$.

  PublisherDOI

• The Geometry of Flux Surfaces with Quasi-Poloidal Symmetry

  arXiv (Cornell University) · 2026-01-20

  preprintOpen accessSenior author

  Quasi-poloidal (QP) magnetic fields have desirable properties for confining plasma: no radial drift of guiding centres (with positive implications for neoclassical transport), zero Pfirsch-Schlüter current, a lower level of damping for poloidal flows, leading to reduced anomalous transport, and possible stability benefits. Despite their attractive properties, QP fields are not amenable to the near-axis expansion, a major theoretical tool for understanding toroidal fields. In this paper, we provide a novel framework for defining and understanding QP flux surfaces. This framework relies on a simplification that transforms the task of finding a quasi-poloidal flux surface from a 3D problem to a 2D problem. This simplification also applies to asymmetric magnetic mirrors with desirable properties. We sketch how this 2D problem can form the basis of an efficient optimisation problem for finding QP flux surfaces. We leverage this 2D problem for theoretical understanding: for instance, we identify one class of QP flux surfaces that are naturally flat mirrors (Velasco et al. 2023). The reduced model is validated against numerically optimised QP equilibria. We further utilise the reduced model to explain the prevalence of cusps, high mirror ratios, and narrow pinch points in these numerical equilibria.

  PublisherDOI

• The Geometry of Flux Surfaces with Quasi-Poloidal Symmetry

  ArXiv.org · 2026-01-20

  articleOpen accessSenior author

  Quasi-poloidal (QP) magnetic fields have desirable properties for confining plasma: no radial drift of guiding centres (with positive implications for neoclassical transport), zero Pfirsch-Schlüter current, a lower level of damping for poloidal flows, leading to reduced anomalous transport, and possible stability benefits. Despite their attractive properties, QP fields are not amenable to the near-axis expansion, a major theoretical tool for understanding toroidal fields. In this paper, we provide a novel framework for defining and understanding QP flux surfaces. This framework relies on a simplification that transforms the task of finding a quasi-poloidal flux surface from a 3D problem to a 2D problem. This simplification also applies to asymmetric magnetic mirrors with desirable properties. We sketch how this 2D problem can form the basis of an efficient optimisation problem for finding QP flux surfaces. We leverage this 2D problem for theoretical understanding: for instance, we identify one class of QP flux surfaces that are naturally flat mirrors (Velasco et al. 2023). The reduced model is validated against numerically optimised QP equilibria. We further utilise the reduced model to explain the prevalence of cusps, high mirror ratios, and narrow pinch points in these numerical equilibria.

  PublisherOA PDF

• Non-thermal particle acceleration in multi-species kinetic plasmas: universal power-law distribution functions and temperature inversion in the solar corona

  ArXiv.org · 2026-01-06

  articleOpen accessSenior author

  Non-thermal power-law distribution functions are ubiquitous in astrophysical, space, and laboratory kinetic plasmas, but their origin remains unclear. A related puzzle is the temperature inversion of the solar corona. We show that these phenomena are deeply connected by developing a self-consistent quasilinear theory for electromagnetically driven, unmagnetized kinetic plasmas. The theory yields a multi-species Fokker-Planck equation with drive-induced diffusion from direct acceleration by broad-band turbulent or narrow-band wave-like fields, indirect acceleration by excited waves, and Balescu-Lenard diffusion/drag from Debye-scale fluctuations and Coulomb collisions. For a super-Debye turbulent electric-field spectrum, $|{\bf E}_{\bf k}|^2\propto k^{-α}$, electrons and ions relax toward a universal $f(v)\propto v^{-5}$, or $N(E)\propto E^{-2}$, attractor, equivalent to the high-energy tail of a $κ=1.5$ distribution, when $α\ge5$. This universality follows from Debye screening: large-scale fields accelerate unscreened fast particles but not screened slow ones. For shallower spectra, $α<5$, the tail scales as $v^{-α}$; incomplete relaxation and anisotropy also break universality. Anisotropic wave drives yield branch- and spectrum-dependent exponents. Because collisions cannot decelerate suprathermal particles, the tails resist Maxwellianization. In the solar atmosphere, such tails may be generated by chromospheric convection or nanoflares despite collisional and radiative losses. Direct wave heating energizes electrons through Landau-resonant interactions with whistler and electron-cyclotron waves, while ions may be accelerated by turbulent ambipolar fields. Resulting $κ\simeq1.5$--$3$ distributions produce an abrupt upper-chromosphere/lower-corona transition and velocity-filtration-driven inverted profiles, yielding coronal temperatures $\sim10^6\,{\rm K}$.

  PublisherOA PDF

• Scale-dependent alignment in compressible magnetohydrodynamic turbulence

  ArXiv.org · 2025-04-22

  preprintOpen accessSenior author

  Using $10,\!080^3$ grid simulations, we analyze scale-dependent alignment in driven, compressible, no net-flux magnetohydrodynamic turbulence. The plasma self-organizes into localized, strongly aligned regions. Alignment spans all primitive variables and their curls. Contrary to incompressible theory, velocity-magnetic alignment scales as $θ(λ) \sim λ^{1/8}$, where $λ$ is the scale, suggesting a distinct three-dimensional eddy anisotropy and a much higher critical transition scale toward a reconnection-mediated cascade.

  PublisherOA PDFDOI

• Computation of magnetohydrodynamic equilibria with Voigt regularization

  Physics of Plasmas · 2025-06-01 · 3 citations

  articleOpen accessSenior author

  This work presents the first numerical investigation of using Voigt regularization as a method for obtaining magnetohydrodynamic (MHD) equilibria without the assumption of nested magnetic flux surfaces. Voigt regularization modifies the MHD dynamics by introducing additional terms that vanish in the infinite-time limit, allowing for magnetic reconnection and the formation of magnetic islands, which can overlap and produce field-line chaos. The utility of this approach is demonstrated through numerical solutions of two-dimensional ideal and resistive test problems. Our results show that Voigt regularization can significantly accelerate the convergence to solutions in resistive MHD problems while also highlighting challenges in applying the method to ideal MHD systems. This research opens up new possibilities for developing more efficient and robust MHD equilibrium solvers, which could contribute to the design and optimization of future fusion devices.

  PublisherDOI

• Global stellarator coil optimization with quadratic constraints and objectives

  Nuclear Fusion · 2025-01-09 · 3 citations

  articleOpen accessSenior author

  Abstract Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mrow> <mml:msup> <mml:mn>10</mml:mn> <mml:mn>2</mml:mn> </mml:msup> <mml:mo>×</mml:mo> </mml:mrow> </mml:math> faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples.

  PublisherDOI

• Periodic Korteweg-de Vries soliton potentials generate quasisymmetric magnetic field strength in a finite plasma-beta equilibrium

  ArXiv.org · 2025-07-09

  preprintOpen accessSenior author

  Quasisymmetry (QS) is a hidden symmetry of the magnetic field strength, B, that enables effective confinement of charged particles in a fully three-dimensional (3D) toroidal plasma equilibrium. Such equilibria are typically modeled by the ideal magnetohydrostatic (MHS) equation. The nonlinear, overdetermined nature of the quasisymmetric MHS equations severely complicates our understanding of the interplay between 3D shaping, equilibrium properties such as pressure and rotational transform, and B. Progress has been made through expansions near the magnetic axis; however, a more comprehensive theory is desirable. Using a combination of analysis and regression on a large dataset of numerically optimized quasisymmetric stellarators, we demonstrate that there is a hidden lower dimensionality of B on a magnetic flux surface with connections to the theory of periodic solitons. We show that $B$ on a flux surface is determined by three or at most four flux functions, each of which is a critical value of the derivative of B along the field line. While being consistent with the near-axis models, our results are global and hold even on the last closed flux surface.

  PublisherOA PDFDOI

• Density fluctuation-Mach number scaling in compressible, high plasma beta turbulence: in-situ space observations and high-Reynolds number simulations

  ArXiv.org · 2025-02-13

  preprintOpen accessSenior author

  The file "mms_int-beta-Mt-dror-r_high-beta.dat" contains the MMS intervals used. The file contains columns of start time, end time, ion beta, turbulent mach number, normalized density fluctuation, and density (in units of cm^-3). The file "mass_density_velocity_disp_10k.csv" contains the median, 16th percentile and 84th percentile of the turbulent mach number and normalized density fluctuation values obtained from the MHD simulation.

  PublisherOA PDFDOI

• The spectrum of magnetized turbulence in the interstellar medium

  Nature Astronomy · 2025-05-13 · 10 citations

  articleOpen accessSenior author
  PublisherOA PDFDOI

Recent grants

• Collaborative Research: Frameworks: A Software Ecosystem for Plasma Science and Space Weather Applications

  NSF · $2.0M · 2022–2027

• SHINE: Instability of Large-Scale Current Sheets, Plasmoid Formation, and Particle Acceleration in Coronal Plasmas

  NSF · $158k · 2013–2015

• Stability of Thin Current Sheets in the Earth's Magnetotail: Theory, Simulations, and Observations

  NSF · $330k · 2004–2008

• Collaborative Research: Integration of Extended Magnetohydrodynamic (MHD) and Kinetic Effects in Global Magnetosphere Models

  NSF · $1.3M · 2013–2019

• GEM: Can Global Magnetohydrodynamics (MHD) Codes Model Magnetic Reconnection in the High Lundquist Number Limit?

  NSF · $300k · 2008–2011

Frequent coauthors

• Yi-Min Huang

  Princeton University

  254 shared

• K. Germaschewski

  University of New Hampshire

  228 shared

• W. Fox

  Princeton University

  142 shared

• Ammar Hakim

  Bandung Institute of Technology

  124 shared

• C. S. Ng

  106 shared

• Jonathan Ng

  University of Maryland, College Park

  101 shared

• N. Bessho

  University of Maryland, College Park

  98 shared

• W. Daughton

  90 shared

Labs

• Amitava Bhattacharjee LaboratoryPI