Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surprisal”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Zonal flow screening in negative triangularity tokamaks

Here, this paper presents a comparative study of zonal flow screening in positive and negative triangularity tokamaks. Neoclassical screening sets the strength of zonal flow shear for fixed drive and damping. Orbit calculations show that the radial excursions of trapped particles are reduced in negative triangularity tokamaks, as compared to positive triangularity. Yet surprisingly, the neoclassical dielectric susceptibility actually increases with decreasing triangularity, such that the susceptibility is higher for negative triangularity than for positive triangularity. This is because the reduction in neoclassical polarization by shrinking the banana width is offset by the increase in neoclassical polarization by the enhancement of trapped fraction for negative triangularity. As a result, the zonal flow screening length is actually enhanced for negative triangularity, as compared to positive triangularity. Hence, the zonal flow residual is smaller for negative triangularity than for positive triangularity. Results from gyrokinetic simulations support the analytic calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A reduced model for the ITER divertor based on SOLPS solutions for ITER Q = 10 baseline conditions: A. identifying options for the control parameters*

Edge codes such as SOLPS coupled to neutral codes such as EIRENE have become so comprehensive and sophisticated that they now constitute, in effect, 'code-experiments' that, as for actual experiments, can benefit from interpretation using simple models and conceptual frameworks, i.e. reduced models. The first task is the identification of options for the reduced model control parameters that are best suited for control of the action of the divertor, i.e. for control of target power loading and sputter-erosion, primarily. A strong correlation between the electron temperature at the divertor target, T e,t , and the neutral deuterium D 2 density at the target, n D2,t , flux-tube resolved, has recently been reported for a number of code studies including SOLPS-4.3 modeling of a set of ~50 ITER baseline cases: Q DT = 10, q 95 = 3, P SOL = 100 MW, metallic walls, and Ne seeding. This part A of the present study reports new results for largely the same ITER cases, confirming the strong correlations reported earlier between local values of T e,t , and (i) n D2,t , and (ii) normalized volumetric losses of power and pressure in the divertor. Strong correlations have now also been found, and are reported here for the first time, between T e,t and all of the divertor target quantities of practical interest. A physical explanation for this surprising result has not been fully identified; nevertheless it has encouraging implications for reduced modeling of the ITER divertor. For such ITER conditions, (i) the global Ne injection rate, Inj Ne [Ne/s], and (ii) the electron temperature at the location on the target where the peak power deposition occurs, T e,t @q ⊥,pk (eV), are found to be promising reduced model control parameters. In the companion report, part B, a reduced model for the ITER divertor is developed and described in detail, based on reversed-direction 2 point modelling, Rev2PM. The input to the reduced model is a value of the variable pair (T et @q ⟂,pk ,Inj Ne ) and the output are values of the various target as well as divertor-entrance quantities of practical interest, e.g. q ⊥,pk , n e,Xpt (the electron density at the poloidal location of the X-point), etc. In part B the reduced model is quantitatively characterized using one half of the code cases; it is then used to successfully predict (replicate) the code values of e.g. n e,Xpt for the other half of the cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A reduced model for the ITER divertor based on SOLPS solutions for ITER Q = 10 baseline conditions: B. A reduced model based on reversed-direction two point modeling

Edge codes such as SOLPS coupled to neutral codes such as EIRENE have become so comprehensive and sophisticated that they now constitute, in effect, 'code-experiments' that, as for actual experiments, can benefit from interpretation using simple models and conceptual frameworks, i.e. reduced models. The first task is the identification of options for the reduced model control parameters that are best suited for control of the action of the divertor, i.e. for control of target power loading and sputter–erosion, primarily. A strong correlation between the electron temperature at the divertor target, T e,t , and the neutral deuterium D 2 density at the target, n D2,t , flux-tube resolved, has recently been reported for a number of code studies including SOLPS-4.3 modeling of a set of ~50 ITER baseline cases: Q DT = 10, q 95 = 3, P SOL = 100 MW, metallic walls, and Ne seeding (Pitts et al 2019). Part A of the present study reports new results for largely the same ITER cases, confirming the strong correlation reported earlier between local values of T e,t , and (i) n D2,t , and (ii) normalized volumetric losses of power and pressure in the divertor. Strong correlations have now also been found, and are reported here for the first time, between T e,t and all of the divertor target quantities of practical interest. A physical explanation for this surprising result has not as yet been fully identified; nevertheless it has encouraging implications for reduced modeling of the ITER divertor. For such ITER conditions, (i) the global Ne injection rate, Inj Ne (Ne s –1 ), and (ii) the electron temperature at the location on the target where the peak power deposition occurs, T e,t @q ⊥,pk (eV), are found to be promising reduced model control parameters. In this part B, a reduced model for the ITER divertor is developed and described in detail, based on reversed-direction two point modeling, Rev2PM. The input to the reduced model is a value of the variable pair (T e,t @q ⊥,pk , Inj Ne )for a chosen case and the output are values of the various target as well as divertor-entrance quantities of practical interest, e.g. q ⊥,pk , the electron density at the X-point, n e,Xpt , etc. The reduced model was quantitatively characterized using one half of the code cases; it was then used to successfully predict (replicate) the code values of e.g. n e,Xpt for the other half.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Staged cooling of a fusion-grade plasma in a tokamak thermal quench

