Engineering Papers⌕ Search

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.

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

At least 631 records · Page 35

Finite density QCD equation of state: Critical point and lattice-based 𝑇′ expansion

Here, we present a novel construction of the QCD equation of state (EoS) at finite baryon density. Our work combines a recently proposed resummation scheme for lattice QCD results with the universal critical behavior at the QCD critical point. This allows us to obtain a family of equations of state in the range 0 ≤ 𝜇 𝐵 ≤ 700 MeV and 25 MeV ≤ 𝑇 ≤ 800 MeV, which match lattice QCD results near 𝜇 𝐵 =0 while featuring a critical point in the 3D Ising model universality class. The position of the critical point can be chosen within the range accessible to beam-energy scan heavy-ion collision experiments. The strength of the singularity and the shape of the critical region are parametrized using a standard parameter set. We impose stability and causality constraints and discuss the available ranges of critical point parameter choices, finding that they extend beyond earlier parametric QCD EoS proposals. We present thermodynamic observables, including baryon density, pressure, entropy density, energy density, baryon susceptibility and speed of sound, that cover a wide range in the QCD phase diagram relevant for experimental exploration.

Astronomy & Astrophysics↗

Renormalization of asymmetric staple-shaped Wilson-line operators in lattice and continuum perturbation theory

In this work, we study the renormalization of nonlocal quark bilinear operators containing an asymmetric staple-shaped Wilson line at the one-loop level in both lattice and continuum perturbation theory. These operators enter the first-principle calculation of transverse momentum-dependent parton distribution functions (TMDPDFs) in lattice QCD using the formulation of large momentum effective theory. We provide appropriate RI ′ -type conditions that address the power and logarithmic divergences, as well as the mixing among staple operators of different Dirac structures, using a number of different possible projectors. A variant of RI ′ , including calculations of rectangular Wilson loops, which cancel the pinch-pole singularities of the staple operators at infinite length and reduce residual power divergences, is also employed. We calculate at one-loop order the conversion matrix, which relates the quasi-TMDPDFs in the RI ′ -type schemes to the reference scheme MS ¯ for arbitrary values of the renormalization momentum scale and of the dimensions of the staple. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Waiting for inflation: A new initial state for the Universe

We propose a cosmological lingering phase for the initial state prior to inflation which would help address the singularity problem of inflation. The Universe begins with a constant (Hagedorn) temperature and then transitions into an inflationary universe while preserving the null energy condition (NEC). We first consider the phase space of positive spatial curvature models within general relativity and with matter sources that respect the NEC. Depending on the duration of the postlingering inflation, these models can produce a small amount of observable spatial curvature in the cosmic microwave background. We also discuss how lingering can arise with or without spatial curvature in theories of quantum gravity when considering the thermodynamic scaling of particles and its impact on the early Universe. The string theory dilaton is essential to the dynamics. There are many open questions that remain. Published by the American Physical Society 2024

Melcher, Brandon (ORCID:0000000280849682)↗

Breakdown of collinear factorization in the exclusive photoproduction of a π 0 γ pair with large invariant mass

We study the exclusive photoproduction of a π 0 γ pair with large invariant mass M γ π 2 , which is sensitive to the exchange of either two quarks or two gluons in the t channel. In this paper, we show that the process involving two-gluon exchanges does factorize in the Bjorken limit at the leading twist. This can be explicitly demonstrated by the fact that there exist diagrams, which contribute at the leading twist, for which are , due to the pinching of the contour integration of the plus and minus component of the Glauber gluon momentum. For the same reason, π 0 -nucleon scattering to two photons also suffers from the same issue. On the other hand, we stress that there are no issues with respect to collinear factorization for the quark channels. By considering an analysis of all potential reduced diagrams of leading pinch-singular surfaces, we argue that the quark channel is safe from Glauber pinches, and therefore, a collinear factorization in that case follows through without any problems. This means that processes where gluon exchanges are forbidden, such as the exclusive photoproduction of π ± γ and ρ 0 , ± γ , are unaffected by the factorization breaking effects we point out in this paper. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantum real-time evolution using tensor renormalization group methods

We introduce an approach for approximate real-time evolution of quantum systems using tensor renormalization group (TRG) methods originally developed for imaginary time. We use higher-order TRG to generate a coarse-grained time evolution operator for a 1+1⁢D transverse Ising model with a longitudinal field. We show that the standard tensor norm used for the singular value decomposition-based truncation is degenerate and propose an alternate method to discriminate. We show that it is effective and efficient in evolving Gaussian wave packets for one and two particles in the disordered phase, while ordered phase behavior is more challenging to capture. We compare our algorithm with local simulators for universal quantum computers and discuss possible benchmarking in the near future.

lattice gauge theory↗

Generalized symmetries of nonsupersymmetric orbifolds

