Search NASA⌕ Search

SEARCH · Search NASA

Results for “order”

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.

At least 217 records · Page 12

1st- and 2nd-order motion and texture resolution in central and peripheral vision

STIMULI. The 1st-order stimuli are moving sine gratings. The 2nd-order stimuli are fields of static visual texture, whose contrasts are modulated by moving sine gratings. Neither the spatial slant (orientation) nor the direction of motion of these 2nd-order (microbalanced) stimuli can be detected by a Fourier analysis; they are invisible to Reichardt and motion-energy detectors. METHOD. For these dynamic stimuli, when presented both centrally and in an annular window extending from 8 to 10 deg in eccentricity, we measured the highest spatial frequency for which discrimination between +/- 45 deg texture slants and discrimination between opposite directions of motion were each possible. RESULTS. For sufficiently low spatial frequencies, slant and direction can be discriminated in both central and peripheral vision, for both 1st- and for 2nd-order stimuli. For both 1st- and 2nd-order stimuli, at both retinal locations, slant discrimination is possible at higher spatial frequencies than direction discrimination. For both 1st- and 2nd-order stimuli, motion resolution decreases 2-3 times more rapidly with eccentricity than does texture resolution. CONCLUSIONS. (1) 1st- and 2nd-order motion scale similarly with eccentricity. (2) 1st- and 2nd-order texture scale similarly with eccentricity. (3) The central/peripheral resolution fall-off is 2-3 times greater for motion than for texture.

Pattern Recognition, Visual/physiology↗

Universal Decoder for PPM of any Order

A recently developed algorithm for demodulation and decoding of a pulse-position- modulation (PPM) signal is suitable as a basis for designing a single hardware decoding apparatus to be capable of handling any PPM order. Hence, this algorithm offers advantages of greater flexibility and lower cost, in comparison with prior such algorithms, which necessitate the use of a distinct hardware implementation for each PPM order. In addition, in comparison with the prior algorithms, the present algorithm entails less complexity in decoding at large orders. An unavoidably lengthy presentation of background information, including definitions of terms, is prerequisite to a meaningful summary of this development. As an aid to understanding, the figure illustrates the relevant processes of coding, modulation, propagation, demodulation, and decoding. An M-ary PPM signal has M time slots per symbol period. A pulse (signifying 1) is transmitted during one of the time slots; no pulse (signifying 0) is transmitted during the other time slots. The information intended to be conveyed from the transmitting end to the receiving end of a radio or optical communication channel is a K-bit vector u. This vector is encoded by an (N,K) binary error-correcting code, producing an N-bit vector a. In turn, the vector a is subdivided into blocks of m = log2(M) bits and each such block is mapped to an M-ary PPM symbol. The resultant coding/modulation scheme can be regarded as equivalent to a nonlinear binary code. The binary vector of PPM symbols, x is transmitted over a Poisson channel, such that there is obtained, at the receiver, a Poisson-distributed photon count characterized by a mean background count nb during no-pulse time slots and a mean signal-plus-background count of ns+nb during a pulse time slot. In the receiver, demodulation of the signal is effected in an iterative soft decoding process that involves consideration of relationships among photon counts and conditional likelihoods of m-bit vectors of coded bits. Inasmuch as the likelihoods of all the m-bit vectors of coded bits mapping to the same PPM symbol are correlated, the best performance is obtained when the joint mbit conditional likelihoods are utilized. Unfortunately, the complexity of decoding, measured in the number of operations per bit, grows exponentially with m, and can thus become prohibitively expensive for large PPM orders. For a system required to handle multiple PPM orders, the cost is even higher because it is necessary to have separate decoding hardware for each order. This concludes the prerequisite background information. In the present algorithm, the decoding process as described above is modified by, among other things, introduction of an lbit marginalizer sub-algorithm. The term "l-bit marginalizer" signifies that instead of m-bit conditional likelihoods, the decoder computes l-bit conditional likelihoods, where l is fixed. Fixing l, regardless of the value of m, makes it possible to use a single hardware implementation for any PPM order. One could minimize the decoding complexity and obtain an especially simple design by fixing l at 1, but this would entail some loss of performance. An intermediate solution is to fix l at some value, greater than 1, that may be less than or greater than m. This solution makes it possible to obtain the desired flexibility to handle any PPM order while compromising between complexity and loss of performance.

