Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized functions”

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 793 records · Page 44

Nuclear-level effective theory of 𝜇→𝑒 conversion: Formalism and applications

Over the next decade new 𝜇→𝑒 conversion searches at Fermilab (Mu2e) and J-PARC (COMET, DeeMe) are expected to advance limits on charged lepton flavor violation (CLFV) by more than four orders of magnitude. Here, by considering the consequence of 𝑃 and 𝐶⁢𝑃 on elastic 𝜇→𝑒 conversion and the structure of possible charge and current densities, we show that rates are governed by six nuclear responses and a single scale, 𝑞/𝑚 𝑁 , where 𝑞 ≈ 𝑚 𝜇 is the momentum transferred from the leptons to the nucleus. To relate this result to microscopic formulations of CLFV, we construct in nonrelativistic effective theory (NRET) the CLFV nucleon-level interaction, pointing out the relevance of the dimensionless scales 𝑦=($^{𝑞⁢𝑏}_2$) 2 ⁢ > |$\overrightarrow{𝑣}_N$⁢|>|$\overrightarrow{𝑣}_μ$⁢|>|$\overrightarrow{𝑣}_T$| , where 𝑏 is the nuclear size, $\overrightarrow{𝑣}_N$ and $\overrightarrow{𝑣}_μ$ are the nucleon and muon intrinsic velocities, and $\overrightarrow{𝑣}_T$ is the target recoil velocity. We discuss previous work, noting the lack of a systematic treatment of the various small parameters. Because the parameter 𝑦 is not small, a proper calculation of 𝜇→𝑒 conversion requires a full multipole expansion of the nuclear response functions, an apparently daunting task with Coulomb-distorted electron partial waves. We demonstrate that the multipole expansion can be carried out to high precision by introducing a simplifying local momentum 𝑞 eff for the electron. Previous work has been limited to simple charge or spin interactions, thereby treating the nucleus effectively as a point particle. We show that such formulations are not compatible with the general form of the 𝜇→𝑒 conversion rate, failing to generate three of the six allowed nuclear response functions. The inclusion of the nucleon velocity $\overrightarrow{𝑣}_N$ yields an NRET with 16 operators and a rate of the general form. Consequently, in the current discovery era for CLFV, it provides the most sensible starting point for experimental analysis, defining what can and cannot be determined about CLFV from the highly exclusive process of 𝜇→𝑒 conversion. Finally, we expand the NRET operator basis to account for the effects of $\overrightarrow{𝑣}_μ$, associated with the muon's lower component, generating corrections to the CLFV coefficients of the point-nucleus response functions. Using advanced shell-model methods, we compute 𝜇→𝑒 conversion rates for a series of experimental targets, deriving bounds on the coefficients of the CLFV operators. These calculations are the first to include a general basis of CLFV operators, full evaluation of the associated nuclear response functions, and an accurate treatment of electron and muon Coulomb effects. We discuss target selection as an experimental “knob” that can be turned to probe the microscopic origins of CLFV. We describe two types of coherence that enhance certain CLFV operators and selection rules that blind elastic 𝜇→𝑒 conversion to others. We discuss the matching of the NRET onto higher level effective field theories, such as those constructed at the light quark level, noting opportunities to build on existing work in direct detection of dark matter. We discuss the relation of 𝜇→𝑒 conversion to 𝜇→𝑒+𝛾 and 𝜇→3⁢𝑒, showing how MEG II and Mu3e results will complement those of Mu2e and COMET. Finally we describe a accompanying script—in Mathematica and Python versions—that can be used to compute 𝜇→𝑒 conversion rates in various nuclear targets for the full set of NRET operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Formation energy puzzle in intermetallic alloys: Random phase approximation fails to predict accurate formation energies

