SEARCH · Engineering Papers

Results for “dependence”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Distributionally robust facility location problem under decision-dependent stochastic demand

While the traditional facility location problem considers exogenous demand, in some applications, locations of facilities could affect the willingness of customers to use certain types of services, e.g., carsharing, and therefore they also affect realizations of random demand. Moreover, a decision maker may not know the exact distribution of such endogenous demand and how it is affected by location choices. In this paper, we consider a distributionally robust facility location problem, in which we interpret the moments of stochastic demand as functions of facility-location decisions. We reformulate a two-stage decision-dependent distributionally robust optimization model as a monolithic formulation, and then derive exact mixed-integer linear programming reformulation as well as valid inequalities when the means and variances of demand are piecewise linear functions of location solutions. We conduct extensive computational studies, in which we compare our model with a decision-dependent deterministic model, as well as stochastic programming and distributionally robust models without the decision-dependent assumption. Here, the results show superior performance of our approach with remarkable improvement in profit and quality of service under various settings, in addition to computational speed-ups given by formulation enhancements. These results draw attention to the need of considering the impact of location decisions on customer demand within this strategic-level planning problem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Spectrally-dependent positive and negative photo-plastic effects in CdZnTe and CdZnTeSe

We studied the Vickers microhardness HV0.025 of CdZnTe and CdZnTeSe samples upon illumination with quasi-monochromatic light with wavelengths ranging between 750 and 1540 nm. We observed that light with a wavelength close to or above the bandgap of the samples results in a positive photo-plastic effect (hardening), and light in the spectral range well below the band gap can lead both to positive or negative effects depending on the parameters of deep levels in the bandgap. We performed measurements of spectrally dependent photoconductivity and compared them with HV0.025. Based on the analysis, we propose a complex qualitative model explaining the spectral dependence of HV0.025 by the charging and neutralization of deep levels by optical transitions between the valence and conduction bands and between the bands and deep levels in the band gap. Tuning of parameters of deep levels could be used to control the sign and magnitude of the photo-plastic effect.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Pressure-dependent kinetics of cyclopentene consumption

Recent experiments have uncovered deficiencies in the ability of current kinetic models for cyclopentane combustion to predict species profiles as a function of O2 concentration. One likely source of these discrepancies flows from the limited availability of relevant kinetic data for key functionalized intermediates, requiring the use of distant analogies to estimate their rates. This work fills this gap for cyclopentene, the most prominent C5 intermediate in low-temperature cyclopentane combustion. Ab initio transition-state-theory based master equation methods are used to calculate pressure-dependent rate constants across four potential energy surfaces (C5H8 + OH, C5H8 + HO2, C5H7, and C5H7 + O2). Potential energy surfaces are characterized at the CCSD(T)-F12/cc-pVDZ-F12//revDSD-PBEP86-D3BJ/def2-TZVPP level of theory. Pressure-dependent rate constants and branching fractions are determined via 1D master equation calculations. The results are used to analyze the competition for radicals between cyclopentene and cyclopentane in the intermediate temperature region, and to identify the dominant pathways flowing from these initiation channels. On the OH initiation surface, we find good agreement with recently published experimental rate constants at high temperature, and we elucidate the complex temperature dependence below 800 K, where 𝜋𝜋 addition becomes prominent. On the HO2 initiation surface, we find an interesting alternative route to bicyclic ether formation via well-skipping over the HO2 adduct well. On the cyclopentenyl radical surface, the major product is cyclopentadiene plus H. On the cyclopentenyl + O2 surface, we find that resonance stabilization largely prevents the allylic radical from undergoing oxidation, whereas the alkylic radical has a somewhat elevated propensity for oxidation pathways relative to cyclopentyl.

ab initio calculations↗

DOE OSTI2025

The Temperature Dependence of the Langmuir Adsorption Model for a Single-Site Metal–Organic Framework

