Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Satellite Metallicity Enhancement I: Suppressed Star Formation, Stellar Mass Loss, and Enriched Inflow of DESI and EAGLE Galaxies around Massive Clusters

Environmental effects are a primary driver of elevated gas-phase metallicities in galaxies around massive clusters, but the underlying physical mechanisms for this satellite metallicity enhancement (SME) are still unclear. Using the Dark Energy Spectroscopic Instrument (DESI) Data Release 1, we present the first measurement of the average SME as a function of projected cluster-centric distance. The resulting profile reveals three distinct regimes: a steep decline from the cluster center, a plateau near the cluster boundary, and an extended downturn across several cluster radii. Remarkably, the complex shape and amplitude of this observed SME profile are successfully reproduced in the EAGLE cosmological simulation. Drawing insights from EAGLE, we develop a novel satellite chemical evolution model to decompose the observed SME into physical contributions from suppressed star formation, stellar mass loss, and enriched gas inflow. Our analysis shows that continuous accretion of enriched intracluster medium dominates the SME plateau within the cluster virial radius, while mass loss and quenching jointly drive the rapid metallicity decline in the cluster core. Our method disentangles the impacts of three environmental processes on galactic chemical enrichment in the cosmic web, providing a powerful framework for understanding cluster galaxy evolution with current and future spectroscopic surveys.

