Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized spherical 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Identification of crystal plasticity model parameters by multi-objective optimization integrating microstructural evolution and mechanical data

Crystal plasticity models evolve a polycrystalline yield surface using meso-scale descriptions of deformation mechanisms. The activation of deformation mechanisms is governed by crystallography and a set of model parameters, which are typically calibrated through the fitting of mechanical data such as stress–strain curves and elastic lattice strains. Microstructural data such as phase fractions and texture evolution are used for verifying crystal plasticity parameters. In this study, we use a multi-objective genetic algorithm to identify hardening parameters from flow stress curves with an option to incorporate texture into the optimization approach. Robust, generalized objective functions are developed and used to identify sets of parameters pertaining to dislocation density-based hardening laws in visco-plastic and elasto-plastic self-consistent (VPSC and EPSC) homogenization models. First, the parameters are identified for pure Nb directly from texture using an objective function based on generalized spherical harmonics. Since texture evolution is driven by the relative contribution of active slip systems, the parameters governing the evolution of slip resistance ratios can be recovered from fitting discrete textures at a series of strains. Next, a comprehensive set of load reversal data for dual phase (DP) 780 steel is used to fit a hardening law and a back-stress law in EPSC. Finally, parameters pertaining to a complex hardening law for the evolution of slip and twinning in pure α-Ti are identified. Remarkably, using texture as an objective in combination with stress–strain objectives constrains the model of Ti to fully reproduce not only stress–strain and texture evolution but also hierarchical twinning measurements as a function of initial grain size and texture. Furthermore, given an appropriate model fit to representative experimental texture evolution, underlying twin volume fractions contributing to texture evolution can be predicted.

42 ENGINEERING↗
DOE OSTI2022

Analysis of core asymmetries in inertial confinement fusion implosions using three-dimensional hot-spot reconstruction

Three-dimensional effects play a crucial role during the hot-spot formation in inertial confinement fusion (ICF) implosions. A data analysis technique for 3D hot-spot reconstruction from experimental observables has been developed to characterize the effects of low modes on 3D hot-spot formations. In nuclear measurements, the effective flow direction, governed by the maximum eigenvalue in the velocity variance of apparent ion temperatures, has been found to agree with the measured hot-spot flows for implosions dominated by mode ℓ = 1. Asymmetries in areal-density (ρR) measurements were found to be characterized by a unique cosine variation along the hot-spot flow axis. In x-ray images, a 3D hot-spot x-ray emission tomography method was developed to reconstruct the 3D hot-spot plasma emissivity using a generalized spherical-harmonic Gaussian function. The gradient-descent algorithm was used to optimize the mapping between the projections from the 3D hot-spot emission model and the measured x-ray images along multiple views. Furthermore, this work establishes a platform to analyze 3D low-mode core asymmetries in ICF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Higher spin partition functions via the quasinormal mode method in de Sitter quantum gravity

In this note we compute the 1-loop partition function of spin- s s fields on Euclidean de Sitter space S^{2n+1} S 2 n + 1 using the quasinormal mode method. Instead of computing the quasinormal mode frequencies from scratch, we use the analytic continuation prescription L_{\text{AdS}}\to iL_{\text{dS}} L AdS → i L dS , appearing in the dS/CFT correspondence, and Wick rotate the normal mode frequencies of fields on thermal \text{AdS}_{2n+1} AdS 2 n + 1 into the quasinormal mode frequencies of fields on de Sitter space. We compare the quasinormal mode and heat kernel methods of calculating 1-loop determinants, finding exact agreement, and furthermore explicitly relate these methods via a sum over the conformal dimension. We discuss how the Wick rotation of normal modes on thermal \text{AdS}_{2n+1} AdS 2 n + 1 can be generalized to calculating 1-loop partition functions on the thermal spherical quotients S^{2n+1}/\mathbb{Z}_{p} S 2 n + 1 / ℤ p . We further show that the quasinormal mode frequencies encode the group theoretic structure of the spherical spacetimes in question, analogous to the analysis made for thermal AdS in [1-3] .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Nuclear charge densities in spherical and deformed nuclei: Toward precise calculations of charge radii

