Polyatomic Classical Models Overview

• Polyatomic classical models are statistical and kinetic frameworks that incorporate internal molecular degrees of freedom, such as rotation and vibration, to simulate gas and material behavior.
• They extend conventional Boltzmann and BGK models by introducing multisite interactions and moment hierarchies (e.g., ET15) to accurately capture non-equilibrium processes and energy transfer.
• These models enable efficient numerical schemes and simulations by addressing multi-timescale relaxation, guiding applications in computational fluid dynamics, material science, and statistical mechanics.

A polyatomic classical model refers to any classical statistical or kinetic model that explicitly incorporates the internal structure and degrees of freedom of polyatomic molecules—typically rotational, vibrational, or multiplicity of atomic sites—within the framework of classical physics. These models are foundational in statistical mechanics, rarefied gas dynamics, materials theory, and computational chemistry, and are essential to capture the thermophysical behavior, relaxation phenomena, and collective effects unique to polyatomic systems. Prominent examples include polyatomic extensions of Boltzmann and BGK kinetic models, moment hierarchies for gas dynamics (e.g., ET15), classical many-body potentials for condensed matter, and molecular models for simulations such as water in atomistic force fields.

1. Structural Foundation: Kinetic and Statistical Models for Polyatomic Gases

Classical kinetic models for polyatomic gases describe the phase-space evolution of a molecular distribution function $f(t,x,v,I,\ldots)$ where $v$ is translational velocity and $I$ encodes continuous or discrete internal energy variables associated with rotations, vibrations, or electronic states (Pirner, 2018). The governing equation is typically a Boltzmann-type or relaxation equation: $\partial_t f + v \cdot \nabla_x f = Q[f, f].$ Here, $Q$ is a binary-collision operator that, for polyatomics, must properly account for both collisional transfer of translational and internal energy. For example, in the polyatomic BGK model, collision frequency and BGK-type relaxation are partitioned between translational and internal relaxation scales, and equilibrium is described by a Maxwellian with appropriate dependence on both $v$ and $I$ (Pirner, 2018).

In addition to Boltzmann-type models, lattice and mixture models in statistical thermodynamics treat polyatomicity by incorporating multisite occupancy, internal partition functions, and specific interaction energies—as in the quasi-chemical approximation for polyatomic mixtures (Dávila et al., 2016).

2. Extended Thermodynamics: Hierarchies of Moments and Closure Schemes

A central theme in polyatomic classical modeling is developing macroscopic equations by taking moments of the kinetic equation. For polyatomic gases, this produces a hierarchy indexed by products of translational velocities and internal energies, leading to systems with, e.g., 15 fields (ET15 model), which encompass mass, momentum, energy, heat flux, stress tensor, dynamic pressure, and high-order couplings between translational and internal modes (Arima et al., 2020, Arima et al., 2022).

The closure of these moment systems is achieved using the maximum entropy principle (MEP), yielding explicit expressions for the non-equilibrium distribution as a local Maxwellian perturbed by Lagrange multipliers conjugate to the macroscopic nonequilibrium variables (e.g., dynamic pressure, heat flux, high-order energy mixing) (Arima et al., 2020).

A distinct feature of ET15 and similar models is the presence of new moments, such as $H_{llmm}$, which track the non-equilibrium correlation between translational and internal energies and enable correct multi-scale relaxation dynamics (e.g., separation of translational and internal temperature equilibration scales) (Arima et al., 2020).

3. Relaxation Mechanisms: BGK and ES-BGK Models for Polyatomic Systems

The BGK family of kinetic models—particularly the ellipsoidal-statistical BGK (ES-BGK) extensions—provide tractable relaxation approximations to the Boltzmann collision operator. In polyatomic versions, the collision term relaxes $f$ toward a Maxwellian that appropriately incorporates internal energy distribution, with tunable rates for translational and internal (rotational, vibrational) degrees of freedom (Park et al., 2017, Pirner, 2018, Dauvois et al., 2020). Parameters such as the internal-to-translational relaxation number (or $\theta$) and Prandtl correction ($\nu$) enable simultaneous control of transport coefficients such as viscosity, bulk viscosity, and thermal conductivity (Kolluru et al., 2022, Dauvois et al., 2020).

