Search NASA⌕ Search

SEARCH · Search NASA

Results for “Singularity”

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 OSTI2023

Local Decomposition of Hexahedral Singular Nodes into Singular Curves

Hexahedral (hex) meshing is a long studied topic in geometry processing with many challenging associated problems. Hex meshes vary from structured to unstructured depending on application or domain of interest. Fully structured meshes require that all interior mesh edges be adjacent to four hexes each. Edges failing this criteria are singular and indicate an unstructured hex mesh. Singular edges join together into singular curves that either form closed cycles, end on the mesh boundary, or end at a singular node, a complex junction of more than two singular curves. Hex meshes with more complex singular nodes tend to have more distorted elements and smaller scaled Jacobian values. In this work, we study the topology of singular nodes. We show that all eight of the most common singular nodes are decomposable into just singular curves. We further show that all singular nodes, regardless of edge valence, are locally decomposable. Finally we demonstrate these decompositions on hex meshes, thereby decreasing their distortion and converting all singular nodes into singular curves. In conclusion, with this decomposition, the enigmatic complexity of 3D singular nodes becomes effectively 2D.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Polar Scaling Technique for the Regularization of Strongly Singular and Strongly Near-Singular Helmholtz Surface Integrals Evaluated Over 2-D Domains

The numerical integration of expressions containing strong singularities or strong near-singularities has long been a challenging problem in the electromagnetics community. Much attention has been paid to this problem, as strong $1/R^{{2}}$ singularities routinely appear when implementing electromagnetic simulation techniques like the method of moments (MoM). To date, several techniques, from singularity extraction to singularity cancellation (SC), have been employed to deal with problems that require the evaluation of 2-D strongly singular integrals. However, no single technique has been proposed that can deal with both strong singularities and strong near-singularities in a fully numerical manner for arbitrary 2-D domains. Moreover, it has been claimed that the Helmholtz-type strongly singular integral found in the MoM is convergent in a principal value sense, but this convergence value has yet to be proven mathematically. In this work, we will conduct the convergence proof and introduce a “polar scaling” change of variables method that may be used to evaluate Helmholtz integrals with both strong and weak singularities/near-singularities. The technique is fully numerical and can in principle be applied to any planar or curved polygon and any nonsingular basis function. We will also provide numerical results showing useful convergence behavior for integrals involving both exact and near-singularities.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

The structure of quantum singularities on a Cauchy horizon

Spacetime singularities pose a long-standing puzzle in quantum gravity. Unlike Schwarzschild, a generic family of black holes gives rise to a Cauchy horizon on which, even in the Hartle-Hawking state, quantum observables such as ⟨T μν ⟩ — the expectation value of the stress-energy tensor — can diverge, causing a breakdown of semiclassical gravity. Because they are diagnosed within quantum field theory (QFT) on a smooth background, these singularities may provide a better-controlled version of the spacetime singularity problem, and merit further study. Here, I highlight a mildness puzzle of Cauchy horizon singularities: the ⟨T μν ⟩ singularity is significantly milder than expected from symmetry and dimensional analysis. I address the puzzle in a simple spacetime ${\mathcal{W}}_P$, which arises universally near all black hole Cauchy horizons: the past of a codimension-two spacelike plane in flat spacetime. Specifically, I propose an extremely broad QFT construction in which, roughly speaking, Cauchy horizon singularities originate from operator insertions in the causal complement of the spacetime. The construction reproduces well-known outer horizon singularities (e.g., in the Boulware state), and remarkably, when applied to ${\mathcal{W}}_P$, gives rise to a universal mild singularity structure for robust singularities, ones whose leading singular behavior is state-independent. I make non-trivial predictions for all black hole Cauchy horizon singularities using this, and discuss extending the results beyond robust singularities and the strict near Cauchy horizon limit.

black holes↗
DOE OSTI2023

Single monkey-saddle singularity of a Fermi surface and its instabilities

