Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter space exploration”

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 505 records · Page 28

3D printing of void-free glass monoliths: rheological and geometric considerations

Direct ink writing (DIW) is a versatile additive manufacturing technique capable of printing 3D architectures from a broad range of materials. The key advantage of DIW is the ability to spatially control both the architecture and the composition of the printed part, which enables the production of components with previously unachievable combinations of functionalities. One emergent application space for DIW has been in additively manufactured glass optics, where a “green body” is patterned layer-by-layer by extruding a silica-based slurry through a microfluidic nozzle and then thermally consolidated to transparent glass. A critical aspect of the process is ensuring that the printed part is completely free of voids that can occur in the interstitial spaces between the extruded filaments. This outcome is governed by the interplay between ink rheology and strains imposed by the printing process and geometric packing of the filaments. Here, we explore the strain- and rate-dependent deformation events in the DIW process to determine conditions that enable the printing of void-free monoliths without sacrificing shape fidelity. We focus on yield stress fluids consisting of fumed silica nanoparticles dispersed in organic solvents at varied volume fractions to control the rheological properties. We investigate the printability of these ink materials as a function of print process variables and devise the appropriate dimensionless parameters that capture the geometry-dependent and rate-dependent effects on the deformation of the extruded ink. Lastly, we use these scaling arguments to construct a phase diagram for void-free monolith printability.

36 MATERIALS SCIENCE↗

First Sterile Neutrino Search at NOvA Experiment Using both Neutrino and Anti-Neutrino Data

Neutrinos are among the most abundant fundamental particles in the universe, and yet they are so elusive that they pass through almost everything without leaving a trace. For decades, researchers have measured neutrino oscillations in great detail, refining our understanding of the theories. Yet, once again, neutrino shows their incredibly surprising nature through some recent experimental results, suggesting that our current understanding might be incomplete. Several experiments have observed an unexpected difference in the observed neutrino behaviour from what the Standard Model predicts, showing an excess in some cases (e.g. LSND, MiniBooNE) and a deficit of neutrino events in others (e.g. SAGE, GALLEX). These anomalies challenge our conventional theories of standard three neutrino flavors and phenomenologically motivate the sterile neutrino hypothesis, which does not interact through the usual weak force. While some results appear consistent with this idea, others have found no supporting evidence (e.g. MicroBooNE), leaving the question unanswered. NOvA is one such long-baseline neutrino experiment with its Near Detector (ND) at Fermilab and a Far Detector (FD) at 810 km in Minnesota. It uses the NuMI beam as its primary source of neutrinos for oscillation studies. Its two-detector setup, extended baseline, and ability to conduct combined $\nu$-$\bar\nu$ analyses position it to explore these phenomena in new ways. We use the joint $\nu_{\mu}$ and neutral current (NC) disappearance channels to probe active-sterile mixing in the 3+1 model. This analysis utilizes a data sample of 27$\times$10$^{20}$ protons on target (POT) in neutrino mode and 12.5$\times$10$^{20}$ protons on target (POT) in antineutrino mode. We found no evidence of sterile neutrino at 90$\%$ C.L. We present the limits in $\sin^{2}\theta_{24}$ vs $\Delta m^{2}_{41}$, $\sin^{2}\theta_{34}$ vs $\Delta m^{2}_{41}$, and $\sin^{2}2\theta_{\mu\tau}$ vs $\Delta m^{2}_{41}$ parameter spaces. The limits provided by this analysis are world leading in most of the parameter spaces.

