SEARCH · Engineering Papers

Results for “density matrix methods”

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

SourceResearch connectionFrom yearThrough year

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

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2024

R-matrix calculations for opacities: I. Methodology and computations

Abstract An extended version of the R -matrix methodology is presented for calculation of radiative parameters for improved plasma opacities. Contrast and comparisons with existing methods primarily relying on the distorted wave approximation are discussed to verify accuracy and resolve outstanding issues, particularly with reference to the opacity project (OP). Among the improvements incorporated are: (i) large-scale Breit–Pauli R -matrix calculations for complex atomic systems including fine structure, (ii) convergent close coupling wave function expansions for the ( e + ion) system to compute oscillator strengths and photoionization cross sections, (iii) open and closed shell iron ions of interest in astrophysics and experiments, (iv) a treatment for plasma broadening of autoionizing resonances as function of energy-temperature-density dependent cross sections, (v) a ‘top-up’ procedure to compare convergence with R -matrix calculations for highly excited levels, and (vi) spectroscopic identification of resonances and bound ( e + ion) levels. The present R -matrix monochromatic opacity spectra are fundamentally different from OP and lead to enhanced Rosseland and Planck mean opacities. An outline of the work reported in other papers in this series and those in progress is presented. Based on the present re-examination of the OP work, opacities of heavy elements might require revisions in high temperature-density plasma sources.

Pradhan, A. K. (ORCID:0000000187753643)↗

DOE OSTI2025

Holographic View of Mixed-State Symmetry-Protected Topological Phases in Open Quantum Systems

We establish a holographic duality between d -dimensional mixed-state symmetry-protected topological (mSPT) phases and ( d + 1 ) -dimensional subsystem symmetry-protected topological (SSPT) states. Specifically, we show that the reduced density matrix of the boundary layer of a ( d + 1 ) -dimensional SSPT state with subsystem symmetry S and global symmetry G corresponds to a d -dimensional mSPT phase with strong S and weak G symmetries. Conversely, we demonstrate that the wave function of an SSPT state can be constructed by replicating the density matrix of the corresponding lower-dimensional mSPT phase. This mapping links the density matrix in lower dimensions to the entanglement properties of higher-dimensional wave functions, providing an approach for analyzing nonlinear quantities and quantum information metrics in mixed-state systems. Our duality offers a new perspective for studying intrinsic mSPT phases that are unique to open quantum systems, without pure-state analogues. We show that strange correlators and twisted Rényi- N correlators can diagnose these nontrivial phases and we explore their connection to strange correlators in pure-state SSPT phases. Furthermore, we discuss several implications of this holographic duality, including a method for preparing intrinsic mSPT phases through the duality. Published by the American Physical Society 2025

Sun, Shijun (ORCID:0000000168880030)↗

DOE OSTI2020

Magnetic properties of γ-Fe 2 O 3 nanoparticles in a porous SiO 2 shell for drug delivery

A method is presented for synthesizing core–shell nanoparticles with a magnetic core and a porous shell suitable for drug delivery and other medical applications. The core contains multiple γ -Fe 2 O 3 nanoparticles (~15 nm) enclosed in a SiO 2 (~100–200 nm) matrix using either methyl (denoted TMOS- γ -Fe 2 O 3 ) or ethyl (TEOS- γ -Fe 2 O 3 ) template groups. Low-temperature Mössbauer spectroscopy showed that the magnetic nanoparticles have the maghemite structure, γ -Fe 2 O 3 , with all the vacancies in the octahedral sites. Saturation magnetization measurements revealed that the density of γ -Fe 2 O 3 was greater in the TMOS- γ -Fe 2 O 3 nanoparticles than TEOS- γ -Fe 2 O 3 nanoparticles, presumably because of the smaller methyl group. Magnetization measurements showed that the blocking temperature is around room temperature for the TMOS- γ -Fe 2 O 3 and around 250 K for the TEOS- γ -Fe 2 O 3 . Three dimensional topography analysis shows clearly that the magnetic nanoparticles are not only at the surface but have penetrated deep in the silica to form the core–shell structure.