We performed density functional calculations to estimate the formation energies of intermetallic alloys. We used two semilocal approximations, the generalized gradient approximation (GGA) by Perdew-Burke- Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) meta-GGA. In addition, we utilized two nonlocal DFT functionals, the hybrid HSE06, and the state-of-the-art random phase approximation (RPA). The nonlocal functionals such as HSE06 and RPA yield accurate formation energies of binary alloys with completely-filled d-band metals, where semilocal functionals underperform. The accuracy at the nonlocal functionals is greatly reduced when a partially-filled d-band metal is present in an alloy, while PBE-GGA outperforms in these cases. We show that the accurate prediction of formation energies by any DFT method depends on its ability to predict the accurate electronic properties, e.g., valence d-band contribution to the density of states (DOS). The SCAN meta-GGA often corrects the PBE-DOS, however, it does not provide accurate formation energies compared to PBE. This is assumed to be due to the lack of proper error cancellation that should be expected due to the similar bulk nature of both alloys and their constituents, which may improve with the modification of meta-GGA ingredients. RPA yields too negative formation energies of alloys with partially-filled d-band metals. RPA results can be corrected by restoring the exchange-correlation kernel, thereby improving the short-range electron-electron correlation in metallic densities.

36 MATERIALS SCIENCE↗

Tidal Love numbers and Green’s functions in black hole spacetimes

Tidal interactions play a crucial role in deciphering gravitational wave signals emitted by the coalescence of binary systems. They are usually quantified by a set of complex coefficients which include tidal Love numbers, describing the conservative response to an external perturbation. In the static case, these are found to vanish exactly for asymptotically flat black holes in general relativity in four space-time dimensions, and recently they have been generalized to dynamical interactions. In the context of response theory, the retarded Green’s function provides the complete description of the behavior of dynamical systems. Here, in this work we investigate the relation between Love numbers and Green’s functions, and highlight the relevance of radiation reaction effects to their connection. As a special case, we discuss Banados-Teitelboim-Zanelli black holes, where the absence of radiative modes allows us to make a direct link between them.

classical black holes↗

Statistical characteristics of laser interference and its effect on IR optoelectronic observation systems

We consider statistical characteristics of laser interference recorded in the focal region of an IR optoelectronic observation system under field and out-of-field illumination by laser radiation. The applicability of the analytical description of experimental histograms of the laser interference signal distribution at the photodetector array output using the probability density approximations in the form of gamma distribution and the Gaussian function is substantiated. It is shown that the specific form of these functions is determined both by the mean value (mathematical expectation) of the interference signal and by the characteristic parameter M, which depends on the statistical properties of the laser interference. (laser applications and other issues in quantum electronics)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamical edge modes and entanglement in Maxwell theory

Abstract Previous work on black hole partition functions and entanglement entropy suggests the existence of “edge” degrees of freedom living on the (stretched) horizon. We identify a local and “shrinkable” boundary condition on the stretched horizon that gives rise to such degrees of freedom. They can be interpreted as the Goldstone bosons of gauge transformations supported on the boundary, with the electric field component normal to the boundary as their symplectic conjugate. Applying the covariant phase space formalism for manifolds with boundary, we show that both the symplectic form and Hamiltonian exhibit a bulk-edge split. We then show that the thermal edge partition function is that of a codimension-two ghost compact scalar living on the horizon. In the context of a de Sitter static patch, this agrees with the edge partition functions found by Anninos et al. in arbitrary dimensions. It also yields a 4D entanglement entropy consistent with the conformal anomaly. Generalizing to Proca theory, we find that the prescription of Donnelly and Wall reproduces existing results for its edge partition function, while its classical phase space does not exhibit a bulk-edge split.

Physics↗

Optimal function estimation with photonic quantum sensor networks

The problem of optimally measuring an analytic function of unknown local parameters each linearly coupled to a qubit sensor is well understood, with applications ranging from field interpolation to noise characterization. Here we resolve a number of open questions that arise when extending this framework to Mach-Zehnder interferometers and quadrature displacement sensing. In particular, we derive lower bounds on the achievable mean square error in estimating a linear function of either local phase shifts or quadrature displacements. In the case of local phase shifts, these results prove, and somewhat generalize, a conjecture by Proctor []. For quadrature displacements, we extend proofs of lower bounds to the case of arbitrary linear functions. We provide optimal protocols achieving these bounds up to small (multiplicative) constants and describe an algebraic approach to deriving new optimal protocols, possibly subject to additional constraints. Using this approach, we prove necessary conditions for the amount of entanglement needed for any optimal protocol for both local phase and displacement sensing. Published by the American Physical Society 2024

