Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interactive model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Origin of Ferromagnetic Exchange Coupling in Donor–Acceptor Biradical Analogues of Charge-Separated Excited States

A new donor–acceptor biradical complex, Tp Cum,Me Zn(SQ-VD) (Tp Cum,Me Zn + = zinc(II) hydro-tris(3-cumenyl-5-methylpyrazolyl)borate complex cation; SQ = orthosemiquinone; VD = oxoverdazyl), which is a ground-state analogue of a charge-separated excited state, has been synthesized and structurally characterized. Here, the magnetic exchange interaction between the S = 1/2 SQ and the S = 1/2 VD within the SQ-VD biradical ligand is observed to be ferromagnetic, with $J_{SQ-VD}$ = +77 cm –1 (H = -2$J_{SQ-VD}\bar{S}_{SQ}·\bar{S}_{VD}$) determined from an analysis of the variable-temperature magnetic susceptibility data. The pairwise biradical exchange interaction in Tp Cum,Me Zn(SQ-VD) can be compared with that of the related donor–acceptor biradical complex Tp Cum,Me Zn(SQ-NN) (NN = nitronyl nitroxide, S = 1/2), where $J_{SQ-NN}$ ≅ +550 cm –1 . This represents a dramatic reduction in the biradical exchange by a factor of ~7, despite the isolobal nature of the VD and NN acceptor radical SOMOs. Computations assessing the magnitude of the exchange were performed using a broken-symmetry density functional theory (DFT) approach. These computations are in good agreement with those computed at the CASSCF NEVPT2 level, which also reveals an S = 1 triplet ground state as observed in the magnetic susceptibility measurements. A combination of electronic absorption spectroscopy and CASSCF computations has been used to elucidate the electronic origin of the large difference in the magnitude of the biradical exchange coupling between Tp Cum,Me Zn(SQ-VD) and Tp Cum,Me Zn(SQ-NN). A Valence Bond Configuration Interaction (VBCI) model was previously employed to highlight the importance of mixing an SQ SOMO → NNLUMO charge transfer configuration into the electronic ground state to facilitate the stabilization of the high-spin triplet (S = 1) ground state in Tp Cum,Me Zn(SQ-NN). Here, CASSCF computations confirm the importance of mixing the pendant radical (e.g., VD, NN) LUMO (VD LUMO and NN LUMO ) with the SOMO of the SQ radical (SQ SOMO ) for stabilizing the triplet, in addition to spin polarization and charge transfer contributions to the exchange. An important electronic structure difference between Tp Cum,Me Zn(SQ-VD) and Tp Cum,Me Zn(SQ-NN), which leads to their different exchange couplings, is the reduced admixture of excited states that promote ferromagnetic exchange into the Tp Cum,Me Zn(SQ-VD) ground state, and the intrinsically weaker mixing between the VD LUMO and the SQ SOMO compared to that observed for Tp Cum,Me Zn(SQ-NN), where this orbital mixing is significant. The results of this comparative study contribute to a greater understanding of biradical exchange interactions, which are important to our understanding of excited-state singlet–triplet energy gaps, electron delocalization, and the generation of electron spin polarization in both the ground and excited states of (bpy)Pt(CAT-radical) complexes.

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

Flavoured warped axion

We consider a 5D extension of the DFSZ axion model that addresses both the axion quality and fermion mass hierarchy problems, and predicts flavour-dependent, off-diagonal axion-fermion couplings. The axion is part of a 5D complex scalar field charged under a U(1)PQ symmetry that is spontaneously broken in the bulk, and is insensitive to explicit PQ breaking on the UV boundary. Bulk Standard Model fermions interact with two Higgs doublets that can be localized on the UV boundary or propagate in the bulk to explain the fermion masses and mixings. When the Higgs doublets are localized on the UV boundary, they induce flavour diagonal couplings between the fermions and the axion. However, when the Higgs doublets propagate in the bulk, the overlap of the axion and fermion profiles generates flavour off-diagonal couplings. The effective scale of these off-diagonal couplings in both the quark and lepton sectors can be as small as 10 11 GeV, and therefore will be probed in future precision flavour experiments.