Physics↗

DOE OSTI2020

Spontaneous Raman–LIF–CO–OH measurements of species concentration in turbulent spray flames

This paper presents new measurements of species concentrations, temperature and mixture fraction in selected regions of a turbulent ethanol spray flame. The line-Raman–LIF–COsingle bondOH setup developed at the Sandia's Combustion Research Facility is utilised to probe regions of a spray flame where laser breakdown of liquid droplets is avoided and the remaining interferences can be corrected. The spray flame is stabilised on the piloted Sydney needle spray burner, where axial translation of the liquid injecting needle in the air-blast stream can transition the spray from dilute to dense. The solution to obtaining successful measurements is found to be multifaceted and includes: the appropriate selection of flame conditions; high sensitivity of the Raman detection system permitting reduced laser energies; development of a pre-processing algorithm to reject strong droplet interferences; and application of the hybrid matrix inversion method combined with wavelet denoising to account for interference corrections and noise at the very low signal levels obtained. Unique and necessary for the successful measurements reported in this paper, a pre-processing algorithm is outlined that removes data points corrupted with strong interferences from droplets. These interferences arise from a range of sources, but the most intense are due to the laser interaction with surrounding mist or liquid fragments, such that measurements near the jet centreline are corrupted and hence discarded. Reliable measurements of mixture fraction, temperature obtained from the sum of the species number densities, and species mole fractions are reported for regions in the flames sufficiently far from the centreline. The paper demonstrates the feasibility of the judicious use of Raman scattering in turbulent spray flames, the results of which will be extremely useful for validating numerical simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Structure and Sulfur: Tuning the Viscoelastic and Surface Properties of Natural Keratin Fibers

Natural keratin fibers, such as wool, possess a complex hierarchical structure that governs their mechanical properties and surface energy. However, the extent to which these characteristics are influenced by combined contributions of structural variations (e.g., fiber diameter, intermediate filament (IF) packing) and chemical composition (e.g., disulfide bond density) remains poorly understood. In this study, we investigate wool fibers from five sheep breeds (Merino, Polwarth, Cheviot, Eider, and Devon) to elucidate how these factors influence viscoelasticity and surface interactions. Using a multimodal approach integrating interfacial and bulk characterization methods, including inverse gas chromatography (IGC), atomic force microscopy-infrared spectroscopy (AFM-IR), X-ray photoelectron spectroscopy (XPS), uniaxial tensile testing, and synchrotron small-angle X-ray scattering (SAXS), we show that the nanometer-thick 18-methyleicosanoic acid (18-MEA) layer is consistently present across all wool types and plays a key role in governing hydrophobicity and surface heterogeneity. A controlled isothermal treatment at 200 °C, designed to cleave disulfide bonds, results in a nearly 40% reduction in specific surface area across all fiber types, accompanied by a significant decrease in tensile strength and 80% reduction in elongation at break for Merino and Devon wool, but limited influence on the mechanical properties of Eider fibers. Furthermore, rate-dependent tensile testing within the elastic regime reveals distinct viscoelastic responses among the fiber types, suggesting that the sulfur-rich protein matrix surrounding IFs and its structure contribute actively to stress partitioning. Altogether, when combined with conclusions from SAXS measurements of IF spacing, our work offers compelling insights into the role of the keratin-associated protein (KAP) matrix in shaping wool fiber mechanics. Differences in mechanical behavior among wool types, despite similar IF spacing or sulfur content, highlight the importance of matrix composition and cross-linking density, suggesting that the molecular architecture of the KAP network may be a dominant factor in determining fiber performance.

X-ray scattering↗

DOE OSTI2021

Tunable plasmon-enhanced second-order optical nonlinearity in transition metal dichalcogenide nanotriangles