The single-site Langmuir adsorption model, also known as the Langmuir isotherm equation, is one of the simplest possible descriptions of adsorption phenomena, and yet finds widespread applicability across a range of disciplines. In its simplest form, it is deployed to treat adsorption equilibria at constant temperature (i.e., along isotherms); however, at the heart of its derivation is a more general class of models that each incorporate an explicit temperature dependence, subject to assumptions about the spatial/translational degrees of freedom of the adsorbed species. In this work, measurements of the temperature dependence of supercritical adsorption of H 2 on a single-site metal-organic framework (MOF) are presented, and fitted using a range of Langmuir models with distinct treatments of degrees of freedom in the adsorbed phase. Surprisingly, all of the models can be used to adequately represent the measured data, despite yielding significantly different values for binding energy and the temperature dependence of the isosteric enthalpy of adsorption (i.e., the isosteric heat, qst). Furthermore, a critical finding is that the mean-temperature isosteric enthalpy of adsorption remains consistent across all models within experimental error, highlighting its reliability for evaluating adsorption thermodynamics.

08 HYDROGEN↗

DOE OSTI2022

Importance of imposing gauge invariance in time-dependent density functional theory calculations with meta-generalized gradient approximations

It has been known for more than a decade that the gauge variance of the kinetic energy density τ leads to additional terms in the magnetic orbital rotation Hessian used in linear-response time-dependent density functional theory (TDDFT), affecting excitation energies obtained with τ-dependent exchange–correlation functionals. While previous investigations found that a correction scheme based on the paramagnetic current density has a small effect on benchmark results, we report more pronounced effects here, in particular, for the popular M06-2X functional and for some other meta-generalized gradient approximations (mGGAs). In the first part of this communication, this is shown by a reassessment of a set of five Ni(II) complexes for which a previous benchmark study that did not impose gauge invariance has found surprisingly large errors for excitation energies obtained with M06-2X. These errors are more than halved by restoring gauge invariance. The variable importance of imposing gauge invariance for different mGGA-based functionals can be rationalized by the derivative of the mGGA exchange energy integrand with respect to τ. In the second part, a large set of valence excitations in small main-group molecules is analyzed. For M06-2X, several selected n → π* and π→π$^{*}_{⊥}$ excitations are heavily gauge-dependent with average changes of –0.17 and –0.28 eV, respectively, while π→π$^{*}_{∥}$ excitations are marginally affected (–0.04 eV). Similar patterns, but of the opposite signs, are found for SCAN0. Here, the results suggest that reevaluation of previous gauge variant TDDFT results based on M06-2X and other mGGA functionals is warranted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

First-principles wave-vector- and frequency-dependent exchange-correlation kernel for jellium at all densities

Here we propose a spatially and temporally nonlocal exchange correlation (XC) kernel for the spin-unpolarized fluid phase of ground-state jellium for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact XC kernel to accurately describe the correlation energies of bulk jellium and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel are maintained in the revised MCP07 (rMCP07) kernel, while others are added. The revision $f^{rMCP07}_{XC}$ (q, ω) differs from MCP07 only for nonzero frequencies ω. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wave vector q tends to zero, the kernel has a -4$πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A Distributionally Robust Resilience Enhancement Strategy for Distribution Networks Considering Decision-Dependent Contingencies

When performing the resilience enhancement for distribution networks, there are two obstacles to reliably model the uncertain contingencies: 1) decision-dependent uncertainty (DDU) due to various line hardening decisions, and 2) distributional ambiguity due to limited outage information during extreme weather events (EWEs). Here, to address these two challenges, this paper develops scenario-wise decision-dependent ambiguity sets (SWDD-ASs), where the DDU and distributional ambiguity inherent in EWE-induced contingencies are simultaneously captured for each possible EWE scenario. Then, a two-stage tri-level decision-dependent distributionally robust resilient enhancement (DD-DRRE) model is formulated, whose outputs include the optimal line hardening, distributed generation (DG) allocation, and proactive network reconfiguration strategy under the worst-case distributions in SWDD-ASs. Subsequently, the DD-DRRE model is equivalently recast to a mixed-integer linear programming (MILP)-based master problem and multiple scenario-wise subproblems, facilitating the adoption of a customized column-and-constraint generation (C&CG) algorithm. Finally, case studies demonstrate a remarkable improvement in the out-of-sample performance of our model, compared to its prevailing stochastic and robust counterparts. Moreover, the potential values of incorporating the ambiguity and distributional information are quantitatively estimated, providing a useful reference for planners with different budgets and risk-aversion levels.