Global Symmetries↗
NASA NTRS2019

Optimal Control within the Context of Multidisciplinary Design, Analysis, and Optimization

Multidisciplinary design, analysis and optimization involves modeling the interactions of complex systems across a variety of disciplines. The optimization of such systems can be a computationally expensive exercise with multiple levels of nested nonlinear solvers running under an optimizer.The application of optimal control in project development often involves performing trajectory optimization for fixed vehicle designs or parametric sweeps across some key vehicle properties.This information is then relayed to the subsystem design teams who update their designs and relay some bulk characteristics back to the trajectory optimization procedure.This iteration is then repeated until the design closes.However, with increasing interest in more tightly coupled systems, such as electric and hybrid-electric aircraft propulsion and boundary layer ingestion, this process is prone to ignore subtle coupling between vehicle subsystem designs and vehicle operation on a given mission.Integrating trajectory optimization into a tightly coupled multidisciplinary design procedure can be computationally prohibitive, depending on the complexity of the subsystem analyses and the optimal control technique applied.To address these issues a new optimal control software tool, Dymos, has been developed.Dymos is built upon NASA's OpenMDAO software and can leverage its capabilities to efficiently compute gradients for the optimization and optimize complex models in parallel on distributed memory systems.This report provides some explanation into the numerical methods employed in Dymos and provides several use cases that demonstrate its performance on traditional optimal control problems and improvements ino techniques have been used extensively in recent decades to solve a variety of optimal control problems, typically in the form of aerospace vehicle trajectory optimization.

pseudospectral↗
DOE OSTI2026

Quantum Theory of Surface Lattice Resonances

The collective interactions of nanoparticles arranged in periodic structures give rise to high‐ in‐plane diffractive modes known as surface lattice resonances. Although these resonances and their broader implications have been extensively studied within the framework of classical electrodynamics and linear response theory, a quantum optical theory capable of describing the dynamics of these structures, especially in the presence of material nonlinearities beyond ad hoc few‐mode approximations, is largely missing. To this end, we consider a lattice of metallic nanoparticles coupled to the electromagnetic field and derive the quantum input–output relations within the electric dipole approximation. As applications, we analyze coupling between the nanoparticle array and external quantum emitters, and show how the formalism extends to molecular optomechanics, where the high ‐factors of SLRs enable coupling to collective vibrational modes. We further consider arrays composed of saturable excitonic emitters, demonstrating how emitter nonlinearities can be used to switch the SLR condition between electronic transitions. Using a perturbative approach that accounts for population dynamics, we show how these effects can be probed in pump–probe experiments and give rise to nonlinear phase‐matching phenomena. Our work provides a microscopic framework for modeling SLRs interacting with quantum emitters without phenomenological descriptions of the electromagnetic environment.

molecular optomechanics↗
DOE OSTI2021

Constructing tensor network influence functionals for general quantum dynamics

Here, we describe an iterative formalism to compute influence functionals that describe the general quantum dynamics of a subsystem beyond the assumption of linear coupling to a quadratic bath. We use a space-time tensor network representation of the influence functional and investigate its approximability in terms of its bond dimension and time-like entanglement in the tensor network description. We study two numerical models, the spin-boson model and a model of interacting hard-core bosons in a 1D harmonic trap. We find that the influence functional and the intermediates involved in its construction can be efficiently approximated by low bond dimension tensor networks in certain dynamical regimes, which allows the quantum dynamics to be accurately computed for longer times than with direct time evolution methods. However, as one iteratively integrates out the bath, the correlations in the influence functional can first increase before decreasing, indicating that the final compressibility of the influence functional is achieved via non-trivial cancellation.

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

Mechanism of endosomal escape by p H -responsive nucleic-acid vectors