The development of nanomaterials with a large nonlinear susceptibility is essential for nonlinear nanophotonics. Here we show that transition metal dichalcogenide (TMD) nanotriangles have a large effective second-order susceptibility [χ (2) ] at midinfrared to near-infrared frequencies owing to their broken centrosymmetry. χ (2) is calculated within the density-matrix formalism that accounts for dissipation and screening. χ (2) peaks in the vicinity of both two-photon resonances (specified by the geometry) and plasmon resonances (tunable via the carrier density). Aligning the resonances yields the values of χ (2) as high as 10 –6 m/V. These findings underscore the potential of TMD nanotriangles for nonlinear nanophotonics, particularly second-harmonic generation.

0-dimensional systems↗

DOE OSTI2021

Quantifying and Controlling Entanglement in the Quantum Magnet Cs 2 CoCl 4

The lack of methods to experimentally detect and quantify entanglement in quantum matter impedes our ability to identify materials hosting highly entangled phases, such as quantum spin liquids. We thus investigate the feasibility of using inelastic neutron scattering (INS) to implement a model-independent measurement protocol for entanglement based on three entanglement witnesses: one-tangle, two-tangle, and quantum Fisher information (QFI). We perform high-resolution INS measurements on Cs 2 CoCl 4 , a close realization of the S = 1/2 transverse-field XXZ spin chain, where we can control entanglement using the magnetic field, and compare with density-matrix renormalization group calculations for validation. The three witnesses allow us to infer entanglement properties and make deductions about the quantum state in the material. Overall, we find QFI to be a particularly robust experimental probe of entanglement, whereas the one and two-tangles require more careful analysis. Our results lay the foundation for a general entanglement detection protocol for quantum spin systems.

1-dimensional spin chains↗

DOE OSTI2024

Quantum chaos on edge

Recently, the physics of many-body quantum chaotic systems close to their ground states has come under intensified scrutiny. Such studies are motivated by the emergence of model systems exhibiting chaotic fluctuations throughout the entire spectrum [the Sachdev-Ye-Kitaev (SYK) model being a renowned representative] as well as by the physics of holographic principles, which likewise unfold close to ground states. Interpreting the edge of the spectrum as a quantum critical point, here we combine a wide range of analytical and numerical methods to the identification and comprehensive description of two different universality classes: the near edge physics of “sparse” and the near edge of “dense” chaotic systems. The distinction lies in the ratio between the number of a system's random parameters and its Hilbert space dimension, which is exponentially small or algebraically small in the sparse and dense case, respectively. Notable representatives of the two classes are generic chaotic many-body models (sparse) and invariant random matrix ensembles or chaotic gravitational systems (dense). While the two families share identical spectral correlations at energy scales comparable to the level spacing, the density of states and its fluctuations near the edge are different. Considering the SYK model as a representative of the sparse class, we apply a combination of field theory and exact diagonalization to a detailed discussion of its edge spectrum. Conversely, Jackiw-Teitelboim gravity is our reference model for the dense class, where an analysis of the gravitational path integral and random matrix theory reveal universal differences to the sparse class, whose implications for the construction of holographic principles we discuss. Published by the American Physical Society 2024

Altland, Alexander (ORCID:0000000229914805)↗

DOE OSTI2022

Revisiting the Hail Radar Reflectivity–Kinetic Energy Flux Relation by Combining T-Matrix and Discrete Dipole Approximation Calculations to Size Distribution Observations