We determine generalized symmetries for 4D theories engineered via type II strings on nonsupersymmetric orbifold backgrounds ℝ 3,1 × ℝ 6 /Γ. Probe branes detect generalized symmetries via the adjacency matrix for fermionic degrees of freedom in an associated quiver gauge theory. In situations where the tachyons are sequestered away from the boundary 𝑆 5 /Γ, this exactly matches the result extracted from singular homology. In situations with an unsequestered tachyon which stretches out to the boundary, the presence of tachyonic pulses partitions up the space into several distinct sectors, and the net contribution again matches with the answer expected via quiver methods. For IIA backgrounds, the presence of a localized closed string tachyon leads to transitions in the spectrum of states, generalized symmetries, higher-group symmetries, as well as the level matrix of the associated symmetry topological field theory (SymTFT). For IIB backgrounds with a stack of spacetime filling probe D3-branes, the onset of a radiatively generated potential leads to similar considerations involving scale dependent transitions in the symmetries of the theory, including structures such as duality defects/interfaces.

compactification↗

Virasoro completeness relation and the inverse Shapovalov form

In this work, we introduce an explicit expression for the inverse of the symmetric bilinear form of Virasoro Verma modules, the so-called Shapovalov form, in terms of singular vector operators and their conformal dimensions. Our proposed expression also determines the resolution of the identity for Verma modules of the Virasoro algebra, and can be thus employed in the computation of Virasoro conformal blocks via the sewing procedure. Published by the American Physical Society 2025

Fortin, Jean-François (ORCID:0000000278135015)↗

Ergodic and nonergodic many-body dynamics in strongly nonlinear lattices

The study of nonlinear oscillator chains in classical many-body dynamics has a storied history going back to the seminal work of Fermi et al. [Los Alamos Scientific Laboratory Report No. LA-1940, 1955 (unpublished)]. Here, we introduce a family of such systems which consist of chains of N harmonically coupled particles with the nonlinearity introduced by confining the motion of each individual particle to a box or stadium with hard walls. The stadia are arranged on a one-dimensional lattice but they individually do not have to be one dimensional, thus permitting the introduction of chaos already at the lattice scale. For the most part we study the case where the motion is entirely one dimensional. We find that the system exhibits a mixed phase space for any finite value of N . Computations of Lyapunov spectra at randomly picked phase space locations and a direct comparison between Hamiltonian evolution and phase space averages indicate that the regular regions of phase space are not significant at large system sizes. While the continuum limit of our model is itself a singular limit of the integrable sinh Gordon theory, we do not see any evidence for the kind of nonergodicity famously seen in the work of Fermi et al. Finally, we examine the chain with particles confined to two-dimensional stadia where the individual stadium is already chaotic and find a much more chaotic phase space at small system sizes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Two-way coupled particle-turbulence interaction: Effect of numerics and resolution on fluid and particle statistics

Euler-Lagrange point-particle simulation has emerged as a premier methodology for studying dispersed particle-laden flows. This method's popularity stems from its ability to resolve fine-scale fluid structures while also tracking individual particles at reduced cost using an appropriate particle acceleration model. However, the point-particle model has known convergence issues in that refinement of the fluid grid can lead to changes in the predicted statistics. The reasons for non-convergence are twofold: the point-particle two-way coupling force in the Navier-Stokes equations requires a numerical regularization and, without careful implementation, yields a singular force on the fluid with grid refinement. The second factor that yields grid-dependent statistics is that the point-particle force model in general depends on the undisturbed fluid velocity. When the undisturbed fluid velocity is not robustly modeled in a grid-insensitive way, the calculated force for both particles and fluid will be grid-dependent, contaminating their respective statistics. While the first issue regarding regularizing the point-particle source term has received attention in the literature, the consequences of robustly modeling the undisturbed velocity in the context of grid refinement of turbulence has received little attention. Here, we consider decaying homogeneous isotropic turbulence laden with particles at different Stokes numbers. For a given Stokes number, we systematically refine the grid and demonstrate that explicitly modeling the undisturbed fluid velocity yields relative grid insensitivity for the energy of the particle and fluid phases, as well as acceleration of the particles. We also demonstrate that an appropriately defined dissipation rate is also grid-insensitive when an undisturbed fluid velocity correction is used. In contrast, when the undisturbed fluid velocity is modeled using the conventional approach of interpolating the local fluid velocity to the particle location, we show this procedure yields divergent statistics with grid refinement. In particular, we show that higher-order interpolation of the fluid velocity in two-way coupled problems is worse than lower-order interpolation, in the absence of a correction procedure to estimate the undisturbed fluid velocity. We also examine velocity derivative statistics of the fluid phase and demonstrate that these statistics are not in general convergent even when the undisturbed fluid velocity is explicitly modeled. Collectively, the observations in this paper are used to present a philosophy on the types of questions which are answerable with point-particle methods.

