Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forcing fields”

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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Utilizing Novel Field and Data Exploration Methods to Explore Hot Moments in High-Frequency Soil Nitrous Oxide Emissions Data: Opportunities and Challenges

Soil nitrous oxide (N 2 O) emissions are an important driver of climate change and are a major mechanism of labile nitrogen (N) loss from terrestrial ecosystems. Evidence increasingly suggests that locations on the landscape that experience biogeochemical fluxes disproportionate to the surrounding matrix (hot spots) and time periods that show disproportionately high fluxes relative to the background (hot moments) strongly influence landscape-scale soil N 2 O emissions. However, substantial uncertainties remain regarding how to measure and model where and when these extreme soil N 2 O fluxes occur. High-frequency datasets of soil N 2 O fluxes are newly possible due to advancements in field-ready instrumentation that uses cavity ring-down spectroscopy (CRDS). Here, we outline the opportunities and challenges that are provided by the deployment of this field-based instrumentation and the collection of high-frequency soil N 2 O flux datasets. While there are substantial challenges associated with automated CRDS systems, there are also opportunities to utilize these near-continuous data to constrain our understanding of dynamics of the terrestrial N cycle across space and time. Finally, we propose future research directions exploring the influence of hot moments of N 2 O emissions on the N cycle, particularly considering the gaps surrounding how global change forces are likely to alter N dynamics in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Arctic Aerosol Sources and Mixing States Field Campaign Report

The Arctic is warming at a faster rate than anywhere else on Earth, with rapidly shrinking sea ice extent transforming the region. Depending on chemical and physical properties, atmospheric aerosols directly scatter and/or absorb radiation, serve as cloud droplet and/or ice crystal nuclei, and/or reduce the reflectiveness of the snow surface, thereby altering the atmospheric energy budget. There is a wide spread in the magnitude of simulated arctic aerosol radiative forcing, and significant differences in aerosol concentration levels and seasonal cycles often exist between models and observations. Increasing local natural and anthropogenic emissions are significant, with uncertain climate impacts due to complex feedbacks. Model evaluations, however, are limited by the dearth of arctic aerosol observations available and an inadequate understanding of arctic aerosol processes. The majority of previous arctic aerosol observations have been made through intensive spring/summer field campaigns, with few intensive measurement studies focused on the fall-winter transition, a period when freeze-up is occurring later and thinning sea ice is resulting in wintertime ice fracturing. Aerosol monitoring at arctic coastal stations has provided knowledge of long-term seasonal trends in aerosol concentrations. The completed ARM field campaign addresses observational and knowledge gaps through detailed aerosol size and chemical composition measurements during the fall-winter transition in the coastal Arctic and through the entire annual cycle in the central Arctic. The observations are improving our understanding of the sources and processes controlling the aerosol population in the rapidly changing Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Topological Spin Textures in an Insulating van der Waals Ferromagnet

Generation and control of topological spin textures constitutes one of the most exciting challenges of modern spintronics given their potential applications in information storage technologies. Of particular interest are magnetic insulators, which due to low damping, absence of Joule heating and reduced dissipation can provide energy-efficient spin-textures platform. Here, it is demonstrated that the interplay between sample thickness, external magnetic fields, and optical excitations can generate a prolific paramount of spin textures, and their coexistence in insulating CrBr 3 van der Waals (vdW) ferromagnets. Using high-resolution magnetic force microscopy and large-scale micromagnetic simulation methods, the existence of a large region in T-B phase diagram is demonstrated where different stripe domains, skyrmion crystals, and magnetic domains exist and can be intrinsically selected or transformed to each-other via a phase-switch mechanism. Lorentz transmission electron microscopy unveils the mixed chirality of the magnetic textures that are of Bloch-type at given conditions but can be further manipulated into Néel-type or hybrid-type via thickness-engineering. The topological phase transformation between the different magnetic objects can be further inspected by standard photoluminescence optical probes resolved by circular polarization indicative of an existence of exciton-skyrmion coupling mechanism. The findings identify vdW magnetic insulators as a promising framework of materials for the manipulation and generation of highly ordered skyrmion lattices relevant for device integration at the atomic level.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dark kinetic heating of exoplanets and brown dwarfs