Moision, Bruce E.↗

Discovery of Charge Order in the Transition Metal Dichalcogenide Fe x NbS 2

The Fe intercalated transition metal dichalcogenide (TMD), Fe 1/3 NbS 2 , exhibits remarkable resistance switching properties and highly tunable spin ordering phases due to magnetic defects. We conduct synchrotron x-ray scattering measurements on both underintercalated ( x = 0.32 ) and overintercalated ( x = 0.35 ) samples. We discover a new charge order phase in the overintercalated sample, where the excess Fe atoms lead to a zigzag antiferromagnetic order. The agreement between the charge and magnetic ordering temperatures, as well as their intensity relationship, suggests a strong magnetoelastic coupling as the mechanism for the charge ordering. In conclusion, our results reveal the first example of a charge order phase among the intercalated TMD family and demonstrate the ability to stabilize charge modulation by introducing electronic correlations, where the charge order is absent in bulk 2H - NbS 2 compared to other pristine TMDs.

36 MATERIALS SCIENCE↗

High-Fidelity Large-Signal Order Reduction Approach for Composite Load Model

With the increasing penetration of electronic loads and distributed energy resources, conventional load models cannot capture their dynamics. Therefore, a new comprehensive composite load model is developed by Western Electricity Coordinating Council (WECC). However, this model is a complex high-order non-linear system with multi-time-scale property, which poses challenges on power system studies with heavy computational burden. In order to reduce the model complexity, the authors firstly develop a large-signal order reduction (LSOR) method using singular perturbation theory. In this method, the fast dynamics are integrated into the slow ones to preserve transient characteristics of the former. Then, accuracy assessment conditions are proposed and embedded into the LSOR to improve and guarantee the accuracy of reduced-order model. Finally, the reduced-order WECC composite load model is derived by using the proposed algorithm. Overall, simulation results show that the reduced-order large-signal model significantly alleviates the computational burden while maintaining similar dynamic responses as the original composite load model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Pressure-induced cation and vacancy disorder-order transition and near-zero τ f in deficient hexagonal perovskite Ba 8 ZnTa 6 O 24 dielectrics

Pressure applications can enable the tuning of atomic/defect ordering and provide access to new functional materials. Here, in this study, we report that pressure-induced structural transformation featuring disorder-order transition of both cations and vacancies in the 8-layer deficient hexagonal perovskite tantalate dielectrics Ba 8 ZnTa 6 O 24 , which transformed the structure from twin to shift and remarkably lowered the temperature coefficient of resonant frequency τ f down to near zero (∼0.56 ppm/°C) from 38 ppm/°C for the twinned precursor. The atomic scale STEM-HAADF and EDS results confirm the ordering of Zn in the Ta host at the nanometer scale in the shifted material featuring well-ordered Ba 8 ZnTa 6 O 24 slabs intergrown with Ba 3 ZnTa 2 O 9 and Ba 5 Ta 4 O 15 monolayers and anti-phase grain boundaries as planar defects. The pressure-induced twin-shift structural transformation of Ba 8 ZnTa 6 O 24 features the rare constant concentration of the hexagonal stacked layers, which is allowed by the vacancy ordering at the central layers of face-shared octahedral (FSO) trimers avoiding the FSO B-B repulsion, and remarkably the faster cationic ordering kinetics compared with the 2:1 ordered complex perovskites. Although the inclusion of numerous planar defects and the oxidizable atomic defects led to significant p-type conduction and inhomogeneous electrical microstructures, resulting in an extraordinarily high extrinsic dielectric loss for the high-pressure shifted Ba 8 ZnTa 6 O 24 pellet, the intrinsically near-zero τ f could make the shifted Ba 8 ZnTa 6 O 24 perovskite an ideal microwave dielectric resonator candidate if the defects could be eliminated.

high pressure↗

Modeling wave propagation in elastic solids via high-order accurate implicit-mesh discontinuous Galerkin methods