Fermi surfaces can undergo sharp transitions under smooth changes of parameters. Such transitions can have a topological character, as is the case when a higher-order singularity, one that requires cubic or higher-order terms to describe the electronic dispersion near the singularity, develops at the transition. When time-reversal and inversion symmetries are present, odd singularities can only appear in pairs within the Brillouin zone. In this case, the combination of the enhanced density of states that accompanies these singularities and the nesting between the pairs of singularities leads to interaction-driven instabilities. Here, we present examples of single n = 3 (monkey-saddle) singularities when time-reversal and inversion symmetries are broken. We then turn to the question of what instabilities are possible when the singularities are isolated. For spinful electrons, we find that the inclusion of repulsive interactions destroys any isolated monkey-saddle singularity present in the noninteracting spectrum by developing Stoner or Lifshitz instabilities. In contrast, for spinless electrons and at the mean-field level, we show that an isolated monkey-saddle singularity can be stabilized in the presence of short-range repulsive interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Linking the singularities of cosmological correlators

Much of the structure of cosmological correlators is controlled by their singularities, which in turn are fixed in terms of flat-space scattering amplitudes. An important challenge is to interpolate between the singular limits to determine the full correlators at arbitrary kinematics. This is particularly relevant because the singularities of correlators are not directly observable, but can only be accessed by analytic continuation. In this paper, we study rational correlators — including those of gauge fields, gravitons, and the inflaton — whose only singularities at tree level are poles and whose behavior away from these poles is strongly constrained by unitarity and locality. We describe how unitarity translates into a set of cutting rules that consistent correlators must satisfy, and explain how this can be used to bootstrap correlators given information about their singularities. We also derive recursion relations that allow the iterative construction of more complicated correlators from simpler building blocks. In flat space, all energy singularities are simple poles, so that the combination of unitarity constraints and recursion relations provides an efficient way to bootstrap the full correlators. In many cases, these flat-space correlators can then be transformed into their more complex de Sitter counterparts. As an example of this procedure, we derive the correlator associated to graviton Compton scattering in de Sitter space, though the methods are much more widely applicable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Direct geometric probe of singularities in band structure

A quantum system’s energy landscape may have points where multiple energy surfaces are degenerate and that exhibit singular geometry of the wave function manifold, with major consequences for the system’s properties. Ultracold atoms in optical lattices have been used to indirectly characterize such points in the band structure. Here we measured the non-Abelian transformation produced by transport directly through the singularities. We accelerated atoms along a quasi-momentum trajectory that enters, turns, and then exits the singularities at linear and quadratic band-touching points of a honeycomb lattice. Measurements after transport identified the topological winding numbers of these singularities to be 1 and 2, respectively. Our work introduces a distinct method for probing singularities that enables the study of non-Dirac singularities in ultracold-atom quantum simulators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Spatially quasi-periodic bifurcations from periodic traveling water waves and a method for detecting bifurcations using signed singular values

We present a method of detecting bifurcations by locating zeros of a signed version of the smallest singular value of the Jacobian. This enables the use of quadratically convergent root-bracketing techniques or Chebyshev interpolation to locate bifurcation points. Only positive singular values have to be computed, though the method relies on the existence of an analytic or smooth singular value decomposition (SVD). The sign of the determinant of the Jacobian, computed as part of the bidiagonal reduction in the SVD algorithm, eliminates slope discontinuities at the zeros of the smallest singular value. We use the method to search for spatially quasi-periodic traveling water waves that bifurcate from large-amplitude periodic waves. The water wave equations are formulated in a conformal mapping framework to facilitate the computation of the quasi-periodic Dirichlet-Neumann operator. We find examples of pure gravity waves with zero surface tension and overhanging gravity-capillary waves. In both cases, the waves have two spatial quasi-periods whose ratio is irrational. We follow the secondary branches via numerical continuation beyond the realm of linearization about solutions on the primary branch to obtain traveling water waves that extend over the real line with no two crests or troughs of exactly the same shape. The pure gravity wave problem is of relevance to ocean waves, where capillary effects can be neglected. Such waves can only exist through secondary bifurcation as they do not persist to zero amplitude. The gravity-capillary wave problem demonstrates the effectiveness of using the signed smallest singular value as a test function for multi-parameter bifurcation problems. This test function becomes mesh independent once the mesh is fine enough.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Altermagnetism from coincident Van Hove singularities: application to κ -Cl