Dark kinetic heating of neutron stars has been previously studied as a promising dark matter detection avenue. Kinetic heating occurs when dark matter is sped up to relativistic speeds in the gravitational well of high-escape velocity objects, and deposits kinetic energy after becoming captured by the object, thereby increasing its temperature. We show that dark kinetic heating can be significant even in objects with low-escape velocities, such as exoplanets and brown dwarfs, increasing the discovery potential of such searches. This can occur if there is a long-range dark force, creating a “dark escape velocity”, leading to heating rates substantially larger than those expected from neutron stars. We consequently set constraints on dark sector parameters using Wide-field Infrared Survey Explorer and JWST data on Super-Jupiter WISE 0855-0714, and map out future sensitivity to the dark matter scattering cross section below 10 −40 cm 2 . We compare dark kinetic heating rates of other lower escape velocity objects such as the Earth, Sun, and white dwarfs, finding complementary kinetic heating signals are possible depending on particle physics parameters.

Astrophysics↗
DOE OSTI2021

On the in-plane vibrations and electromechanical resonance characteristics of non-uniformly polarized rectangular piezoelectric wafers: Selective mode-type excitation and specific mode enhancement

Here, we investigate the in-plane vibrations and electromechanical resonance characteristics of non-uniformly polarized rectangular piezoelectric wafers. Non-uniform polarization is represented as a non-uniform electromechanical coupling coefficient using a polarization function. Governing equations are derived for the forced in-plane vibrations of a thin wafer under the assumption of generalized plane stress. The effect of non-uniform polarization is explicitly obtained in the forcing terms of the governing equations. These equations are then recast into a variational weak form that is then solved using the finite element method to obtain the displacement fields for different modes. The electromechanical response of a non-uniformly polarized piezoelectric wafer is derived in terms of the out-of-plane displacement profile on the surface of the wafer. Using the derived analytical expression, a necessary and sufficient condition for the presence/absence of a vibrational mode in the electromechanical impedance spectrum is obtained. Based on this condition, criteria for selective mode-type excitation and specific mode enhancement of vibrational modes in the electromechanical impedance spectrum are postulated. Selective mode-type excitation of in-plane extensional and shear modes is demonstrated for a square wafer and that of in-plane bending modes is demonstrated for a rectangular wafer. Specific mode enhancement is demonstrated for both square and rectangular wafers. In addition, it is also demonstrated how the criteria can be used to suppress specific vibrational modes in the electromechanical impedance spectrum. The proposed methodology of using non-uniformly polarized piezoelectric wafers finds application in the design of single element transducers with multi-frequency operation, frequency-tuned receivers/sensors, acoustic holograms, designing acoustic beams of prescribed shape/lobes, and other non-traditional applications such as information storage.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Surface Dynamics of Charge Transport in LaAlO 3 /SrTiO 3 with Time-Resolved Kelvin Probe Force Microscopy

Time-resolved Kelvin probe force microscopy (tr-KPFM) was used to map surface potential dynamics in LAO/STO under its critical thickness, where it acts as an insulator, as polarizing and depolarizing electric fields were applied to planar electrodes in ambient conditions. Two time constants, for polarization and relaxation (i.e., depolarization) each, were derived and were shown to be independent of the direction of electrical field, revealing the reproducibility and recovery of the surfaces. Ambient aging or recent electrical history did not change the surface kinetics as determined by tr-KPFM. During polarization, the faster time constant (~3–6 s) is likely associated with the transport of H+ ions through the full length of the channel. The slower time constant (~40–50 s), on the other hand, is observed across only half the channel width near the higher-potential–more positive–electrode, which is likely associated with an electrochemical reaction such as the dissociation of water. This slower time constant drops to ~10–20 s toward the middle of the channel which could represent the transport of OH– ions. Furthermore, during relaxation, the fast time constant is consistent with polarization and is correlated to the same mechanism while the slower time constant is higher than it was during polarization and may represent the reabsorption of moisture from the atmosphere.

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