Abstract In tokamak disruptions where the magnetic connection length becomes comparable to or even shorter than the plasma mean-free-path, parallel transport can dominate the energy loss and the thermal quench of the core plasma goes through four phases (stages) that have distinct temperature ranges and durations. The main temperature drop occurs while the core plasma remains nearly collisionless, with the parallel electron temperature T e ∥ dropping in time t as T e ∥ ∝ t − 2 and a cooling time that scales with the ion sound wave transit time over the length of the open magnetic field line. These surprising physics scalings are the result of effective suppression of parallel electron thermal conduction in an otherwise bounded, quasineutral, and collisionless plasma, which is different from what are known to date on electron thermal conduction along the magnetic field in a nearly collisionless and quasineural plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Plasma performance and operational space with an RMP-ELM suppressed edge

Abstract The operational space and global performance of plasmas with edge-localized modes (ELMs) suppressed by resonant magnetic perturbations (RMPs) are surveyed by comparing AUG, DIII-D, EAST, and KSTAR stationary operating points. RMP-ELM suppression is achieved over a range of plasma currents, toroidal fields, and RMP toroidal mode numbers. Consistent operational windows in edge safety factor are found across devices, while windows in plasma shaping parameters are distinct. Accessed pedestal parameters reveal a quantitatively similar pedestal-top density limit for RMP-ELM suppression in all devices of just over 3 × 10 19 m −3 . This is surprising given the wide variance of many engineering parameters and edge collisionalities, and poses a challenge to extrapolation of the regime. Wide ranges in input power, confinement time, and stored energy are observed, with the achieved triple product found to scale like the product of current, field, and radius. Observed energy confinement scaling with engineering parameters for RMP-ELM suppressed plasmas are presented and compared with expectations from established H and L-mode scalings, including treatment of uncertainty analysis. Different scaling exponents for individual engineering parameters are found as compared to the established scalings. However, extrapolation to next-step tokamaks ITER and SPARC find overall consistency within uncertainties with the established scalings, finding no obvious performance penalty when extrapolating from the assembled multi-device RMP-ELM suppressed database. Overall this work identifies common physics for RMP-ELM suppression and highlights the need to pursue this no-ELM regime at higher magnetic field and different plasma physical size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simultaneous access to high normalized density, current, pressure, and confinement in strongly-shaped diverted negative triangularity plasmas

Strongly-shaped diverted negative triangularity (NT) plasmas in the DIII-D tokamak demonstrate simultaneous access to high normalized density, current, pressure, and confinement. NT plasmas are shown to exist across an expansive parameter space compatible with high fusion power production, revealing surprisingly good core stability properties that compare favorably to conventional positive triangularity plasmas in DIII-D. Non-dimensionalizing the key parameters, expanded operating spaces featuring edge safety factors below 3, normalized betas above 3, Greenwald density fractions above 1, and high-confinement mode (H-mode) confinement qualities above 1 are observed, even simultaneously, and all with a robustly stable edge free from deleterious edge-localized mode instabilities. Scaling of the confinement time with engineering parameters reveals at least a linear dependence on plasma current although with significant power degradation, both in excess of expected H-mode scalings. These results increase confidence that NT plasmas are a viable approach to realize fusion power and open directions for future detailed study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Understanding the L-H isotope effect at the DIII-D tokamak and advancements in synthetic turbulence diagnostics