Here, a high-order accurate implicit-mesh discontinuous Galerkin framework for wave propagation in single-phase and bi-phase solids is presented. The framework belongs to the embedded-boundary techniques and its novelty regards the spatial discretization, which enables boundary and interface conditions to be enforced with high-order accuracy on curved embedded geometries. High-order accuracy is achieved via high-order quadrature rules for implicitly-defined domains and boundaries, whilst a cell-merging strategy addresses the presence of small cut cells. The framework is used to discretize the governing equations of elastodynamics, written using a first-order hyperbolic momentum-strain formulation, and an exact Riemann solver is employed to compute the numerical flux at the interface between dissimilar materials with general anisotropic properties. The space-discretized equations are then advanced in time using explicit high-order Runge–Kutta algorithms. Several two- and three-dimensional numerical tests including dynamic adaptive mesh refinement are presented to demonstrate the high-order accuracy and the capability of the method in the elastodynamic analysis of single- and bi-phases solids containing complex geometries.

42 ENGINEERING↗

Cost-efficient finite-volume high-order schemes for compressible magnetohydrodynamics

We present an efficient dimension-by-dimension finite-volume method which solves the adiabatic magnetohydrodynamics equations at high discretization order, using the constrained-transport approach on Cartesian grids. Results are presented up to tenth order of accuracy. The algorithmic architecture of this method is very close to that of commonly employed second-order schemes: it requires only one reconstructed value per face for each computational cell, independently of the scheme's order. This property is highly beneficial for the numerical efficiency. It results from reusing the required values already available in neighboring grid cells, in contrast to standard algorithms that require a number of reconstructions and evaluations which increases with the scheme's order of accuracy. At a given resolution, these high-order schemes present significantly less numerical dissipation than commonly employed lower-order approaches. Thus, results of comparable accuracy are achievable at a substantially coarser resolution, yielding overall performance gains. We also present a way to include physical dissipative terms: viscosity, magnetic diffusivity and cooling functions, respecting the finite-volume and constrained-transport frameworks. Benefits of this method are shown through applications in turbulent flows.

97 MATHEMATICS AND COMPUTING↗

Exceptional enhancement of mechanical properties in high-entropy alloys via thermodynamically guided local chemical ordering

Understanding the local chemical ordering propensity in random solid solutions, and tailoring its strength, can guide the design and discovery of complex, paradigm-shifting multicomponent alloys. First, we present a simple thermodynamic framework, based solely on binary enthalpies of mixing, to select optimal alloying elements to control the nature and extent of chemical ordering in high-entropy alloys (HEAs). Next, we couple high-resolution electron microscopy, atom probe tomography, hybrid Monte-Carlo, special quasirandom structures, and density functional theory calculations to demonstrate how controlled additions of Al and Ti and subsequent annealing drive chemical ordering in nearly random equiatomic face-centered cubic CoFeNi solid solution. We establish that short-range ordered domains, the precursors of long-range ordered precipitates, inform mechanical properties. Specifically, a progressively increasing local order boosts the tensile yield strengths of the parent CoFeNi alloy by a factor of four while also substantially improving ductility, which breaks the so-called strength–ductility paradox. Finally, we validate the generality of our approach by predicting and demonstrating that controlled additions of Al, which has large negative enthalpies of mixing with the constituent elements of another nearly random body-centered cubic refractory NbTaTi HEA, also introduces chemical ordering and enhances mechanical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Superstructures and magnetic order in heavily Cu-substituted ( Fe 1 – x Cu x ) 1 + y Te