Enhancing the Solubility of Semiconducting Polymers in Eco-Friendly Solvents with Carbohydrate-Containing Side Chains

Semiconducting polymers are at the forefront of next-generation organic electronics due to their robust mechanical and optoelectronic properties. However, their extended π-conjugation often leads to materials with low solubilities in common organic solvents, thus requiring processing in high-boiling-point and toxic halogenated solvents to generate thin-film devices. To address this environmental concern, a natural product-inspired side-chain engineering approach was used to incorporate galactose-containing moieties into semiconducting polymers toward improved processability in greener solvents. Novel isoindigo-based polymers with different ratios of galactose-containing side chains were synthesized to improve the solubilities of the organic semiconductors in alcohol-based solvents. The addition of carbohydrate-containing side chains to π-conjugated polymers was found to considerably impact the intermolecular aggregation of the materials and their microstructures in the solid state as confirmed by atomic force microscopy and grazing-incidence wide-angle X-ray scattering. The charge transport characteristics of the new semiconductors were evaluated by the fabrication of organic field-effect transistors prepared from both toxic halogenated and greener alcohol-based solvents. Importantly, the incorporation of carbohydrate-containing side chains was shown to have very little detrimental impact on the electronic properties of the polymer when processed from green solvents.

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

Incorporation of GTR (generation–transport–recombination) in semiconductor simulations

With the emergence of phase change memory, where the devices experience extreme thermal gradients (~100 K/nm) during transitions between low and high resistive states, the study of thermoelectric effects at small scales becomes particularly relevant. We had earlier observed asymmetric melting of self-heated nano-crystalline silicon micro-wires, where current densities of ~10 7 A/cm 2 were forced through the wires by 1 μs, ~30 V pulses. The extreme asymmetry can be explained by the generation of considerable amount of minority carriers, transport under the electric field, and recombination downstream, a heat transfer process we termed as generation–transport–recombination, which is in opposite direction of the electronic-convective heat carried by the majority carriers. Here, we present a full semiconductor physics treatment of this carrier-lattice heat transport mechanism and the contribution of the minority carriers on the evolution of the melt–solid interface, which can be applied to various high-temperature electronic devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Nano-infrared imaging of epitaxial graphene on SiC revealing doping and thickness inhomogeneities

Here, we report on a nano-infrared (IR) imaging and spectroscopy study of epitaxial graphene on silicon carbide (SiC) by using scattering-type scanning near-field optical microscopy (s-SNOM). With nano-IR imaging, we reveal in real space microscopic domains with distinct IR contrasts. By analyzing the nano-IR, atomic force microscopy, and scanning tunneling microscopy imaging data, we conclude that the imaged domains correspond to single-layer graphene, bilayer graphene (BLG), and higher-doped BLG. With nano-IR spectroscopy, we find that graphene can screen the SiC phonon resonance, and the screening is stronger at more conductive sample regions. Our work offers insights into the rich surface properties of epitaxial graphene and demonstrates s-SNOM as an efficient and effective tool in characterizing graphene and possibly other two-dimensional materials.

2D materials↗
DOE OSTI2020

Mixed-flow design for microfluidic printing of two-component polymer semiconductor systems