decision-dependent uncertainty↗

DOE OSTI2025

Machine-Learning-Based Multiscale Methods for 3D Modelling of Granular Materials by Incorporating History-Dependent State Variables

Over the past decades, the prevalence of machine learning (ML) methods has made the development of ML-based constitutive models for granular materials undoubtedly a popular subject. Numerous studies have been made to feature the loading path or history-dependent stress-strain response of granular media using neural networks. In this work, a novel finite element method (FEM)–ML multiscale approach was developed by incorporating internal variables to improve the simulation accuracy of 3D history-dependent granular materials for the first time. To this end, a surrogate constitutive model based on the single-step-based multi-layer perceptron (MLP) neural network was used to replace representative volume element (RVE) simulations conducted by the discrete element method (DEM) in the multiscale FEM–DEM approach. Although the prediction principle of the MLP aligns with the FEM algorithm, artificially added internal variables are required to differentiate the loading history. To address this issue, history variables associated with the Frobenius norm are proposed to be fed into the MLP coupled with the strain tensor to extract the history-dependent behaviour of granular assemblies. The developed FEM–ML approach was demonstrated in 3D conventional triaxial compression (CTC) simulations. Compared to the multiscale FEM–DEM approach, the proposed FEM–ML method exhibits a significantly improved computational efficiency.

granular materials↗

DOE OSTI2020

Polarization Dependent Excitation and High Harmonic Generation from Intense Mid-IR Laser Pulses in ZnO

The generation of high order harmonics from femtosecond mid-IR laser pulses in ZnO has shown great potential to reveal new insight into the ultrafast electron dynamics on a few femtosecond timescale. In this work we report on the experimental investigation of photoluminescence and high-order harmonic generation (HHG) in a ZnO single crystal and polycrystalline thin film irradiated with intense femtosecond mid-IR laser pulses. The ellipticity dependence of the HHG process is experimentally studied up to the 17th harmonic order for various driving laser wavelengths in the spectral range 3–4 µm. Interband Zener tunneling is found to exhibit a significant excitation efficiency drop for circularly polarized strong-field pump pulses. For higher harmonics with energies larger than the bandgap, the measured ellipticity dependence can be quantitatively described by numerical simulations based on the density matrix equations. The ellipticity dependence of the below and above ZnO band gap harmonics as a function of the laser wavelength provides an efficient method for distinguishing the dominant HHG mechanism for different harmonic orders.

tunneling excitation↗

DOE OSTI2020

Spirometric traits show quantile-dependent heritability, which may contribute to their gene-environment interactions with smoking and pollution