Successful intracellular delivery of nucleic acids (NAs) hinges on many factors, one of them being NAs’ efficacious escape from endosomes. As competent NA vectors, pH-responsive gemini surfactants (GSs) might achieve high efficacy by facilitating endosomal escape. However, how the GSs assist the escape remains debated as many proposed mechanisms still lack experimental support, which hinders replication and further improvement of the efficient delivery. Here, via UV, fluorescence spectroscopy, and small-angle neutron scattering (SANS), we examined a pH-responsive GS's and a pH-unresponsive GS's capabilities to compact DNA and withstand binding competition, and their interactions with model endosomal and lysosomal membranes, at varied pHs. Acidification-driven enhancement of DNA-compaction capability and of stability against binding competition were found specific to the pH-responsive GS. Alongside the pH-responsive GS's structural perturbation to the membranes as observed with SANS, the features suggest that pH-responsive GSs facilitate endosomal escape by releasing excess GS molecules from DNA-GS complexes upon acidification in endosome maturation, with the released GS molecules disrupting endosomal and lysosomal membranes and thereby assisting the escape. Here a general design principle for NA vectors is proposed on the basis of this experimental finding.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

PSCF+: An Extended and Improved Open-Source Software Package for Polymer Self-Consistent Field Calculations

This work introduces PSCF+, a freely available open‑source software tool for modeling how block copolymers self‑assemble into ordered nanoscale structures. The software allows researchers to study a wide range of polymer models and interaction types, making it possible to directly compare theoretical predictions with molecular and simulation‑based studies in a consistent way. PSCF+ incorporates several recently developed numerical methods that significantly reduce computer memory requirements and speed up calculations, enabling efficient use of modern GPU hardware. Here, it also provides automated tools to calculate free‑energy landscapes and phase boundaries, which are essential for understanding and predicting material behavior. Together, these advances make PSCF+ a powerful and flexible platform for studying polymer self‑assembly and support the design of new polymeric materials for scientific and technological applications.

Chemical calculations↗
NASA NTRS2007

Runtime and Pressurization Analyses of Propellant Tanks

Multi-element unstructured CFD has been utilized at NASA SSC to carry out analyses of propellant tank systems in different modes of operation. The three regimes of interest at SSC include (a) tank chill down (b) tank pressurization and (c) runtime propellant draw-down and purge. While tank chill down is an important event that is best addressed with long time-scale heat transfer calculations, CFD can play a critical role in the tank pressurization and runtime modes of operation. In these situations, problems with contamination of the propellant by inclusion of the pressurant gas from the ullage causes a deterioration of the quality of the propellant delivered to the test article. CFD can be used to help quantify the mixing and propellant degradation. During tank pressurization under some circumstances, rapid mixing of relatively warm pressurant gas with cryogenic propellant can lead to rapid densification of the gas and loss of pressure in the tank. This phenomenon can cause serious problems during testing because of the resulting decrease in propellant flow rate. With proper physical models implemented, CFD can model the coupling between the propellant and pressurant including heat transfer and phase change effects and accurately capture the complex physics in the evolving flowfields. This holds the promise of allowing the specification of operational conditions and procedures that could minimize the undesirable mixing and heat transfer inherent in propellant tank operation. It should be noted that traditional CFD modeling is inadequate for such simulations because the fluids in the tank are in a range of different sub-critical and supercritical states and elaborate phase change and mixing rules have to be developed to accurately model the interaction between the ullage gas and the propellant. We show a typical run-time simulation of a spherical propellant tank, containing RP-1 in this case, being pressurized with room-temperature nitrogen at 540 R. Nitrogen, shown in blue on the right-hand side of the figures, enters the tank from the diffuser at the top of the figures and impinges on the RP-1, shown in red, while the propellant is being continuously drained at the rate of 1050 lbs/sec through a pipe at the bottom of the tank. The sequence of frames in Figure 1 shows the resultant velocity fields and mixing between nitrogen and RP-1 in a cross-section of the tank at different times. A vortex is seen to form in the incoming nitrogen stream that tends to entrain propellant, mixing it with the pressurant gas. The RP-1 mass fraction contours in Figure 1 are also indicative of the level of mixing and contamination of the propellant. The simulation is used to track the propagation of the pure propellant front as it is drawn toward the exit with the evolution of the mixing processes in the tank. The CFD simulation modeled a total of 10 seconds of run time. As is seen from Figure 1d, after 5.65 seconds the propellant front is nearing the drain pipe, especially near the center of the tank. Behind this pure propellant front is a mixed fluid of compromised quality that would require the test to end when it reaches the exit pipe. Such unsteady simulations provide an estimate of the time that a high-quality propellant supply to the test article can be guaranteed at the modeled mass flow rate. In the final paper, we will discuss simulations of the LOX and propellant tanks at NASA SSC being pressurized by an inert ullage. Detailed comparisons will be made between the CFD simulations and lower order models as well as with test data. Conditions leading to cryo collapse in the tank will also be identified.