Bringewatt, Jacob↗

Energetics of silicon in the bulk and near surfaces of tungsten: a first-principles study

Abstract Siliconization of the tokamak walls is a candidate method to improve plasma confinement in fusion tokamaks containing tungsten plasma facing components (W PFCs). To understand the interactions of silicon (Si) with W, the Si behavior in bulk W, and near three low-index W surfaces ((100), (110) and (111)) has been investigated using first-principles density functional theory. In bulk W, Si interstitial atoms have a low solution ability and high mobility, and Si atoms can be strongly trapped by W vacancies. The interaction between two Si adatoms is responsible for the stability of adatom superstructures on W surfaces, consistent with previous experimental observation (Tsong and Casanova 1981 Phys. Rev. Lett. 47 113). Although the coverage dependence of Si adsorption and diffusion energetics on surfaces is related to surface orientation, the W(110) surface has lower Si adsorption affinity and higher Si diffusivity than either the W(111) or W(100) surfaces. The most stable Si adatom superstructure on W surfaces is: square c(2 × 2) pattern on W(100) covered with 0.5 ML Si; rectangular c(4 × 2) pattern on W(110) with 0.25 ML Si; and rhombus p(1 × 1) pattern on W(111) with 1 ML Si. The coverage dependence of Si mobility on/toward W surfaces is generally related to the stability of the Si superstructures as a function of coverage on each surface. Interestingly, Si adatoms prefer to transport below the surface and into W subsurface by an exchange mechanism with W atoms, indicating the likelihood of epitaxial growth of W silicide layers on W surfaces during the operation of W PFCs.

Yang, L. (ORCID:0000000322166071)↗

Microscopic Imprints of Learned Solutions in Tunable Networks

In physical networks trained using supervised learning, physical parameters are adjusted to produce desired responses to inputs. An example is an electrical contrastive local learning network of nodes connected by edges that adjust their conductances during training. When an edge conductance changes, it upsets the current balance of every node. In response, physics adjusts the node voltages to minimize the dissipated power. Learning in these systems is therefore a coupled double-optimization process, in which the network descends both a cost landscape in the high-dimensional space of edge conductances and a physical landscape—the power dissipation—in the high-dimensional space of node voltages. Because of this coupling, the physical landscape of a trained network contains information about the learned task. Here, we derive a structure-function relation for trained tunable networks and demonstrate that all the physical information relevant to the trained input-output relation can be captured by a tuning susceptibility, an experimentally measurable quantity. We supplement our theoretical results with simulations to show that the tuning susceptibility is correlated with functional importance and that we can extract physical insight into how the system performs the task from the conductances of highly susceptible edges. Our analysis is general and can be applied directly to mechanical networks, such as networks trained for protein-inspired function such as allostery.

36 MATERIALS SCIENCE↗

Functional characterization of diverse type I-F CRISPR-associated transposons

CRISPR-Cas systems generally provide adaptive immunity in prokaryotes through RNA-guided degradation of foreign genetic elements like bacteriophages and plasmids. Recently, however, transposon-encoded and nuclease-deficient CRISPR-Cas systems were characterized and shown to be co-opted by Tn7-like transposons for CRISPR RNA-guided DNA transposition. As a genome engineering tool, these CRISPR-Cas systems and their associated transposon proteins can be deployed for programmable, site-specific integration of sizable cargo DNA, circumventing the need for DNA cleavage and homology-directed repair involving endogenous repair machinery. Here, we selected a diverse set of type I-F3 CRISPR-associated transposon systems derived from Gammaproteobacteria, predicted all components essential for transposition activity, and deployed them for functionality testing within Escherichia coli. Our results demonstrate that these systems possess a significant range of integration efficiencies with regards to temperature, transposon size, and flexible PAM requirements. Additionally, our findings support the categorization of these systems into functional compatibility groups for efficient and orthogonal RNA-guided DNA integration. This work expands the CRISPR-based toolbox with new CRISPR RNA-guided DNA integrases that can be applied to complex and extensive genome engineering efforts.

59 BASIC BIOLOGICAL SCIENCES↗

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Performant implementation of the atomic cluster expansion (PACE) and application to copper and silicon

