An official website of the United States government
We have developed a Liouville partial-differential-equation (PDE)-based method for computing complex-valued eikonals in real phase space in the multivalued sense in attenuating media with frequency-independent qualify factors, where the new method computes the real and imaginary parts of the complex-valued eikonal in two steps by solving Liouville equations in real phase space. Because the earth is composed of attenuating materials, seismic waves usually attenuate so that seismic data processing calls for properly treating the resulting energy losses and phase distortions of wave propagation. In the regime of high-frequency asymptotics, the complex-valued eikonal is one essential ingredient for describing wave propagation in attenuating media because this unique quantity summarizes two wave properties into one function: Its real part describes the wave kinematics and its imaginary part captures the effects of phase dispersion and amplitude attenuation. Because some popular ordinary-differential-equation (ODE)-based ray-tracing methods for computing complex-valued eikonals in real space distribute the eikonal function irregularly in real space, we are motivated to develop PDE-based Eulerian methods for computing such complex-valued eikonals in real space on regular meshes. Therefore, we solved novel paraxial Liouville PDEs in real phase space so that we can compute the real and imaginary parts of the complex-valued eikonal in the multivalued sense on regular meshes. We call the resulting method the Liouville PDE method for complex-valued multivalued eikonals in attenuating media; moreover, this new method provides a unified framework for Eulerianizing several popular approximate real-space ray-tracing methods for complex-valued eikonals, such as viscoacoustic ray tracing, real viscoelastic ray tracing, and real elastic ray tracing. In addition, we also provide Liouville PDE formulations for computing multivalued ray amplitudes in a weakly viscoacoustic medium. Numerical examples, including a synthetic gas-cloud model, demonstrate that our methods yield highly accurate complex-valued eikonals in the multivalued sense.
more »
« less
Exploring the possibility of a complex-valued non-Gaussianity measure for quantum states of light
We consider a quantity that is the differential relative entropy between a generic Wigner function and a Gaussian one. We prove that said quantity is minimized with respect to its Gaussian argument, if both Wigner functions in the argument of the Wigner differential entropy have the same first and second moments, i.e., if the Gaussian argument is the Gaussian associate of the other, generic Wigner function. Therefore, we introduce the differential relative entropy between any Wigner function and its Gaussian associate and we examine its potential as a non-Gaussianity measure. The proposed, phase-space based non-Gaussianity measure is complex-valued, with its imaginary part possessing the physical meaning of the Wigner function’s negative volume. At the same time, the real part of this measure provides an extra layer of information, rendering the complex-valued quantity a measure of non-Gaussianity, instead of a quantity pertaining only to the negativity of the Wigner function. We prove that the measure (both the real and imaginary parts) is faithful, invariant under Gaussian unitary operations, and find a sufficient condition on its monotonic behavior under Gaussian channels. We provide numerical results supporting the aforesaid condition. In addition, we examine the measure’s usefulness to non-Gaussian quantum state engineering with partial measurements.
more »
« less
- Award ID(s):
- 2122337
- PAR ID:
- 10545421
- Publisher / Repository:
- AIP Publishing
- Date Published:
- Journal Name:
- APL Quantum
- Volume:
- 1
- Issue:
- 3
- ISSN:
- 2835-0103
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Discrete and continuous frames can be considered as positive operator-valued measures (POVMs) that have integral representations using rank-one operators. However, not every POVM has an integral representation. One goal of this paper is to examine the POVMs that have finite-rank integral representations. More precisely, we present a necessary and sufficient condition under which a positive operator-valued measure $$F: \Omega \to B(H)$$ has an integral representation of the form $$F(E) =\sum_{k=1}^{m} \int_{E}\, G_{k}(\omega)\otimes G_{k}(\omega) d\mu(\omega)$$ for some weakly measurable maps $$G_{k} \ (1\leq k\leq m) $$ from a measurable space $$\Omega$$ to a Hilbert space $$\mathcal{H}$$ and some positive measure $$\mu$$ on $$\Omega$$. Similar characterizations are also obtained for projection-valued measures. As special consequences of our characterization we settle negatively a problem of Ehler and Okoudjou about probability frame representations of probability POVMs, and prove that an integral representable probability POVM can be dilated to a integral representable projection-valued measure if and only if the corresponding measure is purely atomic.more » « less
-
In this paper, we propose a controller that stabilizes a holonomic robot with single-integrator dynamics to a target position in a bounded domain, while preventing collisions with convex obstacles. We assume that the robot can measure its own position and heading in a global coordinate frame, as well as its relative position vector to the closest point on each obstacle in its sensing range. The robot has no information about the locations and shapes of the obstacles. We define regions around the boundaries of the obstacles and the domain within which the robot can sense these boundaries, and we associate each region with a virtual potential field that we call a local navigation-like function (NLF), which is only a function of the robot’s position and its distance from the corresponding boundary. We also define an NLF for the remaining free space of the domain, and we identify the critical points of the NLFs. Then, we propose a switching control law that drives the robot along the negative gradient of the NLF for the obstacle that is currently closest, or the NLF for the remaining free space if no obstacle is detected. We derive a conservative upper bound on the tunable parameter of the NLFs that guarantees the absence of locally stable equilibrium points, which can trap the robot, if the obstacles’ boundaries satisfy a minimum curvature condition. We also analyze the convergence and collision avoidance properties of the switching control law and, using a Lyapunov argument, prove that the robot safely navigates around the obstacles and converges asymptotically to the target position. We validate our analytical results for domains with different obstacle configurations by implementing the controller in both numerical simulations and physical experiments with a nonholonomic mobile robot.more » « less
-
Seasonal snowpack is an important predictor of the water resources available in the following spring and early-summer melt season. Total basin snow water equivalent (SWE) estimation usually requires a form of statistical analysis that is implicitly built upon the Gaussian framework. However, it is important to characterize the non-Gaussian properties of snow distribution for accurate large-scale SWE estimation based on remotely sensed or sparse ground-based observations. This study quantified non-Gaussianity using sample negentropy; the Kullback–Leibler divergence from the Gaussian distribution for field-observed snow depth data from the North Slope, Alaska; and three representative SWE distributions in the western USA from the Airborne Snow Observatory (ASO). Snowdrifts around lakeshore cliffs and deep gullies can bring moderate non-Gaussianity in the open, lowland tundra of North Slope, Alaska, while the ASO dataset suggests that subalpine forests may effectively suppress the non-Gaussianity of snow distribution. Thus, non-Gaussianity is found in areas with partial snow cover and wind-induced snowdrifts around topographic breaks on slopes and on other steep terrain features. The snowpacks may be considered weakly Gaussian in coastal regions with open tundra in Alaska and alpine and subalpine terrains in the western USA if the land is completely covered by snow. The wind-induced snowdrift effect can potentially be partitioned from the observed snow spatial distribution guided by its Gaussianity.more » « less
-
Abstract We consider the problem of computing the partition function $$\sum _x e^{f(x)}$$ , where $$f: \{-1, 1\}^n \longrightarrow {\mathbb R}$$ is a quadratic or cubic polynomial on the Boolean cube $$\{-1, 1\}^n$$ . In the case of a quadratic polynomial f , we show that the partition function can be approximated within relative error $$0 < \epsilon < 1$$ in quasi-polynomial $$n^{O(\ln n - \ln \epsilon )}$$ time if the Lipschitz constant of the non-linear part of f with respect to the $$\ell ^1$$ metric on the Boolean cube does not exceed $$1-\delta $$ , for any $$\delta>0$$ , fixed in advance. For a cubic polynomial f , we get the same result under a somewhat stronger condition. We apply the method of polynomial interpolation, for which we prove that $$\sum _x e^{\tilde {f}(x)} \ne 0$$ for complex-valued polynomials $$\tilde {f}$$ in a neighborhood of a real-valued f satisfying the above mentioned conditions. The bounds are asymptotically optimal. Results on the zero-free region are interpreted as the absence of a phase transition in the Lee–Yang sense in the corresponding Ising model. The novel feature of the bounds is that they control the total interaction of each vertex but not every single interaction of sets of vertices.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0219011