Realizing two-dimensional (2D) altermagnets is important for spintronics applications. Here we propose a microscopic template for stabilizing 2D altermagnetism through Van Hove singularities that are coincident in both energy and momentum. These coincident Van Hove singularities are a generic consequence of non-symmorphic symmetries in nine 2D space groups. Due to nontrivial symmetry properties of the Hamiltonian, these coincident Van Hove singularities allow new hopping interactions between the Van Hove singularities that do not appear in analogous Van Hove singularity based patch models for cuprates and graphene. We show these new interactions can give rise to various weak coupling, and BCS-based instabilities, including altermagnetism, nematicity, inter-band d-wave superconductivity, and orbital altermagnetic order. We apply our results to quasi-2D organic κ -Cl in which altermagnetism is known to appear.

electronic properties and materials↗
DOE OSTI2025

A Class of Sparse Johnson–Lindenstrauss Transforms and Analysis of their Extreme Singular Values

The Johnson–Lindenstrauss (JL) lemma is a powerful tool for dimensionality reduction in modern algorithm design. The lemma states that any set of high-dimensional points in a Euclidean space can be projected into lower dimensions while approximately preserving pairwise Euclidean distances. Random matrices satisfying this lemma are called JL transforms (JLTs). Inspired by existing $s$-hashing JLTs with exactly $s$ nonzero elements on each column, the present work introduces an ensemble of sparse matrices encompassing so-called $s$-hashing-like matrices whose expected number of nonzero elements on each column is $s$. The independence of the sub-Gaussian entries of these matrices and the knowledge of their exact distribution play an important role in their analyses. Using properties of independent sub-Gaussian random variables, these matrices are demonstrated to be JLTs, and their smallest nontrivial singular values and largest singular values are estimated nonasymptotically using a technique from geometric functional analysis. As the dimensions of the matrix grow to infinity, these singular values are proved to converge almost surely to fixed quantities (by using the universal Bai–Yin law) and in distribution to the Gaussian orthogonal ensemble Tracy–Widom law after proper rescalings. Understanding the behaviors of extreme singular values is important in general because they are often used to define a measure of stability of matrix algorithms. For example, JLTs were recently used in derivative-free optimization algorithmic frameworks to select random subspaces in which are constructed random models or poll directions to achieve scalability, and hence estimating their smallest singular value in particular helps determine the dimension of these subspaces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Quantum fate of the Choptuik naked singularity

Classical critical collapse provides a dynamical route from smooth initial data to a naked singularity, representing a sharper violation of predictability than ordinary black hole singularities. We argue that this distinction is erased by quantum backreaction. Building on the semiclassical interior analysis, where quantum self-energy of the collapsing matter generates a universal growing mode and a finite mass gap, we study the exterior naked singularity region that determines global visibility in the Einstein-scalar system. We analyze controlled exterior models in both $2+1$ and $3+1$ dimensions. In the former, smooth matching and physical boundary conditions analytically select a vacuum polarization state, whose backreaction cloaks the classically naked region by a quantum trapped branch. In the latter, numerical horizon tracing shows that near a quantum-shifted threshold the exterior develops finite-mass marginally trapped surfaces rather than a zero-mass naked endpoint. These results suggest a global quantum picture in which the Choptuik naked singularity shares the fate of an ordinary black hole singularity: quantum effects push the putative Cauchy horizon behind a quantum-generated horizon, thereby reducing the loss of predictability to the standard black hole evaporation problem.

FOS: Physical sciences↗
DOE OSTI2023

Structure of pressure-gradient-driven current singularity in ideal magnetohydrodynamic equilibrium