42 ENGINEERING↗

Hydrodynamics of two-dimensional compressible fluid with broken parity: Variational principle and free surface dynamics in the absence of dissipation

In this paper, we consider an isotropic compressible nondissipative fluid with broken parity subject to free surface boundary conditions in two spatial dimensions. The hydrodynamic equations describing the bulk dynamics of the fluid and the free surface boundary conditions depend explicitly on the parity-breaking nondissipative odd viscosity term. We construct an effective action which gives both bulk hydrodynamic equations and free surface boundary conditions. The free surface boundary conditions require an additional boundary term in the action which resembles a 1+1D chiral boson field coupled to the background geometry. We solve the linearized hydrodynamic equations for the deep water case and derive the dispersion of chiral surface waves. We show that in the long-wavelength limit the flow profile exhibits an oscillating vortical boundary layer near the free surface. The layer thickness is controlled by the ratio between the odd viscosity (ν o ) and the sound velocity (c s ), δ ~ ν o /c s . In the incompressible limit, c s → ∞, the vortical boundary layer becomes singular with the vorticity within the layer diverging as ω ~ c s . The boundary layer is formed by odd viscosity coupling the divergence of velocity ∇ • v to vorticity ∇ × v. It results in nontrivial chiral free surface dynamics even in the absence of external forces. The structure of the odd-viscosity-induced boundary layer is very different from the conventional free surface boundary layer associated with dissipative shear viscosity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Deep learning for in situ data compression of large turbulent flow simulations

As the size of turbulent flow simulations continues to grow, in situ data compression is becoming increasingly important for visualization, analysis, and restart checkpointing. For these applications, single-pass compression techniques with low computational and communication overhead are crucial. In this paper we present a deep-learning approach to in situ compression using an autoencoder architecture that is customized for three-dimensional turbulent flows and is well suited for contemporary heterogeneous computing resources. The autoencoder is compared against a recently introduced randomized single-pass singular value decomposition (SVD) for three different canonical turbulent flows: decaying homogeneous isotropic turbulence, a Taylor-Green vortex, and turbulent channel flow. Our proposed fully convolutional autoencoder architecture compresses turbulent flow snapshots by a factor of 64 with a single pass, allows for arbitrarily sized input fields, is cheaper to compute than the randomized single-pass SVD for typical simulation sizes, performs well on unseen flow configurations, and has been made publicly available. The results reported here show that the autoencoder dramatically outperforms a randomized single-pass SVD with similar compression ratio and yields comparable performance to a higher-rank decomposition with an order of magnitude less compression in regard to preserving a number of important statistical quantities such as turbulent kinetic energy, enstrophy, and Reynolds stresses.

97 MATHEMATICS AND COMPUTING↗

Learning unknown physics of non-Newtonian fluids

Here, we use physics-informed neural networks (PINNs) to learn viscosity models of two non-Newtonian systems (polymer melts and suspensions of particles) using only velocity measurements. For synthetic velocity data generated with the power-law viscosity model, the PINN-inferred viscosity model agrees with the analytical model for shear rates with large absolute values but deviates for shear rates near zero where the analytical model has an unphysical singularity. Once the viscosity model is learned the PINN method can solve the momentum conservation equation using only the boundary conditions.

42 ENGINEERING↗

Global Phase Diagram of a Spin-Orbital Kondo Impurity Model and the Suppression of Fermi-Liquid Scale

Many correlated metallic materials are described by Landau Fermi-liquid theory at low energies, but for Hund metals the Fermi-liquid coherence scale $T_{FL}$ is found to be surprisingly small. In this Letter, we study the simplest impurity model relevant for Hund metals, the three-channel spin-orbital Kondo model, using the numerical renormalization group (NRG) method and compute its global phase diagram. In this framework, $T_{FL}$ becomes arbitrarily small close to two new quantum critical points that we identify by tuning the spin or spin-orbital Kondo couplings into the ferromagnetic regimes. We find quantum phase transitions to a singular Fermi-liquid or a novel non-Fermi-liquid phase. The new non-Fermi-liquid phase shows frustrated behavior involving alternating overscreenings in spin and orbital sectors, with universal power laws in the spin (ω -1/5 ), orbital (ω 1/ 5 ) and spin-orbital (ω 1 ) dynamical susceptibilities. These power laws, and the NRG eigenlevel spectra, can be fully understood using conformal field theory arguments, which also clarify the nature of the non-Fermi-liquid phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Intrinsic Anomalous Hall Conductivity in a Nonuniform Electric Field