Abstract It is determined that while heat flux differences between hydrogen and deuterium isotope experiments result from natural differences in carbon impurity content at DIII-D, it is not the origin of the low to high confinement mode (L-H) transition isotope effect. More specifically, a two times larger edge radial electric field in hydrogen compared to deuterium is uncovered and believed to play an important role. The origin of this radial electric field difference is determined to have two possible origins: differences in poloidal rotation and turbulent Reynolds stress in the closed field line region, and increased outer strike point temperatures and space potentials on open field lines. Experimental observations from both profile and turbulence diagnostics are supported by nonlinear gyrokinetic simulations using the code CGYRO. Simulations illustrated heat transport isotope effects in the plasma edge and shear layer resulting from differences in impurity content, electron non-adiabaticity, and main ion mass dependent E × B shear stabilization. Turbulence prediction comparisons from flux-matched CGYRO simulations to experimental measurements including electron temperature, density and velocity fluctuations are found to be in good agreement with available data. A dedicated DIII-D experiment in hydrogen was performed to seed more carbon than naturally occurring, to match deuterium experiments, and possibly reduce the L-H power threshold based on gyro-kinetic predictions. To our surprise, while ion temperature gradient (ITG) turbulence was stabilized, nodiscernible change in L-H power threshold were observed in these special hydrogen experiments. In particular, it is noticed that the edge radial electric field and Reynolds stress were observed as nearly unchanging in the presence of ITG stabilization. These experimental data have enabled a more comprehensive picture of the multitude of isotope effects at play in fusion experiments, and the important potential connection between the confined and unconfined plasma regions in regulating L-H transition dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simultaneous enhancement of tritium burn efficiency and fusion power with low-tritium spin-polarized fuel

This study demonstrates that using spin-polarized deuterium-tritium (D-T) fuel with more deuterium than tritium can increase tritium burn efficiency (TBE) by at least an order of magnitude without compromising fusion power output, compared to unpolarized fuel. Although previous studies show that a low tritium fraction can enhance TBE, this strategy resulted in reduced fusion power density. The surprising improvement in TBE at fixed power reported here is due to the TBE increasing nonlinearly with decreasing tritium fraction but the fusion power density increasing roughly linearly with D-T cross section. A study is performed for an ARC-like tokamak producing 481 MW of fusion power with unpolarized 53:47 D-T fuel, finding the minimum startup tritium inventory (I startup,min ) is 0.69 kg. By spin-polarizing half of the fuel and using a 60:40 D-T mix, I startup,min is reduced to 0.08 kg, and fully spin-polarizing the fuel with a 63:37 D-T mix further reduces I startup,min to 0.03 kg. Some ARC-like scenarios are predicted to achieve plasma ignition with relatively modest spin polarization. These findings indicate that, with advancements in helium divertor pumping efficiency, TBE values of approximately 10%–40% could be achieved using low-tritium-fraction and spin-polarized fuel with minimal power loss. This would dramatically lower tritium startup inventory requirements and reduce the amount of on-site tritium. More generally than just for spin-polarized fuels, increased plasma performance can be used to increase TBE. This strongly motivates the development of spin-polarized fuels and low-tritium-fraction operation for burning plasmas.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Advancing Tritium Self-Sufficiency in Fusion Power Plants: Insights from the BABY Experiment

Abstract In the pursuit of fusion power, achieving tritium self-sufficiency stands as a pivotal challenge.
Tritium breeding within molten salts is a critical aspect of next-generation fusion reactors, yet experimental measurements of \gls{tbr} have remained elusive.
Here we present the results of the \gls{baby} experiment, which represents a pioneering effort in tritium research by utilizing high-energy (\SI{14}{\mega\electronvolt}) neutron irradiation of molten salts, a departure from conventional low-energy neutron approaches. 
Using a small-scale (\SI{100}{\milli\litre}) molten salt tritium breeding setup, we not only simulated, but also directly measured a \gls{tbr} \textcolor{red}{(\num{3.57e-4})}. 
This innovative approach provides crucial experimental validation, offering insights unattainable through simulation alone.
Moreover, our findings reveal a surprising outcome: tritium was predominantly collected as HT, contrary to the expected TF. 
This underscores the complexity of tritium behavior in molten salts, highlighting the need for further investigation.
This work lays the foundation for a more sophisticated experimental setup, including increasing the volume of the breeder, enhancing neutron detection, and refining tritium collection systems.
Such improvements are crucial for advancing our understanding of fusion reactor feasibility and paving the way for future experiments.