Background: “Quantile-dependent expressivity” refers to a genetic effect that is dependent upon whether the phenotype (e.g., spirometric data) is high or low relative to its population distribution. Forced vital capacity (FVC), forced expiratory volume in 1 second (FEV 1 ), and the FEV 1 /FVC ratio are moderately heritable spirometric traits. The aim of the analyses is to test whether their heritability (h 2 ) is constant over all quantiles of their distribution. Methods: Quantile regression was applied to the mean age, sex, height and smoking-adjusted spirometric data over multiple visits in 9,993 offspring-parent pairs and 1,930 sibships from the Framingham Heart Study to obtain robust estimates of offspring-parent (β OP ), offspring-midparent (β OM ), and full-sib regression slopes (β FS ). Nonparametric significance levels were obtained from 1,000 bootstrap samples. β OP s were used as simple indicators of quantile-specific heritability (i.e.,h 2 = 2β OP /(1+r spouse ), where r spouse was the correlation between spouses). Results: β OP ± standard error (SE) decreased by 0.0009 ± 0.0003 (P = 0.003) with every one-percent increment in the population distribution of FEV 1 /FVC, i.e., β OP ± SE were: 0.182 ± 0.031, 0.152 ± 0.015; 0.136 ± 0.011; 0.121 ± 0.013; and 0.099 ± 0.013 at the 10th, 25th, 50th, 75th, and 90th percentiles of the FEV 1 /FVC distribution, respectively. These correspond to h 2 ± SEs of 0.350 ± 0.060 at the 10th, 0.292 ± 0.029 at the 25th, 0.262 ± 0.020 at the 50th, 0.234 ± 0.025 at the 75th, and 0.191 ± 0.025 at the 90th percentiles of the FEV 1 /FVC ratio. Maximum mid-expiratory flow (MMEF) h 2 ± SEs increased 0.0025 ± 0.0007 (P = 0.0004) with every one-percent increment in its distribution, i.e.: 0.467 ± 0.046, 0.467 ± 0.033, 0.554 ± 0.038, 0.615 ± 0.042, and 0.675 ± 0.060 at the 10th, 25th, 50th, 75th, and 90th percentiles of its distribution. This was due to forced expiratory flow at 75% of FVC (FEF75%), whose quantile-specific h 2 increased an average of 0.0042 ± 0.0008 for every one-percent increment in its distribution. It is speculated that previously reported gene-environment interactions may be partially attributable to quantile-specific h 2 , i.e., greater heritability in individuals with lower FEV 1 /FVC due to smoking or airborne particles exposure vs. nonsmoking, unexposed individuals. Conclusion: Heritabilities of FEV 1 /FVC, MMEF, and FEF75% from quantile-regression of offspring-parent and sibling spirometric data suggest their quantile-dependent expressivity.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

A study on the scale dependence of mixing indices for Eulerian multiphase models

Abstract Mixing can vary based on the scale at which the system is observed, and a mixing index that can capture the features at different length scales is desirable. In this article, we analyze the scale dependence of the mixing indices developed for Eulerian multiphase models. Relevant length scales are distinguished by filtering solid fraction fields. The scale‐dependence study is first done on manufactured fields of solid fraction to assess the performance of the mixing indices. The study is extended to a two‐dimensional CFD simulation of the segregation of a bidisperse gas–solid mixture. The local mixing index performs well in capturing the spatial variation of mixing at different scales. The scale dependence of two global mixing indices is considered in the study, where the state of mixing is defined based on statistical measures. We demonstrate that the choice of measures influences the sensitivity of mixing indices to mixing at different scales.

Nagawkar, Barlev R.↗

DOE OSTI2023

Hydrogen Evolution on Electrode‐Supported Pt n Clusters: Ensemble of Hydride States Governs the Size Dependent Reactivity

Abstract We report the size‐dependent activity and stability of supported Pt 1,4,7,8 for electrocatalytic hydrogen evolution reaction, and show that clusters outperform polycrystalline Pt in activity, with size‐dependent stability. To understand the size effects, we use DFT calculations to study the structural fluxionality under varying potentials. We show that the clusters can reshape under H coverage and populate an ensemble of states with diverse stoichiometry, structure, and thus reactivity. Both experiment and theory suggest that electrocatalytic species are hydridic states of the clusters (≈2 H/Pt). An ensemble‐based kinetic model reproduces the experimental activity trend and reveals the role of metastable states. The stability trend is rationalized by chemical bonding analysis. Our joint study demonstrates the potential‐ and adsorbate‐coverage‐dependent fluxionality of subnano clusters of different sizes and offers a systematic modeling strategy to tackle the complexities.

Zhang, Zisheng↗

DOE OSTI2022

Density Dependence in Three Snake River Sockeye Salmon Nursery Lakes in Central Idaho

