Thermodynamics of the gray atmosphere. II - Unattainability of detailed balancing.
Unattainability of detailed balancing in thermodynamic equilibrium of gray atmosphere
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Unattainability of detailed balancing in thermodynamic equilibrium of gray atmosphere
It is demonstrated numerically that detailed balance indeed is well satisfied in the multi-state impact parameter formulation of proton-hydrogen atom collisions. It is also demonstrated numerically that detailed balance is a much more sensitive test of computational accuracy than is the charge conservation test usually employed in such impact parameter calculations. It is suggested that detailed balance be used as a check on all impact parameter calculations. That even agreement with detailed balance may not be a sufficient test for accuracy of the results is pointed out.
The detailed balance equations are solved in the ground state terms of 37 ions of C, N, O, Ne, Mg, Si, S and Fe; atomic data for 235 transitions of these ions are tabulated, and 14 line ratios of q = 2,4 ions and eight line ratios of q = 3 ions are graphed. Forbidden emission lines of these ions are in the far and near UV, visible, and near and far IR regions of the spectrum. In addition, detailed calculations of the relative populations of the levels of the ground state are presented as a function of temperature and density.
A far-wing theory in which the validity of the detailed balance principle is maintained in each step of the derivation is presented. The role of the total density matrix including the initial correlations is analyzed rigorously. By factoring out the rapidly varying terms in the complex-time development operator in the interaction representation, better approximate expressions can be obtained. As a result, the spectral density can be expressed in terms of the line-coupling functions in which two coupled lines are arranged symmetrically and whose frequency detunings are omega - 1/2(omega(sub ji) + omega (sub j'i'). Using the approximate values omega - omega(sub ji) results in expressions that do not satisfy the detailed balance principle. However, this principle remains satisfied for the symmetrized spectral density in which not only the coupled lines are arranged symmetrically, but also the initial and final states belonging to the same lines are arranged symmetrically as well.
Traditional inelastic neutron scattering (INS) characterizes excitations—such as phonons and magnons—in condensed matter systems at thermodynamic equilibrium. However, the most intriguing and puzzling many-body effects in open quantum systems often emerge from dissipative dynamics that are inherently out of equilibrium. Here, we use a combination of laser pumping and INS to experimentally observe long-lived nonequilibrium magnons in a two-dimensional (2D) square-lattice Heisenberg antiferromagnet. These nonequilibrium magnons manifest themselves as a violation of detailed balance in the dynamic structure factor and reach steady states under periodic driving, analogous to nonequilibrium steady states in driven dissipative systems. Furthermore, we show that the violation of detailed balance reflects the quantum-mechanical nature of the underlying dynamical system, where out-of-time-ordered correlations of creation and annihilation operators do not satisfy commutation relations. The in operando INS technique developed here provides a new approach to studying nonequilibrium magnons in prototypical 2D quantum magnets and can be extended to other systems, including one-dimensional spin chains and topological many-body spin systems, where nonequilibrium effects are widespread and rich in discovery potential.
Lindblad dynamics and other open-system dynamics provide a promising path towards efficient Gibbs sampling on quantum computers. In these proposals, the Lindbladian is obtained via an algorithmic construction akin to designing an artificial thermostat in classical Monte Carlo or molecular dynamics methods, rather than being treated as an approximation to weakly coupled system-bath unitary dynamics. Recently, Chen, Kastoryano, and Gilyén (arXiv:2311.09207) introduced the first efficiently implementable Lindbladian satisfying the Kubo–Martin–Schwinger (KMS) detailed balance condition, which ensures that the Gibbs state is a fixed point of the dynamics and is applicable to non-commuting Hamiltonians. This Gibbs sampler uses a continuously parameterized set of jump operators, and the energy resolution required for implementing each jump operator depends only logarithmically on the precision and the mixing time. In this work, we build upon the structural characterization of KMS detailed balanced Lindbladians by Fagnola and Umanità, and develop a family of efficient quantum Gibbs samplers using a finite set of jump operators (the number can be as few as one), akin to the classical Markov chain-based sampling algorithm. Compared to the existing works, our quantum Gibbs samplers have a comparable quantum simulation cost but with greater design flexibility and a much simpler implementation and error analysis. Moreover, it encompasses the construction of Chen, Kastoryano, and Gilyén as a special instance.
The slow microstructural evolution of materials often plays a key role in determining material properties. When the unit steps of the evolution process are slow, direct simulation approaches such as molecular dynamics become prohibitive and Kinetic Monte-Carlo (kMC) algorithms, where the state-to-state evolution of the system is represented in terms of a continuous-time Markov chain, are instead frequently relied upon to efficiently predict long-time evolution. The accuracy of kMC simulations however relies on the complete and accurate knowledge of reaction pathways and corresponding kinetics. This requirement becomes extremely stringent in complex systems such as concentrated alloys where the astronomical number of local atomic configurations makes the a priori tabulation of all possible transitions impractical. Machine learning models of transition kinetics have been used to mitigate this problem by enabling the efficient on-the-fly prediction of kinetic parameters. While conventional KMC methods based on transition state theory naturally yield reversible dynamics that exactly obey the detailed balance criterion, providing strong guarantees on the properties of the stationary distribution, many recently-proposed ML-based approaches to barrier predictions provide no such guarantees. In this study, we derive conditions under which physics-informed ML architectures exactly enforce the detailed balance condition by construction, even when relying on non-extensive descriptions of states in terms of local environments around mobile defects. In conclusion, using the diffusion of a vacancy in a concentrated alloy as an example, we show that such ML architectures also exhibit superior performance in terms of prediction accuracy, demonstrating that the imposition of physical constraints can facilitate the accurate learning of barriers at no increase in computational cost.
Two general empirical formalisms are presented for the spectral density which take into account the deviations from the Lorentz line shape in the wing regions of resonance lines. These formalisms satisfy the detailed balance requirement. Empirical line shape functions, which are essential to provide the continuum absorption at different temperatures in various frequency regions for atmospheric transmission codes, can be obtained by fitting to experimental data.
Populations for the first three bound states and the continuum of hydrogen are determined for an isothermal hydrostatic atmosphere at 20,000 K. The atmosphere is treated as optically thin in the Balmer and Paschen continua and illuminated by continuum radiation at these wavelengths with prescribed radiation temperatures. The atmosphere is optically thick in the 2-1, 3-1, 3-2 and c-1 transitions. Three stages of approximation are treated: (1) radiative detailed balance in the 2-1, 3-1 and 3-2 transitions, (2) radiative detailed balance in the 3-1 and 3-2 transitions, and (3) all transitions out of detailed balance. The solution of this problem is nontrivial and presents sufficient difficulty to have caused the failure of at least one rather standard technique. The problem is thus a good archetype against which new methods or new implementations of old methods may be tested.
The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).
Abstract The efficiency potential for single‐junction photovoltaics (PV) is described by the detailed balance model, which requires the elimination of nonradiative recombination and perfect minority carrier collection. Improvements in GaAs, Si, and perovskite PV follow this model. It might be more complex for CdTe, a leading thin‐film PV technology. While lifetime, passivation, and doping goals for 25% efficient CdTe solar cells are largely reached, voltage is ≈20% below the detailed balance limit. Why is that? In Se‐alloyed CdSe x Te 1‐x (Se is required for >20% efficiency) additional losses can occur due to electrostatic and bandgap fluctuations and due to electronic trap states. To understand mechanisms limiting CdSeTe solar cell performance and to suggest improvements, carrier dynamics, and transport in CdSe x Te 1‐x with variation in Se composition and as doping is analyzed. It is shown that trapping, likely due to anion‐site defects and their complexes, is correlated with low charge carrier mobility of 0.1–0.6 cm 2 (Vs) −1 . Even with 1000 ns charge carrier lifetimes, carrier diffusion length is less than the absorber thickness, reducing efficiency to ≈23%. Device simulations are used to analyze the performance of CdSe x Te 1‐x solar cells; thermodynamic models are not sufficient for absorbers with electronic disorder and trapping.
With shifting environmental trends, many Earth system elements may be poised to undergo critical transitions or ‘tipping’. Reliable anticipation of these tipping elements is vital to inform policy decisions. Many of the current methods for tipping point detection are based on loss of resilience or ‘critical slowdown’ of the system as it approaches a tipping point. However, these methods are prone to false alarms; the detected slowdown may be an artifact of nonstationary noise unrelated to tipping behavior. Here, we explore the efficacy of early warning signs based on a nonequilibrium thermodynamics framework. The model-free detection method relies on the increased intrinsic time-irreversibility due to detailed balance breaking, preceding the onset of tipping or instabilities. We demonstrate that these EWSs are effective for tipping point detection and robust against false alarms due to nonstationary noise, using idealized models for two key elements of the Earth system that are prone to tipping: the Atlantic Meridional Overturning Circulation and Arctic sea-ice loss.
Chirality-induced spin selectivity (CISS) describes how chiral molecules and materials generate spin polarization even at thermal equilibrium. This observation has challenged established principles of microscopic reversibility and Onsager reciprocity. We resolve this paradox by formulating a pseudo-Hermitian quantum framework that separates thermodynamic equilibrium from time-reversal symmetry. Within this approach, structural chirality and electron correlations, irrespective of their microscopic origin, are sufficient to produce CISS observables. Chirality enters through a non-local metric η that couples spin and spatial motion, leading to real spectra, unitary evolution, and thermodynamic consistency. The framework predicts a chirality-induced spin magnetic ordering characterized by a spin-displacement order, which reconciles equilibrium spin polarization with detailed balance and explains the persistence of CISS in materials composed of light elements. We derive generalized Onsager-Casimir relations that respect the observed CISS symmetry, i.e., parity $(\mathscr{P})$-odd and time-reversal $(\mathscr{T})$-odd, but exhibiting $\mathscr{PT}$-even symmetry. This approach establishes a coherent foundation for equilibrium CISS and provides a route to link chemical chirality with measurable spin-to-charge conversion effects.
Thermoradiative (TR) cells convert heat to work through emission of thermal radiation. Multijunction thermoradiative cells have received little research interest due to the apparent overlap with energy harvesting limits. Through detailed balance formalism, the present study models both single- and multi-diode TR devices, establishing their performance limits and identifying key factors that influence efficiency. For single-junction TR cells, we derive a relationship for the optimal bandgap and show that higher emitter temperatures increase efficiency. When the receiver is at absolute zero temperature, multijunction TR cells have little advantage over the single-junction cell in terms of the maximum output power density. However, for higher receiver temperatures, the multijunction configuration demonstrates significant improvements in both power and efficiency through the optimization of the chemical potential and bandgap of each diode in the complete ensemble. A 500 K emitter and 300 K receiver TR system can achieve a 21.25% increase in efficiency and a 10% improvement in the power density with multijunction architecture compared to a single junction through chemical potential optimization. These findings suggest that multijunction TR cells offer a promising approach to advancing low-grade heat recovery technologies for efficient heat-to-work conversion.
The suggestion by Dungey (1988) that the gyroviscosity associated with gradients of the off-diagonal elements of the electron pressure tensor can balance a reconnection electric field along a magnetic X line in a collisionless plasma is investigated. The detailed balance of forces in the vicinity of an X line is evaluated using a two-dimensional magnetic field model and a simple model for particle motion. The results show that the gyroviscosity can indeed provide the force required to balance a reconnection electric field in that region, so that neither collisions nor wave turbulence are necessary for reconnection. The results also show that there should not be a significant increase in current from electron acceleration very near an X line. Reasonable numerical estimates are obtained for conditions expected in the vicinity of the distant X line in the geomagnetic tail.
A balance ring (18) which is shrunk fit within each disk (12) of a rotor is selectively ground for detail balance. A plurality of openings (20) through the outer edge of the balance ring receive weights during the assembly balance of the rotor. A snap ring (42) retains the weights within the openings.
A balance ring (18) which is shrunk fit within each disk (12) of a rotor is selectively ground for detail balance. A plurality of openings (20) through the outer edge of the balance ring receive weights during the asssembly balance of the rotor. A snap ring (42) retains the weights within the openings.
A detailed study was made of face and annular seals under conditions where boiling, i.e., phase change of the leaking fluid, occurs within the seal. Many seals operate in this mode because of flashing due to pressure drop and/or heat input from frictional heating. Some of the distinctive behavior characteristics of two phase seals are discussed, particularly their axial stability. The main conclusions are that seals with two phase flow may be unstable if improperly balanced. Detailed theoretical analyses of low (laminar) and high (turbulent) leakage seals are presented along with computer codes, parametric studies, and in particular a simplified PC based code that allows for rapid performance prediction: calculations of stiffness coefficients, temperature and pressure distributions, and leakage rates for parallel and coned face seals. A simplified combined computer code for the performance prediction over the laminar and turbulent ranges of a two phase flow is described and documented. The analyses, results, and computer codes are summarized.