The atomic cluster expansion is a general polynomial expansion of the atomic energy in multi-atom basis functions. Here we implement the atomic cluster expansion in the performant C++ code that is suitable for use in large-scale atomistic simulations. We briefly review the atomic cluster expansion and give detailed expressions for energies and forces as well as efficient algorithms for their evaluation. We demonstrate that the atomic cluster expansion as implemented in shifts a previously established Pareto front for machine learning interatomic potentials toward faster and more accurate calculations. Moreover, general purpose parameterizations are presented for copper and silicon and evaluated in detail. We show that the Cu and Si potentials significantly improve on the best available potentials for highly accurate large-scale atomistic simulations.

36 MATERIALS SCIENCE↗

Performant implementation of the atomic cluster expansion

The atomic cluster expansion is a general polynomial expansion of the atomic energy in multi-atom basis functions. Here we implement the atomic cluster expansion in the performant C++ code PACE that is suitable for use in large scale atomistic simulations. We briefly review the atomic cluster expansion and give detailed expressions for energies and forces as well as efficient algorithms for their evaluation. We demonstrate that the atomic cluster expansion as implemented in PACE shifts a previously established Pareto front for machine learning interatomic potentials towards faster and more accurate calculations. Moreover, general purpose parameterizations are presented for copper and silicon and evaluated in detail. We show that the new Cu and Si potentials significantly improve on the best available potentials for highly accurate large-scale atomistic simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗

State Preparation of Antisymmetrized Geminal Power on a Quantum Computer without Number Projection

The antisymmetrized geminal power (AGP) is equivalent to the number projected Bardeen–Cooper–Schrieffer (PBCS) wave function. It is also an elementary symmetric polynomial (ESP) state. We generalize previous research on deterministically implementing the Dicke state to a state preparation algorithm for an ESP state, or equivalently AGP, on a quantum computer. Our method is deterministic and has polynomial cost, and it does not rely on number symmetry breaking and restoration. We also show that our circuit is equivalent to a disentangled unitary paired coupled cluster operator and a layer of unitary Jastrow operator acting on a single Slater determinant. Here, the method presented herein highlights the ability of disentangled unitary coupled cluster to capture nontrivial entanglement properties that are hardly accessible with traditional Hartree–Fock based electronic structure methods.

97 MATHEMATICS AND COMPUTING↗

Generalized quasiharmonic approximation via space group irreducible derivatives

The quasiharmonic approximation (QHA) is the simplest nontrivial approximation for interacting phonons under constant pressure, bringing the effects of anharmonicity into temperature-dependent observables. Nonetheless, the QHA is often implemented with additional approximations due to the complexity of computing phonons under arbitrary strains, and the generalized QHA, which employs constant stress boundary conditions, has not been completely developed. In this work we formulate the generalized QHA, providing a practical algorithm for computing the strain state and other observables as a function of temperature and true stress. We circumvent the complexity of computing phonons under arbitrary strains by employing irreducible second-order displacement derivatives of the Born-Oppenheimer potential and their strain dependence, which are efficiently and precisely computed using the lone irreducible derivative approach. We formulate two complementary strain parametrizations: a discretized strain grid interpolation and a Taylor series expansion in symmetrized strain. We illustrate our approach by evaluating the temperature and pressure dependence of select elastic constants and the thermal expansion in thoria (ThO 2 ) using density functional theory with three exchange-correlation functionals. The QHA results are compared to our measurements of the elastic constant tensor using time-domain Brillouin scattering and inelastic neutron scattering. Our irreducible derivative approach simplifies the implementation of the generalized QHA, which will facilitate reproducible, data-driven applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dissociative electron attachment to amide bond containing molecules: N-ethylformamide and N-ethylacetamide

To advance our quest to understand the role of low energy electrons in biomolecular systems, we performed investigations on dissociative electron attachment (DEA) to gas-phase N-ethylformamide (NEF) and N-ethylacetamide (NEA) molecules. Both molecules contain the amide bond, which is the linkage between two consecutive amino acid residues in proteins. Thus, their electron-induced dissociation can imitate the resonant behavior of the DEA process in more complex biostructures. Our experimental results indicate that in these two molecules, the dissociation of the amide bond results in a double resonant structure with peaks at ~5 eV and 9 eV. We also determined the energy position of resonant states for several negative ions, i.e., the other dissociation products from NEF and NEA. Our predictions of dissociation channels were supported by density functional theory calculations of the corresponding threshold energies. Our results and those previously reported for small amides and peptides imply the fundamental nature for breakage of the amide bond through the DEA process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Page curve for reflected entropy

