Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling constant”

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 865 records · Page 48

Dark matter and dark radiation from the early universe with a modulus coupled to the PQMSSM

The supersymmetrized DFSZ axion model is especially compelling in that it contains 1. the SUSY solution to the gauge hierarchy problem, 2. the Peccei-Quinn (PQ) solution to the strong CP problem and 3. the Kim-Nilles solution to the SUSY μ problem. In a string setting, where a discrete R-symmetry ($Z$$^{R}_{24}$ for example) may emerge from the compactification process, a high-quality accidental axion (accion) can emerge from the accidental, approximate remnant global U(1) PQ symmetry where the decay constant f a is linked to the SUSY breaking scale, and is within the cosmological sweet zone. In this setup, one also expects the presence of stringy remnant moduli fields Φ i . Here, we consider the situation of a single light modulus Φ coupled to the PQMSSM in the early universe, with mixed axion plus higgsino-like WIMP dark matter. We evaluate dark matter and dark radiation production via nine coupled Boltzmann equations and assess the severity of the cosmological moduli problem (CMP) along with dark matter and dark radiation production rates. We find that typically the light modulus mass should be m Φ ≳ 10 4 TeV to avoid the moduli-induced dark matter overproduction problem. If one is able to (anthropically) tune the modulus field amplitude, we find a value of Φ 0 ≲ 10 –7 m P would be required to solve the overall CMP.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electronically-coupled redox centers in trimetallic cobalt complexes

Synthesis and isolation of molecular building blocks of metal–organic frameworks (MOFs) can provide unique opportunities for characterization that would otherwise be inaccessible due to the heterogeneous nature of MOFs. Herein, we report a series of trinuclear cobalt complexes incorporating dithiolene ligands, triphenylene-2,3,6,7,10,11-hexathiolate (THT) (1 3+ ), and benzene hexathiolate (BHT) (2 3+ ), with 1,1,1,-tris(diphenylphosphinomethyl)ethane (triphos) employed as the capping ligand. Single crystal X-ray analyses of 1 3+ and 2 3+ display three five-coordinate cobalt centers bound to the triphos and dithiolene ligands in a distorted square pyramidal geometry. Cyclic voltammetry studies of 1 3+ and 2 3+ reveal three redox features associated with the formation of mixed valence states due to the sequential reduction of the redox-active metal centers (Co III/II ). Using this electrochemical data, the comproportionality values were determined for 1 and 2 (log K c = 1.4 and 1.5 for 1, and 4.7 and 5.8 for 2), suggesting strong resonance-stabilized coupling of the metal centers, with stronger electronic coupling observed for complex 2 compared to that for complex 1. Cyclic voltammetry studies were also performed in solvents of varying polarity, whereupon the difference in the standard potentials (ΔE 1/2 ) for 1 and 2 was found to shift as a function of the polarity of the solvent, indicating a negative correlation between the dielectric constant of the electrochemical medium and the stability of the mixed valence species. Spectroelectrochemical studies of in situ generated multi-valent (MV) states of complexes 1 and 2 display characteristic NIR intervalence charge transfer (IVCT) bands, and analysis of the IVCT transitions for complex 2 suggests a weakly coupled class II multi-valent species and relatively large electronic coupling factors (1700 cm –1 for the first multi-valent state of 2 2+ , and 1400 and 4000 cm –1 for the second multi-valent state of 2 + ). Here, density functional theory (DFT) calculations indicate a significant deviation in relative energies of the frontier orbitals of complexes 1 3+ , 2 3+ , and 3 + that contrasts those calculated for the analogous trinuclear cobalt dithiolene complexes employing pentamethylcyclopentadienyl (Cp*) as the capping ligand (Co 3 Cp* 3 THT and Co 3 Cp* 3 BHT, respectively), and may be a result of the cationic nature of complexes 1 3+ , 2 3+ , and 3 + .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic commensuration of a spin moiré superlattice in a layered magnetic semimetal