Precise measurements of atomic transitions affected by electron-nucleus hyperfine interactions offer sensitivity to explore basic properties of the atomic nucleus and study fundamental symmetries, including the search for new physics beyond the standard model of particle physics. Specifically, such measurements, augmented by atomic and nuclear calculations, will enable extraction of the higher-order radial moments of the charge-density distribution in spherical and deformed nuclei. The new data impose higher precision requirements on a theoretical description. The nuclear charge density is composed of the proton point distribution folded with the nucleonic charge distributions. The latter induce subtle relativistic corrections due to the coupling of nucleon magnetic moments with the nuclear spin-orbit density. Additional corrections come from the effect of center-of-mass projection. We assess the precision of nuclear charge density calculations by studying the behavior of relativistic and center-of-mass motion corrections to the second and fourth charge radial moments. Special attention has been paid to the magnetic spin-orbit density associated with the local variations of the spin-orbit current. The calculations for semimagic and open-shell nuclei are performed in the framework of self-consistent mean-field theory using quantified energy density functionals and density-dependent pairing forces. We used the general expression for the spin-orbit form factor that is valid for spherical and deformed nuclei. We studied the impact of various correction terms on the charge radii, fourth radial moments, diffraction radii, and surface thickness of spherical and deformed nuclei. The spin-orbit corrections to charge radial moments and surface thickness show strong shell fluctuations which can make an appreciable effect when aiming at high-precision predictions of isotopic shifts. The inclusion of relativistic and center-of-mass corrections impacts the quality of energy density functionals optimized to charge radii data. To establish reliable constraints on the existence of new forces from isotope shift measurements, precise calculations of nuclear charge densities of deformed nuclei are needed. The proper inclusion of the spin-orbit charge density and other correction terms is essential when aiming at extraction of subtle effects which become particularly visible in isotopic trends. It is also important when developing high-quality nuclear energy density functionals optimized using heterogeneous datasets involving absolute charge radii, differential charge radii, and charge form factor properties deduced from electron-scattering data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Analytical expression of a finite, long, conical canted-cosine-theta coil for particle collider interaction regions

Magnets in the accelerator interaction region (IR) present significant challenges because of high field requirements and limited available space. Conical-shaped magnets offer advantages in these environments by allowing closer placement to the interaction point while maintaining clearance from synchrotron radiation. Interestingly, numerical studies have shown that conical canted-cosine-theta (CCT) designs produce a constant field distribution along the axial direction in the IR quadrupoles for the Electron-Ion Collider (EIC) at Brookhaven National Laboratory. However, the field harmonics generated by conical CCT windings are not yet fully understood. This paper presents an analytical approach to describe the magnetic field produced by a conical surface current and proposes a method for designing conical CCT magnets for accelerator applications. First, we begin with a surface current sheet having a general cosine-theta distribution in spherical coordinates and solve the vector potential using the Green’s function. The magnetic fields generated by the conical current sheet are expressed using associated Legendre polynomials. These results are then related to circular field harmonics and integral field harmonics for designing a coil that produces a pure multipole field. Next, a single layer of the conical CCT winding path is produced based on the cosine-theta current distribution. Finally, the magnetic field quality of dipole and quadrupole conical CCT coils with multiple layers is verified using the Biot-Savart law.

Yang, Ye↗
DOE OSTI2019

Topological order in matrix Ising models

We study a family of models for an N_1 \times N_2 N 1 × N 2 matrix worth of Ising spins S_{aB} S a B . In the large N_i N i limit we show that the spins soften, so that the partition function is described by a bosonic matrix integral with a single ‘spherical’ constraint. In this way we generalize the results of to a wide class of Ising Hamiltonians with O(N_1,\mathbb{Z})\times O(N_2,\mathbb{Z}) O ( N 1 , ℤ ) × O ( N 2 , ℤ ) symmetry. The models can undergo topological large N N phase transitions in which the thermal expectation value of the distribution of singular values of the matrix S_{aB} S a B becomes disconnected. This topological transition competes with low temperature glassy and magnetically ordered phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

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

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

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Some aspects of the five-dimensional Lovelock black hole spacetime: Strong homotopy retract, perihelion precession and quasistationary levels