We study the reflected entropy $S_{R}$ in the West Coast Model, a toy model of black hole evaporation consisting of JT gravity coupled to end-of-the-world branes. We demonstrate the validity of the holographic duality relating it to the entanglement wedge cross section away from phase transitions. Further, we analyze the important non-perturbative effects that smooth out the discontinuity in the phase transition. By performing the gravitational path integral, we obtain the reflected entanglement spectrum analytically. The spectrum takes a simple form consisting of superselection sectors, which we interpret as a direct sum of geometries, a disconnected one and a connected one involving a closed universe. We find that area fluctuations of $\textit{O}$$\sqrt{G_N}$ spread out the $S_{R}$ phase transition in the canonical ensemble, analogous to the entanglement entropy phase transition. We also consider a Renyi generalization of the reflected entropy and show that the location of the phase transition varies as a function of the Renyi parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Effect of the Prior and the Experimental Design on the Inference of the Precision Matrix in Gaussian Chain Graph Models

Here, we investigate whether (and how) experimental design could aid in the estimation of the precision matrix in a Gaussian chain graph model, especially the interplay between the design, the effect of the experiment and prior knowledge about the effect. Estimation of the precision matrix is a fundamental task to infer biological graphical structures like microbial networks. We compare the marginal posterior precision of the precision matrix under four priors: flat, conjugate Normal-Wishart, Normal-MGIG and a general independent. Under the flat and conjugate priors, the Laplace-approximated posterior precision is not a function of the design matrix rendering useless any efforts to find an optimal experimental design to infer the precision matrix. In contrast, the Normal-MGIG and general independent priors do allow for the search of optimal experimental designs, yet there is a sharp upper bound on the information that can be extracted from a given experiment. We confirm our theoretical findings via a simulation study comparing (i) the KL divergence between prior and posterior and (ii) the Stein’s loss difference of MAPs between random and no experiment. Our findings provide practical advice for domain scientists conducting experiments to better infer the precision matrix as a representation of a biological network.

54 ENVIRONMENTAL SCIENCES↗

Flexible Service Contracting for Risk Management within Integrated Transmission and Distribution Systems

The general objective of our project has been to investigate the ability of Independent Distribution System Operators (IDSOs), functioning as linkage agents for Integrated Transmission and Distribution (ITD) systems, to facilitate the flexible availability and usage of reserve in support of ITD system operations. This general objective is in accordance with Order 2222 of the U.S. Federal Energy Regulatory Commission (FERC), titled “Participation of Distributed Energy Resource Aggregations in Markets Operated by Regional Transmission Organizations and Independent System Operators”. The Final Rule for FERC Order 2222 was issued on September 17, 2020. The primary contribution of our project is that we have formulated an innovative energy management Transactive Energy Design (TES) approach for ITD systems that provides promising support for our project objective. Specifically, we have developed a new type of IDSO-managed TES design for distribution system operations, as well as new types of swing contracts permitting IDSOs to participate in transmission system operations as providers of reserve harnessed from distribution system resources in return for appropriate compensation. Together, these design elements constitute a scalable market-based approach facilitating efficient reserve procurement for ITD system operations from a fuller range of power resources. The efficacy of our approach has been demonstrated by means of detailed conceptual analyses as well as test case simulations. To implement the latter, we have developed the ITD TES Platform V2.0, a computational platform that permits the modeling and software implementation of ITD system operations over time. The Electric Reliability Council of Texas (ERCOT) energy region has been used as the empirical anchor for the development of this platform. Our conceptual and test-case work has been reported in a Wiley/IEEE Press book, two refereed book chapters, and seven refereed journal articles. The key components comprising the ITD TES Platform V2.0 have been released as documented open-source software at online GitHub repositories.

24 POWER TRANSMISSION AND DISTRIBUTION↗