Electromagnetic vortex lines riding atop null solutions of the Maxwell equations

New method of introducing vortex lines of the electromagnetic field is outlined. The vortex lines arise when a complex Riemann-Silberstein vector $({\bm E} + i{\bm B})/\sqrt{2}$ is multiplied by a complex scalar function $\phi$. Such a multiplication may lead to new solutions of the Maxwell equations only when the electromagnetic field is null, i.e. when both relativistic invariants vanish. In general, zeroes of the $\phi$ function give rise to electromagnetic vortices. The description of these vortices benefits from the ideas of Penrose, Robinson and Trautman developed in general relativity.Comment: NATO Workshop on Singular Optics 2003 To appear in Journal of Optics

A second derivative SQP method: theoretical issues

Sequential quadratic programming (SQP) methods form a class of highly efficient algorithms for solving nonlinearly constrained optimization problems. Although second derivative information may often be calculated, there is little practical theory that justifies exact-Hessian SQP methods. In particular, the resulting quadratic programming (QP) subproblems are often nonconvex, and thus finding their global solutions may be computationally nonviable. This paper presents a second-derivative SQP method based on quadratic subproblems that are either convex, and thus may be solved efficiently, or need not be solved globally. Additionally, an explicit descent-constraint is imposed on certain QP subproblems, which “guides” the iterates through areas in which nonconvexity is a concern. Global convergence of the resulting algorithm is established

A second derivative SQP method: local convergence

In [19], we gave global convergence results for a second-derivative SQP method for minimizing the exact ℓ1-merit function for a fixed value of the penalty parameter. To establish this result, we used the properties of the so-called Cauchy step, which was itself computed from the so-called predictor step. In addition, we allowed for the computation of a variety of (optional) SQP steps that were intended to improve the efficiency of the algorithm. \ud \ud Although we established global convergence of the algorithm, we did not discuss certain aspects that are critical when developing software capable of solving general optimization problems. In particular, we must have strategies for updating the penalty parameter and better techniques for defining the positive-definite matrix Bk used in computing the predictor step. In this paper we address both of these issues. We consider two techniques for defining the positive-definite matrix Bk—a simple diagonal approximation and a more sophisticated limited-memory BFGS update. We also analyze a strategy for updating the penalty paramter based on approximately minimizing the ℓ1-penalty function over a sequence of increasing values of the penalty parameter.\ud \ud Algorithms based on exact penalty functions have certain desirable properties. To be practical, however, these algorithms must be guaranteed to avoid the so-called Maratos effect. We show that a nonmonotone varient of our algorithm avoids this phenomenon and, therefore, results in asymptotically superlinear local convergence; this is verified by preliminary numerical results on the Hock and Shittkowski test set

A second-derivative trust-region SQP method with a "trust-region-free" predictor step

In (NAR 08/18 and 08/21, Oxford University Computing Laboratory, 2008) we introduced a second-derivative SQP method (S2QP) for solving nonlinear nonconvex optimization problems. We proved that the method is globally convergent and locally superlinearly convergent under standard assumptions. A critical component of the algorithm is the so-called predictor step, which is computed from a strictly convex quadratic program with a trust-region constraint. This step is essential for proving global convergence, but its propensity to identify the optimal active set is Paramount for recovering fast local convergence. Thus the global and local efficiency of the method is intimately coupled with the quality of the predictor step.\ud \ud In this paper we study the effects of removing the trust-region constraint from the computation of the predictor step; this is reasonable since the resulting problem is still strictly convex and thus well-defined. Although this is an interesting theoretical question, our motivation is based on practicality. Our preliminary numerical experience with S2QP indicates that the trust-region constraint occasionally degrades the quality of the predictor step and diminishes its ability to correctly identify the optimal active set. Moreover, removal of the trust-region constraint allows for re-use of the predictor step over a sequence of failed iterations thus reducing computation. We show that the modified algorithm remains globally convergent and preserves local superlinear convergence provided a nonmonotone strategy is incorporated

Metric Perturbation Approach to Gravitational Waves in Isotropic Cosmologies

Gravitational waves in isotropic cosmologies were recently studied using the gauge-invariant approach of Ellis-Bruni. We now construct the linearised metric perturbations of the background Robertson-Walker space-time which reproduce the results obtained in that study. The analysis carried out here also facilitates an easy comparison with Bardeen.Comment: 29 pages, Latex file, accepted for publication in Physical Review

Designing a Green Roof for Ireland

A model is presented for the gravity-driven flow of rainwater descending through the soil layer of a green roof, treated as a porous medium on a at permeable surface representing an efficient drainage layer. A fully saturated zone is shown to occur. It is typically a thin layer, relative to the total soil thickness, and lies at the bottom of the soil layer. This provides a bottom boundary condition for the partially saturated upper zone. It is shown that after the onset of rainfall, well-defined fronts of water can descend through the soil layer. Also the rainwater flow is relatively quick compared with the moisture uptake by the roots of the plants in the roof. In a separate model the exchanges of water are described between the (smaller-scale) porous granules of soil, the roots and the rainwater in the inter-granule pores

Shear-Free Gravitational Waves in an Anisotropic Universe

We study gravitational waves propagating through an anisotropic Bianchi I dust-filled universe (containing the Einstein-de-Sitter universe as a special case). The waves are modeled as small perturbations of this background cosmological model and we choose a family of null hypersurfaces in this space-time to act as the histories of the wavefronts of the radiation. We find that the perturbations we generate can describe pure gravitational radiation if and only if the null hypersurfaces are shear-free. We calculate the gauge-invariant small perturbations explicitly in this case. How these differ from the corresponding perturbations when the background space-time is isotropic is clearly exhibited.Comment: 32 pages, accepted for publication in Physical Review