Highlights: • Some topological aspects are explored. • The perihelion precession is calculated. • The general Heun functions are applied to obtain analytical solutions. • The quasibound states and the radial wave eigenfunctions are obtained. In this work we explore some mathematical physics aspects of the spherically symmetric Lovelock black hole in high dimensions. Intended for this aim, we thoroughly consider the metric corresponding to the five-dimensional Lovelock black hole spacetime. We construct the strong retractions by the geodesic equations on the background under consideration. As a result, from the topological point of view, we construct the theory of strong homotopy retract, which will allow us, in principle, to better understand some of its suitable applications on astrophysics and cosmology, in particular, in the analysis of the spacetime singularities. We find the solutions of the equation of motion for both radial and angular coordinates, and then we describe the outer (“exterior”) and lower (“interior”) apparent horizons. Indeed, the outer apparent horizon is the last surface from which the light waves could still escape from the black hole. Thus, it is meaningful to analyze some physical phenomena related to quantum particles propagating outside the exterior apparent horizon, in particular, we discuss the quasistationary levels of scalar fields and their radial wave functions, which are given in terms of the general Heun functions. We also calculate the perihelion precession in this background.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Hyperheavy spherical and toroidal nuclei: The role of shell structure

The properties of toroidal hyperheavy even-even nuclei and the role of toroidal shell structure are extensively studied within covariant density functional theory. The general trends in the evolution of toroidal shapes in the Z ≈ 130–180 region of the nuclear chart are established for the first time. These nuclei are stable with respect to breathing deformations. The most compact fat toroidal nuclei are located in the Z ≈ 136, N ≈ 206 region of the nuclear chart, but thin toroidal nuclei become dominant with increasing proton number and toward proton and neutron drip lines. The roles of toroidal shell structure, its regularity, supershell structure, and shell gaps as well as the role of different groups of the pairs of the orbitals in its formation are investigated in detail. The lowest in energy solutions at axial symmetry are characterized either by large shell gaps or low density of the single-particle states in the vicinity of the Fermi level in at least one of the subsystems (proton or neutron). Related quantum shell effects are expected to act against the instabilities in breathing and sausage deformations for these subsystems. The investigation with a large set of covariant energy density functionals reveals that substantial proton Z = 154 and 186 and neutron N = 228, 308, and 406 spherical shell gaps exist in all functionals. The nuclei in the vicinity of the combination of these particle numbers form the islands of stability of spherical hyperheavy nuclei. Furthermore, the study suggests that the N = 210 toroidal shell gap plays a substantial role in the stabilization of fat toroidal nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Electric and magnetic black holes in a new nonlinear electrodynamics model

Highlights: • We introduce a new NED model which is comparable with the BI model in the weak-field limit. • We find an electric black hole solution in the context of Einstein’s gravity minimally coupled with the new NED model. • We study the thermal stability of the electric black hole solution. • Modified Smarr’s formula consistent with the first law of black hole thermodynamics is obtained. A new nonlinear electrodynamics (NED) model is introduced in the form of a nonpolynomial Lagrangian which admits static spherical and also plane wave solutions in a flat space. Upon coupling with gravity the electric field is finite and comparable with the Born–Infeld counterpart. The electric and magnetic black hole solutions in the Einstein’s gravity coupled with this NED model are presented. The solutions give both asymptotically and in the weak field limit Reissner–Nordström (RN) black hole and unlike the other known models our electric solution is expressed in terms of elementary functions in a closed form. We study the first law and derive the modified Smarr’s formula for the electric type extension of our model. Considerable rich structure, especially thermodynamic ones, ranging from first to the second order phase transitions are added to the RN black hole of linear electrodynamics with this NED model. Having the exact solution for the metric function at our disposal we investigate the stability of the electric black hole from both the thermodynamical and causal points of view.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Full-potential KKR within the removed-sphere method: A practical and accurate solution to the Poisson equation

Here an efficient and accurate generalization of the removed-sphere method (RSM) to solve the Poisson equation for total charge density in a solid with space-filling convex Voronoi polyhedra (VPs) and any symmetry is presented. The generalized RSM avoids the use of multipoles and VP shape functions for cellular integrals, which have associated ill-convergent large, double-internal L sums in spherical-harmonic expansions, so that fast convergence in single-L sums is reached. Our RSM adopts full Ewald formulation to work for all configurations or when symmetry breaking occurs, such as for atomic displacements or elastic constant calculations. The structure-dependent coefficients A L that define RSM can be calculated once for a fixed structure and speed up the whole self-consistent-field procedure. The accuracy and rapid convergence properties are confirmed using two analytic models, including the Coulomb potential and energy. We then implement the full-potential RSM using the Green's function Korringa-Kohn-Rostoker (KKR) method for real applications and compare the results with other first-principle methods and experimental data, showing that they are equally as accurate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Reciprocal swimming at intermediate Reynolds number

