Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tail”

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 217 records · Page 12

Structural basis for ATG9A recruitment to the ULK1 complex in mitophagy initiation

The assembly of the autophagy initiation machinery nucleates autophagosome biogenesis, including in the PINK1- and Parkin-dependent mitophagy pathway implicated in Parkinson’s disease. The structural interaction between the sole transmembrane autophagy protein, autophagy-related protein 9A (ATG9A), and components of the Unc-51–like autophagy activating kinase (ULK1) complex is one of the major missing links needed to complete a structural map of autophagy initiation. We determined the 2.4-Å x-ray crystallographic structure of the ternary structure of ATG9A carboxyl-terminal tail bound to the ATG13:ATG101 Hop1/Rev7/Mad2 (HORMA) dimer, which is part of the ULK1 complex. We term the interacting portion of the extreme carboxyl-terminal part of the ATG9A tail the “HORMA dimer–interacting region” (HDIR). This structure shows that the HDIR binds to the HORMA domain of ATG101 by β sheet complementation such that the ATG9A tail resides in a deep cleft at the ATG13:ATG101 interface. Disruption of this complex in cells impairs damage-induced PINK1/Parkin mitophagy mediated by the cargo receptor NDP52.

59 BASIC BIOLOGICAL SCIENCES↗

m2p (Monomers to Polymers)

The software is Python-based code that takes monomers, represented as SMILES, and "reacts" them together to form long polymer chains, represented as SMILES, based on known polymer chemistries. The software is used in a machine-learning pipeline where a database and machine-learning approaches are used to predict polymer properties base on molecular structure. This software is necessary to represent the polymer chain, instead of simple repeat units, as the software statistically represents different molecular architectures that occur during reactions (e.g. head-to-head, head-to-tail, and tail-to-tail). It is capable of building polymer chains with five different chemistries: vinyl/olefin, polyamide, polycarbonate, polyester, and polyimides. The user is able to specify different monomer pairings and degree of polymerization for different polymer chains.

Wilson, Nolan↗

HL-LHC Beam Dynamics with Hollow Electron Lenses

Each of the two proton beams in the High-Luminosity Large Hadron Collider (HL-LHC) will carry a total energy of 720 MJ. One concern for machine protection is the energy stored in the transverse beam tails, estimated to potentially reach up to 5% of the total stored energy. Several failure scenarios could drive these tails into the collimators, potentially causing damage and therefore severely affecting operational efficiency. Hollow Electron Lenses (HEL) were integrated in the HL-LHC baseline to mitigate this risk by depleting the tails in a controlled way. A hollow-shaped electron beam runs co-axially to the hadron beam over about 3 m, such that halo particles at large amplitudes become unstable, while core particles ideally remain undisturbed. Residual fields from e-beam asymmetries can, however, induce emittance growth of the beam core. Various options for the pulsing of the HEL are considered and are compared using two figures of merit: halo depletion efficiency and core emittance growth. This contribution presents simulations for these two effects with different HEL pulsing modes using the final HL-LHC optics, that was optimized at the location of the lenses.

43 PARTICLE ACCELERATORS↗

Detection and Perception of Sound by Eagles and Surrogate Raptors