We study how the intrinsic anomalous Hall conductivity is modified in two-dimensional crystals with broken time-reversal symmetry due to weak inhomogeneity of the applied electric field. Focusing on a clean noninteracting two-band system without band crossings, we derive the general expression for the Hall conductivity at small finite wave vector q to order q 2 , which governs the Hall response to the second gradient of the electric field. Using the Kubo formula, we show that the answer can be expressed through the Berry curvature, Fubini-Study quantum metric, and the rank-3 symmetric tensor which is related to the quantum geometric connection and physically corresponds to the gauge-invariant part of the third cumulant of the position operator. We further compare our results with the predictions made within the semiclassical approach. By deriving the semiclassical equations of motion, we reproduce the result obtained from the Kubo formula in some limits. We also find, however, that the conventional semiclassical description in terms of the definite position and momentum of the electron is not fully consistent because of singular terms originating from the Heisenberg uncertainty principle. We thus present a clear example of a case when the semiclassical approach inherently suffers from the uncertainty principle, implying that it should be applied to systems in nonuniform fields with extra care.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Monolayer CrCl 3 as an Ideal Test Bed for the Universality Classes of 2D Magnetism

The monolayer halides CrX 3 (X=Cl, Br, I) attract significant attention for realizing 2D magnets with genuine long-range order (LRO), challenging the Mermin-Wagner theorem. In this work, we show that monolayer CrCl 3 has the unique benefit of exhibiting tunable magnetic anisotropy upon applying a compressive strain. This opens the possibility to use CrCl 3 for producing and studying both ferromagnetic and antiferromagnetic 2D Ising-type LRO as well as the Berezinskii-Kosterlitz-Thouless (BKT) regime of 2D magnetism with quasi-LRO. Using state-of-the-art density functional theory, we explain how realistic compressive strain could be used to tune the monolayer's magnetic properties so that it could exhibit any of these phases. Building on large-scale quantum Monte Carlo simulations, we compute the phase diagram of strained CrCl 3 , as well as the magnon spectrum with spin-wave theory. Our results highlight the eminent suitability of monolayer CrCl 3 to achieve very high BKT transition temperatures, around 50 K, due to their singular dependence on the weak easy-plane anisotropy of the material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Charge Density Waves and Electronic Properties of Superconducting Kagome Metals

We report that kagome metals AV 3 Sb 5 (A = K, Rb, and Cs) exhibit intriguing superconductivity below 0.9 ~ 2.5 K, a charge density wave (CDW) transition around 80 ~ 100 K, and $\mathbb{Z}$ 2 topological surface states. The nature of the CDW phase and its relation to superconductivity remains elusive. In this work, we investigate the electronic and structural properties of CDW by first-principles calculations. We reveal an inverse Star of David deformation as the 2 x 2 x 2 CDW ground state of the kagome lattice. The kagome lattice shows softening breathing-phonon modes, indicating the structural instability. However, electrons play an essential role in the CDW transition via Fermi surface nesting and van Hove singularity. The inverse Star of David structure agrees with recent experiments by scanning tunneling microscopy (STM). The CDW phase inherits the nontrivial $\mathbb{Z}$ 2 -type topological band structure. Further, we find that the electron-phonon coupling is too weak to account for the superconductivity T c in all three materials. It implies the existence of unconventional pairing of these kagome metals. Our results provide essential knowledge toward understanding the superconductivity and topology in kagome metals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Integrability Breaking from Backscattering

Herein we analyze the onset of diffusive hydrodynamics in the one-dimensional hard-rod gas subject to stochastic backscattering. While this perturbation breaks integrability and leads to a crossover from ballistic to diffusive transport, it preserves infinitely many conserved quantities corresponding to even moments of the velocity distribution of the gas. In the limit of small noise, we derive the exact expressions for the diffusion and structure factor matrices, and show that they generically have off diagonal components. We find that the particle density structure factor is non-Gaussian and singular near the origin, with a return probability showing logarithmic deviations from diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Vacancy Spectroscopy of Non-Abelian Kitaev Spin Liquids

Spin vacancies in the non-Abelian Kitaev spin liquid are known to harbor Majorana zero modes, potentially enabling topological quantum computing at elevated temperatures. Here, we study the spectroscopic signatures of such Majorana zero modes in a scanning tunneling setup where a non-Abelian Kitaev spin liquid with a finite density of spin vacancies forms a tunneling barrier between a tip and a substrate. In this study, our key result is a well-defined peak close to zero bias voltage in the derivative of the tunneling conductance whose voltage and intensity both increase with the density of vacancies. This “quasi-zero-voltage peak” is identified as the closest analog of the zero-voltage peak observed in topological superconductors that additionally reflects the fractionalized nature of spin-liquid-based Majorana zero modes. We further highlight a single-fermion Van Hove singularity at a higher voltage that reveals the energy scale of the emergent Majorana fermions in the Kitaev spin liquid. Our proposed signatures are within reach of current experiments on the candidate material α-RuCl 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