Most iron-based superconductors exhibit stripe-type magnetism, characterized by the ordering vector Q = ($\frac{1}{2},\frac{1}{2}$). In contrast, Fe 1+y Te, the parent compound of the Fe 1+y Te 1–x Se x superconductors, exhibits double-stripe magnetic order associated with the ordering vector Q = ($\frac{1}{2},0$). Here, we use elastic neutron scattering to investigate heavily Cu-substituted (Fe 1–x Cu x ) 1+y Te compounds and reveal that (1) for x ≳ 0.4, short-range magnetic order emerges around the stripe-type vector at Q = ($\frac{1}{2}$ ± δ, $\frac{1}{2}$ ± δ, $\frac{1}{2}$) with δ ≈ 0.05; (2) the short-range magnetic order is associated with a superstructure modulation at Q = ($\frac{1}{3},\frac{1}{3},\frac{1}{2}$), with the magnetic correlation length shorter than that for the superstructure; and (3) for x ≳ 0.55, we observe an additional intergrown phase with higher Cu content, characterized by a superstructure modulation vector Q = ($\frac{1}{3},\frac{1}{3},0$) and magnetic peaks at Q = ($\frac{2}{3},\frac{1}{3},\frac{1}{2}$)/($\frac{1}{3},\frac{2}{3},\frac{1}{2}$). The positions of superstructure peaks suggest that relative to the tetragonal unit cell of Fe 1+y Te, heavy Cu substitution leads to Fe-Cu orderings that expand the unit cell by $\sqrt{2}$ × 3$\sqrt{2}$ times in the ab plane, corroborated by first-principles calculations that suggest the formation of spin chains and spin ladders. Finally, our findings show that stripe-type magnetism is common in magnetically diluted iron pnictides and chalcogenides, despite the varying associated atomic orderings

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Assessment of Local Observation of Atomic Ordering in Alloys via the Radial Distribution Function: A Computational and Experimental Approach

As a powerful analytical technique, atom probe tomography (APT) has the capacity to acquire the spatial distribution of millions of atoms from a complex sample. However, extracting information at the Ångstrom-scale on atomic ordering remains a challenge due to the limits of the APT experiment and data analysis algorithms. The development of new computational tools enable visualization of the data and aid understanding of the physical phenomena such as disorder of complex crystalline structures. Here, we report progress towards this goal using two steps. We describe a computational approach to evaluate atomic ordering in the crystal structure by generating radial distribution functions (RDF). Atomic ordering is rendered as the Fractional Cumulative Radial Distribution Function (FCRDF) which allows for greater visibility of local compositions at short range in the structure. Further, we accommodate in the analysis additional parameters such as uncertainty in the atomic coordinates and the atomic abundance to ascertain short-range ordering in APT data sets. We applied the FCRDF analysis to synthetic and experimental APT data sets for Ni 3 Al. The ability to observe a signal of atomic ordering consistent with the known L1 2 crystal structure is heavily dependent on spatial uncertainty, irrespective of abundance. Detection of atomic ordering is subject to an upper limit of spatial uncertainty of atoms described with Gaussian distributions with a standard deviation of 1.3 Å. The FCRDF analysis was also applied to the APT data set for a six-component alloy, Al 1.3 CoCrCuFeNi. In this case, we are currently able to visualize elemental segregation at the nanoscale, though unambiguous identification of atomic ordering at the Ångstrom (nearest-neighbor) scale remains a goal.

36 MATERIALS SCIENCE↗

Cascade of vestigial orders in two-component superconductors: Nematic, ferromagnetic, s -wave charge- 4 e , and d -wave charge- 4 e states

Electronically ordered states that break multiple symmetries can melt in multiple stages, similarly to liquid crystals. In a partially melted phase, known as vestigial phase, a bilinear made out of combinations of the multiple components of the primary order parameter condenses. Multicomponent superconductors are thus natural candidates for vestigial order since they break both the U⁡(1)-gauge and also time-reversal or lattice symmetries. Here, we use group theory to classify all possible real-valued and complex-valued bilinears of a generic two-component superconductor on a tetragonal or hexagonal lattice. While the more widely investigated real-valued bilinears correspond to vestigial nematic or ferromagnetic order, the little explored complex-valued bilinears correspond to a vestigial charge-4e condensate, which itself can have an underlying s-wave, d x 2 –y 2 -wave, or d xy -wave symmetry. To properly describe the fluctuating regime of the superconducting Ginzburg-Landau action and thus access these competing vestigial phases, we employ both a large-N and a variational method. We show that while vestigial order can be understood as a weak-coupling effect in the large-N approach, it is akin to a moderate-coupling effect in the variational method. Despite these distinctions, both methods yield similar results in wide regions of the parameter space spanned by the quartic Landau coefficients. Specifically, we find that the nematic and ferromagnetic phases are the leading vestigial instabilities, whereas the various types of charge-4⁢e order are attractive albeit subleading vestigial channels. The only exception is for the hexagonal case, in which the nematic and s-wave charge-4e vestigial states are degenerate. Finally, we discuss the limitations of our approach, as well as the implications of our results for the realization of exotic charge-4⁢e states in material candidates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mechanical analysis of heterogeneous materials with higher-order parameters

