BZPEERSource author record
Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.
Catalog footprint
Published work
We provide a detailed description of the maps associated with spectral interlacing, for rank one perturbations and bordering of symmetric and Hermitian matrices. The arguments rely on standard techniques of nonlinear analysis.
The observed output of an interferometer is the result of interference among the parts of the input light beam traveling along each possible optical path. In complex systems, writing down all these possible optical paths and computing their cumulative effect can become a difficult task. We present an intuitive graph-based method for solving this problem and calculating electric fields within an interferometric setup, classical and quantum. We show how to associate a weighted directed graph to an interferometer and define rules to simplify these associated graphs. Successive application of the rules results in a final graph containing information on the desired field amplitudes. The method is applied to a number of examples in cavity optomechanics and cavity-enhanced interferometers.
Vector fields and line fields, their counterparts without orientations on tangent lines, are familiar objects in the theory of dynamical systems. Among the techniques used in their study, the Morse--Smale decomposition of a (generic) field plays a fundamental role, relating the geometric structure of phase space to a combinatorial object consisting of critical points and separatrices. Such concepts led Forman to a satisfactory theory of discrete vector fields, in close analogy to the continuous case. In this paper, we introduce discrete line fields. Again, our definition is rich enough to provide the counterparts of the basic results in the theory of continuous line fields: a Euler-Poincaré formula, a Morse--Smale decomposition and a topologically consistent cancellation of critical elements, which allows for topological simplification of the original discrete line field.
For large n, take a random n×n permutation matrix and its associated discrete copula Xn. For a,b=0,1,…,n, let yn(na,nb)=n1(Xa,b−nab); define yn:[0,1]2→R by interpolating quadratically on squares of side n1. We prove a Donsker type central limit theorem: nyn approaches a bridged Brownian sheet on the unit square.
According to the Ambrosetti-Prodi theorem, the map F(u)=−Δu−f(u) between appropriate functional spaces is a global fold. Among the hypotheses, the convexity of the function f is required. We show in two different ways that, under mild conditions, convexity is indeed necessary. If f is not convex, there is a point with at least four preimages under F. More, F generically admits cusps among its critical points. We present a larger class of nonlinearities f for which the critical set of F has cusps. The results are true for a class of boundary conditions.
Since the seminal work of Ambrosetti and Prodi, the study of global folds was enriched by geometric concepts and extensions accomodating new examples. We present the advantages of considering fibers, a construction dating to Berger and Podolak's view of the original theorem. A description of folds in terms of properties of fibers gives new perspective to the usual hypotheses in the subject. The text is intended as a guide, outlining arguments and stating results which will be detailed elsewhere.
For Dirichlet boundary conditions on a bounded domain, what happens to the critical set of the Ambrosetti-Prodi operator if the nonlinearity is only a Lipschitz map? It turns out that many properties which hold in the smooth case are preserved, despite of the fact that the operator is not even differentiable at some points. In particular, a global Lyapunov-Schmidt decomposition of great convenience for numerical inversion is still available.
Originally a model for wave propagation on the line, the Toda lattice is a wonderful case study in mechanics and symplectic geometry. In Flaschka's variables, it becomes an evolution given by a Lax pair on the vector space of real, symmetric, tridiagonal matrices. Its very special asymptotic behavior was studied by Moser by introducing norming constants, which play the role of discrete inverse variables in analogy to the solution by inverse scattering of KdV. It is a completely integrable system on the coadjoint orbit of the upper triangular group. Recently, bidiagonal coordinates, which parameterize also non-Jacobi tridiagonal matrices, were used to reduce asymptotic questions to local theory. Larger phase spaces for the Toda lattice lead to the study of isospectral manifolds and different coadjoint orbits. Additionally, the time one map of the associated flow is computed by a familiar algorithm in numerical linear algebra. The text is mostly expositive and self contained, presenting alternative formulations of familiar results and applications to numerical analysis.
A common algorithm for the computation of eigenvalues of real symmetric tridiagonal matrices is the iteration of certain special maps Fσ called shifted QR steps. Such maps preserve spectrum and a natural common domain is TΛ, the manifold of real symmetric tridiagonal matrices conjugate to the diagonal matrix Λ. More precisely, a (generic) shift $s \in \RR$ defines a map Fs:TΛ→TΛ. A strategy $σ: {\cal T}_Λ\to \RR$ specifies the shift to be applied at T so that Fσ(T)=Fσ(T)(T). Good shift strategies should lead to fast deflation: some off-diagonal coordinate tends to zero, allowing for reducing of the problem to submatrices. For topological reasons, continuous shift strategies do not obtain fast deflation; many standard strategies are indeed discontinuous. Practical implementation only gives rise systematically to bottom deflation, convergence to zero of the lowest off-diagonal entry b(T). For most shift strategies, convergence to zero of b(T) is cubic, ∣b(Fσ(T))∣=Θ(∣b(T)∣k) for k=3. The existence of arithmetic progressions in the spectrum of T sometimes implies instead quadratic convergence, k=2. The complete integrability of the Toda lattice and the dynamics at non-smooth points are central to our discussion. The text does not assume knowledge of numerical linear algebra.
We present an algorithm to solve $- \lap u - f(x,u) = g$ with Dirichlet boundary conditions in a bounded domain Ω. The nonlinearities are non-resonant and have finite spectral interaction: no eigenvalue of $-\lap_D$ is an endpoint of $\bar{\partial_2f(Ω,\RR)}$, which in turn only contains a finite number of eigenvalues. The algorithm is based in ideas used by Berger and Podolak to provide a geometric proof of the Ambrosetti-Prodi theorem and advances work by Smiley and Chun for the same problem.
The computation of eigenvalues of real symmetric tridiagonal matrices frequently proceeds by a sequence of QR steps with shifts. We introduce simple shift strategies, functions sigma satisfying natural conditions, taking each n x n matrix T to a real number sigma(T). The strategy specifies the shift to be applied by the QR step at T. Rayleigh and Wilkinson's are examples of simple shift strategies. We show that if sigma is continuous then there exist initial conditions for which deflation does not occur, i.e., subdiagonal entries do not tend to zero. In case of deflation, we consider the rate of convergence to zero of the (n, n-1) entry: for simple shift strategies this is always at least quadratic. If the function sigma is smooth in a suitable region and the spectrum of T does not include three consecutive eigenvalues in arithmetic progression then convergence is cubic. This implies cubic convergence to deflation of Wilkinson's shift for generic spectra. The study of the algorithm near deflation uses tubular coordinates, under which QR steps with shifts are given by a simple formula.
We define cut-and-paste, a construction which, given a quadriculated disk obtains a disjoint union of quadriculated disks of smaller total area. We provide two examples of the use of this procedure as a recursive step. Tilings of a disk Δ receive a parity: we construct a perfect or near-perfect matching of tilings of opposite parities. Let BΔ be the black-to-white adjacency matrix: we factor BΔ=LD~U, where L and U are lower and upper triangular matrices, D~ is obtained from a larger identity matrix by removing rows and columns and all entries of L, D~ and U are equal to 0, 1 or -1.
One of the most widely used methods for eigenvalue computation is the QR iteration with Wilkinson's shift: here the shift s is the eigenvalue of the bottom 2×2 principal minor closest to the corner entry. It has been a long-standing conjecture that the rate of convergence of the algorithm is cubic. In contrast, we show that there exist matrices for which the rate of convergence is strictly quadratic. More precisely, let TX be the 3×3 matrix having only two nonzero entries (TX)12=(TX)21=1 and let TL be the set of real, symmetric tridiagonal matrices with the same spectrum as TX. There exists a neighborhood U⊂TL of TX which is invariant under Wilkinson's shift strategy with the following properties. For T0∈U, the sequence of iterates (Tk) exhibits either strictly quadratic or strictly cubic convergence to zero of the entry (Tk)23. In fact, quadratic convergence occurs exactly when limTk=TX. Let X be the union of such quadratically convergent sequences (Tk): the set X has Hausdorff dimension 1 and is a union of disjoint arcs Xσ meeting at TX, where σ ranges over a Cantor set.
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.