Lin, Yuanye [Shanghai Jiao Tong U.; Shanghai Jiaot↗
DOE OSTI2021

Yang-Lee edge singularity triggered entanglement transition

Here we show that a class of $\mathscr{PT}$ symmetric non-Hermitian Hamiltonians realizing the Yang-Lee edge singularity exhibits an entanglement transition in the long-time steady state evolved under the Hamiltonian. Such a transition is induced by a level crossing triggered by the critical point associated with the Yang-Lee singularity and hence is first order in nature. At the transition, the entanglement entropy of the steady state jumps discontinuously from a volume-law to an area-law scaling. We exemplify this mechanism using a one-dimensional transverse field Ising model with additional imaginary fields, as well as the spin-1 Blume-Capel model and the three-state Potts model. We further make a connection to the forced-measurement induced entanglement transition in a Floquet nonunitary circuit subject to continuous measurements followed by post-selections. Our results demonstrate a new mechanism for entanglement transitions in non-Hermitian systems harboring a critical point.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Diazo Linker Ligand Promotes Flexibility and Induced Fit Binding in a Microporous Copper Coordination Network

Abstract Flexible organic linkers represent an intuitive and effective strategy to design flexible metal–organic materials. We report herein a systematic study concerning the effect of varying the central bond of mixed pyridyl‐benzoate linkers, L, upon the flexibility of three isostructuralkddtopology microporous coordination networks (CNs) of formula ML 2 :X‐kdd‐1‐Cu,1= L = (E)‐4‐(pyridin‐4‐yldiazenyl)benzoate;X‐kdd‐2‐Cu,2= L = (E)‐4‐(2‐(pyridin‐4‐yl)vinyl)benzoate; the previously reportedX‐kdd‐3‐Cu,3= L = 4‐(pyridin‐4‐ylethynyl)benzoate. As revealed by single crystal x‐ray diffraction (SCXRD) and gas sorption studies,X‐kdd‐1‐Cu, exhibited gate‐opening during CO 2 and hydrocarbon (C2 and C8) sorption experiments whereas the other two CNs did not. Insight into these phase transformations was gained from in situ variable‐pressure and variable temperature powder X‐ray diffraction (PXRD), SCXRD, and modeling. Rotation of ligand1around the diazo bond, torsion angle changes between phenyl and carboxylate moieties, and deformation of the Cu‐based rod building blocks enabled activatedX‐kdd‐1‐Cuto form new phases with C8 isomers and CH 2 Cl 2 , CH 2 Cl 2 inducing contraction of the activated phase. Computational studies suggest that1enables flexibility thanks to its lower barrier of deformation versus2or3. This study teaches that diazo moieties could offer a general strategy to enhance the flexibility of CNs.

Chemistry↗
DOE OSTI2020

Squaring the fermion: The threefold way and the fate of zero modes

In this work, we investigate topological properties and classification of mean-field theories of stable bosonic systems. Of the three standard classifying symmetries, only time reversal represents a real symmetry of the many-boson system, while the other two, particle-hole and chiral, are simply constraints that manifest as symmetries of the effective single-particle problem. For gapped systems in arbitrary space dimension, we establish three fundamental no-go theorems that prove the absence of parity switches, symmetry-protected-topological quantum phases, and localized bosonic zero modes under open boundary conditions. We then introduce a squaring, kernel-preserving map connecting noninteracting Hermitian theories of fermions and stable boson systems, which serves as a playground to reveal the role of topology in bosonic phases and their localized midgap boundary modes. Finally, we determine the symmetry classes inherited from the fermionic tenfold-way classification, unveiling an elegant threefold-way topological classification of noninteracting bosons. We illustrate our main findings in one- and two-dimensional bosonic lattice and field-theory models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Resolution of zigzag magnetic correlations in Na-deficient Na x IrO 3 without long-range ordering

The materials search for Kitaev quantum spin liquids led to the discovery of many honeycomb lattice materials. Much attention has been paid to materials without magnetic order down to the lowest temperatures. The newly synthesized Na-deficient Na x IrO 3 has been found to bear no sign of long-range magnetic order above 1 K from physical property measurements. In this paper, we report momentum-and energy-resolved excitation spectra in Na-deficient Na x IrO 3 measured using a resonant inelastic x-ray scattering spectrometer. Orbital excitation spectra show that the octahedral and trigonal crystal field splittings are larger in Na x IrO 3 than in Na 2 IrO 3 . On the other hand, the low-energy spectrum at low temperature shows a wave-vector dispersion and a spectral weight distribution that are similar to those of Na 2 IrO 3 , revealing that the two-dimensional zigzag magnetic correlations in Na x IrO 3 are similar to those in Na 2 IrO 3 in terms of the ordered magnetic moment direction and three dynamically fluctuating zigzag orders. The azimuth angle dependence of the low-energy spectrum corroborates these results. The two-dimensional zigzag magnetic correlations rapidly weaken until 50 K. At high temperatures, the spectral weight distribution of the low-energy excitation resembles that of the pure Kitaev model, indicating that the Kitaev interaction dominates the dynamic magnetic response at high temperature. Here we suggest that the larger crystal field and distortion and the weakened longer-range Heisenberg exchange interactions due to the Na deficiency in Na x IrO 3 contribute to bring Na x IrO 3 away from the zigzag long-range magnetic order phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Enhanced Deuteron Coalescence Probability in Jets

The transverse-momentum ($p_T$) spectra and coalescence parameters $B_2$ of (anti)deuterons are measured in p-p collisions at $\sqrt{s}=13$ $\text{TeV}$ for the first time in and out of jets. In this measurement, the direction of the leading particle with the highest $p_T$ in the event ($p^{lead}_T > 5 GeV/c$) is used as an approximation for the jet axis. The event is consequently divided into three azimuthal regions, and the jet signal is obtained as the difference between the toward region, that contains jet fragmentation products in addition to the underlying event (UE), and the transverse region, which is dominated by the UE. The coalescence parameter in the jet is found to be approximately a factor of 10 larger than that in the underlying event. This experimental observation is consistent with the coalescence picture and can be attributed to the smaller average phase-space distance between nucleons in the jet cone as compared with the underlying event. The results presented in this Letter are compared to predictions from a simple nucleon coalescence model, where the phase-space distributions of nucleons are generated using pythia8 with the Monash 2013 tuning, and to predictions from a deuteron production model based on ordinary nuclear reactions with parametrized energy-dependent cross sections tuned on data. The latter model is implemented in pythia8.3. Both models reproduce the observed large difference between in-jet and out-of-jet coalescence parameters, although the almost flat trend of the $B^{Jet}_2$ is not reproduced by the models, which instead give a decreasing trend.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

High Fidelity Modeling of WELL58-32 DFITS

Between 2017 and 2019, a series of injection tests were performed at pilot vertical well 58-32 at the Utah FORGE site. These injection activities were conducted in three distinct zones, isolated with packers and bridge plugs during the tests. Zone 1 consists of an open-hole section at the toe of the well, whereas Zones 2 and 3 are three-meter long cased and perforated sections. Several injection cycles were conducted in each zone. A detailed description of the geology of the FORGE site, of well 58-32 and of all injection cycles can be found in (Xing et al., 2020). The injection cycles included Diagnostic Fracture Injection Tests (DFITs). DFITs comprise two primary stages: an injection (pump-in) phase, where controlled fluid injection is employed to create a small hydraulic fracture in the surrounding rock, followed by a shut-in phase, where injection ceases, allowing observation of the formation’s pressure response. Utilizing methodologies such as the G-function analysis, DFITs provide valuable insights, including fracture closure pressure and the Instantaneous Shut-in Pressure (ISIP) and can be employed to estimate the minimum horizontal stress.

42 ENGINEERING↗
DOE OSTI2022

Experiments and numerical modeling of surrogate droplet combustion for a 93 octane gasoline calibration fuel

Gasoline is the world’s most prolific liquid transportation fuel for ground transportation vehicles. Yet, a fundamental understanding of its burning characteristics is still lacking. Complicating this matter is the highly multicomponent nature of gasoline (comprised of hundreds of chemical species) and the seasonable variability of its composition due to variations in the crude oil feedstocks used in refining. The need to bring order to this complex problem has led to development of calibration fuels (Tier II in presented study) which on average will provide repeatable metrics that replicate performance of refined gasoline for various engine designs. This study presents results of experiments using the configuration of an isolated droplet to study the combustion dynamics of calibration fuel and three different surrogates (S4, S5, S6) for the calibration fuel. Droplet burning is under microgravity conditions, where gas transport is determined entirely by fuel evaporation in order to promote spherical symmetry in the burning process. The initial droplet diameters ranged between 0.55 mm and 0.60 mm. The droplets were spark-ignited and their burning histories were recorded by digital video cameras to obtain quantitative data. For each fuel, the evolutions of droplet, flame, and soot-shell diameters were measured to calculate the burn rates, flame-droplet diameter ratios (FSR) and soot-droplet diameter ratios. Additionally, a detailed numerical modeling is presented. This simulation includes multicomponent phase equilibrium, unsteady gas and liquid transport, radiative heat transfer, internal liquid diffusion and detailed combustion chemistry of different fuels during their burning process. Experiment results show that the calibration fuel and surrogates have a close burn rates and FSR, and S5 presents a relatively similar results with calibration fuel. Predictions from the numerical simulation using the surrogate are compared with experiment measurements, which showed a good agreement with the experiment results, including the droplet burning rate and FSR.

Guo, Songtao↗
DOE OSTI2020

Sr-Nd-Pb isotope systematics of Australasian tektites: Implications for the nature and composition of target materials and possible volatile loss of Pb

The Australasian tektite (AAT) strewn field is the largest strewn field on the Earth with about ~10–30% coverage, both land and ocean, but a clearly identified source impact crater is absent despite the young age of AAT of ca. 790 ka. A genetic link between the Australasian tektites and their unequivocal parental materials is therefore largely impossible to establish. Nevertheless, the nature of the parental materials and the extent of volatilization can be constrained using the splash form tektites, carrying the chemical signatures of high-temperature processes, and the layered (so-called Muong Nong-type) tektites, which are less chemically homogenized and exceptionally abundant in the AAT field. New high-precision Sr, Nd and Pb isotopic measurements were obtained for a chemically and petrographically well-characterized suite of AAT, which included the Muong Nong-type (MN-AAT) with precisely known field locations in Laos and splash forms (SF-AAT) from different parts of the strewn field. Additionally, optically dark and light zones of the MN-AAT were also separately analyzed. Homogeneous ε Nd values from -11.8 to -11.2, combined with a narrow range of two-stage Nd model ages from 1.67 to 1.72 Ga for the entire AAT suite, point to a well-mixed source, in terms of REE, of the crustal segment from which the sedimentary material for tektites was ultimately derived. The Sr isotopic data largely overlap for SF-AAT and MN-AAT ( 87 Sr/ 86 Sr = 0.71636–0.72021) and indicate Paleozoic to Mesozoic sedimentary parentage. However, late Neogene to early Quaternary re-deposition and formation of a thick silt-sized sedimentary section with vertical stratification is required to comply with 10Be data. Lead isotope systematics documents at least three different components which can perhaps be represented by different mineral phases, such as feldspar, zircon, organic matter adsorbed on young sediments etc., sorted during fluvial transport and final deposition. In addition, the SF-AAT have systematically lower Pb contents than the MN-AAT, and generally show isotopically heavier Pb isotopic ratios. This is theoretically consistent with a preferential volatilization of lighter Pb isotopes during evaporation and considerably larger Pb loss from SF-AAT when compared to MN-AAT. Nevertheless, further experimental work would be necessary to unambiguously distinguish kinetic fractionation from source mixing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

In Situ High-Temperature Synchrotron Diffraction Studies of (Fe,Cr,Al) 3 O 4 Spinels

The modeling of a loss-of-coolant-accident scenario involving nuclear fuels with FeCrAl cladding materials in consideration to replace a Zircaloy requires knowledge of the thermodynamics of oxidized structures. At temperatures higher than 1500 degrees C, oxidation of FeCrAl alloys forms (Fe,Cr,Al) 3 O 4 spinels. In situ high-energy X-ray diffraction in a conical nozzle levitator installed at beamline 6-ID-D of the APS was used to study the structural evolution of the oxides as a function of the temperature. Single-phase (spinel) and multiphase (spinel-corundum-FeAlO 3 ) regions are mapped as a function of the temperature for three different compositions of FeCrAl oxidation products. The thermal expansion coefficients and cation distribution in the spinel structure have been refined. The temperature at which complete melting of the fuel cladding is expected has been determined by the liquidus temperatures of the oxidized products to be between 1657 and 1834 degrees C in a 20% O 2 /Ar atmosphere using the cooling trace method. The liquidus temperature increases with increasing Al and Cr content in the spinel phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A purely kinetic description of the evaporation of water droplets

The process of water evaporation, although deeply studied, does not enjoy a kinetic description that captures known physics and can be integrated with other detailed processes such as drying of catalytic membranes embedded in vapor-fed devices and chemical reactions in aerosol whose volumes are changing dynamically. In this work, we present a simple, three-step kinetic model for water evaporation that is based on theory and validated by using well-established thermodynamic models of droplet size as a function of time, temperature, and relative humidity as well as data from time-resolved measurements of evaporating droplet size. Additionally, the kinetic mechanism for evaporation is a combination of two limiting processes occurring in the highly dynamic liquid-vapor interfacial region: direct first order desorption of a single water molecule and desorption resulting from a local fluctuation, described using third order kinetics. The model reproduces data over a range of relative humidities and temperatures only if the interface that separates bulk water from gas phase water has a finite width, consistent with previous experimental and theoretical studies. The influence of droplet cooling during rapid evaporation on the kinetics is discussed; discrepancies between the various models point to the need for additional experimental data to identify their origin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Evaluation of RMI Ejecta Sourcing from Reconstructed Interfaces in 2D FLAG Flyer-Plate Simulations

This study investigates the performance of RMI ejecta sourcing from a reconstructed interface in two-dimensional FLAG ejecta flyer-plate simulations. Enabling the RMI ejecta model to operate on reconstructed interfaces has been a longstanding challenge in FLAG, limiting its applicability to experiments in which interface reconstruction is not required. Addressing this limitation is therefore an important step toward expanding the model’s predictive capability. The first phase of this effort focused on identifying a flyer-plate/target configuration that produces large ejecta clouds from two distinct shocks using the standard RMI ejecta model. A series of simulations examined the effects of flyer composition, layer thickness, and mesh resolution on shock structure and ejecta formation. VISAR traces were used to evaluate shock timing and identify configurations that produced a clear double-shock response. A copper middle-layer flyer with a doubled tin layer yielded the most favorable shock characteristics and was selected as the baseline configuration. After confirming this configuration performed as expected when running a reference FLAG version, it became the basis for three simulations: (1) a standard two-dimensional ejecta flyer-plate calculation using a reference FLAG version, (2) the same calculation with the reconstructed interface option enabled, and (3) a calculation employing ALE strategies to force the ejecta model to operate on the reconstructed interface. While the baseline simulation produced the expected ejecta response, the reconstructed interface approach did not result in measurable ejecta, indicating that this capability remains under development. These results highlight both the progress made and the remaining challenges in extending the RMI ejecta model to reconstructed interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Supercurrent-Induced Weyl Superconductivity

Here, we show that Weyl superconductivity can be induced by finite supercurrent in noncentrosymmetric spin–orbit-coupled superconductors with line nodes. We introduce a three-dimensional tight-binding model of a tetragonal superconductor in a D+p-wave pairing state with a finite center-of-mass momentum, and elucidate that a line-nodal to point-nodal spectral transition occurs by applying an infinitesimal supercurrent. We also clarify that the higher-order effect in spin–orbit coupling is particularly important for this phenomenon. The point nodes are protected by topologically nontrivial Weyl charges, and therefore gapless arc states appear on the surface of the superconductor. Furthermore, both the positions and the Weyl charges of the point nodes depend on the direction of the current. In addition, a quantized Berry phase defined on high-symmetry planes characterizes the Weyl nodes when the in-plane supercurrent is considered. Our proposition paves a new way for controlling the superconducting gap structures by using an external field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

“Magic” Orientation Angles to Suppress Spin-Driven Hall Currents in Anisotropic 2D Materials with an Ideal Skyrmion Gas

The Hall effect depending on conduction electron spin projection becomes very different for anisotropic 2D crystals. In this case the spin-dependent electron current strongly determined by the orientation angle, θ, of the sample with respect to an applied electric eld. The spin-up and -down components of the direct and Hall charge currents oscillate with the angle 2θ. The direct and Hall components of the current have the structure where there are the angle-independent part, oscillation amplitude, and phase shift. All three quantitates strongly depend on the electron spin projection, electron mass ratio, and skyrmion size. We nd that there are "magic" orientation angles where the spin-up, spin-down, and total Hall currents vanish. There is a great interest in computations based on 2D materials with skyrmions. Such properties can be useful for computer logic operations based on skyrmions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Fluctuating intertwined stripes in the strange metal regime of the Hubbard model

Strongly correlated electron systems host a variety of poorly understood correlations in their high-temperature normal state. Unlike ordered phases defined by order parameters, regions of the normal state are often defined through unconventional properties such as strange metallic transport or spectroscopic pseudogaps. Characterizing the microscopic correlations in the normal state is necessary to elucidate mechanisms that lead to these properties and their connection to ground-state orders. Here, in this paper, we establish the presence of intertwined charge and spin stripes in the strange metal normal state of the Hubbard model using determinant quantum Monte Carlo calculations. The charge and spin density waves constituting the stripes are fluctuating and short ranged; yet they obey a mutual commensurability relation and remain microscopically interlocked, as evidenced through measurements of three-point spin-spin-hole correlation functions. Our findings demonstrate the ability of many-body numerical simulations to unravel the microscopic correlations that define quantum states of matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

On the stability of the open-string QED neutron and dark matter

We study the stability of a hypothetical QED neutron, which consists of a color-singlet system of two d quarks and a u quark interacting with the QED interaction. As a quark cannot be isolated, the intrinsic motion of the three quarks in the lowest-energy state may lie predominantly in 1 + 1 dimensions, as in a d-u-d open string. The attractive d-u and u-d QED interactions may overcome the weaker repulsive d-d QED interaction to bind the three quarks together. We examine the QED neutron in a phenomenological three-body problem in 1 + 1 dimensions with an effective interaction extracted from Schwinger’s exact QED solution in 1 + 1 dimensions. The phenomenological model in a variational calculation yields a stable QED neutron at 44.5 MeV. The analogous QED proton with two u quarks and a d quark has been found to be too repulsive to be stable and does not have a bound or continuum state, onto which the QED neutron can decay via the weak interaction. Consequently, the QED neutron is stable against the weak decay, has a long lifetime, and is in fact a QED dark neutron. It may be produced following the deconfinement-to-confinement phase transition of the quark gluon plasma in high-energy heavy-ion collisions. Because of the long lifetime of the QED dark neutron, self-gravitating assemblies of QED dark neutrons or dark antineutrons may be good candidates for a part of the primordial dark matter produced during the phase transition of the quark gluon plasma in the evolution of the early Universe.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