Dust Acoustic Solitons Tracked Analytically in Inhomogeneous Nonthermal Plasmas

A new theoretical study has worked out an analytical solution for dust acoustic waves travelling through an inhomogeneous dusty plasma whose electrons and ions deviate from thermal equilibrium. Published in the journal Discover Space, the analysis by Himangshu Kalita, Hirak Jyoti Dehingia and Karabi Devi derives a damped Korteweg–de Vries equation governing the nonlinear evolution of the dust acoustic mode and then solves it exactly, producing explicit expressions for the amplitude and width of the resulting solitary structures. The work is framed as the first nonlinear theoretical investigation to combine equilibrium inhomogeneity with Tribeche-type nonthermal electron and ion distributions in the study of dust acoustic solitary wave propagation.

Dusty plasmas are mixtures of conventional plasma constituents—electrons, ions and neutral particles—together with nanometer- to micrometer-sized dust grains. Such systems have been observed in interstellar clouds, interplanetary space, cometary tails, the mesosphere, planetary ring systems and laboratory devices, and they have attracted sustained attention from the plasma physics community over recent decades. Because electrons are far more mobile than ions, dust grains immersed in a plasma generally collect electrons and acquire a negative charge, and this charging modifies the plasma wave spectrum in significant ways. Depending on how mobile the grains are, different wave modes emerge: when dust is treated as immobile while ions move, a dust–ion acoustic wave appears, a modified version of the ordinary ion acoustic wave; when the grains themselves are mobile, a distinct low-frequency mode called the dust acoustic wave arises.

The experimental pedigree of the dust acoustic wave matters for the credibility of any new theoretical treatment. Its existence was confirmed experimentally by Barkan and colleagues in 1995 and 1996 and later by Hayashi in 2001, and a large body of theoretical and experimental work has followed, much of it assuming spatially uniform plasmas. The authors of the new study point out, however, that both astrophysical and laboratory plasmas are often inherently inhomogeneous, with inhomogeneities arising from spatial variations in density, temperature or magnetic field strength. Earlier work has shown that such gradients matter: Singh and Rao demonstrated that dust density inhomogeneity can significantly affect the amplitude of dust acoustic waves, while other groups examined solitary waves in strongly coupled inhomogeneous dusty plasmas and in four-component plasmas containing oppositely charged dust species.

A further limitation of much of the earlier literature, according to the authors, is that investigations of equilibrium inhomogeneity generally assumed electrons and ions obeying Maxwellian or Boltzmann velocity distributions. Yet observations in space and laboratory plasmas indicate that particle populations often deviate substantially from thermal equilibrium, exhibiting nonthermal features such as energetic tails and trapped particle populations. Some recent studies have incorporated nonthermal effects in dusty plasma systems, but the combined influence of equilibrium inhomogeneity and specific nonthermal distributions had, in the authors’ assessment, not been adequately explored. The new model therefore adopts the Tribeche distribution for both electrons and ions, which allows controlled departures from Maxwellian behaviour through two nonthermal parameters, denoted beta-e for electrons and beta-i for ions. In the limit where both parameters tend to zero, the distribution functions reduce to the usual Boltzmann forms, and the parameters are restricted to values that keep the particle densities positive-definite.

The physical model is an unmagnetized, collisionless, one-dimensional dusty plasma composed of negatively charged dust grains, nonthermal electrons and nonthermal ions. The massive dust component is described by fluid equations—a continuity equation and a momentum equation—with the dust treated as cold, so that dust thermal pressure is neglected. The electrostatic potential is governed by Poisson’s equation, and overall charge neutrality at equilibrium ties the equilibrium ion density to the equilibrium electron density and the product of the dust charge number and equilibrium dust density. The equations are normalized using a fully specified electrostatic dusty-plasma convention: lengths are scaled by the electron Debye length, time by the inverse electron plasma frequency, the potential by the electron thermal energy, densities by the equilibrium plasma density, and the dust fluid velocity by the characteristic dust acoustic speed derived from balancing dust inertia against the electrostatic force. Dust charging dynamics and higher-order charge-to-mass corrections are neglected at leading order, on the grounds that their influence is small compared with the dominant electrostatic and inertial terms.

The mathematical core of the paper is the reductive perturbation technique. The authors introduce stretched coordinates containing a small expansion parameter that represents the strength of nonlinearity and dispersion, together with a normalized Mach speed for the dust acoustic wave. The dependent variables—electron, ion and dust densities, dust fluid velocity and electrostatic potential—are expanded in powers of this parameter, with the subscript zero denoting the equilibrium state and successive terms representing first- and second-order perturbations due to weak nonlinearity and inhomogeneity. Substituting the expansions into the governing equations, neglecting powers of the expansion parameter beyond the relevant order and matching coefficients of like powers yields, at first order, the phase velocity of the dust acoustic solitary waves, which depends on the equilibrium densities, the nonthermal parameters, the electron-to-ion temperature ratio and the dust charge number.

