SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Nanoscale Mapping and Defect-Assisted Manipulation of Surface Plasmon Resonances in 2D Bi 2 Te 3 /Sb 2 Te 3 In-Plane Heterostructures

Here, the Bi 2 Te 3 /Sb 2 Te 3 in-plane heterostructure is reported as a low-dimensional tunable chalcogenide well suited as plasmonic building block for the visible-UV spectral range. Electron-driven plasmon excitations of low-dimensional Bi 2 Te 3 /Sb 2 Te 3 are investigated by monochromated electron energy loss spectroscopy spectrum imaging. To resolve the nanoscale spatial distribution of various local plasmonic resonances, singular value decomposition is used to disentangle the spectral data and identify the individual spectral contributions of various corner, edge, and face modes. Furthermore, defect-plasmon interactions are investigated both for nanoscale intrinsic and thermally induced extrinsic polygonal defects (in situ sublimation). Signature of defect-induced red shift ranging from a several hundreds of millielectronvolts to a few electronvolts, broadening of various plasmon response, together with selective enhancement and significant variations in their intensity are detected. This study highlights the presence of a heterointerface and identifies defects as physical tuning pathways to modulate the plasmonic response over a broad spectral range. Finally, the experimental observations are compared qualitatively and validated with numerical simulations using the electron-driven discrete dipole approximation. Low-dimensional Bi 2 Te 3 /Sb 2 Te 3 as a less explored plasmonic system holds great promises as emerging platform for integrated plasmonics. Furthermore, introducing controlled structural defects can open the door for nanoengineering of plasmonic properties in such systems.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Obstacles to constructing de Sitter space in string theory

There have been many attempts to construct de Sitter space-times in string theory. While arguably there have been some successes, this has proven challenging, leading to the de Sitter swampland conjecture: quantum theories of gravity do not admit stable or metastable de Sitter space. Here we explain that, within controlled approximations, one lacks the tools to construct de Sitter space in string theory. Such approximations would require the existence of a set of (arbitrarily) small parameters, subject to severe constraints. But beyond this one also needs an understanding of big-bang and big-crunch singularities that is not currently accessible to standard approximations in string theory. The existence or non-existence of metastable de Sitter space in string theory remains a matter of conjecture.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Rescuing a black hole in the large- q coupled SYK model

In this paper, we develop a general effective theory for two copies of the Sachdev-Ye-Kitaev (SYK) model with a time-dependent bilinear coupling. For a quantum quench problem with an initial state of the thermofield double state, we show how the evolution of the system is described by a complex reparametrization field with a classical Hamiltonian. We study correlation functions in this system and compare the large-q theory with the bulk low energy effective theory. In particular, we study the special case of a “rescued black hole”, which describes how a time-evolved thermofield double state can evolve to the ground state of a coupled SYK model by a carefully tuned time-dependent coupling. In the low energy region, there is a holographic dual interpretation, which is a geometry that crosses over from an eternal black hole to a global AdS2 vacuum. This family of geometries allow us to access the bulk region that would be the black hole interior without the rescue process. By comparing the large-q and low energy theory, we find that even in the low energy region the deviation from the low energy theory cannot be neglected if the rescue process starts late. This provides evidence that the low energy effective theory of the bulk fails near the inner horizon of the black hole. We note the possibility of a connection to a two-dimensional analog of the higher-dimensional black hole singularity.

2D Gravity↗

DOE OSTI2021

Holography abhors visible trapped surfaces

We prove that consistency of the holographic dictionary implies a hallmark prediction of the weak cosmic censorship conjecture: that in classical gravity, trapped surfaces lie behind event horizons. In particular, the existence of a trapped surface implies the existence of an event horizon, and that furthermore this event horizon must be outside of the trapped surface. More precisely, we show that the formation of event horizons outside of a strong gravity region is a direct consequence of causal wedge inclusion, which is required by entanglement wedge reconstruction. We make few assumptions beyond the absence of evaporating singularities in strictly classical gravity. We comment on the implication that spacetimes with naked trapped surfaces do not admit a holographic dual, note a possible application to holographic complexity, and speculate on the dual CFT interpretation of a trapped surface.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Shape-constrained Input Estimation for Efficient Multi-shaker Vibration Testing