Delaporte-Mathurin, Rémi (ORCID:0000000310648882)↗
DOE OSTI2026

Persistence of deuterium and tritium nuclear spin-polarization in presence of high-frequency plasma waves

We present first-principles numerical calculations of the depolarization rate of spin-polarized deuterium and tritium nuclei in realistic tokamak plasmas, driven by resonant interactions with plasma waves. Backed up by first-of-a-kind linear and nonlinear simulations, we find that alpha particle-driven Alfvénic modes cause only negligible depolarization, which is contrary to expectations in prior literature. Other Alfvénic instabilities can in principle degrade polarization, but only under conditions unlikely to be realized on transport timescales. By combining full-orbit particle tracing with a dedicated depolarization solver, we demonstrate that wave-driven depolarization is surprisingly weak in SPARC and ITER-scale devices. These results provide strong evidence that spin-polarized fuel can maintain its polarization long enough to boost fusion reactivity, opening a viable path toward substantially enhanced performance in magnetic confinement fusion power plants.

Cook, James William Setchfield [UK Atomic Energy A↗
DOE OSTI2022

Design of detectors at the electron ion collider with artificial intelligence

Abstract Artificial Intelligence (AI) for design is a relatively new but active area of research across many disciplines. Surprisingly when it comes to designing detectors with AI this is an area at its infancy. The electron ion collider is the ultimate machine to study the strong force. The EIC is a large-scale experiment with an integrated detector that extends for about ±35 meters to include the central, far-forward, and far-backward regions. The design of the central detector is made by multiple sub-detectors, each in principle characterized by a multidimensional design space and multiple design criteria also called objectives. Simulations with Geant4 are typically compute intensive, and the optimization of the detector design may include non-differentiable terms as well as noisy objectives. In this context, AI can offer state of the art solutions to solve complex combinatorial problems in an efficient way. In particular, one of the proto-collaborations, ECCE, has explored during the detector proposal the possibility of using multi-objective optimization to design the tracking system of the EIC detector. This document provides an overview of these techniques and recent progress made during the EIC detector proposal. Future high energy nuclear physics experiments can leverage AI-based strategies to design more efficient detectors by optimizing their performance driven by physics criteria and minimizing costs for their realization.

Instruments & Instrumentation↗
DOE OSTI2021

Comment on ‘Carbon Intensity of corn ethanol in the United States: state of the science’

Abstract In their recent contribution, Scully et al (2021 Environ. Res. Lett. 16 043001) review and revise past life cycle assessments of corn-grain ethanol’s carbon (C) intensity to suggest that a current ‘central best estimate’ is considerably less than all prior estimates. Their conclusion emerges from selection and recombination of sector-specific greenhouse gas emission predictions from disparate studies in a way that disproportionately favors small values and optimistic assumptions without rigorous justification nor empirical support. Their revisions most profoundly reduce predicted land use change (LUC) emissions, for which they propose a central estimate that is roughly half the smallest comparable value they review (figure 1). This LUC estimate represents the midpoint of (a) values retained after filtering the predictions of past studies based on a set of unfounded criteria; and (b) a new estimate they generate for domestic (i.e. U.S.) LUC emissions. The filter the authors apply endorses a singular means of LUC assessment which they assert as the ‘best practice’ despite a recent unacknowledged review (Malins et al 2020 J. Clean. Prod. 258 120716) that shows this method almost certainly underestimates LUC. Moreover, their domestic C intensity estimate surprisingly suggests that cropland expansion newly sequesters soil C, counter to ecological theory and empirical evidence. These issues, among others, prove to grossly underestimate the C intensity of corn-grain ethanol and mischaracterize the state of our science at the risk of perversely affecting policy outcomes.

Spawn-Lee, Seth A. (ORCID:0000000188215345)↗
DOE OSTI2022

Rapidly changing high-latitude seasonality: implications for the 21st century carbon cycle in Alaska

Seasonal variations in high-latitude terrestrial carbon (C) fluxes are predominantly driven by air temperature and radiation. At present, high-latitude net C uptake is largest during the summer. Recent observations and modeling studies have demonstrated that ongoing and projected climate change will increase plant productivity, microbial respiration, and growing season lengths at high-latitudes, but impacts on high-latitude C cycle seasonality (and potential feedbacks to the climate system) remain uncertain. Here we use ecosys, a well-tested and process-rich mechanistic ecosystem model that we evaluate further in this study, to explore how climate warming under an RCP8.5 scenario will shift C cycle seasonality in Alaska throughout the 21st century. The model successfully reproduced recently reported large high-latitude C losses during the fall and winter and yet still predicts a high-latitude C sink, pointing to a resolution of the current conflict between process-model and observation-based estimates of high-latitude C balance. We find that warming will result in surprisingly large changes in net ecosystem exchange (NEE; defined as negative for uptake) seasonality, with spring net C uptake overtaking summer net C uptake by year 2100. This shift is driven by a factor of 3 relaxation of spring temperature limitation to plant productivity that results in earlier C uptake and a corresponding increase in magnitude of spring NEE from –19 to –144 gC m –2 season –1 by the end of the century. Although a similar relaxation of temperature limitation will occur in the fall, radiation limitation during those months will limit increases in C fixation. Additionally, warmer soil temperatures and increased carbon inputs from plants lead to combined fall and winter C losses (163 gC m –2 ) that are larger than summer net uptake (123 gC m –2 season –1 ) by year 2100. However, this increase in microbial activity leads to more rapid N cycling and increased plant N uptake during the fall and winter months that supports large increases in spring NPP. Due to the large increases in spring net C uptake, the high-latitude atmospheric C sink is projected to sustain throughout this century. Our analysis disentangles the effects of key environmental drivers of high-latitude seasonal C balances as climate changes over the 21st century.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Dynamic land use implications of rapidly expanding and evolving wind power deployment

Abstract The expansion of wind power poses distinct and varied geographic challenges to a sustainable energy transition. However, current knowledge of its land use impacts and synergies is limited by reliance on static characterizations that overlook the role of turbine technology and plant design in mediating interactions with the environment. Here, we investigate how wind technology development and innovation have shaped landscape interactions with social and ecological systems within the United States and contribute to evolving land area requirements. This work assesses trends in key land use facets of wind power using a holistic set of metrics to establish an evidence base that researchers, technology designers, land use managers, and policymakers can use in envisioning how future wind-intensive energy systems may be jointly optimized for clean energy, social, and environmental objectives. Since 2000, we find dynamic land occupancy patterns and regional trends that are driven by advancing technology and geographic factors. Though most historical U.S. wind deployment has been confined to the temperate grassland biome in the nation’s interior, regional expansion has implicated diverse land use and cover types. A large percentage of the typical wind plant footprint (∼96% to > 99%) is not directly impacted by permanent physical infrastructure, allowing for multiple uses in the spaces between turbines. Surprisingly, turbines are commonly close to built structures. Moreover, rangeland and cropland have supported 93.4% of deployment, highlighting potential synergies with agricultural lands. Despite broadly decreasing capacity densities, offsetting technology improvements have stabilized power densities. Land use intensity, defined as the ratio of direct land usage to lifetime power generation of wind facilities, has also trended downwards. Although continued deployment on disturbed lands, and in close proximity to existing wind facilities and other infrastructure, could minimize the extent of impacts, ambitious decarbonization trajectories may predispose particular biomes to cumulative effects and risks from regional wind power saturation. Increased land-use and sustainability feedback in technology and plant design will be critical to sustainable management of wind power.

17 WIND ENERGY↗
DOE OSTI2024

A new data-driven map predicts substantial undocumented peatland areas in Amazonia

Tropical peatlands are among the most carbon-dense terrestrial ecosystems yet recorded. Collectively, they comprise a large but highly uncertain reservoir of the global carbon cycle, with wide-ranging estimates of their global area (441 025–1700 000 km 2 ) and below-ground carbon storage (105–288 Pg C). Substantial gaps remain in our understanding of peatland distribution in some key regions, including most of tropical South America. Here we compile 2413 ground reference points in and around Amazonian peatlands and use them alongside a stack of remote sensing products in a random forest model to generate the first field-data-driven model of peatland distribution across the Amazon basin. Our model predicts a total Amazonian peatland extent of 251 015 km 2 (95th percentile confidence interval: 128 671–373 359), greater than that of the Congo basin, but around 30% smaller than a recent model-derived estimate of peatland area across Amazonia. The model performs relatively well against point observations but spatial gaps in the ground reference dataset mean that model uncertainty remains high, particularly in parts of Brazil and Bolivia. For example, we predict significant peatland areas in northern Peru with relatively high confidence, while peatland areas in the Rio Negro basin and adjacent south-western Orinoco basin which have previously been predicted to hold Campinarana or white sand forests, are predicted with greater uncertainty. Similarly, we predict large areas of peatlands in Bolivia, surprisingly given the strong climatic seasonality found over most of the country. Very little field data exists with which to quantitatively assess the accuracy of our map in these regions. Data gaps such as these should be a high priority for new field sampling. This new map can facilitate future research into the vulnerability of peatlands to climate change and anthropogenic impacts, which is likely to vary spatially across the Amazon basin.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Critical percolation on the kagome hypergraph

In this report we study the percolation critical surface of the kagome lattice in which each triangle is allowed an arbitrary connectivity. Using the method of critical polynomials, we find points along this critical surface to high precision. This kagome hypergraph contains many unsolved problems as special cases, including bond percolation on the kagome and (3, 12 2 ) lattices, and site percolation on the hexagonal, or honeycomb, lattice, as well as a single point for which there is an exact solution. We are able to compute enough points along the critical surface to find a very accurate fit, essentially a Taylor series about the exact point, that allows estimations of the critical point of any system that lies on the surface to precision rivaling Monte Carlo and traditional techniques of similar accuracy. We find also that this system sheds light on some of the surprising aspects of the method of critical polynomials, such as why it is so accurate for certain problems, like the kagome and (3, 12 2 ) lattices. The bond percolation critical points of these lattices can be found to 17 and 18 digits, respectively, because they are in close proximity, in a sense that can be made quantitative, to the exact point on the critical surface. We also discuss in detail a parallel implementation of the method which we use here for a few calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The Thomas theorem and the Efimov States within a generalized Bohr model

Using a generalized Bohr model and the hyper-spherical formalism for a three-body system, we derive the universal energy functions. We apply our model to 4He atom where only the Coulomb potential is dominant and find that its binding energy is well reproduced (less than 5.5% difference). Later, we focus on the equal mass three-body systems and derive the Thomas theorem assuming a simple interaction depending on the range of the potential. We discuss the conditions for which an unbound two-body system produces a bound three-body system and apply our model to 4He and triton atoms as well as to the triton nucleus. Using their scattering lengths and effective ranges, we are able to reproduce the two-body or the three-body binding energies (less than 5% difference) with only one parameter fitted. Prediction for excited (Efimov) levels are also given and in particular we demonstrate that for some hyper-angles two equal minima appear which indicate a phase (shape) transition similar to the Landau's theory of phase transition. We suggest that the observed excited levels in two different experiments for the triton nucleus are indeed Efimov levels and there may be more surprises.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Network analysis of memristive device circuits: dynamics, stability and correlations

Abstract Networks with memristive devices are a potential basis for the next generation of computing devices. They are also an important model system for basic science, from modeling nanoscale conductivity to providing insight into the information-processing of neurons. The resistance in a memristive device depends on the history of the applied bias and thus displays a type of memory. The interplay of this memory with the dynamic properties of the network can give rise to new behavior, offering many fascinating theoretical challenges. But methods to analyze general memristive circuits are not well described in the literature. In this paper we develop a general circuit analysis for networks that combine memristive devices alongside resistors, capacitors and inductors and under various types of control. We derive equations of motion for the memory parameters of these circuits and describe the conditions for which a network should display properties characteristic of a resonator system. For the case of a purely memresistive network, we derive Lyapunov functions, which can be used to study the stability of the network dynamics. Surprisingly, analysis of the Lyapunov functions show that these circuits do not always have a stable equilibrium in the case of nonlinear resistance and window functions. The Lyapunov function allows us to study circuit invariances, wherein different circuits give rise to similar equations of motion, which manifest through a gauge freedom and node permutations. Finally, we identify the relation between the graph Laplacian and the operators governing the dynamics of memristor networks operators, and we use these tools to study the correlations between distant memristive devices through the effective resistance.

97 MATHEMATICS AND COMPUTING↗