The retrieval of hail kinetic energy with weather radars or its simulation in numerical models is challenging because of the shape complexity and variable density of hailstones. We combine 3D scans of individual hailstones with measurements of the particle size distributions (PSD) and T-matrix calculations to understand how hail reflectivity Z changes when approximating hailstones as spheroids, as compared to the realistic shapes obtained by 3D scanning technology. Additionally, recent terminal velocity relations are used to compare Z to the hail kinetic energy flux E ˙ . We parameterize the hail backscattering cross sections at L, S, C, and X bands as a function of size between 0.5 and 5.0 cm, matching the range of the observed PSDs. The scattering calculations use the T-matrix method for size parameters below 1.0 and the discrete dipole approximation (DDA) method otherwise. The DDA calculations are done for 48 digital models of realistic hailstones of sizes between 1 and 5 cm. The DDA cross sections are calculated for multiple orientations and averaged assuming a fully random orientation distribution to provide a single value per hailstone. The T-matrix reflectivity assuming solid ice spheres presents negligible differences to DDA results for size parameters below 1.0. Therefore, T matrix was used to fill in the gaps left by the DDA calculations. The results are mapped to the same size bins of the observed PSDs, allowing the calculation of the radar reflectivity. This is then correlated to E ˙ , allowing a potential improvement of past retrieval methods of E ˙ from Z in multiple wavelengths.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Physics-Informed Machine Learning Model for Ceramic Matrix Composite Creep

A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.

ceramic matrix composites↗

DOE OSTI2021

Reactivity of CO 2 with Utica, Marcellus, Barnett, and Eagle Ford Shales and Impact on Permeability

We report that in order to reduce greenhouse gas emissions while recovering hydrocarbons from unconventional shale formations, processes that make use of carbon dioxide to enhance oil recovery while storing carbon dioxide (CO 2 ) should be considered. Here, we examine samples from three shale basins across the United States (Utica and Marcellus Shales in the Appalachian Basin, Barnett Shale in the Bend Arch-Ft. Worth Basin, and Eagle Ford in the Western Gulf Basin) to address the following questions: (1) do changes from reaction with CO 2 and fluids at the micrometer and nanometer scale alter flow pathways and, in turn, impact hydrocarbon production, CO 2 storage, and seal integrity and (2) can CO 2 or fluid reactivity be predicted based on physical or chemical properties of shale formations? Experiments were conducted at 40 °C and 10.3 MPa to characterize the interaction between CO 2 and shale using X-ray diffraction (XRD), carbon and sulfur analysis, in situ Fourier transform infrared spectroscopy (FT-IR), feature relocation scanning electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), mercury (Hg) intrusion porosimetry, and Brunauer-Emmett-Teller (BET) surface area and pore size analysis coupled with density functional theory (DFT) methods. Changes in mechanical, physical, and flow properties of shale cores due to CO 2 exposure were addressed using a New England Research Autolab 1500 and Xenon X-ray computed tomography (CT) scanning. Results showed that CO 2 did not promote significant reactivity with the shale if water was not present; only shales with swelling clays or residual interstitial pore water reacted with dry CO 2 to promote reactivity in shale. When water was added as a reactant, CO 2 formed carbonic acid and reacted with the shale to dissolve carbonate pockets, etched and pitted the shale matrix surfaces, and increased the microporosity and decreased nanoporosity. Porosity and permeability increased appreciably in core shale samples after exposure to CO 2 saturated fluid due to dissolution of carbonate. Shale mechanical properties were not altered. Trends were not observed that could tie CO 2 or fluid reactivity to physical or chemical properties of the shale formations at the basin scale from the samples we examined. However, if the shale contained significant amounts of carbonate and water was available to react with the CO 2 , pore sizes were altered in the matrix and permeability and porosity increased.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2025

Extending Orbital-Optimized Density Functional Theory to L-Edge XPS and Beyond: Spin–Orbit Coupling via Nonorthogonal Quasi-Degenerate Perturbation Theory