Spin moiré superlattices (SMSs) have been proposed as a magnetic analog of crystallographic moiré systems and a source of electron minibands offering vector-field moiré tunability and Berry curvature effects. However, it has proven challenging to realize an SMS in which a large exchange coupling J is transmitted between conduction electrons and localized spins. Furthermore, most systems have carrier mean free paths l mfp shorter than their spin moiré lattice constant a spin , inhibiting miniband formation. Here, we discover that the layered magnetic semimetal EuAg 4 Sb 2 overcomes these challenges by forming an interface with J ~ 100 milli–electron volts transferred between a Eu triangular lattice and anionic Ag 2 Sb bilayers hosting a two-dimensional electron band in the ballistic regime (l mfp >> a spin ). The system realizes an SMS with a spin commensurate with the Fermi momentum, leading to a marked quenching of the transport response from miniband formation. Our findings demonstrate an approach to magnetically engineering moiré superlattices and a potential route to an emergent spin-driven quantum Hall state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Energy-conserving physics for nonhydrostatic dynamics in mass coordinate models

Abstract. Motivated by reducing errors in the energy budget related to enthalpy fluxes within the Energy Exascale Earth System Model (E3SM), we study several physics–dynamics coupling approaches. Using idealized physics, a moist rising bubble test case, and the E3SM's nonhydrostatic dynamical core, we consider unapproximated and approximated thermodynamics applied at constant pressure or constant volume. With the standard dynamics and physics time-split implementation, we describe how the constant-pressure and constant-volume approaches use different mechanisms to transform physics tendencies into dynamical motion and show that only the constant-volume approach is consistent with the underlying equations. Using time step convergence studies, we show that the two approaches both converge but to slightly different solutions. We reproduce the large inconsistencies between the energy flux internal to the model and the energy flux of precipitation when using approximate thermodynamics, which can only be removed by considering variable latent heats, both when computing the latent heating from phase change and when applying this heating to update the temperature. Finally, we show that in the nonhydrostatic case, for physics applied at constant pressure, the general relation that enthalpy is locally conserved no longer holds. In this case, the conserved quantity is enthalpy plus an additional term proportional to the difference between hydrostatic pressure and full pressure.

58 GEOSCIENCES↗

Effect of flow on the acoustic performance of extended reaction lined ducts

A model is developed for the effects of uniform and boundary-layer mean flow on the attenuation and propagation of harmonically excited sound waves in an extended reaction lined cylindrical duct. A duct geometry consisting of an annular outer region of bulk material surrounding an inner cylinder of air is utilized. A numerical solution is obtained for the coupled wave equations governing the motion of the sound in both the inner and annular regions. It is found that the numerically predicted attenuation and propagations constants are in excellent agreement with measured values using Kevlar as the liner material for plane-wave mode (O,O) excitation over a wide range of mean flows and sound frequency. The boundary-layer effects are determined to be unimportant, at least for plane-wave sound. In addition, numerical studies indicate small differences between the use of either the radial velocity or the radial displacement boundary conditions.

Hersh, A. S.↗

A viscoplastic theory applied to copper

A phenomenologically based viscoplastic model is derived for copper. The model is thermodynamically constrained by the condition of material dissipativity. Two internal state variables are considered. The back stress accounts for strain-induced anisotropy, or kinematic hardening. The drag stress accounts for isotropic hardening. Static and dynamic recovery terms are not coupled in either evolutionary equation. The evolution of drag stress depends on static recovery, while the evolution of back stress depends on dynamic recovery. The material constants are determined from isothermal data. Model predictions are compared with experimental data for thermomechanical test conditions. They are in good agreement at the hot end of the loading cycle, but the model overpredicts the stress response at the cold end of the cycle.

Freed, Alan D.↗

Ultra-high temperature stability Joule-Thomson cooler with capability to accomodate pressure variations

A Joule-Thomson cryogenic refrigeration system capable of achieving high temperature stabilities in the presence of varying temperature, atmospheric pressure, and heat load is provided. The Joule-Thomson cryogenic refrigeration system includes a demand flow Joule-Thomson expansion valve disposed in a cryostat of the refrigeration system. The expansion valve has an adjustable orifice that controls the flow of compressed gas therethrough and induces cooling and partial liquefaction of the gas. A recuperative heat exchanger is disposed in the cryostat and coupled to the expansion valve. A thermostatically self-regulating mechanism is disposed in the cryostat and coupled to the J-T expansion valve. The thermostatically self-regulating mechanism automatically adjusts the cross sectional area of the adjustable valve orifice in response to environmental temperature changes and changes in power dissipated at a cold head. A temperature sensing and adjusting mechanism is coupled to a cold head for adjusting the temperature of the cold head in response to the change in heat flow in the cold head. The temperature sensing and adjusting mechanism comprises a temperature sensitive diode, a wound wire heater, and an electrical feedback control circuit coupling the diode to the heater. An absolute pressure relief valve is interposed between the output of the cryostat and an exhaust port for maintaining a constant exhaust temperature in the refrigerating system, independent of the changes in atmospheric pressure.

