Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modified gravity”

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 127 records · Page 7

Dynamics of active liquid interfaces

Controlling interfaces of phase-separating fluid mixtures is key to the creation of diverse functional soft materials. Traditionally, this is accomplished with surface-modifying chemical agents. Using experiment and theory, we studied how mechanical activity shapes soft interfaces that separate an active and a passive fluid. Chaotic flows in the active fluid give rise to giant interfacial fluctuations and noninertial propagating active waves. At high activities, stresses disrupt interface continuity and drive droplet generation, producing an emulsion-like active state composed of finite-sized droplets. When in contact with a solid boundary, active interfaces exhibit nonequilibrium wetting transitions, in which the fluid climbs the wall against gravity. These results demonstrate the promise of mechanically driven interfaces for creating a new class of soft active matter.

Science & Technology - Other Topics↗

Spectrograph stabilization using a single-delay interferometer on the Hale Telescope

We describe a technique for spectrograph stabilization useful when conventional mitigation techniques of vacuum tanks, thermal insulation, and laser frequency comb may be impractical, expensive, heavy, or bulky. This includes spectrographs on airborne platforms or mounted on telescopes where they suffer a changing gravity vector or other drifts. Placing a fixed-delay interferometer in series with a spectrograph forms an externally dispersed interferometer (EDI). This produces a uniform sinusoidal comb multiplying input spectrum, creating (through heterodyning) beats (moiré patterns). In Fourier space for low frequencies up to the comb frequency, the moiré generated signal counter-rotates to ordinary spectra under an unknown disperser wavenumber drift Δx. This generates a large negative feedback signal useful in a conceptual control loop, to converge rapidly to a stable spectrum and yield Δx. A modified EDI data analysis algorithm (“crossfading”) combines frequency-weighted moiré with conventional spectrum to cancel net output spectrum reaction to Δx. Needing only a single-delay, this is a practical improvement over prior crossfading analyses requiring multiple delays. We test crossfading on ThAr data near 4850 cm−1 taken on Hale telescope in an earlier project. In a single pass, we reduce drift 20 times. Using seven iterations, we reduce 0.5 cm−1 (31 km/s Doppler equivalent) drift to 4×10−7 cm−1 (2.5 cm/s). The interferometer delay can wander, because linearity of phase versus wavenumber interpolates science features between bracketing calibrating spectral references. Second, mathematically reversing the heterodyning effect doubles effective spectral resolution without changing disperser slit.