Abstract Snake River Sockeye Salmon Oncorhynchus nerka, listed as an endangered species in 1991, currently inhabit three nursery lakes (Redfish, Pettit, and Alturas lakes) in the Sawtooth Valley, Idaho. Conspecific kokanee (lacustrine Sockeye Salmon) are also present in the lakes. Snake River Sockeye Salmon recovery efforts, initially focused on genetic conservation, are now attempting to rebuild naturally spawning populations using hatchery supplementation. However, in Sockeye Salmon nursery lakes, density dependence is frequently observed when elevated O. nerka abundance leads to declines in zooplankton biomass, body size, and shifts in community composition. In turn, these changes lead to reductions in juvenile O. nerka growth rates, survival, and adult returns. We examined a long-term data set of O. nerka population metrics and associated zooplankton community metrics. We found evidence of density dependence within and among nursery lakes. We detected differences in zooplankton biomass, lengths of preferred zooplankton prey (Daphnia spp. and cyclopoid copepods), parr growth rates, and age-1 smolt size among the three lakes. We found negative relationships between O. nerka density and zooplankton biomass and size. We identified positive relationships between zooplankton biomass and two response variables: smolt size at migration and growth rates of hatchery parr. The relationships were generally similar among lakes. Variable outcomes were a result of differences in O. nerka density (or zooplankton biomass), controlled primarily by the relative proportion of spawning and rearing habitat in each lake. Understanding unique lake habitats, ecological interactions, and the role of density dependence is germane to management of Snake River Sockeye Salmon populations.

Griswold, Robert G.↗

DOE OSTI2025

First‐Order Empirical Interpolation Method for Real‐Time Solution of Parametric Time‐Dependent Nonlinear PDEs

ABSTRACT We present a model reduction approach for the real‐time solution of time‐dependent nonlinear partial differential equations (PDEs) with parametric dependencies. A major challenge in constructing efficient and accurate reduced‐order models for nonlinear PDEs is the efficient treatment of nonlinear terms. We address this by unifying the implementation of hyperreduction methods to deal with nonlinear terms. Furthermore, we introduce a first‐order empirical interpolation method (EIM) to provide an efficient approximation of the nonlinear terms in time‐dependent PDEs. We demonstrate the effectiveness of our approach on the Allen–Cahn equation, which models phase separation, and the Buckley–Leverett equation, which describes two‐phase fluid flow in porous media. Numerical results highlight the accuracy, efficiency, and stability of the proposed method compared with both the Galerkin–Newton approach and hyper‐reduced models using the standard EIM.