Field, Robert E.↗
DOE OSTI2026

Strong Sensitivity of Simulated Biomass Burning Aerosol Transport and Radiative Effects Over the South Atlantic to Carbonaceous Aerosol Aging and Particle Density

Biomass burning aerosol (BBA) impacts climate through aerosol‐cloud‐radiation interactions, but models disagree on the sign and magnitude of BBA radiative effects. We quantify the sensitivity of BBA radiative effects and transport to three BBA‐relevant processes and properties: parameterized oxidative aging of organic aerosol (OA), a combined change to black carbon (BC) density and the method for calculating aerosol refractive index, and reduction in OA density. We evaluate Unified Model simulations against two aircraft campaigns from summer 2017 over the Southeast Atlantic. The model generally performs well, such that discrepancies between the observational data sets may sometimes limit the precision of the evaluation. Our newly developed aging parameterization reproduces observed OA:BC mass ratios well and allows modeled OA:BC to decrease with smoke age, but increases bias in aerosol extinction and changes the BBA radiative effect little (+0.12 W m -2 ). We calculate aerosol refractive index using either a volume‐weighted component average or the Maxwell‐Garnett (MG) mixing assumption, which represents BC as small inclusions in a host material. Compared to MG mixing, the volume‐weighted average refractive index and reduced BC density increase aerosol absorption, substantially increasing the total BBA radiative effect (+2.66 W m -2 ) and amount of BBA transported across the ocean through BC self‐lofting. Reducing OA density to better match literature values changes the total BBA radiative effect by −1.96 W m -2 . Changes to direct radiative effects exceed changes to cloud radiative effects. Our findings emphasize the sensitivity of aerosol radiative effects and transport to these processes and properties, which we suggest could be improved in climate models.

Southeast Atlantic↗
DOE OSTI2023

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the “gold standard” coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller−Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

Chemistry↗
DOE OSTI2021

BETO 2021 Peer Review - WBS 4.1.2.32: Bioeconomy Scenario Analysis

The Bioeconomy Scenario Analysis project uses systems thinking and analysis to assess current and/or prospective techno-economics, research and development, deployment strategies, policy, and market conditions and their impact on the potential development trajectories of the bioenergy industry over time. Results from this project include identification of opportunities and constraints to industrial development, quantification of multiple metrics (energy, economic, environmental) and informing researchers, decision makers, and industry of the steps needed for a sustainable, nationwide biofuels industry. Analyses from this project enable the creation of a bioenergy industry by (1) inciting policy-makers to explore scenarios for nationwide biofuels production, identifying policy actions consistent with pathways for growth; (2) improving industry’s understanding of industry growth potential under different technology and investment conditions, better targeting their development efforts; and (3) providing universities and other interested stakeholders with tools and analyses that can be adapted to meet research and teaching objectives, connecting students with careers that build the industry. One of the many modeling tools used in this project, the Biomass Scenario Model (BSM) is a publicly-available, unique, validated, state-of-the-art, award-winning, fourth-generation model of the domestic biofuels supply chain which explicitly focuses on how and under what conditions biofuel technologies might be deployed to contribute to the U.S. transportation energy sector. We use models like the BSM to examine the implications of policies and incentives as well as their potential side-effects. The BSM uses a system-dynamics simulation to model dynamic interactions and transitions across the supply chain; it tracks the deployment of biofuels given industrial learning and the reaction of the investment community in the context of land availability, projected oil markets, consumer demand for biofuels, and government policies over time. Under expected market conditions, analyses using the BSM suggest that the biofuels industry may require significant external actions in the early years to thrive. Interventions that accelerate the industrial learning process (e.g. operation of pre-commercial and commercial facilities) have been identified as having strong influence in starting the growth of a commercial biofuel industry. Policies which are coordinated across the whole supply chain in BSM foster the growth of the biofuels industry and production of tens of billions of gallons of biofuels may occur under sufficiently favorable conditions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS2021