Abstract Even though heterogeneous porous materials are widely used in a variety of engineering and scientific fields, such as aerospace, energy-storage technology, and bio-engineering, the relationship between effective material properties of porous materials and their underlying morphology is still not fully understood. To contribute to this knowledge gap, this paper adopts a higher-order asymptotic homogenization method to numerically investigate the effect of complex micropore morphology on the effective mechanical properties of a porous system. Specifically, we use the second-order scheme that is an extension of the first-order computational homogenization framework, where a generalized continuum enables us to introduce length scale into the material constitutive law and capture both pore size and pore distribution. Through several numerical case studies with different combinations of porosity, pore shapes, and distributions, we systematically studied the relationship between the underlying morphology and effective mechanical properties. The results highlight the necessity of higher-order homogenization in understanding the mechanical properties and reveal that higher-order parameters are required to capture the role of realistic pore morphologies on effective mechanical properties. Furthermore, for specific pore shapes, higher-order parameters exhibit dominant influence over the first-order continuum.

42 ENGINEERING↗

Atomic-scale probing of short-range order and its impact on electrochemical properties in cation-disordered oxide cathodes

Abstract Chemical short-range-order has been widely noticed to dictate the electrochemical properties of Li-excess cation-disordered rocksalt oxides, a class of cathode based on earth abundant elements for next-generation high-energy-density batteries. Existence of short-range-order is normally evidenced by a diffused intensity pattern in reciprocal space, however, derivation of local atomic arrangements of short-range-order in real space is hardly possible. Here, by a combination of aberration-corrected scanning transmission electron microscopy, electron diffraction, and cluster-expansion Monte Carlo simulations, we reveal the short-range-order is a convolution of three basic types: tetrahedron, octahedron, and cube. We discover that short-range-order directly correlates with Li percolation channels, which correspondingly affects Li transport behavior. We further demonstrate that short-range-order can be effectively manipulated by anion doping or post-synthesis thermal treatment, creating new avenues for tailoring the electrochemical properties. Our results provide fundamental insights for decoding the complex relationship between local chemical ordering and properties of crystalline compounds.

25 ENERGY STORAGE↗

Pressure-induced charge orders and their postulated coupling to magnetism in hexagonal multiferroic LuFe 2 O 4

Hexagonal LuFe 2 O 4 is a promising charge order (CO) driven multiferroic material with high charge and spin-ordering temperatures. The coexisting charge and spin orders on Fe 3+ /Fe 2+ sites result in magnetoelectric behaviors, but the coupling mechanism between the charge and spin orders remains elusive. Here, by tuning external pressure, we reveal three charge-ordered phases with suggested correlation to magnetic orders in LuFe 2 O 4 : (i) a centrosymmetric incommensurate three-dimensional CO with ferrimagnetism, (ii) a non-centrosymmetric incommensurate quasi-two-dimensional CO with ferrimagnetism, and (iii) a centrosymmetric commensurate CO with antiferromagnetism. Experimental in situ single-crystal X-ray diffraction and X-ray magnetic circular dichroism measurements combined with density functional theory calculations suggest that the charge density redistribution caused by pressure-induced compression in the frustrated double-layer [Fe 2 O 4 ] cluster is responsible for the correlated spin-charge phase transitions. The pressure-enhanced effective Coulomb interactions among Fe-Fe bonds drive the frustrated (1/3, 1/3) CO to a less frustrated (1/4, 1/4) CO, which induces the ferrimagnetic to antiferromagnetic transition. Our results not only elucidate the coupling mechanism among charge, spin, and lattice degrees of freedom in LuFe 2 O 4 , but also provide a new way to tune the spin-charge orders in a highly controlled manner.