In Stokes flow, Purcell's scallop theorem forbids objects with time-reversible (reciprocal) swimming strokes from moving. In the presence of inertia, this restriction is eased and reciprocally deforming bodies can swim. A number of recent works have investigated dimer models that swim reciprocally at intermediate Reynolds numbers Re≈1–1000. These show interesting results (e.g. switches of the swim direction as a function of inertia) but the results vary and seem to be case specific. Here, in this work, we introduce a general model and investigate the behaviour of an asymmetric spherical dimer of oscillating length for small-amplitude motion at intermediate Re. In our analysis we make the important distinction between particle and fluid inertia, both of which need to be considered separately. We asymptotically expand the Navier–Stokes equations in the small-amplitude limit to obtain a system of linear partial differential equations. Using a combination of numerical (finite element) and analytical (reciprocal theorem, method of reflections) methods we solve the system to obtain the dimer's swim speed and show that there are two mechanisms that give rise to motion: boundary conditions (an effective slip velocity) and Reynolds stresses. Each mechanism is driven by two classes of sphere–sphere interactions, between one sphere's motion and (1) the oscillating background flow induced by the other's motion, and (2) a geometric asymmetry induced by the other's presence. We can thus unify and explain behaviours observed in other works. Our results show how sensitive, counterintuitive and rich motility is in the parameter space of finite inertia of particles and fluid.

42 ENGINEERING↗
DOE OSTI2024

Galaxy bispectrum in the spherical Fourier-Bessel basis

The bispectrum, the three-point correlation in Fourier space, is a crucial statistic for studying many effects targeted by the next-generation galaxy surveys, such as primordial non-Gaussianity (PNG) and general relativistic (GR) effects on large scales. In this work we develop a formalism for the bispectrum in the spherical Fourier-Bessel (SFB) basis—a natural basis for computing correlation functions on the curved sky, as it diagonalizes the Laplacian operator in spherical coordinates. Working in the SFB basis allows for line-of-sight effects such as redshift space distortions and GR to be accounted for exactly, i.e., without having to resort to perturbative expansions to go beyond the plane-parallel approximation. Only analytic results for the SFB bispectrum exist in the literature given the intensive computations needed. We numerically calculate the SFB bispectrum for the first time, enabled by a few techniques: We implement a template decomposition of the redshift-space kernel Z 2 into Legendre polynomials, and separately treat the PNG and velocity-divergence terms. We derive an identity to integrate a product of three spherical harmonics connected by a Dirac delta function as a simple sum and use it to investigate the limit of a homogeneous and isotropic Universe. Furthermore, we present a formalism for convolving the signal with separable window functions and use a toy spherically symmetric window to demonstrate the computation and give insights into the properties of the observed bispectrum signal. While our implementation remains computationally challenging, it is a step toward a feasible full extraction of information on large scales via a SFB bispectrum analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Systematic Way to Extend the Debye–Hückel Theory beyond Dilute Electrolyte Solutions

An extended Debye–Hückel theory with fourth order gradient term is developed for electrolyte solutions; namely, the electric potential φ(r) of the bulk electrolyte solution can be described by ∇ 2 φ(r) = κ 2 φ(r) + L Q 2 ∇ 4 φ(r), where the parameters κ and L Q are chosen to reproduce the first two roots of the dielectric response function of the bulk solution. Three boundary conditions for solving the electric potential problem are proposed based upon the continuity conditions of involving functions at the dielectric boundary, with which a boundary element method for the electric potential of a solute with a general geometrical shape and charge distribution is derived. Solutions for the electric potential of a spherical ion and a diatomic molecule are found and used to calculate their electrostatic solvation energies. In this work, the validity of the theory is successfully demonstrated when applied to binary as well as multicomponent primitive models of electrolyte solutions.

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

Theory of the effect of external stress on the activated dynamics and transport of dilute penetrants in supercooled liquids and glasses

