Longest Path Search

When are fields bounded over all possible designs?

Posted 2022-04-17

There are a few important questions that popped up while writing many of the bounds present in this paper. One of the main questions that popped up repeatedly was: when are the solutions to the physics equation bounded over all possible designs?

One specific case in which this question is useful is in the field bounds paragraph in pages 18 and 19 of the paper, though more generally this question can also help answer a number of other important results (which we do not mention here).

Unfortunately, this post will end in a bit of a disappointing note: the result given here depends on some condition which is likely not easy to check in practice. On the other hand, it does lead to a suggestive definition I have never seen before of an "elementwise nonexpansive operator." I would be quite curious to see if anyone had any references!

Problem formulation

In this case, we will again focus on the diagonal physics equation. Here, the physics equation is:

(A+diag(θ))z=b, (A + \mathbf{diag}(\theta))z = b,

where θRn\theta \in \mathbf{R}^n are the design parameters (usually the permittivities in many problems) while zRnz \in \mathbf{R}^n are the fields. In general, we are only allowed to choose parameters within a certain range, so we will write this as

1θ1, -1 \le \theta \le 1,

without loss of generality. (In particular, if θ\theta is constrained to lie within any range, we can always rescale the physics equation to make θ\theta lie between 1-1 and 11. For more details on how to do this, see section 1.2 of the paper.)

Eliminating the design variable

Taking this formulation, we can "eliminate" the design parameter. In other words, we will write a number of equations, depending only on the variable zz, such that, when zz satisfies all the equations, there exists some design θ[1,1]n\theta \in [-1, 1]^n which makes the physics equation true.

To do this, note that, we can take the initial physics equation and rearrange it as follows:

Azb=diag(θ)z. Az - b = -\mathbf{diag}(\theta)z.

Taking the elementwise absolute value of both sides gives

Azb=diag(θ)z. |Az - b| = |\mathbf{diag}(\theta)z|.

(Here, we interpret |\cdot| to be elementwise.) Because θ1|\theta| \le 1, we can see that

diag(θ)zz, |\mathbf{diag}(\theta)z| \le |z|,

so

Azbz. |Az - b| \le |z|.

In fact, diag(θ)zz|\mathbf{diag}(\theta)z| \le |z| is true if, and only if, θ1|\theta| \le 1. Meaning the inequality we just derived, depending only on zz, is true if, and only if, there exists some design θ\theta satisfying θ1|\theta| \le 1 that makes the original physics equation true.

The nice part about this equation is that it encapsulates all of the important parts of the problem in a simple-to-reason-about format. (It also suggests some interesting heuristics, but that's for another time!)

When are fields bounded?

Now, finally, to answer the question! At least partially.

The proof technique here is relatively simple: we will

Built with Franklin.jl and Julia.