Multi-shaker vibration testing is gaining interest from structural dynamics test engineers as it can provide a much more accurate match to complicated field vibration responses than traditional single-axis shaker tests. However, the force capabilities of the small modal shakers typically used in multi-shaker vibration tests has limited the achievable response levels. To date, most multi-shaker vibration tests have been performed using a variety of standard, commercially-available control systems. While these control systems are adequate for a wide range of multiple-input/multiple-output tests, their control algorithms have not been tailored for the specific problem of multi-shaker vibration tests: efficiently coordinating the various shakers to work together to achieve a desired response. Here, a new input estimation algorithm is developed and demonstrated using simulations and actual test data. This algorithm, dubbed shape-constrained input estimation, is shown to effectively coordinate multiple shakers using a set of constraint vectors based on the deflection shapes of the test structure. This is accomplished by using the singular vector shapes of the system frequency response matrix, which allows the constraint vectors to automatically change as a function of frequency. Simulation and test results indicate a significant reduction in the input forces required to achieve a desired response. Finally, the results indicate that shape-constrained input estimation is an effective method to achieve higher response levels from limited shaker forces which will enable higher level multi-shaker vibration tests to be performed.

42 ENGINEERING↗

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 OSTI2025

Advantages of the Samarskii-type schemes on the Shishkin mesh

The schemes of the Samarskii type are simple modifications of the upwind scheme. We use them on the Shishkin mesh and discuss their advantages over the upwind scheme when applied to the linear one-dimensional singularly perturbed convection–diffusion problem. One of the advantages is that the Samarskii-type schemes have exact first-order accuracy uniform in the perturbation parameter, as opposed to the upwind scheme which is almost first-order uniformly accurate because its accuracy is diminished by logarithmic factors. Although this is not a new result, we re-emphasize it in the paper. We also demonstrate another advantage, that the Samarskii-type schemes are almost second-order uniformly accurate on the layer component of the solution. Motivated by this fact, we present a further improvement of the numerical method.

Convection–diffusion↗

DOE OSTI2022

Towards sharp error analysis of extended Lagrangian molecular dynamics

The extended Lagrangian molecular dynamics (XLMD) method provides a useful framework for reducing the computational cost of a class of molecular dynamics simulations with constrained latent variables. The XLMD method relaxes the constraints by introducing a fictitious mass ε for the latent variables and solving a set of singularly perturbed ordinary differential equations. While favorable numerical performance of XLMD has been demonstrated in several different contexts in the past decade, mathematical analysis of the method remains scarce. Here, we propose the first error analysis of the XLMD method in the context of a classical polarizable force field model. While the dynamics with respect to the atomic degrees of freedom are general and nonlinear, the key mathematical simplification of the polarizable force field model is that the constraints on the latent variables are given by a linear system of equations. We prove that when the initial value of the latent variables is compatible in a sense that we define, XLMD converges as the fictitious mass ε is made small with $\mathscr{O}$(ε) error for the atomic degrees of freedom and with $\mathscr{O}$($\sqrt{ε}$) error for the latent variables, when the dimension of the latent variable d' is 1. Furthermore, when the initial value of the latent variables is improved to be optimally compatible in a certain sense, we prove that the convergence rate can be improved to $\mathscr{O}$(ε) for the latent variables as well. Numerical results verify that both estimates are sharp not only for d'=1, but also for arbitrary d'. In the setting of general d', we do obtain convergence, but with the non-sharp rate of $\mathscr{O}$($\sqrt{ε}$) for both the atomic and latent variables.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

An efficient quantum circuit for block encoding a pairing Hamiltonian

We present an efficient quantum circuit for block encoding a pairing Hamiltonian often studied in nuclear physics. Our block encoding scheme does not require mapping the creation and annihilation operators to the Pauli operators and representing the Hamiltonian as a linear combination of unitaries. Instead, we show how to encode the Hamiltonian directly using controlled swap operations. We analyze the gate complexity of the block encoding circuit and show that it scales polynomially with respect to the number of qubits required to represent a quantum state associated with the pairing Hamiltonian. We also show how the block encoding circuit can be combined with the quantum singular value transformation to construct an efficient quantum circuit for approximating the density of states of a pairing Hamiltonian. The techniques presented can be extended to encode more general second-quantized Hamiltonians.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Causal horizons, geodesic completeness and stability in slow contraction cosmology

We show that cosmological models with a semi-infinite phase of slow contraction (ekpyrosis) possess a combination of properties that can address several fundamental problems in cosmology, otherwise faced in contracting de Sitter phases or standard big bang expansion. In particular, flat or open slow contraction admits a stable past attractor that asymptotes to Minkowski space and is past geodesically complete, as well as a stable, flat, homogeneous, and isotropic future attractor with negligible Weyl curvature (and, therefore, negligible gravitational entropy). In bouncing cosmologies, this contracting attractor is terminated by a smooth, non-singular bounce that transforms the attractor properties at the end of contraction into the initial conditions for the subsequent expanding phase. Cosmologies incorporating a slow contraction phase have no particle horizon and therefore avoid the causal horizon problem. The past Minkowski attractor also generates an initial spectrum of vacuum-like quantum fluctuations on all wavelengths. Moreover, because the averaged expansion rate along past-directed geodesics is non-positive, models incorporating a semi-infinite phase of slow contraction also evade the Borde–Guth–Vilenkin theorem. By contrast, contracting de Sitter space possesses a finite particle horizon and becomes unstable in the presence of scalar fields, matter, or radiation.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗

