Galactic warps are common in disk galaxies. While often attributed to galaxy--galaxy tides, a non-spherical dark matter (DM) halo has also been proposed as a driver of disk warping. We investigate links among warp morphology, satellite distribution, and large-scale structure using the Sloan Digital Sky Survey catalog of warped disks compiled by Zee et al.\ (2022). Warps are classified into 244 S and 127 U types, hosting 1,373 and 740 satellites, respectively, and are compared to an unwarped control matched in stellar mass, redshift, and local density. As an indirect, population-level proxy for the host halo's shape and orientation, we analyze the stacked spatial distribution of satellites. Warped hosts show a significant anisotropy: an excess at $45^{\circ}<φ<90^{\circ}$ (measured from the host major axis), peaking at $P(φ)\simeq 0.003$, versus nearly isotropic controls. Satellites of S-type warps preferentially align with the nearest cosmic filament, whereas those of U-type warps are more often perpendicular. The incidence of warps increases toward filaments ($r_{\rm fil}<4,{\rm Mpc},h^{-1}$), while the number of satellites around warped hosts remains approximately consta
The field of warp research has been dominated by analytical methods to investigate potential solutions. However, these approaches often favor simple metric forms that facilitate analysis but ultimately limit the range of exploration of novel solutions. So far the proposed solutions have been unphysical, requiring energy condition violations and large energy requirements. To overcome the analytical limitations in warp research, we introduce Warp Factory: a numerical toolkit designed for modeling warp drive spacetimes. By leveraging numerical analysis, Warp Factory enables the examination of general warp drive geometries by evaluating the Einstein field equations and computing energy conditions. Furthermore, this comprehensive toolkit provides the determination of metric scalars and insightful visualizations in both 2D and 3D, offering a deeper understanding of metrics and their corresponding stress-energy tensors. The paper delves into the methodology employed by Warp Factory in evaluating the physicality of warp drive spacetimes and highlights its application in assessing commonly modeled warp drive metrics. By leveraging the capabilities of Warp Factory, we aim to further warp dri
We present the first local simulations of disc breaking/tearing in a warped accretion disc. Warps can arise due to a misalignment between the disc and the rotation axis of the central object, or a misalignment with the orbital plane of a binary (or planetary) companion. Warped discs can break into rings, as found in observations of circumbinary protoplanetary discs and global simulations of tilted discs around spinning black holes. In this work we isolate the mechanism of disc breaking in high-resolution, quasi-2D, local (shearing box), hydrodynamic simulations of a Keplerian disc. We consider the evolution of a free (unforced) warp in the wavelike ($α< H/r$) regime. At large warp amplitudes ($ψ_{\text{max}} \gtrsim 1$) the disc breaks into four rings on timescales of around 20 orbits which are separated by gaps of around $\sim 10H_0$. The warp exhibits a rich tapestry of small-scale dynamics, including horizontal sloshing motions, vertical oscillations or bouncing, warp steepening, and shocks. The shocks act as a source of enhanced dissipation which facilitates gap opening and thus disc breaking. At smaller warp amplitudes $ψ_{\text{max}} \lesssim 1$, for which we also develop
Classical Cepheids (CCs) are important probes for the large-scale warp structure of the Milky Way. Using Gaia DR3 CCs, we establish an optimal time-dependent warp model, where the warp height increases with radius following a power-law, the line of nodes (LONs) exhibit linear twisting with radius, following a leading spiral pattern, and the LONs undergo prograde evolution over time. Structurally, we identify significant warp features in the $5-9$ kpc region of the Galactic disk, where the warp model performs better than the flat model. Beyond 15 kpc, the model with the second Fourier term does not fit the observations well, whereas the model with twisted LONs better matches the data. Kinematically, we derived expressions for the vertical velocities using direct differentiation and then calculated the precession rates for each CC. Our results intuitively indicate a nearly uniform and low warp precession rate of $ω= 4.86 \pm (0.88)_{stat} \pm (2.14)_{sys}$ km s$^{-1}$ kpc$^{-1}$ beyond 12.5 kpc, in agreement with classical kinematic estimates. Based on these findings, we propose a simple yet comprehensive time-dependent warp model, $Z_{w}(t) = 0.00019R^{3.08}\sin(φ- (3.87R-41.79 + 4.
The last few decades of warp drive research have focused on analytic methods to explore warp solutions to Einstein's field equations. These analytic solutions tend to favor simple metric forms which are easier to analyze but limit the space of exploration. In addition, all solutions to date have involved unphysical qualities, such as negative energy, violation of energy conditions, or enormous energy requirements. In an effort to explore the space of physically meaningful warp drives, the Advanced Propulsion Laboratory (APL) at Applied Physics has developed Warp Factory, a toolkit written in MATLAB for numerically analyzing and optimizing warp drive geometries. Warp Factory consists of a series of three primary modules: the solver, the analyzer, and the optimizer. Together, these modules allow users to solve the Einstein field equations, compute energy conditions and scalars, and perturbatively optimize general metrics. Finally, the toolkit offers insightful 2D and 3D visualizations of general metrics and stress-energy tensors. The methods used in Warp Factory, along with their application in evaluating and optimizing common metrics, are discussed. With Warp Factory, APL hopes to a
Following the work of H. Ellis, we study warp drives in the gravitational field of a Schwarzschild black hole. We find that as long as the warp drive crosses the black hole horizon at a subluminal speed, the horizon would be effectively absent inside the warp bubble. Moreover, we discover that the black hole's gravitational field can alleviate the violations of the weak energy condition (WEC) and the null energy condition (NEC) and therefore decrease the amount of negative energy required to sustain a warp drive, which may be instrumental for creating microscopic warp drives in lab experiments. We also consider the thermodynamics of a warp bubble interacting with a black hole and point out some paradoxes that may indicate a gap in our understanding of them from the thermodynamic point of view.
In recent years star formation has been discovered in the Milky Way's warp. These stars formed in the warp (warp stars) must eventually settle into the plane of the disc. We use an $N$-body$+$smooth particle hydrodynamics model of a warped galaxy to study how warp stars settle into the disc. By following warp stars in angular momentum space, we show that they first tilt to partially align with the main disc in a time scale of $\sim1$ Gyr. Then, once differential precession halts this process, they phase mix into an axisymmetric distribution on a time scale of $\sim 6$ Gyr. The warp stars end up contaminating the geometric thick disc. Because the warp in our fiducial simulation is growing, the {\it warp stars} settle to a distribution with a negative vertical age gradient as younger stars settle further from the mid-plane. While vertically extended, warp star orbits are still nearly circular and they are therefore subject to radial migration, with a net movement inwards. As a result warp stars can be found throughout the disc. The density distribution of a given population of warp stars evolves from a torus to an increasingly centrally filled-in density distribution. Therefore we ar
It is commonly accepted that superluminal travel may be used to facilitate time travel. This is a purely special-relativistic argument, using the fact that for observers in two frames of reference, separated by a spacelike interval, the non-causal (spacelike) future of one observer includes part of the causal past of the other. In this paper we provide a concrete realization of this argument in a curved general-relativistic spacetime, using warp drives as the means of faster-than-light travel. By generalizing the usual warp drive metric to allow for a non-unit lapse function, we allow the warp drive to switch between reference frames in a purely geometric way. With an additional modification allowing the warp drive to have compact support, this permits us to glue two warp drives together to construct a closed timelike geodesic, such that a test particle following the geodesics of the two warp drives travels back to its own past. This provides a precise mathematical model for the connection between faster-than-light travel and time travel in general relativity, and the first such model to be explicitly formulated using two warp drives. We also give a detailed discussion of weak ener
We propose a correspondence between the Morris--Thorne wormhole metric and a warp drive metric, which generalizes an earlier result by H. Ellis regarding the Schwarzschild black hole metric and makes it possible to embed a warp drive in a wormhole background. We demonstrate that in order to do that, one needs to also generalize the Natario--Alcubierre definition of warp drive and introduce nonzero intrinsic curvature. However, we also find out that in order to be traversable by a warp drive, the wormhole should have a horizon: in other words, humanly traversable wormholes cannot be traversed by a warp drive, and vice versa. We also discuss possible loopholes in this "no-go" theorem
This paper presents a performance model tailored for warp specialization kernels, focusing on factors such as warp size, tilling size, input matrix size, memory bandwidth, and thread divergence. Our model offers accurate predictions of execution time by leveraging differential equations validated through simulations and experiments. The insights gained from this model not only enhance our understanding of warp specialization techniques but also have practical implications for optimizing GPU-accelerated applications through compiler optimizations, kernel parameter tuning, and algorithm design.
We determine the influence of the Milky Way's warp on the kinematics of stars across the disc, and therefore measure its precession rate and line of nodes under different assumptions. We do this by applying Jeans' first equation to a model of a rigidly precessing warp. The predictions of these models are fit to the average vertical velocities of stars with measured line-of-sight velocities in Gaia DR3 data. We test models in which the warp's line of nodes and precession speed are fixed, and models in which they are allowed to vary linearly with radius. We also test models in which the velocity of stars radially in the disc is included in Jeans' equation. The kinematic data is best fit by models with a line of nodes that is 40 degrees offset from the Sun's Galactic azimuth, significantly leading the line of nodes found from the positions of stars. These models have a warp precession speed of around 13 km/s/kpc in the direction of Galactic rotation, close to other recent estimates. We find that including the velocity of stars radially in the disc in our kinematic model leads to a significantly worse fit to the data, and implausible warp parameters. We conclude that the Milky Way's wa
What happens when a warp bubble has mass? This seemingly innocent question forces one to carefully formalize exactly what one means by a warp bubble, exactly what one means by having the warp bubble "move" with respect to the fixed stars, and forces one to more carefully examine the notion of mass in warp-drive spacetimes. This is the goal of the present article. In this process, we will see that often-made throw-away comments regarding "payloads" are even simpler than commonly assumed, while there are two further, distinct yet subtle ways in which a mass can appear in connection with a warp drive space-time: One, that the warp bubble (not its payload) has the mass; two, that the mass is a background feature in front of which the warp drive moves. For simplicity, we consider generic Natário warp drives with zero-vorticity flow field. The resulting spacetimes are sufficiently simple to allow an exact and fully explicit computation of all of the stress-energy components, and verify that (as expected) the null energy condition (NEC) is violated. Likewise the weak, strong, and dominant energy conditions (WEC, SEC, DEC) are violated. Indeed, this confirms the community's folk wisdom, an
Modern GPUs synchronize threads grouped in a warp at every instruction. These results in improving SIMD efficiency and makes sharing fetch and decode resources possible. The number of threads included in each warp (or warp size) affects divergence, synchronization overhead and the efficiency of memory access coalescing. Small warps reduce the performance penalty associated with branch and memory divergence at the expense of a reduction in memory coalescing. Large warps enhance memory coalescing significantly but also increase branch and memory divergence. Dynamic workload behavior, including branch/memory divergence and coalescing, is an important factor in determining the warp size returning best performance. Optimal warp size can vary from one workload to another or from one program phase to the next. Based on this observation, we propose Dynamic Warp Resizing (DWR). DWR takes innovative microarchitectural steps to adjust warp size during runtime and according to program characteristics. DWR outperforms static warp size decisions, up to 1.7X to 2.28X, while imposing less than 1% area overhead. We investigate various alternative configurations and show that DWR performs better for
Achieving high efficiency in modern photorealistic rendering hinges on using Monte Carlo sampling distributions that closely approximate the illumination integral estimated for every pixel. Samples are typically generated from a set of simple distributions, each targeting a different factor in the integrand, which are combined via multiple importance sampling. The resulting mixture distribution can be far from the actual product of all factors, leading to sub-optimal variance even for direct-illumination estimation. We present a learning-based method that uses normalizing flows to efficiently importance sample illumination product integrals, e.g., the product of environment lighting and material terms. Our sampler composes a flow head warp with an emitter tail warp. The small conditional head warp is represented by a neural spline flow, while the large unconditional tail is discretized per environment map and its evaluation is instant. If the conditioning is low-dimensional, the head warp can be also discretized to achieve even better performance. We demonstrate variance reduction over prior methods on a range of applications comprising complex geometry, materials and illumination.
There are a number of design decisions that impact a GPU's performance. Among such decisions deciding the right warp size can deeply influence the rest of the design. Small warps reduce the performance penalty associated with branch divergence at the expense of a reduction in memory coalescing. Large warps enhance memory coalescing significantly but also increase branch divergence. This leaves designers with two choices: use a small warps and invest in finding new solutions to enhance coalescing or use large warps and address branch divergence employing effective control-flow solutions. In this work our goal is to investigate the answer to this question. We analyze warp size impact on memory coalescing and branch divergence. We use our findings to study two machines: a GPU using small warps but equipped with excellent memory coalescing (SW+) and a GPU using large warps but employing an MIMD engine immune from control-flow costs (LW+). Our evaluations show that building coalescing-enhanced small warp GPUs is a better approach compared to pursuing a control-flow enhanced large warp GPU.
A process for using curvature invariants is applied to evaluate the metrics for the Alcubierre and the Natario warp drives at a constant velocity.Curvature invariants are independent of coordinate bases, so plotting these invariants will be free of coordinate mapping distortions. As a consequence, they provide a novel perspective into complex spacetimes such as warp drives. Warp drives are the theoretical solutions to Einstein's field equations that allow the possibility for faster-than-light (FTL) travel. While their mathematics is well established, the visualisation of such spacetimes is unexplored. This paper uses the methods of computing and plotting the warp drive curvature invariants to reveal these spacetimes. The warp drive parameters of velocity, skin depth and radius are varied individually and then plotted to see each parameter's unique effect on the surrounding curvature. For each warp drive, this research shows a safe harbor and how the shape function forms the warp bubble. The curvature plots for the constant velocity Natario warp drive do not contain a wake or a constant curvature indicating that these are unique features of the accelerating Natario warp drive.
This short technical note presents a formal generalization of the Time Warp Edit Distance (TWED) proposed by Marteau (2009) to arbitrary metric spaces. By viewing both the observation and temporal domains as metric spaces $(X, d)$ and $(T, Δ)$, we define a Generalized TWED (GTWED) that remains a true metric under mild assumptions. We provide self-contained proofs of its metric properties and show that the classical TWED is recovered as a special case when $X = \mathbb{R}^d$, $T \subset \mathbb{R}$, and $g(x) = x$. This note focuses on the theoretical structure of GTWED and its implications for extending elastic distances beyond time series, which enables the use of TWED-like metrics on sequences over arbitrary domains such as symbolic data, manifolds, or embeddings.
"Warp drive" spacetimes are useful as "gedanken-experiments" that force us to confront the foundations of general relativity, and among other things, to precisely formulate the notion of "superluminal" communication. We verify the non-perturbative violation of the classical energy conditions of the Alcubierre and Natario warp drive spacetimes and apply linearized gravity to the weak-field warp drive, testing the energy conditions to first and second order of the non-relativistic warp-bubble velocity. We are primarily interested in a secondary feature of the warp drive that has not previously been remarked upon, if it could be built, the warp drive would be an example of a "reaction-less drive". For both the Alcubierre and Natario warp drives we find that the occurrence of significant energy condition violations is not just a high-speed effect, but that the violations persist even at arbitrarily low speeds. An interesting feature of this construction is that it is now meaningful to place a finite mass spaceship at the center of the warp bubble, and compare the warp field energy with the mass-energy of the spaceship. There is no hope of doing this in Alcubierre's original version of
A novel warp for natural image stitching is proposed that utilizes the property of cylindrical warp and a horizontal pixel selection strategy. The proposed ratio-preserving half-cylindrical warp is a combination of homography and cylindrical warps which guarantees alignment by homography and possesses less projective distortion by cylindrical warp. Unlike previous approaches applying cylindrical warp before homography, we use partition lines to divide the image into different parts and apply homography in the overlapping region while a composition of homography and cylindrical warps in the non-overlapping region. The pixel selection strategy then samples the points in horizontal and reconstructs the image via interpolation to further reduce horizontal distortion by maintaining the ratio as similarity. With applying half-cylindrical warp and horizontal pixel selection, the projective distortion in vertical and horizontal is mitigated simultaneously. Experiments show that our warp is efficient and produces a more natural-looking stitched result than previous methods.
We investigate the dynamical interaction of galactic warps with the surrounding dark matter halo, using analytic perturbation theory. A precessing warp induces a density wake in the collisionless dark matter, which acts back on the original warp, transferring energy and angular momentum between the warp and halo (dynamical friction). In most cases dynamical friction damps the warp, but in unusual circumstances (such as a halo that rotates in the same direction as the precession of the warp, or a warp in the equatorial plane of an axisymmetric prolate halo) friction can excite a warp. The damping/excitation time is usually short compared to the Hubble time for realistic systems. Thus most warps cannot be primordial; they must be maintained by some ongoing excitation mechanism.