Nguyen, Ngoc Cuong [Center for Computational Engin↗

DOE OSTI2026

Backbone Stiffness‐Dependent Photoluminescence of Pendant Fluorophores in Organic Nanoparticles

Fluorescent organic nanoparticles (FoNPs) with backbone stiffness‐dependent photoluminescence were synthesized via microemulsion atom transfer radical polymerization (ATRP) of 2‐(2‐bromoisobutyryloxy)ethyl methacrylate (BiBEM), ethylene glycol dimethacrylate (EGDMA), and methacrylate monomers bearing pendant fluorophores, 1‐pyrenemethyl methacrylate (PyMMA), or 4‐(1,2,2‐triphenylethenyl)benzenemethyl methacrylate (TPEMMA). The crosslinking density precisely tuned the intraparticle rigidity, enabling systematic control over emission mechanisms. Pyrene‐containing FoNPs exhibited a rigidity‐dependent transition from excimer‐dominated to monomer‐dominated fluorescence, whereas TPE‐based FoNPs displayed aggregation‐induced emission (AIE) enhancement as intramolecular motion was restricted. Solvent‐dependent studies revealed that increased polarity and viscosity, particularly in benzyl alcohol and DMSO, promoted cooperative rigidification and enhanced emission intensity through specific polymer–solvent interactions. Furthermore, the retained alkyl bromide chain ends on FoNPs enabled dual roles as initiators and crosslinkers in UV‐induced polymerization of poly(ethylene glycol) acrylates, forming luminescent FoNP–OEG hybrid gels with improved mechanical robustness. This work establishes a versatile platform for integrating tunable optical and mechanical properties into a single polymeric nanoparticle framework, offering new design principles for multifunctional soft materials and providing a platform for future sensing, imaging, and photonic applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Centrality dependence of J /ψ and ψ (2S) production and nuclear modification in p-Pb collisions at $ \sqrt{s_{\mathrm{NN}}} $ = 8 . 16 TeV

The inclusive production of the J/ ψ and ψ (2S) charmonium states is studied as a function of centrality in p-Pb collisions at a centre-of-mass energy per nucleon pair $ \sqrt{s_{\mathrm{NN}}} $ = 8 . 16 TeV at the LHC. The measurement is performed in the dimuon decay channel with the ALICE apparatus in the centre-of-mass rapidity intervals – 4 . 46 < y cms < – 2 . 96 (Pb-going direction) and 2 . 03 < y cms < 3 . 53 (p-going direction), down to zero transverse momentum ( p T ). The J/ ψ and ψ (2S) production cross sections are evaluated as a function of the collision centrality, estimated through the energy deposited in the zero degree calorimeter located in the Pb-going direction. The p T -differential J/ ψ production cross section is measured at backward and forward rapidity for several centrality classes, together with the corresponding average $ \left\langle {p}_{\mathrm{T}}\right\rangle $ and $ \left\langle {p}_{\mathrm{T}}^2\right\rangle $ values. The nuclear effects affecting the production of both charmonium states are studied using the nuclear modification factor. In the p-going direction, a suppression of the production of both charmonium states is observed, which seems to increase from peripheral to central collisions. In the Pb-going direction, however, the centrality dependence is different for the two states: the nuclear modification factor of the J/ ψ increases from below unity in peripheral collisions to above unity in central collisions, while for the ψ (2S) it stays below or consistent with unity for all centralities with no significant centrality dependence. The results are compared with measurements in p-Pb collisions at $ \sqrt{s_{\mathrm{NN}}} $ = 5 . 02 TeV and no significant dependence on the energy of the collision is observed. Finally, the results are compared with theoretical models implementing various nuclear matter effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Non-isometry, state dependence and holography

We establish an equivalence between non-isometry of quantum codes and state dependence of operator reconstruction, and discuss implications of this equivalence for holographic duality. Specifically, we define quantitative measures of non-isometry and state dependence and describe bounds relating these quantities. In the context of holography we show that, assuming known gravitational path integral results for overlaps between semiclassical states, non-isometric bulk-to-boundary maps with a trivial kernel are approximately isometric and bulk reconstruction approximately state-independent. In contrast, non-isometric maps with a non-empty kernel always lead to state-dependent reconstruction. We also show that if a global bulk-to-boundary map is non-isometric, then there exists a region in the bulk which is causally disconnected from the boundary. Finally, we conjecture that, under certain physical assumptions for the definition of the Hilbert space of effective field theory in AdS space, the presence of a global horizon implies a non-isometric global bulk-to-boundary map.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Long- and short-range correlations and their event-scale dependence in high-multiplicity pp collisions at $ \sqrt{s} $ = 13 TeV

Two-particle angular correlations are measured in high-multiplicity proton-proton collisions at √s = 13 TeV by the ALICE Collaboration. The yields of particle pairs at short-(Δη ~ 0) and long-range (1.6 < |Δη| < 1.8) in pseudorapidity are extracted on the near-side (Δ~ 0). They are reported as a function of transverse momentum (p T ) in the range 1 < p T < 4 GeV/c. Furthermore, the event-scale dependence is studied for the first time by requiring the presence of high-p T leading particles or jets for varying p T thresholds. The results demonstrate that the long-range “ridge” yield, possibly related to the collective behavior of the system, is present in events with high-p T processes as well. The magnitudes of the short- and long-range yields are found to grow with the event scale. The results are compared to EPOS LHC and PYTHIA 8 calculations, with and without string-shoving interactions. It is found that while both models describe the qualitative trends in the data, calculations from EPOS LHC show a better quantitative agreement for the p T dependency, while overestimating the event-scale dependency.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