One overarching objective of this program of study was the accumulation of objective, scientifically valid information relating to auditory performance of bald and golden eagles that may be used to guide the development of acoustic alerting/deterrence technologies intended to discourage encroachment into wind energy air spaces. To that end, analyses aimed at the characterization of sensitivity to sound in bald and golden eagles, along with findings in the supra-threshold, dynamic frequency spaces related to response latencies and amplitudes, leads us to conclude that bald, and golden eagles navigate the same basic working auditory space, as in other known and thus far characterized members of the diurnal raptor family. Specifically, bald and golden eagles, along with other raptor species within the group, operate in an auditory space characterized by a frequency band at least four octaves wide and centered on 2 kHz, with an upper frequency limit between 6 and 10 kHz at 80 dB SPL and a lower frequency limit that almost certainly extends below 0.2 kHz. Consequently, we recommend that signal designers use these data as a guideline in efforts to design effective and efficient acoustic alerting/deterrent systems. It is important to note that signal energy broadcast outside of this frequency band at moderate levels will not contribute to the efficacy of a deterrent but will add an unnecessary fraction to the overall acoustic pollution budget. The importance of this consideration is heightened by contemporaneous concerns related to the transmission of noise broadcast by wind energy farms. In addition, based on analyses of data acquired from red-tailed hawks using the same experimental paradigm and data acquisition system, we conclude that auditory function in the red-tailed hawk is sufficiently like that observed in bald and golden eagles to permit its use as a surrogate species. Response waveforms, threshold-frequency curves, and input-output characteristics match those of eagles closely. It should be noted however, that differences in sensitivity and slightly extended high-frequency limits of hearing should be taken into account when extrapolating findings from one species to the others. Although the inclusion of behavioral tests of red-tailed hawks to acoustic stimuli was beyond the scope of this investigation, future efforts to assess response parameters like signal-type preference and habituation rate will further elucidate their suitability to serve as eagle surrogates in behavioral studies; nonetheless, the species in question are well matched with respect to basic auditory performance. A second essential objective of this program of study was the acoustic characterization of a subset of calls comprising the vocal repertoires of bald and golden eagles that may be used to supplement auditory performance findings in the effort to guide the development of acoustic alerting signals. With regard to that objective, the vocal repertoires of both bald and golden eagle species are rich and varied. While similar in spectrographic structure, distinctive differences are also clear. Generally, golden eagles produce some calls with shorter durations, and similar “sounding” calls exhibit distinctively different spectrographic patterns than those of bald eagles. Both species produce calls that contain a wide variety of nonlinear elements that operate to enhance the rich and varied nature of commonly observed vocal products. Comparison of the average power spectra of commonly observed bald and golden eagle calls with threshold-frequency curves leads to the conclusion that call energies fall within the frequency bounds of hearing. Further, the acoustic energy of calls considered in this report tend to fall into overlapping, but different frequency ranges of the acoustic sensitivity curve. This condition may encourage signal designers to vary the frequency content of acoustic deterrence signals in the field. Finally, preliminary observations relating to the tendencies and proclivities of bald eagles to attend to the acoustic landscape lead to the conclusion that eagles monitor their immediate sound environment assiduously. Individuals respond to a variety of natural and synthetic sound signals reliably and, perhaps most relevant in the context of the engineering of acoustic alerting/deterrence technologies, habituation to most sounds considered in this effort was minimal. These preliminary results, while calling for extended behavioral testing, are promising and set the stage for the exportation of behavioral studies into real world scenarios.

17 WIND ENERGY↗

Risk-Adaptive Experimental Design for High-Consequence Systems: LDRD Final Report

Constructing accurate statistical models of critical system responses typically requires an enormous amount of data from physical experiments or numerical simulations. Unfortunately, data generation is often expensive and time consuming. To streamline the data generation process, optimal experimental design determines the 'best' allocation of experiments with respect to a criterion that measures the ability to estimate some important aspect of an assumed statistical model. While optimal design has a vast literature, few researchers have developed design paradigms targeting tail statistics, such as quantiles. In this project, we tailored and extended traditional design paradigms to target distribution tails. Our approach included (i) the development of new optimality criteria to shape the distribution of prediction variances, (ii) the development of novel risk-adapted surrogate models that provably overestimate certain statistics including the probability of exceeding a threshold, and (iii) the asymptotic analysis of regression approaches that target tail statistics such as superquantile regression. To accompany our theoretical contributions, we released implementations of our methods for surrogate modeling and design of experiments in two complementary open source software packages, the ROL/OED Toolkit and PyApprox.

97 MATHEMATICS AND COMPUTING↗

Energy-dependent optimization of the prompt fission neutron spectrum with CGMF

