skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 10:00 PM ET on Friday, December 8 until 2:00 AM ET on Saturday, December 9 due to maintenance. We apologize for the inconvenience.

Title: Some matrix properties preserved by generalized matrix functions
Abstract Generalized matrix functions were first introduced in [J. B. Hawkins and A. Ben-Israel, Linear and Multilinear Algebra, 1(2), 1973, pp. 163-171]. Recently, it has been recognized that these matrix functions arise in a number of applications, and various numerical methods have been proposed for their computation. The exploitation of structural properties, when present, can lead to more efficient and accurate algorithms. The main goal of this paper is to identify structural properties of matrices which are preserved by generalized matrix functions. In cases where a given property is not preserved in general, we provide conditions on the underlying scalar function under which the property of interest will be preserved by the corresponding generalized matrix function.  more » « less
Award ID(s):
Author(s) / Creator(s):
Date Published:
Journal Name:
Special Matrices
Page Range / eLocation ID:
27 to 37
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Beginning with the work of Landau, Pollak and Slepian in the 1960s on time‐band limiting, commuting pairs of integral and differential operators have played a key role in signal processing, random matrix theory, and integrable systems. Previously, such pairs were constructed by ad hoc methods, which essentially worked because a commuting operator of low order could be found by a direct calculation. We describe a general approach to these problems that proves that every pointWof Wilson's infinite dimensional adelic Grassmannian gives rise to an integral operator , acting on for a contour , which reflects a differential operator with rational coefficients in the sense that on a dense subset of . By using analytic methods and methods from integrable systems, we show that the reflected differential operator can be constructed from the Fourier algebra of the associated bispectral function . The exact size of this algebra with respect to a bifiltration is in turn determined using algebro‐geometric methods. Intrinsic properties of four involutions of the adelic Grassmannian naturally lead us to consider the reflecting property above in place of plain commutativity. Furthermore, we prove that the time‐band limited operators of the generalized Laplace transforms with kernels given by the rank one bispectral functions always reflect a differential operator. A 90° rotation argument is used to prove that the time‐band limited operators of the generalized Fourier transforms with kernels admit a commuting differential operator. These methods produce vast collections of integral operators with prolate‐spheroidal properties, associated to the wave functions of all rational solutions of the KP hierarchy vanishing at infinity, introduced by Krichever in the late 1970s.

    more » « less
  2. Commuting integral and differential operators connect the topics of signal processing, random matrix theory, and integrable systems. Previously, the construction of such pairs was based on direct calculation and concerned concrete special cases, leaving behind important families such as the operators associated to the rational solutions of the Korteweg–de Vries (KdV) equation. We prove a general theorem that the integral operator associated to every wave function in the infinite-dimensional adelic Grassmannian G r a d of Wilson always reflects a differential operator (in the sense of Definition 1 below). This intrinsic property is shown to follow from the symmetries of Grassmannians of Kadomtsev–Petviashvili (KP) wave functions, where the direct commutativity property holds for operators associated to wave functions fixed by Wilson’s sign involution but is violated in general. Based on this result, we prove a second main theorem that the integral operators in the computation of the singular values of the truncated generalized Laplace transforms associated to all bispectral wave functions of rank 1 reflect a differential operator. A 9 0 ○ rotation argument is used to prove a third main theorem that the integral operators in the computation of the singular values of the truncated generalized Fourier transforms associated to all such KP wave functions commute with a differential operator. These methods produce vast collections of integral operators with prolate-spheroidal properties, including as special cases the integral operators associated to all rational solutions of the KdV and KP hierarchies considered by [Airault, McKean, and Moser, Commun. Pure Appl. Math. 30, 95–148 (1977)] and [Krichever, Funkcional. Anal. i Priložen. 12, 76–78 (1978)], respectively, in the late 1970s. Many examples are presented. 
    more » « less
  3. Abstract

    Additive manufacturing, no longer reserved exclusively for prototyping components, can create parts with complex geometries and locally tailored properties. For example, multiple homogenous material sources can be used in different regions of a print or be mixed during printing to define properties locally. Additionally, heterogeneous composites provide an opportunity for another level of tuning properties through processing. For example, within particulate-filled polymer matrix composites before curing, the presence of an applied electric and/or magnetic fields can reorient filler particles and form hierarchical structures depending on the fields applied. Control of particle organization is important because effective material properties are highly dependent on the distribution of filler material within composites once cured. While previous work in homogenization and effective medium theories have determined properties based upon ideal analytic distributions of particle orientations and spatial location, this work expands upon these methods generating discrete distributions from quasi-Monte Carlo simulations of the electromagnetic processing event. Results of simulations provide predicted microarchitectures from which effective properties are determined via computational homogenization.

    These particle dynamics simulations account for dielectric and magnetic forces and torques in addition to hydrodynamic forces and hard particle separation. As such, the distributions generated are processing field dependent. The effective properties for a composite represented by this distribution are determined via computational homogenization using finite element analysis (FEA). This provides a path from constituents, through processing parameters to effective material properties. In this work, we use these simulations in conjunction with a multi-objective optimization scheme to resolve the relationships between processing conditions and effective properties, to inform field-assisted additive manufacturing processes.

    The constituent set providing the largest range of properties can be found using optimization techniques applied to the aforementioned simulation framework. This key information provides a recipe for tailoring properties for additive manufacturing design and production. For example, our simulation results show that stiffness for a 10% filler volume fraction can increase by 34% when aligned by an electric field as compared to a randomly distributed composite. The stiffness of this aligned sample is also 29% higher in the direction of the alignment than perpendicular to it, which only differs by 5% from the random case [1]. Understanding this behavior and accurately predicting composite properties is key to producing field processed composites and prints. Material property predictions compare favorably to effective medium theory and experimentation with trends in elastic and magnetic effective properties demonstrating the same anisotropic behavior as a result of applied field processing. This work will address the high computational expense of physics simulation based objective functions by using efficient algorithms and data structures. We will present an optimization framework using nested gradient searches for micro barium hexaferrite particles in a PDMS matrix, optimizing on composite magnetization to determine the volume fraction of filler that will provide the largest range of properties by varying the applied electric and magnetic fields.

    more » « less
  4. Abstract

    Intestinal organoid protocols rely on the use of extracellular scaffolds, typically Matrigel, and upon switching from growth to differentiation promoting media, a symmetry breaking event takes place. During this stage, the first bud like structures analogous to crypts protrude from the central body and differentiation ensues. While organoids provide unparalleled architectural and functional complexity, this sophistication is also responsible for the high variability and lack of reproducibility of uniform crypt‐villus structures. If function follows form in organoids, such structural variability carries potential limitations for translational applications (e.g., drug screening). Consequently, there is interest in developing synthetic biomaterials to direct organoid growth and differentiation. It has been hypothesized that synthetic scaffold softening is necessary for crypt development, and these mechanical requirements raise the question, what compressive forces and subsequent relaxation are necessary for organoid maturation? To that end, allyl sulfide hydrogels are employed as a synthetic extracellular matrix mimic, but with photocleavable bonds that temporally regulate the material's bulk modulus. By varying the extent of matrix softening, it is demonstrated that crypt formation, size, and number per colony are functions of matrix softening. An understanding of the mechanical dependence of crypt architecture is necessary to instruct homogenous, reproducible organoids for clinical applications.

    more » « less
  5. Abstract

    Topology optimization has been proved to be an automatic, efficient and powerful tool for structural designs. In recent years, the focus of structural topology optimization has evolved from mono-scale, single material structural designs to hierarchical multimaterial structural designs. In this research, the multi-material structural design is carried out in a concurrent parametric level set framework so that the structural topologies in the macroscale and the corresponding material properties in mesoscale can be optimized simultaneously. The constructed cardinal basis function (CBF) is utilized to parameterize the level set function. With CBF, the upper and lower bounds of the design variables can be identified explicitly, compared with the trial and error approach when the radial basis function (RBF) is used. In the macroscale, the ‘color’ level set is employed to model the multiple material phases, where different materials are represented using combined level set functions like mixing colors from primary colors. At the end of this optimization, the optimal material properties for different constructing materials will be identified. By using those optimal values as targets, a second structural topology optimization is carried out to determine the exact mesoscale metamaterial structural layout. In both the macroscale and the mesoscale structural topology optimization, an energy functional is utilized to regularize the level set function to be a distance-regularized level set function, where the level set function is maintained as a signed distance function along the design boundary and kept flat elsewhere. The signed distance slopes can ensure a steady and accurate material property interpolation from the level set model to the physical model. The flat surfaces can make it easier for the level set function to penetrate its zero level to create new holes. After obtaining both the macroscale structural layouts and the mesoscale metamaterial layouts, the hierarchical multimaterial structure is finalized via a local-shape-preserving conformal mapping to preserve the designed material properties. Unlike the conventional conformal mapping using the Ricci flow method where only four control points are utilized, in this research, a multi-control-point conformal mapping is utilized to be more flexible and adaptive in handling complex geometries. The conformally mapped multi-material hierarchical structure models can be directly used for additive manufacturing, concluding the entire process of designing, mapping, and manufacturing.

    more » « less