We generalize the self-consistent cooperative hopping theory for a dilute spherical penetrant or tracer activated dynamics in dense metastable hard sphere fluids and glasses to address the effect of external stress, the consequences of which are systematically established as a function of matrix packing fraction and penetrant-to-matrix size ratio. All relaxation processes speed up under stress, but the difference between the penetrant and matrix hopping (alpha relaxation) times decreases significantly with stress corresponding to less time scale decoupling. A dynamic crossover occurs at a critical “slaving onset” stress beyond which the matrix activated hopping relaxation time controls the penetrant hopping time. This characteristic stress increases (decreases) exponentially with packing fraction (size ratio) and can be well below the absolute yield stress of the matrix. Below the slaving onset, the penetrant hopping time is predicted to vary exponentially with stress, differing from the power law dependence of the pure matrix alpha time due to system-specificity of the stress-induced changes in the penetrant local cage and elastic barriers. An exponential growth of the penetrant alpha relaxation time with size ratio under stress is predicted, and at a fixed matrix packing fraction, the exponential relation between penetrant hopping time and stress for different size ratios can be collapsed onto a master curve. Direct connections between the short- and long-time activated penetrant dynamics and between the penetrant (or matrix) alpha relaxation time and matrix thermodynamic dimensionless compressibility are also predicted. The presented results should be testable in future experiments and simulations.

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

A discrete dislocation dynamics study of precipitate bypass mechanisms in nickel-based superalloys

Order strengthening in nickel-based superalloys is associated with the extra stress required for dislocations to bypass the $\gamma'$ precipitates distributed in the $\gamma$ matrix. Depending on the operating conditions and microstructure, a rich variety of bypass mechanism has been identified, with various shearing and Orowan looping processes gradually giving way to climb bypass as the operating conditions change from the low/intermediate temperatures and high stress regime, to the high temperature and low stress regime. When anti phase boundary (APB) shearing and Orowan looping mechanisms operate, the classical picture is that, at for a given volume fraction, the bypass mechanism changes from shearing to looping with increased particle size and within a broad coexistence size window. Another possibility, which is supported by indirect experimental evidence, is that a third “hybrid” transition mechanism may operate. Here, in this paper we use discrete dislocation dynamics (DDD) simulations to study dislocation bypass mechanisms in Ni-based superalloys. We develop a new method to compute generalized stacking fault forces in DDD simulations, based on a concept borrowed from complex analysis and known as the winding number of a closed curve about a point. We use this method to study the mechanisms of bypass of a square lattice of spherical $\gamma'$ precipitates by $a/2\langle110\rangle$ {111} edge dislocations, as a function of the precipitates volume fraction and size. We show that not only the hybrid mechanism is possible, but also that it operates as the transition mechanism between the shearing and looping regimes over a wide range of precipitates volume fraction and radii. Based on our simulation results, we propose a simple model for the strength of this mechanism. We also consider the effects of a $\gamma$/$\gamma'$ lattice misfit on the bypass mechanisms, which we approximate by an additional precipitate stress computed according to Eshelby’s inclusion theory. We show that in the shearing and hybrid looping-shearing regimes, a lattice misfit generally results in an increased bypass stress. For sufficiently high lattice misfit, the critical bypass configuration in attractive dislocation-precipitates interactions changes dramatically, and the bypass stress is controlled by the pinning of the trailing dislocation on the exit side of the precipitates, similar to what has been reported in the high-temperature creep literature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nanoparticle cellular internalization is not required for RNA delivery to mature plant leaves

Rapidly growing interest in the nanoparticle-mediated delivery of DNA and RNA to plants requires a better understanding of how nanoparticles and their cargoes translocate in plant tissues and into plant cells. However, little is known about how the size and shape of nanoparticles influence transport in plants and the delivery efficiency of their cargoes, limiting the development of nanotechnology in plant systems. In this study we employed non-biolistically delivered DNA-modified gold nanoparticles (AuNPs) of various sizes (5-20 nm) and shapes (spheres and rods) to systematically investigate their transport following infiltration into Nicotiana benthamiana leaves. Generally, smaller AuNPs demonstrated more rapid, higher and longer-lasting levels of association with plant cell walls compared with larger AuNPs. We observed internalization of rod-shaped but not spherical AuNPs into plant cells, yet, surprisingly, 10 nm spherical AuNPs functionalized with small-interfering RNA (siRNA) were the most efficient at siRNA delivery and inducing gene silencing in mature plant leaves. These results indicate the importance of nanoparticle size in efficient biomolecule delivery and, counterintuitively, demonstrate that efficient cargo delivery is possible and potentially optimal in the absence of nanoparticle cellular internalization. Additionally, our results highlight nanoparticle features of importance for transport within plant tissues, providing a mechanistic overview of how nanoparticles can be designed to achieve efficacious biocargo delivery for future developments in plant nanobiotechnology.

59 BASIC BIOLOGICAL SCIENCES↗