Throughout the course of FY21, significant effort was put into investigating models within the LANL developed Hauser-Feshbach fission fragment decay code, CGMF, to understand and potentially solve the long-standing challenge of a too-soft prompt fission neutron spectrum, PFNS. Several inputs and models to CGMF were investigated, including the discrete nuclear levels, the optical model potential, level densities, and the fission fragment initial conditions. Some of the global models within CGMF led to a slight hardening of the neutron spectrum—particularly the likely incomplete discrete levels in through which γ-rays decay—but none of the changes where large enough for the tail of the PFNS to reproduce experimental data. A significant hardening of the spectrum tail was observed when the fission fragment initial conditions were optimized based on their sensitivities to the PFNS data for thermal incident neutrons. In this way, the parameters for the CGMF mass and total kinetic energy distributions, along with the spin cutoff factor were adjusted to better reproduce the experimental PFNS measurements. This optimization hardened the tail of the PFNS slightly but led to unphysical mass distributions for the fission fragments before neutron emission. It was clear from the above that we do not expect to be able to produce an evaluation-quality PFNS with CGMF in the near future. Challenges at thermal will persist–and possibly worsen–with increasing incident energy, where more models are needed to completely describe the fission. Basic-science research funding exceeding the amount available and scope of our NCSP funds would be needed to tackle this decade-long challenge impacting many fission-fragment event generator. And, in fact, Amy Lovell won LDRD ECR funding to do so over the next few years. Therefore, we focused in FY22 on extending evaluation capabilities beyond thermal incident neutrons, to take into account the incident energy dependence of the PFNS and fission fragment initial condition distributions in CGMF. We chose to set up the evaluation methodology to perform PFNS evaluations with CGMF across incident-neutron energies, in order to have it readily available for future NCSP evaluations when the PFNS from CGMF has improved. In this report, we outline the evaluation methodology, along with the results of the optimization, including full model calculations with CGMF using the evaluated parameters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

VPS26 Moonlights as an Arrestin-like Adapter for a 7-transmembrane RGS 2 protein in Arabidopsis thaliana

Extracellular signals perceived by 7-transmembrane (7TM)-spanning receptors on the plasma membrane utilize a cytoplasmic adaptor that propagates signaling and initiates a feedback circuit for desensitization by removal of these receptors at the plasma membrane leading to desensitization to that signal. Signal (agonist) binding often evokes phosphorylation at the C-terminal tail of many 7TM G-protein-coupled receptors in animal cells which then recruits a cytoplasmic intermediate adaptor, β-arrestin, that sets in motion clathrin-mediated endocytosis (CME). Some 7TM receptors have agonist-induced phosphorylation and CME that do not involve β-arrestin, therefore it is unclear how these phosphorylated 7TM receptors are internalized. Arrestins, neither α– or β-types, are not encoded in the Arabidopsis genome, yet Arabidopsis cells have a well characterized signal-induced CME of a 7TM protein, Regulator of G Signaling 1 (AtRGS1). We show that a component of the retromer complex, VPS26, moonlights as an arrestin-like adapter. VPS26A and VPS26B form heterodimer instead of homodimer or monomer. flg22, a bacterium-derived signal that alerts plant cells of the presence of a potential pathogen, induces phosphorylation of AtRGS1 within a cluster of serines in its C-terminal tail, dimerization, and its CME. The phosphorylated C-terminal tail peptide binds to VPS26A/B heterodimer with greater affinity than the unphosphorylated peptide. We propose that VPS26 serves as an arrestin-like adaptor in the CME of AtRGS1 dimerization biased for VPS26 signaling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterization and effects of impurities on carbonate quantification in heterogenous matrices