Extreme Ozone Loss Following Nuclear War Results in Enhanced Surface Ultraviolet Radiation

For the first time, we use a modern climate model with interactive chemistry including the effects of aerosols on photolysis rates to simulate the consequences of regional and global scale nuclear wars (injecting 5 and 150 Tg of soot respectively) for the ozone layer and surface ultraviolet (UV) light. For a global nuclear war, heating in the stratosphere, reduced photolysis, and an increase in catalytic loss from the Hox cycle cause a 15 year-long reduction in the ozone column, with a peak loss of 75% globally and 65% in the tropics. This is larger than predictions from the 1980s, which assumed large injections of nitrogen oxides (NOx),but did not include the effects of smoke. NOx from the fireball and the fires provide a small (5%) Increase to the global average ozone loss for the first few years. Initially, soot would shield the surface from UV-B, but UV Index values would become extreme: greater than 35in the tropics for 4years, and greater than 45during the summer in the southern polar regions for 3 years. For a regional war, global column ozone would be reduced by 25% with recovery taking 12 years. This is similar to previous simulations, but with a faster recovery time due to a shorter lifetime for soot in our simulations. In-line photolysis provides process specific action spectra enabling future integration with biogeochemistry models and allows output that quantifies the potential health impacts from changes in surface UV for this and other larger aerosol injections.

nuclear war↗
DOE OSTI2023

Multi-level optimization with the koopman operator for data-driven, domain-aware, and dynamic system security

Cyber-Physical Systems (CPSs) like the power grid are critically important but also increasingly vulnerable; ensuring reliable system operation in the face of disruptions is becoming more and more challenging. Multi-Level Optimization (MLO) is a powerful way to model adversarial interactions, which naturally makes it applicable to studying CPS security. However, MLO typically does not address underlying system dynamics, and incorporating nonlinear dynamics is generally infeasible. In this paper, we show how to combine MLO with the Koopman Operator (KO) to remedy this. The KO maps nonlinear dynamics to a lifted space in which those dynamics are linear, thus making it ideal for use with MLO. Moreover, the structure of the KO also provides convenient ways to incorporate domain knowledge into the data-driven process of learning the KO representation of a given system. Here we then demonstrate the use of MLO-KO on a small example problem taken from the power grid domain, discuss the scalability and computational cost of MLO-KO, and identify future research directions for this work.

42 ENGINEERING↗
DOE OSTI2025

Modeling and Analysis of Clean Energy and Storage Technologies (CRADA Final Report, Project 1)

The goal of this project is to provide Southern Company Services, Inc. ("Participant") with custom scripts that can be used to create an average PV energy production profile, calculate lifetime energy value, calculate capacity value, and calculate the resultant financial metrics considering those value streams. Secondly, a fuel-cell model will be added to the public version of System Advisor Model (SAM). This standalone technology will incorporate PV and battery storage, allowing the participant to model the interaction of these three technologies. By adding this capability to a public version of SAM, a broad audience will be able to consider the system performance and financial benefits of installing a fuel cell as a baseline generator with PV. Thirdly, automated dispatch algorithms will be developed and added to SAM. These algorithms will enable Southern Company to dispatch a DC-connected front-of-the-meter battery system while considering price signals and PV clipping behavior. By adding these capabilities to a public version of SAM, users will be able to consider more complex and realistic ways of dispatching a battery system.