DOE OSTI2020

A holey cavity for single-transducer 3D ultrasound imaging with physical optimization

Within the compressive sensing (CS) framework, one effective way to increase the likelihood of successful signal reconstruction is to employ random processes in the construction of the sensing matrix. This work presents a 3D holey cavity, with diverse frequency modes, to spectrally code, that is, randomize, the ultrasound wave fields. The simulated results show that the use of such a cavity enables imaging simple or complex targets, such as spheres or the letter E, by only a single transceiver—something that is not possible without the use of a coding structure like the cavity. The effect of noise on imaging results and the size of the targets on the first-order Born approximation (BA) are also investigated. Moreover, this study attempts to optimize the cavity, based on a single numerical metric, such as the sum of singular values (SSV) or mutual coherence (MC). Yet, it will be shown that neither of these metrics can consistently compare the norm-one imaging performance between two cavities of different materials or hole sizes. This leaves finding a quantitative metric for these kinds of optimizations an open problem.

42 ENGINEERING↗

DOE OSTI2024

A flexible class of priors for orthonormal matrices with basis function-specific structure

Statistical modeling of high-dimensional matrix-valued data motivates the use of a low-rank representation that simultaneously summarizes key characteristics of the data and enables dimension reduction. Low-rank representations commonly factor the original data into the product of orthonormal basis functions and weights, where each basis function represents an independent feature of the data. However, the basis functions in these factorizations are typically computed using algorithmic methods that cannot quantify uncertainty or account for basis function correlation structure a priori. While there exist Bayesian methods that allow for a common correlation structure across basis functions, empirical examples motivate the need for basis function-specific dependence structure. We propose a prior distribution for orthonormal matrices that can explicitly model basis function-specific structure. The prior is used within a general probabilistic model for singular value decomposition to conduct posterior inference on the basis functions while accounting for measurement error and fixed effects. We discuss how the prior specification can be used for various scenarios and demonstrate favorable model properties through synthetic data examples. Finally, we apply our method to two-meter air temperature data from the Pacific Northwest, enhancing our understanding of the Earth system’s internal variability.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

High-Fidelity Large-Signal Order Reduction Approach for Composite Load Model

With the increasing penetration of electronic loads and distributed energy resources, conventional load models cannot capture their dynamics. Therefore, a new comprehensive composite load model is developed by Western Electricity Coordinating Council (WECC). However, this model is a complex high-order non-linear system with multi-time-scale property, which poses challenges on power system studies with heavy computational burden. In order to reduce the model complexity, the authors firstly develop a large-signal order reduction (LSOR) method using singular perturbation theory. In this method, the fast dynamics are integrated into the slow ones to preserve transient characteristics of the former. Then, accuracy assessment conditions are proposed and embedded into the LSOR to improve and guarantee the accuracy of reduced-order model. Finally, the reduced-order WECC composite load model is derived by using the proposed algorithm. Overall, simulation results show that the reduced-order large-signal model significantly alleviates the computational burden while maintaining similar dynamic responses as the original composite load model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

Operator-level quantum acceleration of non-logconcave sampling

Sampling from probability distributions of the form 𝝈 ∝ e −𝜷V , where V is a continuous potential, is a fundamental task across physics, chemistry, biology, computer science, and statistics. However, when V is nonconvex, the resulting distribution becomes non-logconcave, and classical methods such as Langevin dynamics often exhibit poor performance. We introduce a quantum algorithm that provably accelerates a broad class of continuous-time sampling dynamics. For Langevin dynamics, our method encodes the target Gibbs measure into the amplitudes of aquantum state, identified as the kernel of a block matrix derived from a factorization of the Witten Laplacian operator. This connection enables Gibbs sampling via singular value thresholding and yields up to a quartic quantum speedup over best-knownclassical Langevin-based methods in the non-logconcave setting. Building on this framework, we further develop the first quantum algorithm that accelerates replica exchange Langevin diffusion, a widely used method for sampling from complex, rugged energy landscapes.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Solution Irregularity Remediation for Spatial Discretization Error Estimation for S N Transport Solutions