The rational creation of two-component conjugated polymer systems with high levels of phase purity in each component is challenging but crucial for realizing printed soft-matter electronics. Here, we report a mixed-flow microfluidic printing (MFMP) approach for two-component π -polymer systems that significantly elevates phase purity in bulk-heterojunction solar cells and thin-film transistors. MFMP integrates laminar and extensional flows using a specially microstructured shear blade, designed with fluid flow simulation tools to tune the flow patterns and induce shear, stretch, and pushout effects. This optimizes polymer conformation and semiconducting blend order as assessed by atomic force microscopy (AFM), transmission electron microscopy (TEM), grazing incidence wide-angle X-ray scattering (GIWAXS), resonant soft X-ray scattering (R-SoXS), photovoltaic response, and field effect mobility. For printed all-polymer (poly[(5,6-difluoro-2-octyl-2H-benzotriazole-4,7-diyl)-2,5-thiophenediyl[4,8-bis[5-(2-hexyldecyl)-2-thienyl]benzo[1,2-b:4,5-b′]dithiophene-2,6-diyl]-2,5-thiophenediyl]) [J51]:(poly{[N,N′-bis(2-octyldodecyl)naphthalene-1,4,5,8-bis(dicarboximide)-2,6-diyl]-alt-5,5′-(2,2′-bithiophene)}) [N2200]) solar cells, this approach enhances short-circuit currents and fill factors, with power conversion efficiency increasing from 5.20% for conventional blade coating to 7.80% for MFMP. Moreover, the performance of mixed polymer ambipolar [poly(3-hexylthiophene-2,5-diyl) (P3HT):N2200] and semiconducting:insulating polymer unipolar (N2200:polystyrene) transistors is similarly enhanced, underscoring versatility for two-component π -polymer systems. Mixed-flow designs offer modalities for achieving high-performance organic optoelectronics via innovative printing methodologies.

42 ENGINEERING↗
DOE OSTI2025

Numerical investigation of two-dimensional buoyancy-driven eddies in liquid metal magnetohydrodynamic flows in breeding blankets *

Abstract Lead–lithium flows are key features in the design of tokamak breeding blanket concepts such as the dual-coolant lead–lithium (DCLL). Since they flow under magnetic fields, they are affected by magnetohydrodynamic (MHD) effects. The neutron flux originating in the tokamak plasma heats the breeding blanket channels in a non-uniform manner, inducing buoyancy forces in the liquid metal. Buoyancy may become a source of quasi-two-dimensional (Q2D) turbulence, and the appearance of eddies may affect the transport of heat and tritium across the blanket. Blankets characterized by high-speed liquid metal flows (such as DCLL) will need ceramic insulating walls to reduce the MHD-related pressure drop in the channels. In our simulations, we have used the Q2D model proposed by Sommeria and Moreau (SM82) which is especially suitable for modeling electrically insulating channel flows. On top of that, we have modelled buoyancy forces in the momentum equation using the Oberbeck–Boussinesq approximation. In this work, we include a validation of the implemented Q2D model in buoyancy-driven cases and identify a thermohydraulic configuration that promotes the generation of eddies and the accumulation of tritium. We also show the results provided by our post-processing tool based on the bi-dimensional fast Fourier transform for eddy detection and characterization. We complete our investigation by performing an initial assessment of how eddies can accumulate tritium in the breeding blanket. We conclude by discussing the relevance of preventing accumulation of tritium within the blanket and suggest a possible solution.

Suarez, Daniel (ORCID:0000000251740100)↗
DOE OSTI2022

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background >E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more (> = 10) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. Here this work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Black holes from Fermi ball collapse

Fermi balls are nontopological solitons that can naturally form in an early Universe containing a dark sector with heavy fermions and an attractive interaction mediated by a light scalar field. We compute the Fermi ball mass and radius scaling relations when the potential of the scalar field 𝜑 has a non-negligible quartic coupling 𝜆⁢𝜑4. The resulting Fermi balls reach “saturation” very rapidly, even when their radius is much smaller than the effective Yukawa force range. These objects can therefore grow by mergers or by the accretion of ambient dark fermions until they become so dense that they fall within their Schwarzschild radius and collapse to black holes. This setup, therefore, provides an example of a rather natural and economical dark sector scenario for the formation of primordial black holes.

Hypothetical particle physics models↗
DOE OSTI2025

Why Do Nuclei Stick Together?

Although quantum chromodynamics (QCD) has long been accepted as the underlying theory of the strong nuclear force, connecting the structure and interactions of nuclei to the fundamental parameters of QCD remains challenging. In this talk, I will present new results from lattice QCD and effective field theory that explore how nuclear interactions look in a world where the quarks are much heavier than they are in nature, and discuss lessons that can be learned for nuclear physics as well as searches for physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Ferroelectrics everywhere: Ferroelectricity in magnesium substituted zinc oxide thin films