Erskine, David J [Lawrence Livermore National Labo↗

Convolutional double copy in (anti) de Sitter space

The double copy is a remarkable relationship between gauge theory and gravity that has been explored in a number of contexts, most notably scattering amplitudes and classical solutions. The convolutional double copy provides a straightforward method to bridge the two theories via a precise map for the fields and symmetries at the linearised level. This method has been thoroughly investigated in flat space, offering a comprehensive dictionary both with and without fixing the gauge degrees of freedom. In this paper, we extend this to curved space with an (anti) de Sitter background metric. We work in the temporal gauge, and employ a modified convolution that involves the Mellin transformation in the time direction. As an example, we show that the point-like charge in gauge theory double copies to the (dS-) Schwarzschild black hole solution.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Enhanced, continuous, liquid-liquid extraction and in-situ separation of volatile fatty acids from fermentation broth

In 2018 alone, the US landfilled 35.3 million tons of food waste, about 24% of the total landfilled mass. In addition to the negative impacts landfills have demonstrated on the environment and human health, some states have begun to outlaw or dissuade the disposal of food waste and sewage sludge into landfills altogether. An urgent need has thus been created for the development of digestion processes like anaerobic digestion (AD) and arrested methanogenesis (AM) to convert food waste into valuable chemical products. Unfortunately, the buildup of volatile fatty acids (VFAs) during these processes eventually halts the reaction, and energy efficient methodologies for VFA removal are critical for the operation of fermenters. Additionally, VFAs themselves can serve as valuable chemical precursors, and recently AD processes have been modified to increase VFA production during fermentation. However, even with significant research over the past three decades, the separation of VFAs from the fermenter broth has remained expensive. Moreover, the separation of these VFAs from the fermenter broth may cost up to 50% of the entire process budget, hindering the widespread commercial adoption of AD and AM. Here we present a novel liquid-liquid extraction process termed CLEANS (Continuous Liquid-liquid Extraction And iN-situ Separation) as a highly efficient method for continuously separating VFAs from a real fermentation broth solely under gravity. Our optimized process (using an aqueous broth feed pH of 2.5, tri-noctylamine as an extractant, and a 10:1 ratio of aqueous broth to organic extractant), achieved a VFA distribution constant K D = 44.5 ± 7.9, a single-pass recovery = 81.3 ± 2.5%, and an extraction factor = 8.1 ± 0.3. These KD values are over an order of magnitude higher than what has been previously reported for comparable processes. A high aqueous-to-organic flowrate ratio, enabled for the first time by CLEANS, was found to be particularly crucial for achieving optimal extraction. Our separation process demonstrates excellent reproducibility and potential for scalability. The economic and environmental implications of this work are briefly discussed.

42 ENGINEERING↗

Surface stress tensor and junction conditions on a rotating null horizon

The general form of the surface stress tensor of an infinitesimally thin shell located on a rotating null horizon is derived, when different interior and exterior geometries are joined there. Even though the induced metric on the surface must be the same approached from either side, the first derivatives of the metric need not be. Such discontinuities lead to a Dirac δ-distribution in the Einstein tensor localized on the horizon. For a general stationary axisymmetric geometry the surface stress tensor can be expressed in terms of two geometric invariants that characterize the surface, namely the discontinuities [κ] and [J] of the surface gravity κ and angular momentum density J. The Komar energy and angular momentum are given in coordinates adapted to the Killing symmetries, and the surface contributions to each determined in terms of [κ] and [J] . Guided by these, a simple modification of the original Israel junction conditions is verified directly from the Einstein tensor density to give the correct finite result for the surface stress, when the normal n to the surface is allowed to tend continuously to a null vector. The relation to Israel’s original junction conditions, which fail on null surfaces, is given. The modified junction conditions are suitable to the matching of a rotating “black hole” exterior to any interior geometry joined at the Kerr null horizon surface, even when the surface normal is itself discontinuous and the Barrabès-Israel formalism is also inapplicable. This joining on a rotating null horizon is purely of the matter shell type and does not contain a propagating gravitational shock wave.

79 ASTRONOMY AND ASTROPHYSICS↗

Astrophysical and Theoretical Physics Implications from Multimessenger Neutron Star Observations

Here, the Neutron Star Interior Composition Explorer (NICER) recently measured the mass and equatorial radius of the isolated neutron star PSR J0030+0451. We use these measurements to infer the moment of inertia, the quadrupole moment, and the surface eccentricity of an isolated neutron star for the first time, using relations between these quantities that are insensitive to the unknown equation of state of supranuclear matter. We also use these results to forecast the moment of inertia of neutron star A in the double pulsar binary J0737-3039, a quantity anticipated to be directly measured in the coming decade with radio observations. Combining this information with the measurement of the tidal Love number with LIGO/Virgo observations, we propose and implement the first theory-agnostic and equation-of-state-insensitive test of general relativity. Specializing these constraints to a particular modified theory, we find that consistency with general relativity places the most stringent constraint on gravitational parity violation to date, surpassing all other previously reported bounds by 7 orders of magnitude and opens the path for a future test of general relativity with multimessenger neutron star observations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Phase transitions from the fifth dimension

We study the cosmological transition of 5D warped compactifications, from the high-temperature black-brane phase to the low-temperature Randall-Sundrum I phase. The transition proceeds via percolation of bubbles of IR-brane nucleating from the black-brane horizon. The violent bubble dynamics can be a powerful source of observable stochastic gravitational waves. While bubble nucleation is non-perturbative in 5D gravity, it is amenable to semiclassical treatment in terms of a “bounce” configuration interpolating between the two phases. We demonstrate how such a bounce configuration can be smooth enough to maintain 5D effective field theory control, and how a simple ansatz for it places a rigorous lower-bound on the transition rate in the thin-wall regime, and gives plausible estimates more generally. When applied to the Hierarchy Problem, the minimal Goldberger-Wise stabilization of the warped throat leads to a slow transition with significant supercooling. We demonstrate that a simple generalization of the Goldberger-Wise potential modifies the IR-brane dynamics so that the transition completes more promptly. Supercooling determines the dilution of any (dark) matter abundances generated before the transition, potentially at odds with data, while the prompter transition resolves such tensions. We discuss the impact of the different possibilities on the strength of the gravitational wave signals. Via AdS/CFT duality the warped transition gives a theoretically tractable holographic description of the 4D Composite Higgs (de)confinement transition. Our generalization of the Goldberger-Wise mechanism is dual to, and concretely models, our earlier proposal in which the composite dynamics is governed by separate UV and IR RG fixed points. The smooth 5D bounce configuration we introduce complements the 4D dilaton/radion dominance derivation presented in our earlier work.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Model-dependence of minimal-twist OPEs in d > 2 holographic CFTs

Following recent work on heavy-light correlators in higher-dimensional conformal field theories (CFTs) with a large central charge C T , we clarify the properties of stress tensor composite primary operators of minimal twist, [ T m ], using arguments in both CFT and gravity. We provide an efficient proof that the three-point coupling $\left\langle {\mathcal{O}}_L{\mathcal{O}}_L\left[{T}^m\right]\right\rangle$, where ${\mathcal{O}}_L$ is any light primary operator, is independent of the purely gravitational action. Next, we consider corrections to this coupling due to additional interactions in AdS effective field theory and the corresponding dual CFT. When the CFT contains a non-zero three-point coupling $\left\langle TT{\mathcal{O}}_L\right\rangle$, the three-point coupling $\left\langle {\mathcal{O}}_L{\mathcal{O}}_L\left[{T}^2\right]\right\rangle$ is modified at large C T if $\left\langle TT{\mathcal{O}}_L\right\rangle \sim \sqrt{C_T}$. This scaling is obeyed by the dilaton, by Kaluza-Klein modes of prototypical supergravity compactifications, and by scalars in stress tensor multiplets of supersymmetric CFTs. Quartic derivative interactions involving the graviton and the light probe field dual to ${\mathcal{O}}_L $ can also modify the minimal-twist couplings; these local interactions may be generated by integrating out a spin- ℓ ≥ 2 bulk field at tree level, or any spin ℓ at loop level. These results show how the minimal-twist OPE coefficients can depend on the higher-spin gap scale, even perturbatively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

$T\bar{T}$ in JT Gravity and BF Gauge Theory

JT gravity has a first-order formulation as a two-dimensional BF theory, which can be viewed as the dimensional reduction of the Chern-Simons description of 3d 3 d gravity. We consider {T\overbar{T}} T T ¯ -type deformations of the (0+1) ( 0 + 1 ) -dimensional dual to this 2d 2 d BF theory and interpret the deformation as a modification of the BF theory boundary conditions. The fundamental observables in this deformed BF theory, and in its 3d 3 d Chern-Simons lift, are Wilson lines and loops. In the 3d 3 d Chern-Simons setting, we study modifications to correlators involving boundary-anchored Wilson lines which are induced by a {T\overbar{T}} T T ¯ deformation on the 2d 2 d boundary; results are presented at both the classical level (using modified boundary conditions) and the quantum-mechanical level (using conformal perturbation theory). Finally, we calculate the analogous deformed Wilson line correlators in 2d 2 d BF theory below the Hagedorn temperature where the principal series dominates over the discrete series.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Energy Exascale Earth System Model v2.0

First release of version 2 of the Energy Exascale Earth System Model. The atmosphere component remains EAM. Major changes since version 1 include: all column-physics parameterizations are computed on a separate grid that has approximately half the number of points of the dynamics grid, a new nonhydrostatic dynamical core (running in hydrostatic mode) with semi-Lagrangian tracer transport, CLUBB updated from v1 to v2, a new convective trigger (dCAPE/ULL) based on the dynamic Convective Available Potential Energy (CAPE) (dCAPE) and the Unrestricted Launch Level (ULL) concepts is used in ZM. minimum cloud droplet number changed, gravity wave drag energy conservation fixed and new tunings used, dust emission size distribution changed to emit more coarse dust particles The land component is still ELM. Major changes since version 1 include: using SNICAR-AD for radiation in snow to match the sea-ice model and fixing bugs in snow compaction and water state calculation. The ocean component remains MPAS-ocean. Major change since version 1 include: Redi isopycnal mixing has been updated, tested, and tuned in combination with the Gent-McWilliams parameterization, a sign error was fixed in the 3rd-order flux routines, the mesh used in low-resolution coupled cases was modified and the time steps adjusted, new regionally refined meshes were created, one focused on North America and another on the Southern Ocean, including ice shelf cavities. The sea-ice component remains MPAS-seaice. Major changes since version 1 include: A new heat- and freshwater-conserving coupling of frazil ice, turning off of SSH filtering, addition of SNICAR-AD and snow grain aging. The land-ice component remains MPAS-Albany-landIce (MALI) and is a static ice sheet. There is more out-of-the-box support for cryosphere configurations including the new regionally refined configuration around Antarctica with ice shelf cavities. The river model is MOSART. The half-degree river mesh was redone so it no longer treats Black and Caspian seas as ocean. The coupler remains cpl7/MCT. Major changes since version 1 include handling of ice shelf melt fluxes (heat / freshwater exchange with the ocean), and data icebergs. All components allow regional refinement of their meshes and two separate example of refinement, one in and around North America and one around Antarctica and the Southern Ocean, are provided.

E3SM Project, DOE↗

Engineered Water for Improved Oil Recovery from Fractured Reservoirs (Final Scientific/Technical Report)

The goal of this proposal is to field test the novel technology of engineered water (modifying ionic composition and/or adding surfactants or nanoparticles) injection to improve oil recovery from the carbonate matrix in fractured reservoirs. Carbonate reservoirs tend to be oil-wet/mixed-wet due to positively charged carbonate mineral surfaces and negatively charged acidic/asphaltic components in the oil. Many carbonate reservoirs are also highly fractured. Water injection does not work efficiently in fractured oil-wet reservoirs because water does not imbibe into the matrix. The water composition can be engineered with ionic modification, surfactants and nanoparticles which can change the wettability of carbonate rocks towards more water-wet conditions which can imbibe water into the matrix and push the oil into the fractures. The engineered water can also lower interfacial tension and minimize capillary forces so that gravity can push the oil up and out of the matrix. Wettability alteration also increases water relative permeability which can also improve oil recovery in non-fractured carbonate reservoirs. This technology, if successful, can be a game changer for West Texas carbonates.

02 PETROLEUM↗

Neutron stars with exceptionally light QCD axions

We present a comprehensive study of axion condensed neutron stars that arise in models of an exceptionally light axion that couples to quantum chromodynamics (QCD). These axions solve the strong-charge-parity (𝐶⁢𝑃) problem, but have a mass-squared lighter than that due to QCD by a factor of 𝜖 < 1. Inside dense matter, the axion potential is altered, and much of the matter in neutron stars resides in the axion condensed phase where the strong-𝐶⁢𝑃 parameter 𝜃 = 𝜋 and 𝐶⁢𝑃 remains a good symmetry. In these regions, masses and interactions of nucleons are modified, in turn changing the equation of state (EOS), structure, and phenomenology of the neutron stars. We take the first steps toward the study of the EOS of neutron star matter at 𝜃 = 𝜋 within chiral effective field theory and use relativistic mean field theory to deduce the resulting changes to nuclear matter and the neutron star low-density EOS. We derive constraints on the exceptionally light axion parameter space based on observations of the thermal relaxation of accreting neutron stars, isolated neutron star cooling, and pulsar glitches, excluding the region up to 5 ×10 −7 ≲ 𝜖 ≲ 0.2 for 𝑚 𝑎 ≳ 2 × 10 −9 eV. We comment on potential changes to the neutron star mass-radius relationship, and discuss the possibility of novel, nuclear-density compact objects with 𝜃 = 𝜋 that are stabilized not by gravity but by the axion potential.

axions↗

Diagonal Approximation for Holographic Rényi Entropies

Recently, Dong et al., [A modified cosmic brane proposal for holographic Renyi entropy, J. High Energy Phys. 06 (2024) 120] proposed a modified cosmic brane prescription for computing the Rényi entropy 𝑆𝛼 of a holographic system in the presence of multiple extremal surfaces. This prescription was found by assuming a diagonal approximation, where the Rényi entropy is computed after first measuring the areas of all extremal surfaces. We derive this diagonal approximation for the case of two extremal surfaces and show that it accurately computes Rényi entropies up to 𝑂⁡(log⁡𝐺) corrections. For 𝛼 <1, this allows us to derive the modified cosmic brane prescription, which differs from the original cosmic brane prescription at leading order in 𝐺. For 𝛼 >1, it leads to the original cosmic brane prescription without needing to assume that replica symmetry is unbroken in the bulk.

FOS: Physical sciences↗

Non-invertible global symmetries and completeness of the spectrum

It is widely believed that consistent theories of quantum gravity satisfy two basic kinematic constraints: they are free from any global symmetry, and they contain a complete spectrum of gauge charges. For compact, abelian gauge groups, completeness follows from the absence of a 1-form global symmetry. However, this correspondence breaks down for more general gauge groups, where the breaking of the 1-form symmetry is insufficient to guarantee a complete spectrum. We show that the correspondence may be restored by broadening our notion of symmetry to include non-invertible topological operators, and prove that their absence is sufficient to guarantee a complete spectrum for any compact, possibly disconnected gauge group. In addition, we prove an analogous statement regarding the completeness of twist vortices: codimension-2 objects defined by a discrete holonomy around their worldvolume, such as cosmic strings in four dimensions. We discuss how this correspondence is modified in various, more general contexts, including non-compact gauge groups, Higgsing of gauge theories, and the addition of Chern-Simons terms. Finally, we discuss the implications of our results for the Swampland program, as well as the phenomenological implications of the existence of twist strings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Strong gravity extruding peaks in speed of sound profiles of massive neutron stars

The speed of sound squared (SSS) s 2 in massive neutron stars (NSs) characterizes not only the stiffness of supradense neutron-rich matter within but also equivalently properties of the curved geometry due to the strong-field gravity and matter-geometry coupling. A peaked density or radius profile of s 2 has been predicted for massive NSs using various NS equation of state (EOS) models. However, the nature, cause, location and size of the peak in s 2 profiles are still very EOS model dependent. In this work, we investigate systematically s 2 profiles in massive NSs in a new approach that is independent of the nuclear EOS model and without any presumption about the NS structure and/or composition. In terms of the small quantities (reduced radius, the energy density and pressure scaled by their central values), we perform double-element perturbative expansions in solving perturbatively the scaled Tolman-Oppenheimer-Volkoff (TOV) equations and analyzing s 2 profiles from the Newtonian limit to the general relativistic (GR) case. The GR term in the TOV equations plays a twofold role: it compresses NS matter and modifies the pressure/energy density ratio from small values in Newtonian stars showing no s 2 peak to large ones for massive NSs possessing a peak in their s 2 profiles, and eventually takes away the peak in extremely compact/massive NSs approaching the causality limit. In particular, the peaked behavior in s 2 is expected to emerge near the center of massive NSs like PSR J⁢0740+6620, while a sharp phase transition is unlikely to occur there. Furthermore, these features revealed from our analyses are universal as they are intrinsic properties of the GR stellar structure equations independent of the still very uncertain EOS of supradense neutron-rich matter in NSs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

On a Critical Acceleration Scale of Dark Matter in ΛCDM and Dynamical Dark Energy

Abstract Universal acceleration a 0 emerges in various empirical laws, yet its fundamental nature remains unclear. Using Illustris and Virgo N -body simulations, we focus on the velocity and acceleration fluctuations in collisionless dark matter involving long-range gravity. For comparison, in the kinetic theory of gases, molecules undergo random elastic collisions involving short-range interactions, where only velocity fluctuations are relevant. Hierarchical structure formation proceeds through the merging of smaller halos to form larger halos, which facilitates a continuous energy cascade from small to large halos at a constant rate ε u ≈ −10 −7 m 2 s −3 . Velocity fluctuations involve a critical velocity u c ∝ (1 + z ) −3/4 . Acceleration fluctuations involve a critical acceleration a c ∝ (1 + z ) 3/4 . Two critical quantities are related by the rate of energy cascade ε u ≈ − a c u c /[2(3 π ) 2 ], where factor 3 π is from the angle of incidence during merging. With critical velocity u c on the order of 300 km s −1 at z = 0, the critical acceleration is determined to be a c 0 ≡ a c ( z = 0) ≈ 10 −10 m s −2 , suggesting a c might explain the universal acceleration a 0 ≈ 10 −10 m s −2 in the empirical Tully–Fisher relation or modified Newtonian dynamics. The redshift evolution a c ∝ (1 + z ) 3/4 is in good agreement with Magneticum and EAGLE simulations and in reasonable agreement with limited observations. This suggests a larger a 0 at a higher redshift such that galaxies of fixed mass rotate faster at a higher redshift. Note that for dark energy (DE) density ρ DE 0 ≈ a c 0 2 / G = 1 0 − 10 J m −3 , we postulate an entropic origin of the DE from acceleration fluctuations of dark matter, analogous to the gas pressure from velocity fluctuations. This leads to a dynamical DE coupled to the structure evolution involving a relatively constant DE density followed by a slow weakening phase, suggesting possible deviations from the standard ΛCDM paradigm.

N-body simulations↗

Rapidity-dependent spin decomposition of the nucleon

We revisit the two-dimensional Fourier transform of generalized parton distributions (GPDs) at nonzero skewness. At 𝜂 = 0 it reduces to the standard impact-parameter density, while at 𝜂 ≠ 0 it is an off-forward amplitude that we interpret as a genuine parton–nucleon correlation. Its overall strength (the transverse-plane integral of the density) is fixed by the GPD at the kinematic point 𝑡 =−𝑐 𝜂 =−4⁢𝜂 2 ⁢𝑚$^{2}_{𝑁}$/(1 − 𝜂 2 ) and decreases monotonically with the rapidity gap Δ⁢𝑦 = ln⁡[(1+𝜂)/(1−𝜂)] = 2 artanh⁡(𝜂). This rapidity dependence implies rapidity-modified Ji identities that connect helicity, orbital, and total angular momenta of the correlation in closed form. To quantify these effects, we construct leading-twist quark and gluon GPDs in a string-based conformal framework: conformal moments are parametrized by linear open- and closed-string Regge trajectories with slopes constrained by parton distribution functions (PDFs), hadron/glueball spectroscopy, and form-factor data, and GPDs are reconstructed over the full (𝑥,𝜂,𝑡) domain by Mellin-Barnes inversion with next-to-leading order evolution. We find qualitative agreement (and fair quantitative agreement within quoted uncertainties) for several moments and selected nonsinglet 𝑥-space channels at 𝜇 = 2 GeV when compared with lattice QCD, while we also identify channels with visible tension and discuss likely sources (PDF priors and 𝑡-slope systematics).

Gauge-gravity dualities↗

A modified cosmic brane proposal for holographic Renyi entropy

We propose a new formula for computing holographic Renyi entropies in the presence of multiple extremal surfaces. Our proposal is based on computing the wave function in the basis of fixed-area states and assuming a diagonal approximation for the Renyi entropy. For Renyi index n ≥ 1, our proposal agrees with the existing cosmic brane proposal for holographic Renyi entropy. For n < 1, however, our proposal predicts a new phase with leading order (in Newton’s constant G) corrections to the cosmic brane proposal, even far from entanglement phase transitions and when bulk quantum corrections are unimportant. Recast in terms of optimization over fixed-area states, the difference between the two proposals can be understood to come from the order of optimization: for n < 1, the cosmic brane proposal is a minimax prescription whereas our proposal is a maximin prescription. We demonstrate the presence of such leading order corrections using illustrative examples. In particular, our proposal reproduces existing results in the literature for the PSSY model and high-energy eigenstates, providing a universal explanation for previously found leading order corrections to the n < 1 Renyi entropies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