Abstract Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic (MHD) equilibria with a continuum of nested flux surfaces and a continuous rotational transition. These currents have two components: a surface current (Dirac δ -function in flux surface labeling) that prevents the formation of magnetic islands, and an algebraically divergent Pfirsch–Schlüter current density when a pressure gradient is present across the rational surface. On flux surfaces adjacent to the rational surface, the traditional treatment gives the Pfirsch–Schlüter current density scaling as J ∼ 1 / Δ ι , where Δ ι is the difference of the rotational transform relative to the rational surface. If the distance s between flux surfaces is proportional to Δ ι , the scaling relation J ∼ 1 / Δ ι ∼ 1 / s will lead to a paradox that the Pfirsch–Schlüter current is not integrable. In this work, we investigate this issue by considering the pressure-gradient-driven singular current in the Hahm–Kulsrud–Taylor problem, which is a prototype for singular currents arising from resonant magnetic perturbations. We show that not only the Pfirsch–Schlüter current density but also the diamagnetic current density are divergent as ∼ 1 / Δ ι . However, due to the formation of a Dirac δ -function current sheet at the rational surface, the neighboring flux surfaces are strongly packed with s ∼ ( Δ ι ) 2 . Consequently, the singular current density J ∼ 1 / s , making the total current finite, thus resolving the paradox. Furthermore, the strong packing of flux surfaces causes a steepening of the pressure gradient near the rational surface, with ∇ p ∼ d p / d s ∼ 1 / s . In general non-axisymmetric MHD equilibrium, contrary to Grad’s conjecture that the pressure profile is flat around densely distributed rational surfaces, our result suggests a pressure profile that densely steepens around them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Explicit Monotone Stable Super-Time-stepping Methods for Finite Time Singularities

We explore a novel way to numerically resolve the scaling behavior of finite-time singularities in solutions of nonlinear parabolic PDEs. The Runge–Kutta–Legendre (RKL) and Runge–Kutta–Gegenbauer (RKG) super-time-stepping methods were originally developed for nonlinear complex physics problems with diffusion. These are multistage single step second-order, forward-in-time methods with no implicit solves. The advantage is that the time-step size for stability scales with stage number 𝑠 as $\mathcal{O}$⁡(𝑠 2 ). Many interesting nonlinear PDEs have finite-time singularities, and the presence of diffusion often limits one to using implicit or semi-implicit time-step methods for stability constraints. Finite-time singularities are particularly challenging due to the large range of scales that one desires to resolve, often with adaptive spatial grids and adaptive time steps. Here, in this study, we show two examples of nonlinear PDEs for which the self-similar singularity structure has time and space scales that are resolvable using the RKL and RKG methods, without forcing even smaller time steps. Compared to commonly used implicit numerical methods, we achieve a significantly smaller run time while maintaining comparable accuracy. We also prove numerical monotonicity for both the RKL and RKG methods under their linear stability conditions for the constant coefficient heat equation, in the case of infinite domain and periodic boundary condition, leading to a theoretical guarantee of the superiority of the RKL and RKG methods over traditional super-time-stepping methods, such as the Runge-Kutta-Chebyshev and the orthogonal Runge-Kutta-Chebyshev methods. Code can be found at https://github.com/ZT220501/SRK-Singularity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Structural response reconstruction using a system-equivalent singular vector basis

Here, this paper develops a novel method for reconstructing the full-field response of structural dynamic systems using sparse measurements. The singular value decomposition is applied to a frequency response matrix relating the structural response to physical loads, base motion, or modal loads. The left singular vectors form a non-physical reduced basis that can be used for response reconstruction with far fewer sensors than existing methods. The contributions of the singular vectors to measured response are termed singular-vector loads (SVLs) and are used in a regularized Bayesian framework to generate full-field response estimates and confidence intervals. The reconstruction framework is applicable to the estimation of single data records and power spectral densities from multiple records. Reconstruction is successfully performed in configurations where the number of SVLs to identify is less than, equal to, and greater than the number of sensors used for reconstruction. In a simulation featuring a seismically excited shear structure, SVL reconstruction significantly outperforms modal FRF-based reconstruction and successfully estimates full-field responses with as few as two uniaxial accelerometers. SVL reconstruction is further verified in a simulation featuring an acoustically excited cylinder. Finally, response reconstruction and uncertainty quantification are performed on an experimental structure with three shaker inputs and 27 triaxial accelerometer outputs.