Chaudhary, Shivam [Indian Inst. of Info. Tech. Guw↗

Using Parameter Sweep in WaterTAP to Analyze New Water Treatment Technologies

We describe a powerful and generalized parameter sweep tool in this report that was originally developed to analyze the performance of existing and novel water treatment models being developed in WaterTAP. Since WaterTAP is built upon IDAES and Pyomo, the parameter sweep tool can be used to systematically explore and debug the behavior of most Pyomo and IDAES numerical models. In order to enable meaningful analyses, the parameter sweep tool has been designed with the following features: 1) Model flexibility: The parameter sweep tool does not enforce any restrictions on the types of models that can be used with it. As long as a Pyomo model can be solved and the parameter is active and mutable, the tool only needs functions that describe how to run the model, the sweep parameters, and the output quantities of interest. 2) Flexible sampling: The parameter sweep tool has inbuilt functions to generate samples from a random distribution or a multidimensional Euclidean space. Furthermore, the users have to ability to supply samples generated from a tool of their choice. 3) Multiple sweep types: A user can choose from one of 3 types of parameter sweeps depending on their needs. 4) Detailed outputs: Outputs generated by the parameter sweep tool can be stored in detailed H5 file or user-friendly CSV files for post processing. 5) Parallel computing: The parameter sweep supports shared and distributed memory parallel computing to enable the use of high performance computers (HPC) for large-scale analyses. 6) Modular: The parameter sweep tool is self-contained and can easily be integrated within an outer-loop analysis or as desired by the user. 7) Ease of use: The tool is well documented and a simple sweep can be easily executed by following the online documentation in a few lines of code. We demonstrate the use of the parameter sweep tool on a simple water treatment system from the WaterTAP repository and show its parallel scaling performance on an Apple laptop and NREL's Eagle HPC. The parameter sweep tool is actively being used with models currently being developed within WaterTAP and we expect its use to grow beyond it to other IDAES and Pyomo models.

97 MATHEMATICS AND COMPUTING↗

Dark Matter-Induced Nuclear De-Excitation at SBND with Ab Initio Nuclear Theory

We explore the sensitivity of the Short-Baseline Near Detector (SBND) experiment to light dark matter using MeV-scale electromagnetic activity. Inelastic scattering of dark matter with argon nuclei can excite nuclear states that subsequently de-excite via the emission of MeV-scale photons, producing localized low-energy "blip" signatures in a liquid argon time projection chamber. We perform state-of-the-art ab initio nuclear calculations, including all relevant argon excited states with energies up to 18 MeV, to provide reliable predictions for these signals. After accounting for relevant backgrounds, we find that SBND can probe previously unexplored regions of parameter space for light dark matter.