Mineral carbonation simulates a natural weathering phenomena by breaking down silicates and oxides to form Ca & Mg carbonates. Various mineralization methods have been demonstrated as a potential technique to improve the quality of slags and tailings through neutralization and stabilization of problematic species to yield a product better suited for use in concrete. This study aims to characterize, quantify and analyze carbonates in various carbonated products such as mineralized CaCO 3 , CO 2 mineralized Steel Slags and Mine Tailings. More than 10 samples were analyzed for carbonate measurement and verification from industrial and academic partners that pioneer commercial CO 2 mineralization technologies. The samples were characterized primarily by using X-ray diffraction (XRD), Thermogravimetric Analyses (TGA), and Scanning Electron Microscopy (SEM) to gain insights into CaCO 3 content. A baseline characterization of lab-grade CaCO 3 and MgCO 3 also revealed important considerations for CaCO 3 measurement using TGA alone. The experiments using synthetic lab-grade samples also revealed that the presence of MgCO 3 /MgO can accelerate the decomposition of CaCO 3 and thus can affect measurement parameters. Lab-grade CaCO 3 samples dosed into steel slag and mine tailing also showed significant deviation in their decomposition behavior. These insights are used to inform the development of a standardized protocol for the measurement and verification of carbonate-bearing products.

97 MATHEMATICS AND COMPUTING↗

New Determination of Fundamental Properties of Palomar 5 Using Deep DESI Imaging Data

The legacy imaging surveys for the Dark Energy Spectroscopic Instrument project provides multiple-color photometric data, which are about 2 mag deeper than those from the SDSS. In this study, we redetermine the fundamental properties for an old halo globular cluster of Palomar 5 based on these new imaging data, including structure parameters, stellar population parameters, and luminosity and mass functions. These characteristics, together with its tidal tails, are key for dynamical studies of the cluster and constraining the mass model of the Milky Way. By fitting the King model to the radial surface density profile of Palomar 5, we derive the core radius of [[Formula]], tidal radius of [[Formula]], and concentration parameter of c = 0.78 ± 0.04. We apply a Bayesian analysis method to derive the stellar population properties and get an age of 11.508 ± 0.027 Gyr, metallicity of [Fe/H] = −1.798 ± 0.014, reddening of E(B − V) = 0.0552 ± 0.0005, and distance modulus of (m−M){sub 0} = 16.835±0.006. The main-sequence luminosity and mass functions for both the cluster center and tidal tails are investigated. The luminosity and mass functions at different distances from the cluster center suggest that there is obvious spatial mass segregation. Many faint low-mass stars have been evaporated at the cluster center, and the tidal tails are enhanced by low-mass stars. Both the concentration and relaxation times suggest that Palomar 5 is a totally relaxed system.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electron Energy Partition across Interplanetary Shocks. III. Analysis

An analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within ±2 hr of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 au, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partialto-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF, but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions), while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature; even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗

X-Ray Afterglows from the Gamma-Ray Burst “Large-angle” Emission

In this work, we derive basic analytical results for the timing and decay of the gamma-ray burst (GRB) counterpart and delayed afterglow light curves for a brief emission episode from a relativistic surface endowed with angular structure, consisting of a uniform core of size ${\theta }_{c}$ (Lorentz factor ${{\rm{\Gamma }}}_{c}$ and surface emissivity ${i}_{\nu ^{\prime} }^{{\prime} }$ are angle independent) and an axially symmetric power-law envelope (${\rm{\Gamma }}\sim {\theta }^{-g}$). In this "large-angle emission" model, radiation produced during the prompt emission phase (GRB) at angles $\theta \gt {\theta }_{c}$ arrives at the observer well after the burst (delayed emission). The dynamical time range of the very fast decaying GRB "tail" and of the flat afterglow "plateau" and the morphology of the GRB counterpart/afterglow are all determined by two parameters: the core's parameter ${{\rm{\Gamma }}}_{c}{\theta }_{c}$ and the envelope's Lorentz factor index g, leading to three types of light curves that display three post-GRB phases (type 1: tail, plateau/slow decay, post-plateau/normal decay), two post-GRB phases (type 2: tail and fast decay), or just one (type 3: normal decay). We show how X-ray light-curve features can be used to determine core and envelope dynamical and spectral parameters. Testing of the large-angle emission model is done using the Swift/XRT X-ray emission of two afterglows of type 1 (GRB 060607A, GRB 061121), one of type 2 (GRB 061110A), and one of type 3 (GRB 061007). We find that the X-ray afterglows with plateaus require an envelope Lorentz factor ${\rm{\Gamma }}\sim {\theta }^{-2}$ and a comoving-frame emissivity ${i}_{\nu ^{\prime} }^{{\prime} }\sim {\theta }^{2}$; thus, for a typical afterglow spectrum ${F}_{\nu }\sim {\nu }^{-1}$, the lab-frame energy release is uniform over the emitting surface.