Key physical consequences include the emergence of two (or more) relaxation timescales, explicit exponential convergence to equilibrium with spectral decay rates governed by the slowest process (redistribution among internal and translational modes vs. velocity-space relaxation), and an H-theorem ensuring entropy dissipation (Pirner, 2018, Park et al., 2017). The transition between "fast" ($Z_k\gg \nu_{kk}$) and "slow" ($Z_k\ll \nu_{kk}$) regimes is captured analytically in decay-rate formulas (Pirner, 2018), and models can handle multiple species with cross-collisional dynamics.

4. Many-body and Lattice Models: Classical Potentials and Mixtures

The polyatomic classical model framework extends to many-body, lattice, and fluid models. Bond-counting potentials capture covalent valency effects in elemental melts by enforcing local coordination constraints and assigning potential energy based on the allowed number of bonds per atom (valency), with parameters for bonding shells, bond strengths, and hard-core repulsions (Matityahu et al., 2018). Analytic solutions for these models in 1D yield exact partition functions, equations of state, and reveal phase transitions such as liquid–liquid transitions (LLPT), nucleation kinetics, and defect energetics.

Statistical thermodynamics of adsorbed mixtures—such as rods (k-mers, l-mers) on a lattice—are governed by multisite occupancy constraints, internal partition sums, and local interaction energies, and are analytically tractable in exact 1D or with quasi-chemical approximation in higher $d$ (Dávila et al., 2016). These models predict strong entropic effects due to multi-site exclusion, nontrivial adsorption isotherms, and ordering transitions not captured by simpler (e.g., Bragg-Williams) mean-field approaches.

5. Model Reduction and Numerical Schemes

Reduction techniques simplify polyatomic kinetic models by integrating out internal variables via moment closure or by representing internal energy via a minimal set of auxiliary fields (e.g., a rotational energy density governed by an advection-diffusion-relaxation equation) (Kolluru et al., 2022). Such reduction leads to models with a single distribution function coupled to a handful of macroscopic internal-energy fields and enables efficient numerical schemes (such as lattice Boltzmann implementations) that retain correct conservation, H-theorems, and hydrodynamic limits with independently tunable transport coefficients.

Numerical methods for the kinetic equations—particularly semi-Lagrangian schemes—leverage the relaxation structure to achieve unconditional stability for arbitrary Knudsen numbers, preserve conservation laws, positivity, and reproduce the correct Prandtl number (Boscarino et al., 2020).

6. Molecular and Computational Models: Applications to Polyatomic Fluids and Materials

Classical polyatomic models are pervasive in atomistic simulations. Water models, such as the Optimal Point-charge (OPC) model, are constructed by optimizing multipole moments and Lennard-Jones parameters to best reproduce bulk thermodynamic and dielectric properties, abandoning arbitrary geometric constraints in favor of fitting to key electrostatic observables (Izadi et al., 2014).

In radiative association and reaction dynamics, polyatomic extensions of classical quasiclassical trajectory (QCT) methods treat internal degrees of freedom via semiclassical quantization, compute radiation using the classical Larmor formula, and rely on global potential and dipole surfaces for dynamics in high-dimensional phase space (Szabó et al., 2021).

7. Physical Insights, Regimes, and Applicability

Polyatomic classical models provide a rigorous yet tractable way to capture the essential physics of systems with molecular complexity, enabling the description of multiple relaxation channels, entropy production, non-equilibrium transitions, and material properties across kinetic and hydrodynamic regimes (Pirner, 2018, Arima et al., 2020, Park et al., 2017). Their validity rests on the timescale separation between internal relaxation and macroscopic flow or on the appropriateness of treating collective effects within classical statistics.

Constraints and limitations arise from model assumptions (e.g., rigid vs. flexible molecules, form of collision operators, restriction to nearest-neighbor interactions or specific valency constraints), but the polyatomic classical model paradigm underpins—and unifies—a wide range of contemporary theoretical, numerical, and applied studies in gases, plasmas, materials, and liquid-state physics.