At the higher orders of the expansion, the procedure produces a damped Korteweg–de Vries equation for the first-order potential. The damping is not due to collisions—the plasma is collisionless—but arises from an additional term generated by the background plasma density gradient. The equation carries three coefficients: a nonlinear coefficient, a dispersion coefficient and the damping coefficient associated with the density gradient. All three depend on the equilibrium plasma parameters, and the nonlinear and dispersion coefficients depend on both the density inhomogeneity and the dust charge. The authors emphasize that this damped KdV equation describes the nonlinear propagation of dust acoustic solitary waves in a non-uniform dusty plasma containing negatively charged dust grains, and that in several limiting cases their results reduce, in a piecewise manner, to previously reported special and simplified models, showing agreement with earlier studies.

To solve the equation analytically, the authors apply a nonlinear transformation that absorbs the damping term, reducing the problem to a standard KdV form whose coefficients are functions of the equilibrium plasma parameters. Assuming that the coefficient variations are small and moving to a wave frame, they obtain a sech-squared solitary wave solution, the classic soliton profile, with explicit formulas for the amplitude and width in terms of the transformed nonlinear and dispersive coefficients and the phase velocity in the stretched frame. The solution is valid under the assumption of weak nonlinearity and dispersion, and the authors state that it captures the essential features of localized electrostatic structures in nonthermal dusty plasmas.

The parametric study that follows yields a set of specific trends. The phase velocity increases monotonically with both the electron and ion nonthermal parameters, which the authors attribute to stronger effective pressure response and restoring force; the three-dimensional surface plot shows a steeper gradient at higher ion nonthermality, suggesting ion nonthermality has the comparatively stronger impact in the considered regime. Increasing the equilibrium electron or ion density slightly reduces the phase velocity, attributed to enhanced shielding of the dust charge. The phase velocity rises monotonically with equilibrium dust density for all values of the temperature ratio and increases almost linearly with the dust charge number, while higher electron temperature relative to ion temperature lowers the wave speed. The soliton amplitude behaves non-monotonically with respect to the nonthermal parameters: it initially increases, reaches a maximum at an intermediate value of each parameter and then decreases, reflecting a competition between nonlinear steepening, enhanced by nonthermal populations, and dispersive effects that become dominant at larger nonthermality. The amplitude increases with equilibrium electron density but decreases rapidly with increasing dust density and dust charge number, which the authors attribute to enhanced dust inertia and charge weakening nonlinear steepening. The soliton width decreases monotonically with both nonthermal parameters, with increasing dust density—sharply so near zero dust density—and with increasing dust charge number, the last showing an asymptotic decay at high charge. Weak equilibrium inhomogeneity enters through spatial variations in the equilibrium densities, modulating the amplitude and width gradually as the wave propagates, although the principal trends remain governed by electron and ion nonthermality.

The authors argue that the theoretical predictions can be examined under controlled laboratory conditions, since dusty plasma experiments allow the generation of plasmas containing negatively charged dust grains with adjustable spatial and dynamical properties. They outline possible strategies: employing Q-machines or RF discharge devices to produce well-defined dusty plasma environments, varying density, temperature or dust concentration to probe their influence on wave behaviour, and using high-speed imaging diagnostics to monitor the formation and propagation of solitary structures for comparison with theory, noting that previous experimental investigations have demonstrated such measurements are feasible. They also sketch future extensions, including variable dust charge, external magnetic fields, collisional dissipation, higher-dimensional geometries, ion viscosity, relativistic or quantum corrections and more generalized nonthermal distributions, as well as applications to space environments such as interstellar molecular clouds, cometary tails and planetary ring systems. For now, the study remains a theoretical contribution: no datasets were generated or analysed, and the analytical solution, with its predicted optimal nonthermal regime for maximum soliton amplitude and its saturating width behaviour at high dust charge, awaits direct experimental test against contemporary observations in complex plasmas.

Subject of Research: Analytical solutions of dust acoustic solitary waves governed by a damped KdV equation in inhomogeneous dusty plasmas with nonthermal electrons and ions

Article Title: Analytical solution of dust acoustic waves governed by the damped KdV equation in inhomogeneous nonthermal plasmas

Article References: Kalita, H., Dehingia, H. J., & Devi, K. (2026). Analytical solution of dust acoustic waves governed by the damped KdV equation in inhomogeneous nonthermal plasmas. Discover Space, 130(1), Article 12. https://doi.org/10.1007/s11038-026-09591-5

Image Credits: AI Generated

DOI: 10.1007/s11038-026-09591-5

Keywords: dusty plasma, dust acoustic waves, solitons, damped KdV equation, nonthermal distributions, plasma inhomogeneity, reductive perturbation method, dust charging, phase velocity, space plasmas, laboratory plasmas, theoretical physics