We demonstrate ferroelectricity in Mg-substituted ZnO thin films with the wurtzite structure. Zn 1–x Mg x O films are grown by dual-cathode reactive magnetron sputtering on (111)-Pt // (0001)-Al 2 O 3 substrates at temperatures ranging from 26 to 200 °C for compositions spanning from x = 0 to x = 0.37. X-ray diffraction indicates a decrease in the c-lattice parameter and an increase in the a-lattice parameter with increasing Mg content, resulting in a nearly constant c/a axial ratio of 1.595 over this composition range. Transmission electron microscopy studies show abrupt interfaces between Zn 1–x Mg x O films and the Pt electrode. When prepared at pO 2 = 0.025, film surfaces are populated by abnormally oriented grains as measured by atomic force microscopy for Mg concentrations >29%. Raising pO 2 to 0.25 eliminates the misoriented grains. Optical measurements show increasing bandgap values with increasing Mg content. When prepared on a 200 °C substrate, films display ferroelectric switching with remanent polarizations exceeding 100 μC cm –2 and coercive fields below 3 MVcm –1 when the Mg content is between ~30% and ~7%. Substrate temperature can be lowered to ambient conditions, and when doing so, capacitor stacks show only minor sacrifices to crystal orientation and nearly identical remanent polarization values; however, coercive fields drop below 2 MV/cm. Using ambient temperature deposition, we demonstrate ferroelectric capacitor stacks integrated directly with polymer substrate surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Bactrian effect: multiple resonances and light Dirac dark matter

The possibility of light dark matter (DM) annihilating through a dark photon (DP) which kinetically mixes (KM) with the Standard Model (SM) hypercharge field is a very attractive scenario. For DM in the interesting mass range below ~1 GeV, it is well known that bounds from the CMB provide a very strong model building constraint forcing the DM annihilation cross section to be roughly 3 orders of magnitude below that needed to reproduce the observed relic density. Under most circumstances this removes the possibility of an s-wave annihilation process for DM in this mass range as would be the case, e.g., if the DM were a Dirac fermion. In an extra-dimensional setup explored previously, it was found that the s-channel exchange of multiple gauge bosons could simultaneously encompass a suppressed annihilation cross section during the CMB era while also producing a sufficiently large annihilation rate during freeze-out to recover the DM relic density. In this paper, we analyze more globally the necessary requirements for this mechanism to work successfully and then realize them within the context of a simple model with two ‘dark’ gauge bosons having masses of a similar magnitude and whose contributions to the annihilation amplitude destructively interfere. We show that if the DM mass threshold lies appropriately in the saddle region of this destructive interference between the two resonance humps it then becomes possible to satisfy these requirements simultaneously provided several ancillary conditions are met. The multiple constraints on the parameter space of this setup are then explored in detail to identify the phenomenologically successful regions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Affleck-Dine cogenesis of baryon and dark matter

We propose a mechanism for cogenesis of baryon and dark matter (DM) in the universe via the Affleck-Dine (AD) route. An AD field which breaks the lepton number symmetry, leads to the generation of lepton asymmetry by virtue of its cosmic evolution, which then gets transferred into lepton and dark sectors. While the lepton asymmetry gets converted into baryon asymmetry via sphalerons, the dark sector asymmetry leads to the final DM abundance with the symmetric part being annihilated away due to resonantly enhanced annihilation, which we choose to be provided by a gauged B – L portal. Stringent constraints from DM direct detection forces DM and B – L gauge boson masses to be light, in the few GeV ballpark. While a large portion of the model parameter space is already ruled out, the remaining parameter space is within sensitivity of laboratory as well as cosmology based experiments. The AD field also plays the role of inflaton with the required dynamics by virtue of its non-minimal coupling to gravity, consistent with observations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