79 ASTRONOMY AND ASTROPHYSICS↗

Substructure at High Speed. I. Inferring the Escape Velocity in the Presence of Kinematic Substructure

Abstract The local escape velocity provides valuable inputs to the mass profile of the galaxy, and requires understanding the tail of the stellar speed distribution. Following Leonard & Tremaine, various works have since modeled the tail of the stellar speed distribution as ∝ ( v esc − v ) k , where v esc is the escape velocity, and k is the slope of the distribution. In such studies, however, these two parameters were found to be largely degenerate and often a narrow prior is imposed on k in order to constrain v esc . Furthermore, the validity of the power-law form can breakdown in the presence of multiple kinematic substructures or other mis-modeled features in the data. In this paper, we introduce a strategy that for the first time takes into account the presence of kinematic substructure. We model the tail of the velocity distribution as a sum of multiple power laws as a way of introducing a more flexible fitting framework. Using mock data and data from FIRE simulations of Milky Way-like galaxies, we show the robustness of this method in the presence of kinematic structure that is similar to the recently discovered Gaia Sausage. In a companion paper, we present the new measurement of the escape velocity and subsequently the mass of the Milky Way using Gaia eDR3 data.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurements of a 2.1 MeV H$^-$ beam with an Allison scanner

Transverse 2D phase space distribution of a 2.1 MeV, 5 mA H$^-$ beam is measured at the PIPIT test accelerator at Fermilab with an Allison scanner. The paper describes the design, calibration, and performance of the scanner as well as the main results of the beam measurements. Analyses of the recorded phase portraits are performed primarily in action-phase coordinates; the stability of the action under linear optics makes it easier to compare measurements taken with different beamline conditions, e.g. in various locations. The intensity of a single measured point (\pixel") is proportional to the phase density in the corresponding portion of the beam. When the Twiss parameters are calculated using only the high-phase density part of the beam, the pixel intensity in the beam core is found to be decreasing exponentially with action and to be phase-independent. Outside of the core, the intensities decrease with action at a significantly slower rate than in the core. This `tail' comprises 10-30% of the beam, with 0.1% of the total measured intensity extending beyond the action 10-20 times larger than the rms emittance. The transition from the core to the tail is accompanied by the appearance of a strong phase dependence, which is characterized in action-phase coordinates by two `branches' extending beyond the core. A set of selected measurements shows, in part, that there is no measurable emittance dilution along the beam line in the main portion of the beam; the beam parameters are practically constant over a 0.5 ms pulse; and scraping in various parts of the beam line is an effective way to decrease the transverse tails by removing the branches.

43 PARTICLE ACCELERATORS↗

Novel deformation function creating or destroying any number of even kink solutions

Here, we present a one-parameter family of deformation functions f(Φ) which have novel properties. Firstly, the deformation function is its own inverse. We show that a class of potentials remains invariant under this deformation. Further, when applied to a certain class of kink bearing potentials, one obtains potentials which are not bounded from below. Besides, we show that there is a wide class of potentials having power law kink tails which are connected by this deformation but the corresponding kink tails of the two potentials have different asymptotic behavior. Finally, we show that when this deformation function is applied to an appropriate one-parameter family of potentials having two kink solutions, it creates new potentials with an arbitrary even number (2n) of kink solutions. Conversely, by starting from this potential with 2n kink solutions the deformation function annihilates 2n-2 kinks and we get back the potential which we started with having only two kink solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effects of Dead‐End Fractures on Non‐Fickian Transport in Three‐Dimensional Discrete Fracture Networks

Abstract Understanding mechanistic causes of non‐Fickian transport in fractured media is important for many hydrogeologic processes and subsurface applications. This study elucidates the effects of dead‐end fractures on non‐Fickian transport in three‐dimensional (3D) fracture networks. Although dead‐end fractures have been identified as low‐velocity regions that could delay solute transport, the direct relation between dead‐end fractures and non‐Fickian transport has been elusive. We systematically generate a large number of 3D discrete fracture networks with different fracture length distributions and fracture densities. We then identify dead‐end fractures using a novel graph‐based method. The effect of dead‐end fractures on solute residence time maximizes at the critical fracture density of the percolation threshold, leading to strong late‐time tailing. As fracture density increases beyond the percolation threshold, the network connectivity increases, and dead‐end fractures diminish. Consequently, the increase in network connectivity leads to a reduction in the degree of late‐time tailing. We also show that dead‐end fractures can inform about main transport paths, such as the mean tortuosity of particle trajectories. This study advances our mechanistic understanding of solute transport in 3D fracture networks.

58 GEOSCIENCES↗

Quantifying Distribution System Resilience From Utility Data: Large Event Risk and Benefits of Investments

We focus on blackouts in electric distribution systems that have a large cost to customers. To quantify resilience to these events, we show how to calculate risk metrics from the historical outage data routinely collected by utilities' outage management systems. Risk is defined using a customer cost exceedance curve. The exceedance curve has a heavy tail that implies large fluctuations in large blackout costs, and this makes estimating the mean large cost in the usual way impractical. To avoid this problem, we use new resilience metrics describing the large event risk; these metrics are the probability of a large cost event, the annual log cost resilience index, and the average of the logarithm of the cost of large-cost events or the slope magnitude of the tail on a log–log exceedance curve. Resilience can be improved by planned investments to upgrade system components or speed up restoration. The benefits that these investments would have had if they had been made in the past can be quantified by “rerunning history” with the effects of the investment included, and then recalculating the large event risk to find the improvement in resilience. An example using utility data shows a 2% reduction in the probability of a large cost event due to 10% wind hardening and 6%–7% reduction due to 10% faster restoration in two different areas of a distribution utility. This new data-driven approach to quantify resilience and resilience investments is realistic and much easier to apply than complicated approaches based on modeling all the phases of resilience. Moreover, an appeal to improvements to past lived experience may well be persuasive to customers and regulators in making the case for resilience investments.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Stochastic Transport Model for the Cumulative Number of Fissions and Deposited Fission Energy

The stochastic theory of neutron transport is extended to describe the cumulative distribution of fission numbers and deposited fission energy in a multiplying assembly. Solutions for the probability distributions are obtained using analytical approximations and Monte Carlo simulation in lumped geometry and in symmetric homogeneous and heterogeneous spheres. The results show the development of a power-law tail in the steady state fission number and deposited energy distributions when the medium is critical, independent of the fission neutron multiplicity distribution and domain heterogeneity. In contrast, the asymptotic decay is faster than exponential in subcritical media due to rapid chain extinction and in supercritical media due to the increasing probability of chain divergence. Here, a formal asymptotic analysis of the problem in lumped geometry with an arbitrary fission neutron multiplicity confirms the existence of power-law tails at critical.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Explicit kink solutions in several one-parameter families of higher-order field theory models

Here, we present several one-parameter families of higher-order field theory models some of which admit explicit kink solutions with an exponential tail while others admit explicit kink solutions with a power-law tail. Various properties of these families of kink solutions are examined in detail. Further, by applying the recent Manton formalism, we provide estimates for the kink–kink (KK) and antikink–kink (AK-K) acceleration and hence the ratio of the corresponding AK-K and KK forces.

97 MATHEMATICS AND COMPUTING↗