The discrete ordinates linear Boltzmann transport equation is typically solved in its spatially discretized form, incurring spatial discretization error. Quantification of this error for purposes such as adaptive mesh refinement or error analysis requires an a posteriori estimator, which utilizes the numerical solution to the spatially discretized equation to compute an estimate. Because the quality of the numerical solution informs the error estimate, irregularities, present in the true solution for any realistic problem configuration, tend to cause the largest deviation in the error estimate vis-a-vis the true error. In this paper, an analytical partial singular characteristic tracking (pSCT) procedure for reducing the estimator’s error is implemented within our novel residual source estimator for a zeroth-order discontinuous Galerkin scheme, at the additional cost of a single inner iteration. Here, a metric-based evaluation of the pSCT scheme versus the standard residual source estimator is performed over the parameter range of a Method of Manufactured Solutions test suite. The pSCT scheme generates near-ideal accuracy in the estimate in problems where the dominant source of the estimator’s error is the solution irregularity, namely, problems where the true solution is discontinuous and problems where the true solution’s first derivative is discontinuous and the scattering ratio is low. In problems where the scattering ratio is high and the true solution is discontinuous in the first derivative, the error in the scattering source, which is not converged by the pSCT scheme, is greater than the error incurred due to the irregularity. Ultimately, a pSCT scheme is judged to be useful for error estimation in problems where the computational cost of the scheme is justified. In the presence of many irregularities, such a scheme may be intractable for general use, but in benchmarks, as an analytical tool, or in problems that have nondissipative discontinuities, the scheme may prove invaluable.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Competing orders at higher-order Van Hove points

Van Hove points are special points in the energy dispersion, where the density of states exhibits analytic singularities. When a Van Hove point is close to the Fermi level, tendencies towards density wave orders, Pomeranchuk orders, and superconductivity can all be enhanced, often in more than one channel, leading to a competition between different orders and unconventional ground states. In this study we consider the effects from higher-order Van Hove points, around which the dispersion is flatter than near a conventional Van Hove point, and the density of states has a power-law divergence. We argue that such points are present in intercalated graphene and other materials. We use an effective low-energy model for electrons near higher-order Van Hove points and analyze the competition between different ordering tendencies using an unbiased renormalization-group approach. For purely repulsive interactions, we find that two key competitors are ferromagnetism and chiral superconductivity. For a small attractive interaction, we find an unconventional spin Pomeranchuk order, in wich the spin oder parameter winds around the Fermi surface. The supermetal state, predicted for a single higher-order Van Hove point, is an unstable fixed point in our case.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Correlated metals and unconventional superconductivity in rhombohedral trilayer graphene: A renormalization group analysis

Here, motivated by recent experimental observations of correlated metallic phases and superconductivity in rhombohedral trilayer graphene (RTG), we perform an unbiased study of electronic ordering instabilities in hole-doped RTG. Specifically, we focus on electronic states energetically proximate to Van Hove singularities (VHSs), where a large density of states promotes different interaction-induced symmetry-breaking electronic orders. To resolve the Fermi surface near VHSs, we construct a fermionic hot-spot model and demonstrate that a perpendicular electric field can tune different nesting structures of the Fermi surface. Subsequently, we apply a renormalization group analysis to describe the low-energy phase diagrams of our model under both short-range repulsive interactions as well as realistic (long-range) Coulomb interactions. Our analysis shows instabilities towards either intervalley coherent metallic phases or superconducting phases. The dominant pairing channel depends crucially on the nature of Fermi surface nesting—repulsive Coulomb interaction favors spin-singlet d-wave pairing for relatively small displacement field and spin-singlet i-wave pairing for larger displacement field. We argue that the phase diagram of RTG can be well-understood by modeling the realistic Coulomb interaction as the sum of repulsive density-density interaction and ferromagnetic spin-triplet intervalley coherence (IVC) Hund's coupling, while phonon-mediated electronic interactions have a negligible effect on this system, in sharp contrast to twisted graphene multilayers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Large critical fields in superconducting Ti 4 ⁢Ir 2⁢ O from spin-orbit coupling

The recently synthesized 𝜂-carbide-type superconductors exhibit large critical fields. A notable example is Ti 4 ⁢Ir2⁢ O, for which the upper critical field strongly violates the Pauli paramagnetic limit, a behavior that is unusual for cubic materials that preserve inversion symmetry. Here, by combining density functional theory (DFT) and analytic modeling, we provide an explanation for this enhanced Pauli limiting field. We show that the nonsymmorphic Fd⁢-3 ⁢m symmetry implies that the electronic states near the X points exhibit strong spin-orbit coupling (SOC), which leads to a vanishing effective 𝑔 factor and enables the enhanced Pauli limiting field. Furthermore, our DFT results reveal a Van Hove singularity peak near the X points, accounting for ∼65% of the total density of states (DOS), occurring near the chemical potential. We propose that the strong SOC and enhanced DOS in the vicinity of the X points provide the origin of the observed enhancement of the critical field. This leads to a prediction that the magnetic field will lead to a strongly momentum-dependent gap suppression. As a result, the gap due to electronic states away from (near to) the X points will be rapidly (slowly) suppressed by fields.

Critical field↗