42 ENGINEERING↗
DOE OSTI2023

A fast, accurate and easy to implement Kapur–Rokhlin quadrature scheme for singular integrals in axisymmetric geometries

Many applications in magnetic confinement fusion require the efficient calculation of surface integrals with singular integrands. The singularity subtraction approaches typically used to handle such singularities are complicated to implement and low-order accurate. In contrast, we demonstrate that the Kapur–Rokhlin quadrature scheme is well-suited for the logarithmically singular integrals encountered for a toroidally axisymmetric confinement system, is easy to implement and is high-order accurate. As an illustration, we show how to apply this quadrature scheme for the efficient and accurate calculation of the normal component of the magnetic field due to the plasma current on the plasma boundary, via the virtual-casing principle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Spin density wave and van Hove singularity in the kagome metal CeTi 3 Bi 4

Kagome metals with van Hove singularities near the Fermi level can host intriguing quantum phenomena such as chiral loop currents, electronic nematicity, and unconventional superconductivity. However, to our best knowledge, unconventional magnetic states driven by van Hove singularities–like spin-density waves–have not been observed experimentally in kagome metals. Here, we report the magnetic and electronic structure of the layered kagome metal CeTi 3 Bi 4 , where Ti kagome electronic structure interacts with a magnetic sublattice of Ce 3+ J eff = 1/2 moments. Neutron diffraction reveals an incommensurate spin-density wave ground state of the Ce 3+ moments, coexisting with commensurate antiferromagnetic order across most of the temperature-field phase diagram. The commensurate component is preferentially suppressed by thermal fluctuations and magnetic field, yielding a rich phase diagram involving an intermediate single-Q spin-density wave phase. First-principles calculations and angle-resolved photoemission spectroscopy identify van Hove singularities near the Fermi level, with the observed magnetic propagation vectors connecting their high density of states, strongly suggesting a van Hove singularity-assisted spin-density wave. These findings establish kagome metals LnTi 3 Bi 4 as a model platform where the characteristic electronic structure of the kagome lattice plays a pivotal role in magnetic order.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Superfluid weight cross-over and critical temperature enhancement in singular flat bands

Nonanalytic Bloch eigenstates at isolated band degeneracy points exhibit singular behavior in the quantum metric. Here, a description of superfluid weight for zero-energy flat bands in proximity to other high-energy bands is presented, where they together form a singular band gap system. When the singular band gap closes, the geometric and conventional contributions to the superfluid weight as a function of the superconducting gap exhibit different cross-over behaviors. The scaling behavior of superfluid weight with the band gap is studied in detail, and the effect on the Berezinskii–Kosterlitz–Thouless transition temperature is explored. It is found that tuning the singular band gap provides a unique mechanism for enhancing the supercurrent and critical temperature of two-dimensional superconductors.

Science & Technology - Other Topics↗
DOE OSTI2022

Code verification for practically singular equations

We report the method-of-moments implementation of the electric-field integral equation (EFIE) yields many code-verification challenges due to the various sources of numerical error and their possible interactions. Matters are further complicated by singular integrals, which arise from the presence of a Green's function. To address these singular integrals, an approach is presented in wherein both the solution and Green's function are manufactured. Because the arising equations are poorly conditioned, they are reformulated as a set of constraints for an optimization problem that selects the solution closest to the manufactured solution. In this paper, we demonstrate how, for such practically singular systems of equations, computing the truncation error by inserting the exact solution into the discretized equations cannot detect certain orders of coding errors. On the other hand, the discretization error from the optimal solution is a more sensitive metric that can detect orders less than those of the expected convergence rate.

97 MATHEMATICS AND COMPUTING↗