Quantum mechanical calculations of core electron binding energies (CEBEs) are relevant to interpreting X-ray photoelectron spectroscopy (XPS). Orbital-optimized density functional theory (OO-DFT) accurately predicts K-edge CEBEs but is challenged by the presence of significant spin–orbit coupling (SOC) at L- and higher edges involving inner-shell orbitals with nonzero angular momentum. Here, to extend OO-DFT to L-edges and higher, our method utilizes scalar-relativistic, spin-restricted open-shell OO-DFT to construct a minimal, quasi-degenerate basis of core-hole states corresponding to a chosen inner-shell (e.g., ionizing all six possible 2p spin orbitals). Nonorthogonal configuration interaction (NOCI) is then used to obtain the matrix elements of the full Hamiltonian including SOC in this quasi-degenerate model space of determinants. Using a screened 1-electron SOC operator parametrized with the Dirac-Coulomb-Breit (DCB) Hamiltonian results in doublet splitting (DS) values for third row elements that are nearly in quantitative agreement with experiment. The resulting NOCI eigenvalues are shifted by the average of the (scalar) OO-DFT CEBEs to yield CEBEs (split by SOC) corrected for dynamic correlation. Comparing calculations on gas phase molecules with experimental results establishes that NO-QDPT with the SCAN functional (NO-QDPT/SCAN), using the DCB screened 1-electron SOC operator is accurate to about 0.2 eV for L-edge CEBEs of molecules containing third row atoms. However, this NO-QDPT approach becomes less accurate for fourth-row elements starting in the middle of the 3d transition metal series, with errors increasing as atomic number increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Locally purified maximally mixed states at scale: Entanglement pruning and symmetries

Locally Purified Density Operators (LPDOs) are state-of-the-art tensor network ansatze candidates that efficiently represent mixed quantum states at scale. However, given their non-uniqueness, their representational complexity is generally sub-optimal in practical computations. Here, in this work we perform a comprehensive numerical and analytical analysis and resolve this issue in the experimentally relevant limit where noise depolarizes the density operator into a maximally mixed state. To resolve the sub-optimality issue, we analyze two numerical tools, one analytic method, and detail the relations between them. The numerical tools used are fidelity-preserving truncations and isometric gauge transformations leveraging Riemannian optimizations over entropic objective functions. In addition, by invoking the injectivity and symmetry constraints of the maximally mixed LPDO, we also present analytical closed-form expressions for the disentangler and discuss their relation to numerical optimizers. Further, away from the maximally mixed state, our simulations highlight how the truncation threshold smoothly interpolate, as a function of depolarization, between established matrix product results and our new results. Our work shows how, by minimizing the resources required to represent key states of practical interest in experiment, the efficiency of tensor network algorithms can be substantially increased. This paves the path for uncovering tensor network’s fundamental scalability limits and latent potential in representing the wide locus of mixed quantum states that are accessible on near-term quantum devices.