Bard, Steven↗

Gravity and the geoid in the Nepal Himalaya

Materials within the Himalaya are rising due to convergence between India and Asia. If the rate of erosion is comparable to the rate of uplift the mean surface elevation will remain constant. Any slight imbalance in these two processes will lead to growth or attrition of the Himalaya. The process of uplift of materials within the Himalaya coupled with surface erosion is similar to the advance of a glacier into a region of melting. If the melting rate exceeds the rate of downhill motion of the glacier then the terminus of the glacier will receed up-valley despite the downhill motion of the bulk of the glacier. Thus although buried rocks, minerals and surface control points in the Himalaya are undoubtably rising, the growth or collapse of the Himalaya depends on the erosion rate which is invisible to geodetic measurements. Erosion rates are currently estimated from suspended sediment loads in rivers in the Himalaya. These typically underestimate the real erosion rate since bed-load is not measured during times of heavy flood, and it is difficult to integrate widely varying suspended load measurements over many years. An alternative way to measure erosion rate is to measure the rate of change of gravity in a region of uplift. If a control point moves vertically it should be accompanied by a reduction in gravity as the point moves away from the Earth's center of mass. There is a difference in the change of gravity between uplift with and without erosion corresponding to the difference between the free-air gradient and the gradient in the acceleration due to gravity caused by a corresponding thickness of rock. Essentially gravity should change precisely in accord with a change in elevation of the point in a free-air gradient if erosion equals uplift rate. We were funded by NASA to undertake a measurement of absolute gravity simultaneously with measurements of GPS height within the Himalaya. Since both absolute gravity and time are known in an absolute sense to 1 part in 10(exp 10) it is possible to estimate gravity with a precision of 0.1 mu gal. Known systematic errors reduce the measurement to an absolute uncertainty of 6 mu gal. The free air gradient at the point of measurement is typically about 3 mu gals/cm. At Simikot where our experiment was conducted we determined a vertical gravity gradient of 4.4 mu gals/cm.

Bilham, Roger↗

Development of Superconducting Transition Edge Sensors Based on Electron-Phonon Decoupling

We have successfully fabricated a superconducting transition edge sensor (TES), bolometer that centers on the use of electron-phonon decoupling (EPD) for thermal isolation. We have selected a design approach that separates the two functions of far-infrared and THz radiative power absorption and temperature measurement, allowing separate optimization of the performance of each element. We have integrated molybdenum/gold (Mo/Au) bilayer TES and ion assisted thermally evaporated (IAE) bismuth (Bi) films as radiation absorber coupled to a low-loss microstripline from niobium (Nb) ground plane to a twin-slot antenna structure. The thermal conductance (G) and the time constant for the different geometry device have been measured. For one such device, the measured G is 1.16 x 10(exp -10) W/K (plus or minus 0.61 x 10(exp- 10) W/K) at 60 mK, which corresponds to noise equivalent power (NEP) = 1.65 X 10(exp -18)W/vHz and time constant of approximately 5 microseconds.

Jethava, Nikhil↗

Implicit and coupled fluid plasma solver with adaptive Cartesian mesh and its applications to non-equilibrium gas discharges

In this work, we present a new fluid plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (nonlinear, implicit) scheme for non-equilibrium gas discharge plasma. The electrons and ions are described using drift-diffusion approximation coupled to Poisson equation for the electric field. The electron-energy transport equation is solved to account for electron thermal conductivity, Joule heating, and energy loss of electrons in collisions with neutral species. The rate of electron-induced ionization is a function of electron temperature and could also depend on electron density (important for plasma stratification). The ion and gas temperature are kept constant. The transport equations are discretized using a non-isothermal Scharfetter-Gummel scheme to resolve possible large temperature gradients in the sheaths. We demonstrate the new solver for simulations of direct current (DC) and radiofrequency (RF) discharges. The implicit treatment of the coupled equations allows using large time steps. The full-Newton method (FNM) enables fast nonlinear convergence at each time step, offering significantly improved simulation efficiency. We discuss the selection of time steps for solving different plasma problems. The new solver enables solving several problems we could not solve before with existing software: two- and three-dimensional structures of the entire DC discharges including cathode and anode regions, electric field reversals and double-layer formation, the normal cathode spot and an anode ring, moving striations in diffuse and constricted DC discharges, and standing striations in RF discharges. The developed FNM-ACM technique offers many benefits for tackling the disparity of gas discharge plasma systems' time scales and nonlinearity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of the Impedance Resonance of Piezoelectric Multi-Fiber Composite Stacks