36 MATERIALS SCIENCE↗

Charge order landscape and competition with superconductivity in kagome metals

In the kagome metals AV 3 Sb 5 (A = K, Rb, Cs), three-dimensional charge order is the primary instability that sets the stage for other collective orders to emerge, including unidirectional stripe order, orbital flux order, electronic nematicity and superconductivity. In this report we use high-resolution angle-resolved photoemission spectroscopy to determine the microscopic structure of three-dimensional charge order in AV 3 Sb 5 and its interplay with superconductivity. Our approach is based on identifying an unusual splitting of kagome bands induced by three-dimensional charge order, which provides a sensitive way to refine the spatial charge patterns in neighbouring kagome planes. Here we found a marked dependence of the three-dimensional charge order structure on composition and doping. The observed difference between CsV 3 Sb 5 and the other compounds potentially underpins the double-dome superconductivity in CsV 3 (Sb,Sn) 5 and the suppression of T c in KV 3 Sb 5 and RbV 3 Sb 5 . Our results provide fresh insights into the rich phase diagram of AV 3 Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermodynamic evidence for a two-component superconducting order parameter in Sr 2 RuO 4

Sr 2 RuO 4 has stood as the leading candidate for a spin-triplet superconductor for 26 years. However, recent NMR experiments have cast doubt on this candidacy and it is difficult to find a theory of superconductivity that is consistent with all experiments. Here, the order parameter symmetry for this material therefore remains an open question. Symmetry-based experiments are needed that can rule out broad classes of possible superconducting order parameters. Here, we use resonant ultrasound spectroscopy to measure the entire symmetry-resolved elastic tensor of Sr 2 RuO 4 through the superconducting transition. We observe a thermodynamic discontinuity in the shear elastic modulus c 66 , which implies that the superconducting order parameter has two components. A two-component p-wave order parameter, such as p x + ip y , naturally satisfies this requirement. As this order parameter appears to have been precluded by recent NMR experiments, we suggest that two other two-component order parameters, namely {d xz ,d yz } and {d x 2 –y 2 ,g xy (x 2 –y 2 ) }, are now the prime candidates for the order parameter of Sr 2 RuO 4 .

36 MATERIALS SCIENCE↗

Structure of Pb(Fe 2/3 W 1/3 )O 3 single crystals with partial cation order

Despite intensive studies on the complex perovskite Pb(Fe 2/3 W 1/3 )O 3 (PFWO) relaxor, understanding the exact nature of its multifunctional properties has remained a challenge for decades. In this work we report a comprehensive structural study of the PFWO single crystals using a combination of synchrotron X-ray diffraction and high-resolution electron microscopy. The set of {h + ½, k + ½, l + ½} superlattice reflections was observed for the first time based on single-crystal synchrotron X-ray experiments (100–450 K) and transmission electron microscopy investigations, which indicates some kind of B-cation ordering in PFWO which had been thought to be totally disordered. It was found that (1) the crystal structure of PFWO should be described by a partly ordered cubic perovskite (i.e. Fm – 3m ), (2) the weak ferromagnetic properties and excess magnetic moment of PFWO can be understood based on non-random distribution of Fe cations between the 4a and 4b sites, and (3) the Pb displacement disorder is present in this material and the cations are probably displaced along the <100> directions. The X-ray diffraction results of this investigation show that partial cation ordering indeed exists in PFWO, which makes it necessary to revisit the generally accepted interpretations of the results obtained up to date. In agreement with X-ray diffraction study the main results of TEM study include: (1) a long range order that can be described with the Fm – 3m symmetry is reliably detected, (2) the coherence length of that long range order is in the order of 1–2 nm and (3) no remarkable chemical inhomogeneity is found in the tested PFWO crystal, excluding the possibility of a compositional ordering arising from substitutional defects in the perovskite structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose–Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh–Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born–Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick’s theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn–Luttinger-type nonconvergence.

74 ATOMIC AND MOLECULAR PHYSICS↗