Gangapuram, Amit Jamadagni [Oak Ridge National Lab↗

DOE OSTI2024

Noise reduction of stochastic density functional theory for metals

Density Functional Theory (DFT) has become a cornerstone in the modeling of metals. However, accurately simulating metals, particularly under extreme conditions, presents two significant challenges. First, simulating complex metallic systems at low electron temperatures is difficult due to their highly delocalized density matrix. Second, modeling metallic warm-dense materials at very high electron temperatures is challenging because it requires the computation of a large number of partially occupied orbitals. This study demonstrates that both challenges can be effectively addressed using the latest advances in linear-scaling stochastic DFT methodologies. Despite the inherent introduction of noise into all computed properties by stochastic DFT, this research evaluates the efficacy of various noise reduction techniques under different thermal conditions. Our observations indicate that the effectiveness of noise reduction strategies varies significantly with the electron temperature. Furthermore, we provide evidence that the computational cost of stochastic DFT methods scales linearly with system size for metal systems, regardless of the electron temperature regime.

Chemistry↗

DOE OSTI2020

Electronic energy gap closure and metal-insulator transition in dense liquid hydrogen

Here, using quantum Monte Carlo (QMC) calculations, we investigate the insulator-metal transition observed in liquid hydrogen at high pressure. Below the critical temperature of the transition from the molecular to the atomic liquid, the fundamental electronic gap closure occurs abruptly, with a small discontinuity reflecting the weak first-order transition in the thermodynamic equation of state. Above the critical temperature, molecular dissociation sets in while the gap is still open. When the gap closes, the decay of the off-diagonal reduced density matrix shows that the liquid enters a gapless, but localized, phase: there is a crossover between the insulating and the metallic liquids. Compared to different density functional theory (DFT) functionals, our QMC calculations provide larger values for the fundamental gap and the electronic density of states close to the band edges, indicating that optical properties from DFT potentially benefit from error cancellations.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Potential absence of observed π2 linear-chain structures in 14O via 10C(α,α) resonances scattering

Background: The preference for light nuclear systems to coagulate into -particle clusters has been well-studied. The possibility of a linear chain configuration of -particles would allow for a new way to study this phenomenon. Purpose: A rotational band of states in C has been claimed showing a linear chain structure. The mirror system, O, has been studied here to examine how this linear chain structure is affected by replacing the valence neutrons with protons. Method: A beam of C was incident into a chamber filled with He:CO gas with the tracks recorded inside the TexAT Time Projection Chamber and the recoil -particles detected by a silicon detector array to measure the cross section. Results: The experimental cross section was compared with previous studies and fit using R-Matrix theory with the previously-observed O states being transformed to the C using mirror symmetry. The measured cross section does not replicate the claimed states, with the predicted cross section exceeding that observed at several energies and angles. Conclusion: A series of possibilities are highlighted with the most likely being that the originally-seen C states did not constitute a rotational band with a potentially incorrect spin assignment due to the limitations of the angular correlation method with non-zero spin particles. The work highlights the difficulties in measuring broad resonances corresponding to a linear chain state in a high level density.

Bishop, J.↗

DOE OSTI2025

Impact of tracking resolution on 𝜙-meson spin alignment measurements

Measurements of global spin alignment of vector mesons in relativistic heavy-ion collisions can provide unique insights into spin-orbit interactions and vector meson dynamics in the quark-gluon plasma produced in those collisions. The global spin alignment is measured by the 00th coefficient of the spin density matrix, 𝜌 00 , via the polar angle (𝜃*) of the decay-kaon momentum in the parent rest frame with respect to the direction of the orbital angular momentum of the collision. Such measurements are affected by the angular and momentum resolutions of the reconstructed tracks in the experiment. Such effects are nontrivial because of kinematic complications caused by the boost to the parent rest frame, and could be important given that the global spin alignment signal is weak. In this paper, we investigate the effects of experimental tracking resolutions on measurements of the 𝜙 (1020) meson 𝜌 00 . We study these effects for two methods of 𝜌 00 measurements, the conventional method analyzing the 𝜙⁢−meson yield versus cos 2 ⁡𝜃* and the invariant mass (𝑚 inv ) method utilizing ⟨cos 2 ⁡𝜃*⟩ versus 𝑚 inv . Using typical resolution values from experiments, we find that the effect of track resolution on 𝜌 00 is small, well within typical measurement uncertainties.

Relativistic heavy-ion collisions↗

DOE OSTI2026

Pareto-optimal target definition for multi-axis random vibration testing

In random vibration testing with multiple control channels, existing control laws require specification of a complete spectral density matrix at each control frequency. Spectral density matrices include autospectral densities on the diagonal and cross-spectral densities on the off-diagonal. In practice, the off-diagonal terms are often unknown, and recent vibration testing research has focused on fixing the diagonal and specifying the off-diagonal to minimize the required control energy, subject to a constraint that the target matrix is positive semidefinite. This paper shows that, even with a fixed diagonal, off-diagonal terms strongly affect control residuals. This overlooked effect occurs in both square and rectangular systems. By jointly considering input energy and control residuals, open-loop inputs are derived directly from the diagonal without specifying the off-diagonal terms. Vibration targets that can be used in closed-loop control are then derived using the optimal inputs, with positive semidefinite constraints applied during the derivation. The result is a set of Pareto-optimal control solutions. For each solution in the set, any other possible solution produces greater control error, greater input energy, or both. A balanced solution is selected automatically, though others can be chosen based on test needs. Simulations and experiments show that the proposed method outperforms state-of-the-art energy-minimizing approaches, achieving significant reductions in both control error and input energy.

Autospectral density↗