Multi-Fiber CompositesTM (MFC's) produced by Smart Materials Corp behave essentially like thin planar stacks where each piezoelectric layer is composed of a multitude of fibers. We investigate the suitability of using previously published inversion techniques for the impedance resonances of monolithic co-fired piezoelectric stacks to the MFCTM to determine the complex material constants from the impedance data. The impedance equations examined in this paper are those based on the derivation. The utility of resonance techniques to invert the impedance data to determine the small signal complex material constants are presented for a series of MFC's. The technique was applied to actuators with different geometries and the real coefficients were determined to be similar within changes of the boundary conditions due to change of geometry. The scatter in the imaginary coefficient was found to be larger. The technique was also applied to the same actuator type but manufactured in different batches with some design changes in the non active portion of the actuator and differences in the dielectric and the electromechanical coupling between the two batches were easily measureable. It is interesting to note that strain predicted by small signal impedance analysis is much lower than high field stains. Since the model is based on material properties rather than circuit constants, it could be used for the direct evaluation of specific aging or degradation mechanisms in the actuator as well as batch sorting and adjustment of manufacturing processes.

Multi-fiber composites↗

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. Here, a vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Magnetoelectric effect arising from a field-induced pseudo Jahn-Teller distortion in a rare-earth magnet

Magnetoelectric materials are attractive for several applications, including actuators, switches, and magnetic field sensors. Typical mechanisms for achieving a strong magnetoelectric coupling are rooted in transition metal magnetism. In sharp contrast, here we identify CsEr(MoO 4 ) 2 as a magnetoelectric material without magnetic transition metal ions, thus ensuring that the Er ions play a key role in achieving this interesting property. Our detailed study includes measurements of the structural, magnetic, and magnetoelectric properties of this material. Bulk characterization and neutron powder diffraction show no evidence for structural phase transitions down to 0.3 K and therefore CsEr(MoO 4 ) 2 maintains the room temperature P2/c space group over a wide temperature range without external magnetic field. These same measurements also identify collinear antiferromagnetic ordering of the Er 3+ moments below T N =0.87K. Complementary dielectric constant and pyroelectric current measurements reveal that a ferroelectric phase with a maximum polarization P~0.6 nC/cm 2 emerges when applying a modest external magnetic field, which indicates that this material has a strong magnetoelectric coupling. Finally, we argue that the magnetoelectric coupling in this system arises from a pseudo Jahn-Teller distortion induced by the magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Iron redox behavior and oxygen reduction activity of Fe-N-C electrocatalysts in different electrolytes

The iron redox behavior and oxygen reduction reaction (ORR) activity of Fe-N-C ORR electrocatalysts synthesized by a variety of techniques were investigated as a function of the identity of the electrolyte anion (bisulfate/sulfate or perchlorate) at a constant pH. In situ X-ray absorption spectroscopy data support the assignment of the redox peaks in the voltammograms to the Fe 3+ /Fe 2+ redox couple. It was found that for a given Fe-N-C catalyst, there is a correlation between the Fe redox couple peak potential and the ORR activity in perchloric acid electrolyte, but not in sulfuric acid electrolyte. While a higher Fe redox couple potential (≥ 110 mV higher) was observed in perchloric acid electrolyte, a higher ORR activity was obtained in sulfuric acid electrolyte. Here, the higher ORR activity observed in sulfuric acid than perchloric acid was correlated with the higher peak current and larger faradaic charge for the Fe redox couple. A study of the Fe redox behavior using a cavity microelectrode, eliminating the impact of ionomer, showed that the interaction of H 2 SO 4 with Fe-N-C is stronger than that of HClO 4 and that Fe redox in both electrolytes is a reversible surface electrochemical reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Weak-Bonding Elements Lead to High Thermoelectric Performance in BaSnS 3 and SrSnS 3 : A First-Principles Study

SnS2, an earth-abundant and ecofriendly material, is limited as a thermoelectric material because of the high lattice thermal conductivity κ L and low carrier mobility μ. By introducing weak-bonding elements Ba or Sr into the SnS 2 framework, we discovered two SnS 2 -based materials BaSnS3 and SrSnS3 with the calculated low κL values of 0.15 and 0.17 W m -1 K -1 , respectively, along the a-axis. The low group velocity and high lattice anharmonicity originating from the weakened and distorted Sn–S bonding network are found in both systems. Moreover, the vibrations of Ba and Sr induce low-lying optical phonons, which strongly couple with the acoustic phonons and strengthen the phonon scattering rates. Compared to SnS 2 , both compounds present lower single-band effective masses, smaller deformation potential constants, and better band convergence, which enhance μ with an insignificantly reduced effective mass. By solving the linearized Boltzmann transport equation with a nonempirical carrier lifetime, we predict excellent ZT values of 2.89 and 2.77 along the a-axis at 900 K in BaSnS 3 and SrSnS 3 , respectively. Further phase diagram calculations of Ba 1–x Sr x SnS 3 solid solutions propose a new compound, Ba 0.5 Sr 0.5 SnS 3 , with an even higher ZT of 3.0. Our work analyzes explicitly how weak-bonding elements enhance μ and suppress κL simultaneously in SnS 2 -analogous systems with a series of compounds nominated as potential high-performance thermoelectric materials.

36 MATERIALS SCIENCE↗

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗

Low-cost calculation and analysis of 2D IR spectra of model diiron trinitrosyl complexes in the NO stretch region with vibrational perturbation theory

Two-dimensional infrared spectroscopy offers unique capabilities for probing vibrational coupling in complex metal–ligand systems. In this paper, we combine two-dimensional infrared spectroscopy with vibrational perturbation theory to investigate vibrational coupling in a diiron trinitrosyl complex across three stable redox states. Although these systems are challenging for electronic structure methods, we demonstrate that key features of experimental 2D IR spectra can be accurately reproduced using reduced-dimensional anharmonic calculations with a small harmonic frequency scaling. Analysis reveals that N–O stretching modes maintain high locality across all redox states, with coupling patterns that directly reflect variations in Fe–N bond strength. Using curvilinear coordinate analysis, we demonstrate these differences result from systematic changes in cubic anharmonic force constants rather than mode delocalization. Our results establish N–O stretches as sensitive probes of metal–ligand bonding strength, expanding the toolkit for studying biologically relevant nitrosyl complexes.

Moran, Hayden A. [Texas A&M University, College St↗

Nonperturbative and perturbative dynamics of a light QCD axion: Dark matter and the strong 𝐶⁢𝑃 problem

Considerable theoretical efforts have gone into expanding the reach of the quantum chromodynamics (QCD) axion beyond its canonical mass–decay-constant relation. The 𝑍 𝒩 QCD axion model reduces the QCD axion mass naturally, by invoking a discrete 𝑍 𝒩 symmetry through which the axion field is coupled to 𝒩 copies of the Standard Model. Before the QCD phase transition at temperature 𝑇 QCD , the 𝑍 𝒩 potential has a minimum at misalignment angle 𝜃 = 𝜋. At 𝑇 QCD , 𝜃 = 𝜋 becomes a maximum; the axion potential becomes exponentially suppressed and develops 𝒩 minima—only one of which actually solves the strong 𝐶⁢𝑃 problem. Before 𝑇 QCD , 𝜃 relaxes toward 𝜋. After 𝑇 QCD , the axion field starts from around the hilltop and may have sufficient kinetic energy to overcome the newly suppressed potential barriers. Such a field evolution leads to nonperturbative effects via the self-interactions near the hilltop, which can cause the exponential growth of fluctuations and backreaction on the coherent motion. This behavior can influence the relic density of the field and the minimum in which it settles. We conduct the first lattice simulations of the 𝑍 𝒩 QCD axion using 𝒞osmoℒattice to accurately calculate dark matter abundances and find nonperturbative dynamics reduce the abundance by up to a factor of two. We furthermore find that the probability of solving the strong 𝐶⁢𝑃 problem tends to diverge considerably from the naïve expectation of 1/𝒩.

Axions↗