Dutta, Bhaskar [Bucharest, IFIN-HH; Texas A-M] (OR↗

Coordination and control – limits in standard representations of multi-reservoir operations in hydrological modeling

Abstract. Major multi-reservoir cascades represent a primary mechanism for dealing with hydrologic variability and extremes within institutionally complex river basins worldwide. These coordinated management processes fundamentally reshape water balance dynamics. Yet, multi-reservoir coordination processes have been largely ignored in the increasingly sophisticated representations of reservoir operations within large-scale hydrological models. The aim of this paper is twofold, namely (i) to provide evidence that the common modeling practice of parameterizing each reservoir in a cascade independently from the others is a significant approximation and (ii) to demonstrate potential unintended consequences of this independence approximation when simulating the dynamics of hydrological extremes in complex reservoir cascades. We explore these questions using the Water Balance Model, which features detailed representations of the human infrastructure coupled to the natural processes that shape water balance dynamics. It is applied to the Upper Snake River basin in the western US and its heavily regulated multi-reservoir cascade. We employ a time-varying sensitivity analysis that utilizes the method of Morris factor screening to explicitly track how the dominant release rule parameters evolve both along the cascade and in time according to seasonal high- and low-flow events. This enables us to address aim (i) by demonstrating how the progressive and cumulative dominance of upstream releases significantly dampens the ability of downstream reservoir rules' parameters to influence flow conditions. We address aim (ii) by comparing simulation results with observed reservoir operations during critical low-flow and high-flow events in the basin. Our time-varying parameter sensitivity analysis with the method of Morris clarifies how independent single-reservoir parameterizations and their tacit assumption of independence leads to reservoir release behaviors that generate artificial water shortages and flooding, whereas the observed coordinated cascade operations avoided these outcomes for the same events. To further explore the role of (non-)coordination in the large deviations from the observed operations, we use an offline multi-reservoir water balance model in which adding basic coordination mechanisms drawn from the observed emergency operations is sufficient to correct the deficiencies of the independently parameterized reservoir rules from the hydrological model. These results demonstrate the importance of understanding the state–space context in which reservoir releases occur and where operational coordination plays a crucial role in avoiding or mitigating water-related extremes. Understanding how major infrastructure is coordinated and controlled in major river basins is essential for properly assessing future flood and drought hazards in a changing world.

54 ENVIRONMENTAL SCIENCES↗

Exploring regular and turbulent flow states in active nematic channel flow via Exact Coherent Structures and their invariant manifolds

This work is a unified study of stable and unstable steady states of 2D active nematic channel flow using the framework of Exact Coherent Structures (ECS). ECS are stationary, periodic, quasiperiodic, or traveling wave solutions of the governing equations that, together with their invariant manifolds, organize the dynamics of nonlinear continuum systems. We extend our earlier work on ECS in the preturbulent regime by performing a comprehensive study of stable and unstable ECS for a wide range of activity values spanning the preturbulent and turbulent regimes. In the weakly turbulent regime, we compute more than 200 unstable ECS that co-exist at a single set of parameters, and uncover the role of symmetries in organizing the phase space geometry. We provide conclusive numerical evidence that in the preturbulent regime, generic trajectories shadow a series of unstable ECS before settling onto an attractor. Lastly, our studies hint at shadowing of quasiperiodic type ECS in the turbulent regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

An interoperable implementation of collective‐variable based enhanced sampling methods in extended phase space within the OpenMM package

Collective variable (CV)-based enhanced sampling techniques are widely used today for accelerating barrier-crossing events in molecular simulations. A class of these methods, which includes temperature accelerated molecular dynamics (TAMD)/driven-adiabatic free energy dynamics (d-AFED), unified free energy dynamics (UFED), and temperature accelerated sliced sampling (TASS), uses an extended variable formalism to achieve quick exploration of conformational space. These techniques are powerful, as they enhance the sampling of a large number of CVs simultaneously compared to other techniques. Extended variables are kept at a much higher temperature than the physical temperature by ensuring adiabatic separation between the extended and physical subsystems and employing rigorous thermostatting. Here, in this work, we present a computational platform to perform extended phase space enhanced sampling simulations using the open-source molecular dynamics engine OpenMM. The implementation allows users to have interoperability of sampling techniques, as well as employ state-of-the-art thermostats and multiple time-stepping. This work also presents protocols for determining the critical parameters and procedures for reconstructing high-dimensional free energy surfaces. As a demonstration, we present simulation results on the high dimensional conformational landscapes of the alanine tripeptide in vacuo, tetra-N-methylglycine (tetra-sarcosine) peptoid in implicit solvent, and the Trp-cage mini protein in explicit water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Explainable AI for Multivariate Time Series Pattern Exploration: Latent Space Visual Analytics With Temporal Fusion Transformer and Variational Autoencoders in Power Grid Event Diagnosis

Detecting and analyzing complex patterns in multivariate time-series data is crucial for decision-making in urban and environmental system operations. However, challenges arise from the high dimensionality, intricate complexity, and interconnected nature of complex patterns, which hinder the understanding of their underlying physical processes. Existing AI methods often face limitations in interpretability, computational efficiency, and scalability, reducing their applicability in real-world scenarios. This paper proposes a novel visual analytics framework that integrates two generative AI models, Temporal Fusion Transformer (TFT) and Variational Autoencoders (VAEs), to reduce complex patterns into lower-dimensional latent spaces and visualize them in 2D using dimensionality reduction techniques such as PCA, t-SNE, and UMAP with DBSCAN. These visualizations, presented through coordinated and interactive views and tailored glyphs, enable intuitive exploration of complex multivariate temporal patterns, identifying patterns’ similarities and uncover their potential correlations for a better interpretability of the AI outputs. The framework is demonstrated through a case study on power grid signal data, where it identifies multi-label grid event signatures, including faults and anomalies with diverse root causes. Additionally, novel metrics and visualizations are introduced to validate the models and assess the performance, efficiency, and consistency of latent maps generated by VAE, which have been utilized in prior studies for latent space cartography and used as a benchmark in this study, and the emerging TFT architecture under various configurations. These analyses provide actionable insights for model parameter tuning and reliability improvements. Comparative results highlight that TFT achieves shorter run times and superior scalability to diverse time-series data shapes compared to VAE. This work advances fault diagnosis in multivariate time series, fostering explainable AI to support critical system operations.

Explainable AI↗

Observable proton decay in flipped SU(5)

We explore proton decay in a class of realistic supersymmetric flipped SU(5) models supplemented by a U(1) R symmetry which plays an essential role in implementing hybrid inflation. Two distinct neutrino mass models, based on inverse seesaw and type I seesaw, are identified, with the latter arising from the breaking of U(1) R by nonrenormalizable superpotential terms. Depending on the neutrino mass model an appropriate set of intermediate scale color triplets from the Higgs superfields play a key role in proton decay channels that include p → (e + , μ + ) π 0 , p → (e + , μ + ) K 0 , p → $\bar v$π + , and p → $\bar v$K + . We identify regions of the parameter space that yield proton lifetime estimates which are testable at Hyper-Kamiokande and other next generation experiments. We discuss how gauge coupling unification in the presence of intermediate scale particles is realized, and a Z 4 symmetry is utilized to show how such intermediate scales can arise in flipped SU(5). Finally, we compare our predictions for proton decay with previous work based on SU(5) and flipped SU(5).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Various Low-Symmetry Phases in High- χ and Conformationally Asymmetric PDMS- b -PTFEA Copolymers

Since the discovery of the Frank–Kasper (FK) phase a decade ago, single-component block copolymer systems have offered unique opportunities to explore the structural and phase-transition properties of topologically close-packed FK phases composed of point-particle motifs. Here, we present the formation of complex low-symmetry phases and their thermodynamic stability using high-χ and high-ε linear diblock copolymers, where χ and ε are the Flory–Huggins interaction parameter and conformational asymmetry between the two blocks, respectively. In this work, a series of polydimethylsiloxane-b-poly(2,2,2-trifluoroethyl acrylate) (PDMS-b-PTFEA) copolymers are prepared by tuning the volume fraction ($f_{\text{PTFEA}}$) of a minor PTFEA block. The packing structures of PDMS-b-PTFEAs are mapped in the space of temperature versus $f_{\text{PTFEA}}$ using X-ray scattering measurements. In addition to the packing structures reported earlier, we identify another A15 phase at higher $f_{\text{PTFEA}}$ close to a hexagonally packed cylinder (HEX) phase and establish the phase sequence of HEX–A15–σ–C14–BCC with decreasing $f_{\text{PTFEA}}$. Theoretical calculation of relative free-energy densities confirms the experimental phase sequence of the composition-dependent sphere-packing structures. Interestingly, the domain sizes of three-dimensional (3D) sphere-packing structures are nearly temperature-invariant compared to that of a two-dimensional HEX structure, and we attribute this dimensional stability to the strong enthalpic effects in the formation of 3D packing structures and the geometrical characteristics of spherical domains. Our results demonstrate that the interaction parameter χ, together with the known role of conformational asymmetry ε, serves as the primary thermodynamic parameter for selecting and stabilizing various low-symmetry packing structures in sphere-forming linear diblock copolymer systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Zero-truncated Poisson regression for sparse multiway count data corrupted by false zeros

Abstract We propose a novel statistical inference methodology for multiway count data that is corrupted by false zeros that are indistinguishable from true zero counts. Our approach consists of zero-truncating the Poisson distribution to neglect all zero values. This simple truncated approach dispenses with the need to distinguish between true and false zero counts and reduces the amount of data to be processed. Inference is accomplished via tensor completion that imposes low-rank tensor structure on the Poisson parameter space. Our main result shows that an $N$-way rank-$R$ parametric tensor $\boldsymbol{\mathscr{M}}\in (0,\infty )^{I\times \cdots \times I}$ generating Poisson observations can be accurately estimated by zero-truncated Poisson regression from approximately $IR^2\log _2^2(I)$ non-zero counts under the nonnegative canonical polyadic decomposition. Our result also quantifies the error made by zero-truncating the Poisson distribution when the parameter is uniformly bounded from below. Therefore, under a low-rank multiparameter model, we propose an implementable approach guaranteed to achieve accurate regression in under-determined scenarios with substantial corruption by false zeros. Several numerical experiments are presented to explore the theoretical results.

97 MATHEMATICS AND COMPUTING↗

Cascade of vestigial orders in two-component superconductors: Nematic, ferromagnetic, s -wave charge- 4 e , and d -wave charge- 4 e states

Electronically ordered states that break multiple symmetries can melt in multiple stages, similarly to liquid crystals. In a partially melted phase, known as vestigial phase, a bilinear made out of combinations of the multiple components of the primary order parameter condenses. Multicomponent superconductors are thus natural candidates for vestigial order since they break both the U⁡(1)-gauge and also time-reversal or lattice symmetries. Here, we use group theory to classify all possible real-valued and complex-valued bilinears of a generic two-component superconductor on a tetragonal or hexagonal lattice. While the more widely investigated real-valued bilinears correspond to vestigial nematic or ferromagnetic order, the little explored complex-valued bilinears correspond to a vestigial charge-4e condensate, which itself can have an underlying s-wave, d x 2 –y 2 -wave, or d xy -wave symmetry. To properly describe the fluctuating regime of the superconducting Ginzburg-Landau action and thus access these competing vestigial phases, we employ both a large-N and a variational method. We show that while vestigial order can be understood as a weak-coupling effect in the large-N approach, it is akin to a moderate-coupling effect in the variational method. Despite these distinctions, both methods yield similar results in wide regions of the parameter space spanned by the quartic Landau coefficients. Specifically, we find that the nematic and ferromagnetic phases are the leading vestigial instabilities, whereas the various types of charge-4⁢e order are attractive albeit subleading vestigial channels. The only exception is for the hexagonal case, in which the nematic and s-wave charge-4e vestigial states are degenerate. Finally, we discuss the limitations of our approach, as well as the implications of our results for the realization of exotic charge-4⁢e states in material candidates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evaluation of GlassNet for physics-informed machine learning of glass stability and glass-forming ability

Glassy materials form the basis of many modern applications, including nuclear waste immobilization, touch-screen displays, and optical fibers, and also hold great potential for future medical and environmental applications. However, their structural complexity and large composition space make design and optimization challenging for certain applications. Of particular importance for glass processing and design is an estimate of a given composition's glass-forming ability (GFA). However, there remain many open questions regarding the underlying physical mechanisms of glass formation, especially in oxide glasses. It is apparent that a proxy for GFA would be highly useful in glass processing and design, but identifying such a surrogate property has proven itself to be difficult. While glass stability (GS) parameters have historically been used as a GFA surrogate, recent research has demonstrated that most of these parameters are not accurate predictors of the GFA of oxide glasses. Here, in this work, we explore the application of an open-source pre-trained neural network model, GlassNet, that can predict the characteristic temperatures necessary to compute GS with reasonable performance and assess the feasibility of using these physics-informed machine learning (PIML)-predicted GS parameters to estimate GFA. In doing so, we track the uncertainties at each step of the computation—from the original ML prediction errors to the compounding of errors during GS estimation, and finally to the final estimation of GFA. While GlassNet exhibits reasonable accuracy on all individual properties, we observe a large compounding of error in the combination of these individual predictions for the PIML prediction of GS, finding that random forest models offer similar accuracy to GlassNet. We also break down the performance of GlassNet on different glass families and find that the error in GS prediction is correlated with the error in crystallization peak temperature prediction. Lastly, we utilize this finding to assess the relationship between top-performing GS parameters and GFA for two ternary glass systems: sodium borosilicate and sodium iron phosphate glasses. We conclude that to obtain true ML predictive capability of GFA, significantly more data needs to be collected.

36 MATERIALS SCIENCE↗

Selective Deposition and Fusion of AISI 316L: An Additive Manufacturing Process for Space Environments via Direct Ink Writing and Laser Processing

Unlocking the potential of additive manufacturing (AM) for space exploration hinges on overcoming key challenges, notably the ability to manufacture or repair parts on-site during exploration missions with consideration of quality, feedstock utilization, and challenges involved in microgravity environments. While there are multiple efforts to investigate the use of existing metal AM processes such as powder bed fusion (PBF), directed energy deposition (DED), and filament-based material extrusion, each process comes with a different set of challenges in space environments. Here, in this work, we introduce a new AM method that integrates the benefits of direct ink writing (DIW) to selectively deposit metallic pastes with laser-based processing to locally debind and subsequently melt and fuse metal powder, layer by layer, enabling the manufacturing of AISI 316L samples with densities exceeding 99.0%. The impact of process parameters on single-track dimensions, surface morphology, and porosity was characterized. The efficacy of laser debinding was assessed via secondary-ion mass spectrometry, permitting the carbon content to be estimated at 0.0152%, which is safely below the acceptable limit (0.03 wt%) for AISI 316L.

36 MATERIALS SCIENCE↗

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for γ astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of γ cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of γ absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift z=2 and to constrain or detect γ halos up to intergalactic-magnetic-field strengths of at least 0.3 pG . Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from γ astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of γ cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Neural networks for modeling and control of particle accelerators

Charged particle accelerators support a wide variety of scientific, industrial, and medical applications. They range in scale and complexity from systems with just a few components for beam acceleration and manipulation, to large scientific user facilities that span many kilometers and have hundreds-to-thousands of individually-controllable components. Specific operational requirements must be met by adjusting the many controllable variables of the accelerator. Meeting these requirements can be challenging, both in terms of the ability to achieve specific beam quality metrics in a reliable fashion and in terms of the time needed to set up and maintain the optimal operating conditions. One avenue toward addressing this challenge is to incorporate techniques from the fields of machine learning (ML) and artificial intelligence (AI) into the way particle accelerators are modeled and controlled. While many promising approaches within AI/ML could be used for particle accelerators, this dissertation focuses on approaches based on neural networks. Neural networks are particularly well-suited to modeling, control, and diagnostic analysis of nonlinear systems, as well as systems with large parameter spaces. They are also very appealing for their ability to process high-dimensional data types, such as images and time series (both of which are ubiquitous in particle accelerators). In this work, key studies that demonstrated the potential utility of modern neural network-based approaches to modeling and control of particle accelerators are presented. The context for this work is important: at the start of this work in 2012, there was little interest in AI/ML in the particle accelerator community, and many of the advances in neural networks and deep learning that enabled its present success had not yet been made at that time. As such, this work was both an exploration of possible application areas and a generator of proof-of-concept demonstrations in these areas.

43 PARTICLE ACCELERATORS↗

The Juno mission as a probe of long-range new physics

Orbits of celestial objects, especially the geocentric and heliocentric ones, have been well explored to constrain new long-range forces beyond the Standard Model (SM), often referred to as fifth forces. In this paper, for the first time, we apply the motion of a spacecraft around Jupiter to probe fifth forces that don’t violate the equivalence principle. The spacecraft is the Juno orbiter, and ten of its early orbits already allow a precise determination of the Jovian gravitational field. We use the shift in the precession angle as a proxy to test non-gravitational interactions between Juno and Jupiter. Requiring that the contribution from the fifth force does not exceed the uncertainty of the precession shift inferred from data, we find that a new parameter space with the mass of the fifth-force mediator around 10 −14 eV is excluded at 95% C.L.

new light particles↗

Entanglement rearrangement in self-consistent nuclear structure calculations

Background: Entanglement plays a central role in a diverse array of increasingly important research areas, including quantum computation, simulation, measurement, sensing, and communication. Extensive suites of investigations have been performed to better understand entanglement in atomic and molecular quantum many-body systems, while the exploration of entanglement in the structure of nuclei and their reactions is presently in its infancy. Purpose: The goal of this work is to begin investigating the entanglement properties of nuclei from first-principles nuclear many-body calculations. We attempt to identify common features and emergent structures of entanglement that could ultimately lead to new and natural many-body schemes. With an eye toward quantum accelerators in future hybrid-supercomputers, criteria for partitioning nuclear many-body calculations into quantum and classical components may provide advantages in future large-scale computations. Along the way we look for explanations of the relative success of phenomenological models such as the nuclear shell model, and for better ways to match to low-energy nuclear effective field theories and lattice QCD calculations to nuclear many-body techniques that are based upon entanglement. Method: We explore the entanglement between single-particle states in 4 He and 6 He. The patterns of entanglement emerging from different single-particle bases are compared, and possible links with the convergence of observables are explored, in particular, ground-state energies. The nuclear wave functions are obtained by performing active-space no-core configuration-interaction calculations using a two-body nucleon-nucleon interaction derived from chiral effective field theory. Entanglement measures within single-particle bases exhibiting different degrees of complexity are determined, in particular, harmonic oscillator (HO), Hartree-Fock (HF), natural (NAT) and variational natural (VNAT) bases. Specifically, single-orbital entanglement entropy, two-orbital mutual information, and negativity are studied. Results: The entanglement structures in 4 He and 6 He are found to be more localized within NAT and VNAT bases than within a HO basis for the optimal HO parameters we have worked with. In particular a core-valence structure clearly emerges from the full no-core calculation of 6 He. Here, the two-nucleon mutual information shows that the VNAT basis, which typically exhibits good convergence properties, effectively decouples the active and inactive spaces. Conclusions: Measures of one- and two-nucleon entanglement are found to be useful in analyzing the structure of nuclear wave functions, in particular the efficacy of basis states, and may provide useful metrics toward developing more efficient schemes for ab initio computations of the structure and reactions of nuclei, and quantum many-body systems more generally.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