14 SOLAR ENERGY↗
NASA NTRS1986

Interactions between stationary planetary waves in the stratosphere

The interactions between stationary planetary waves 1 and 2 are investigated using a truncated quasi-geostrophic model in a midlatitude beta channel. The model includes eight interacting horizontal modes with three zonal wave numbers and two meridional scales in a quasi-geostrophic beta channel 40 deg wide; all wave-wave interactions in the model are a result of dissipation. The procedures for solving the model equations are examined. The theory of dissipation-induced interactions is discussed in terms of interactions between weak waves under severe truncation. The behavior of small amplitude waves 1 and 2 in basic states is discussed; it is observed that wave 1 is amplified in a wave 2 basic state, and a wave 2 disturbance in a wave 1 basic state is affected by the relative phases of the waves. The dynamics of these wave interactions are analyzed using potential enstrophy budgets of the waves. The relevance of the waves to the atmosphere is studied. It is noted that the wave behavior is a combination of the linear results with an amplification of wave 1 and a strong dependence on the relative phases of waves 1 and 2 in the stratosphere.

Robinson, W. A.↗
NASA NTRS1983

Solar wind flow about the terrestrial planets. II - Comparison with gas dynamic theory and implications for solar-planetary interactions

Bow shock models are employed to determine and compare the solar wind flows about Venus, Mars, and the earth. Initially, gas dynamic theory predictions are compared with the large data base available for the near-earth environment. The observed shape and location of the magnetopause proved sufficient for predicting the average dayside bow shock position to within 2% error. Use of the same gas dynamics theory produced highly disparate results for the flows past Venus and Mars. A variation of 510-1000 km altitude (high) for the solar wind-obstacle interface was found for Mars, thus requiring the presence of an effective magnetic moment of within 0.6 of 1.4 x 10 to the 22nd G/cu cm in the Martian magnetosphere, consistent with Viking data. The Venus bow wave was calculated to be closer to the planet than possible with a purely ionospheric interaction, indicating that solar wind-neutral atmosphere interactions in the lower ionosheath must necessarily be included in a gas dynamic modeling of the Venus obstacle to the solar wind.

Slavin, J. A.↗
DOE OSTI2023

Reversal of spin-polarization near the Fermi level of the Rashba semiconductor BiTeCl

Abstract Spin–orbit coupling forms the physical basis for quantum materials with non-trivial topology and potential spintronics applications. The Rashba interaction is a textbook model of spin–orbit interactions, with charge carriers undergoing linear, isotropic spin-splitting in momentum space. Recently, non-centrosymmetric semiconductors in the family BiTe X ( X = Cl, Br, I) have been identified as exemplary Rashba materials due to the strong splitting of their bulk bands, yet a detailed investigation of their spin textures, and their relationships to local crystal symmetry, is currently lacking. We perform high-efficiency spin-resolved photoemission spectroscopy to directly image the spin texture of surface states of BiTeCl, and we find dramatic deviations from idealized behavior, including a reversal of the spin-polarization near the Fermi level. We show that this behavior can be described by higher-order contributions to the canonical Rashba model with the surface states localized to individual trilayers of the crystal. Due to the prominence of these effects near the Fermi level, they should have a strong impact on the spin-dependent transport of carriers.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Proton-neutron entanglement in the nuclear shell model

Abstract We compute the proton-neutron entanglement entropy in the interacting nuclear shell model for a variety of nuclides and interactions. Some results make intuitive sense, for example, that the shell structure, as governed by single-particle and monopole energies, strongly affects the energetically available space and thus the entanglement entropy. We also find a surprising result: that the entanglement entropy at low excitation energy tends to decrease for nuclides when N ≠ Z . While we provide evidence this arises from the physical nuclear force by contrasting with random two-body interactions which shows no such decrease, the exact mechanism is unclear. Nonetheless, the low entanglement suggests that in models of neutron-rich nuclides, the coupling between protons and neutrons may be less computationally demanding than one might otherwise expect.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