Mathematics Colloquium
The UL Lafayette Mathematics Colloquium is an ongoing seminar series that features speakers from other universities and from our department. The topics cover all areas of mathematics and statistics. We try to schedule an interesting mix of topics ranging from very applied to more abstract in nature. These lectures are open to all UL Lafayette students, faculty and community members for the purpose of fostering continued discussion and networking in the various areas of mathematics. Please contact Paul L. Salceanu with questions or suggestions about the colloquium series.
Our colloquia are normally held on Thursday at 3:30 p.m. in room 208 of Maxim Doucet Hall. Refreshments are served at 3:15 in room 208. To accommodate outside speakers, the colloquium is occasionally held on a different day of the week, e.g., Tuesday at 3:30, instead of Thursday.
Remember, our colloquium is open to the public and everyone who is interested is encouraged to attend.
Fall 2017 Schedule

24 August 2017
Mathematical modeling of interactive wild and paratransgenic mosquitoes
Jia Li
Department of Mathematical Sciences
University of Alabama in HuntsvilleWe formulate homogeneous and stagestructured models for the interactive wild mosquitoes and paratransgenic mosquitoes carrying geneticallymodified bacteria which are resistant to malaria transmissions. We establish conditions for the existence and stability of equilibria for the model systems and provide numerical examples to demonstrate our findings. We then investigate how the bacteria uptake rate of wild mosquitoes and the bacteria loads of paratransgenic mosquitoes carrying affect their interactive dynamics.

21 September 2017
Computed Flow and Fluorescence Over the Ocular Surface
Longfei Li
Mathematics Department
UL LafayetteFluorescein is perhaps the most commonly used substance to visualize tear film thickness and dynamics; better understanding of this process aids understanding of dry eye syndrome which afflicts millions of people. We study a mathematical model for tear film flow, evaporation, solutal transport and fluorescence over the exposed ocular surface during the interblink. Transport of the fluorescein ion by fluid flow in the tear film affects the intensity of fluorescence via changes in concentration and tear film thickness. Evaporation causes increased osmolarity and potential irritation over the ocular surface; it also alters fluorescein concentration and thus fluorescence. Using thinning rates from in vivo measurements together with thin film equations for flow and transport of multiple solutes, we compute dynamic results for tear film quantities of interest. We compare our computed fluorescent intensity distributions with in vivo observations. A number of experimental features are recovered by the model.

28 September 2017
Discretetime Structured Model for Malaria Transmission with constant releasing sterile mosquitoes
Yang Li
Mathematics Department
UL LafayetteTo incorporate the interactive mosquitoes into malaria transmissions, we formulate susceptibleexposedinfectiverecovered (SEIR) compartmental discretetime models, which are of high dimensions, and then include the interactive mosquito models into these disease models. We derive formulas for the reproductive number R_0 of infection for the malaria models with or without sterile mosquitoes and explore the existence of endemic Fixed points as well. We then study the impact of sterile mosquitoes releases on the disease transmissions by investigating the effects of varying the releases of sterile mosquitoes. We use numerical simulations to verify our results for all cases and finally give brief discussions of our findings.

3 October 2017 (TUESDAY)
Z^dOdometers : a very interesting class of free minimal actions on the Cantor set
Thierry Giordano
University of OttawaZodometers form a very rigid class of minimal Cantor systems: any two orbit equivalent Zodometers are conjugate. The situation is totally different in higher dimensions. In this talk, I will review the construction of Z^dodometers and show that their first group of cohomology is a complete algebraic invariant of conjugation. Examples of orbit equivalent Z^2odometers which are not conjugate, but orbit equivalent will be presented! (This work is joint with I. Putnam and C. Skau.)

12 October 2017
Diagonalizability and the Pythagorean Theorem
Jireh Loreaux
Southern Illinois University, EdwardsvilleThere is a natural interpretation of the Pythagorean Theorem in terms of elementary linear algebra. While this interpretation is straightforward, it leads to a much more general statement in operator theory with surprising consequences. In this talk, we explore the relationship between these more general statements and the diagonalizability of matrices by means of a unitary which is a small (HilbertSchmidt) perturbation of the identity.

2 November 2017
Hyperrectangular Tolerance and Prediction Regions for Setting Multivariate Reference Regions in Laboratory Medicine
Derek S. Young
University of KentuckyReference regions are widely used in clinical chemistry and laboratory medicine to interpret the results of biochemical or physiological tests of patients. There are wellestablished methods in the literature for reference limits for univariate measurements, however, only limited methods are available for the construction of multivariate reference regions. This is because traditional multivariate statistical regions (e.g., confidence, prediction, and tolerance regions) are not constructed based on a hyperrectangular geometry. We address this problem by developing multivariate hyperrectangular nonparametric tolerance regions for setting the reference regions. Our approach utilizes statistical data depth to determine which points to trim and then the extremes of the trimmed dataset are used as the faces of the hyperrectangular region. We also specify the number of points to trim based on previouslyestablished asymptotic results. Extensive coverage results show the favorable performance of our algorithm provided a minimum sample size criterion is met. Our procedure is used to obtain reference regions for addressing two important clinical problems: (1) characterizing insulinlike growth factor concentrations in the serum of adults and (2) assessing kidney function in adolescents. This is a joint work with Thomas Mathew (UMBC).

9 November 2017
Algebraic Ktheory and polynomial maps
Saul Glasman
University of MinnesotaAssuming no prior background, I'll give a brief introduction to algebraic Ktheory, an elusive invariant of rings with profound connections to algebraic topology and number theory. Many properties of Ktheory are best understood via /categorified algebra,/ where algebraic constructions are performed at the level of categories rather than individual objects. Armed with this doctrine, and a very concrete version of the theory of polynomial functors, we'll see that such functors give rise to polynomial maps on zeroth Kgroups and thence uncover some of the arcane algebraic structure present on these groups. Insofar as time permits, I'll formally discuss lambdarings and the theory of spectral lambdarings developed in joint work in progress with C. Barwick, A. Mathew and T. Nikolaus.

16 November 2017
How to Write a Research Paper for Publication
T. Wu
Southern Illinois University Edwardsville
Edwardsville, IllinoisTo write a research paper is not an easy task. It is even more difficulty to write a research paper for publication. For a beginning or junior faculty she or he will face some promotion pressure, job security, family pressure, and social problem. In a short time period she or he must at least publish one paper in a prestige journal. To do this she or he must have good English background in reading and writing and thoroughly understanding the subject matter and wide range of related disciplines. This talk is based upon the university requirements, department environment, personal interest, and individual ability. The speaker will split entire teaching career into three stages beginning stage, mature period, and enjoy time. The speaker will provide one example for each time period respectively. (1) The speaker provides a paper that is published in the Journal of Number Theory and that is used of five pages to prove of two theorems while the original proof used thirteen pages. (2) The speaker published a paper in the ACM Transactions on Mathematical Software, the reviewer comments on the paper “this is a very good paper that I ever read in many years.” (3) The third example is a paper that is developed to present at an IEEE annual conference that is completed in one week. Hope all audiences will enjoy this talk.
Spring 2017 Schedule

26 January 2017
Nonstandard Finite Difference Methods for Dynamical Systems in Biology
Hristo Kojouharov
University of Texas at ArlingtonA brief overview of the nonstandard finite difference methods is presented. Next, using the nonstandard discretization approach, a positive and elementary stable numerical method is developed for productivedestructive systems. Finally, a nonstandard finitedifference method for general autonomous dynamical systems is constructed. The proposed numerical methods preserve the positivity of solutions and the local behavior of the corresponding dynamical systems near equilibria; and are also computationally efficient and easy to implement. Applications to select problems in biology are given to demonstrate the performance of the new methods.

31 January 2017 (TUESDAY 3:45)
Variable Selection for discrete spatial data using Penalized Quasilikelihood estimating equations
Abdhi Sarkar
Department of Statistics and Probability
Michigan State UniversityTo study real world applications of discrete data on a geographical domain we still face fundamental issues such as not being able to express the likelihood of correlated multivariate data. We circumvent this by assuming a parametric structure on the moments of a multivariate random variable and use a quasilikelihood approach. In this talk, I propose a method that is able to select relevant variables and estimate their corresponding coefficients simultaneously. Under increasing domain asymptotics after introducing a misspecified working correlation matrix that satisfies a certain mixing condition we show that this estimator possess the” oracle” property as first suggested by (Fan and Li, 2001) for the nonconvex SCAD penalty. Several simulation results and a real data example are provided to illustrate the performance of our proposed estimator.

2 February 2017
Estimation and Inference in High Dimensional Error‐in‐Variables Models and an Application to Microbiome Data
Abhishek Kaul
Biostatistics and Computational Biology
National Institute of Environmental Health Sciences
Research Triangle Park, North CarolinaWe discuss three closely related problems in high dimensional error in variables regression, 1.Additive measurement error in covariates, 2.Missing at random covariates and 3.Precision matrix recovery. We propose a two stage methodology that performs estimation post variable selection in such high dimensional measurement error models. We show that our method provides optimal rates of convergence with only a sub‐block of the bias correction matrix, while also providing a higher computational efficiency in comparison to available methods. We then apply the proposed method to human microbiome data, where we classify observations to geographical locations based on corresponding microbial compositions. Lastly, we provide methods for constructing confidence intervals on target parameters in these high dimensional models, our approach is based on the construction of moment conditions that have an additional orthogonality property with respect to nuisance parameters. All theoretical results are also supported by simulations.

3 February 2017 (FRIDAY 11:00)
Symmetric Gini Covariance and Correlation
Yongli Sang
Department of Mathematics
The University of Mississippi
University, MississippiThe most commonly used measure of dependence is the Pearson correlation. This measure is based on the covariance between two variables, which is optimal for the linear relationship between bivariate normal variables. However, the Pearson correlation performs poorly for variables with heavilytailed or asymmetric distributions, and may be seriously impacted even by a single outlier. As a robust alternative, the Spearman correlation is defined as the covariance between the cumulative distribution functions (or ranks) of two variables, but it may loss the efficiency. Complementing these two measures, the traditional Gini correlations are based on the covariance between one variable and the rank of the other, and hence well balance in efficiency and robustness. The Gini correlations, however, are not symmetric due to different roles of two variables. This asymmetry violates the axioms of correlation measurement. We have proposed a symmetric Ginitype covariance and correlation based on the joint rank function, which takes more dependence information than the marginal rank in the traditional Gini correlations. The properties of the symmetric Gini correlation are fully explored. Theoretical results on efficiency and robustness are obtained. Numerical studies demonstrate that the proposed correlation have satisfactory performance under a variety of situations. The proposed symmetric Gini correlation provides an attractive option for measuring correlation.

14 February 2017 (TUESDAY 3:30)
Degenerate Diffusion in Phase Separations
Shibin Dai
New Mexico State UniversityThe CahnHilliard equation is a widely used phenomenological diffuseinterface model for the simulations of phase separation and microstructure evolution in binary systems. We consider a popular form of the equation with a smooth doublewell potential, and with phasedependent diffusion mobilities. The latter is a feature of many materials systems and makes both theoretical analysis and accurate numerical simulations challenging. In this talk, we discuss three aspects of the equation:
1. theoretical predictions on the coarsening dynamics based on asymptotic analysis;
2. numerical simulations that confirm the theoretical predictions;
3. the existence of weak solutions that potentially accommodate the features predicted by asymptotic analysis and exhibited in numerical simulations. 
15 February 2017 (WEDNESDAY 2:30)
Special presentation: The Role of the Department Head
Seth Oppenheimer
Mississippi State UniversityThis is a special presentation by Prof Oppenheimer who is interviewing for the position of Department Head. Prof Oppenheimer will introduce himself and discuss his vision for the Department; all faculty and graduate students are encouraged to attend. The talk will last about 50 minutes. Please do your best to attend! After the talk, faculty will have an opportunity for an extended questionandanswer period with Dr Oppenheimer.

16 February 2017 (usual time)
On contraction of large perturbations of shock waves
MoonJin Kang
University of Texas at AustinAlthough mathematical understanding on hyperbolic conservation laws has made huge contributions across many fields of science, there remain many important unsolved questions. In particular, a global wellposedness of entropy solutions to the system of conservation laws in a class of large initial datas is completely open even in one space dimension. Recently, we have obtained a contraction (up to shift) of entropy shock waves to the hyperbolic systems in a class of large perturbations satisfying strong trace property. Moreover, concerning viscous systems, we have verified the contraction of large perturbations of viscous shock waves to the isentropic NavierStokes system with degenerate viscosity. Since the contraction of viscous shocks is uniformly in time and independent of viscosity coefficient, based on inviscid limit, we have the contraction (thus, uniqueness) of entropy shocks to the isentropic Euler in a class of large perturbation without any local regularity such as strong trace property. In this talk, I will present this kind of contraction property for entropy inviscid shocks and viscous shocks.

17 February 2017 (FRIDAY 11:00)
Solitary water waves
Miles Wheeler
Courant Institute of Mathematical Sciences
New York UniversityThe water wave equations describe the motion of a fluid (water) bounded above by a free surface. This free surface is subject to constant (atmospheric) pressure, while gravity acts as an external force. Traveling waves which are localized (solitary) and have small amplitude can be described by models such as the Korteweg–de Vries equation. To investigate their largeamplitude cousins, however, it is necessary to work with the full (Euler) equations. In this talk we will use continuation arguments to construct curves of largeamplitude solitary waves. We will also discuss whether the free surface of such a wave is necessarily a graph, and show that the wave speed exceeds the critical value appearing in the Kortewegde Vries approximation.

20 February 2017 (MONDAY 2:30)
Computational Modeling of Multiphase Complex Fluids with Applications
Jia Zhao
University of North Carolina at Chapel HillComplex fluids are ubiquitous in nature and in synthesized materials, such as biofilms, synthetic and biological polymeric solutions. Modeling and simulation of complex fluids has been listed as one of the 21st century mathematical challenges by DARPA, which is therefore of great mathematical and scientific significance. In this talk, I will firstly explain our research motivations by introducing several complex fluids examples, and traditional modeling techniques. Integrating the phase field approach, we then derive hydrodynamic theories for modeling multiphase complex fluid flows. Secondly, I will discuss a general technique for developing second order, linear, unconditionally energy stable numerical schemes solving hydrodynamic models. The numerical strategy is rather general that it can be applied for a host of complex fluids models. All numerical schemes developed are implemented in C2FD, a GPUbased software package developed by our group for highperformance computing/simulations. Finally, I will present several applications in cell biology, materials science and soft matter physics. 3D numerical simulations will be given. The modeling, numerical analysis and highperformance simulation tools are systematic and applicable to a large class of fluid flow problems in science and engineering.

23 February 2017
Efficient numerical schemes for the VlasovMaxwell system in plasma applications
Wei Guo
Michigan State UniversityUnderstanding complex behaviors of plasmas plays an increasingly important role in modern science and engineering. A fundamental model in plasma physics is the VlasovMaxwell system, which is a nonlinear kinetic transport model describing the dynamics of charged particles due to the selfconsistent electromagnetic forces. As predictive simulation tools in studying such a complex system, efficient, reliable and accurate transport schemes are of fundamental significance. The main numerical challenges lie in the high dimensionality, nonlinear coupling, and inherent multiscale nature of the system. In this talk, I will present several numerical methodologies to address these challenges. In this first part, I introduce a sparse grid discontinuous Galerkin (DG) method for solving the Vlasov equation, which is able to not only break the curse of dimensionality via a novel sparse approximation space, leading to remarkable computational savings, but also retain attractive properties of DG methods. In the second part, an asymptotic preserving Maxwell's solver is developed. The scheme is shown to be able to recover the correct asymptotic limit known as the Darwin limit and hence address the scale separation issue arising from plasma simulations. Theoretical and numerical results will be presented to demonstrate the efficiency and efficacy of the proposed schemes.

24 February 2017 (FRIDAY 11:00)
Bounding average quantities in dynamical systems using semidefinite programming
David Goluskin
University of MichiganI will discuss the task of proving bounds on average quantities in dissipative dynamical systems, including time averages in finitedimensional systems and spatiotemporal averages in PDE systems. In the finitedimensional case, I will describe computerassisted methods for computing bounds by constructing nonnegative polynomials with certain properties, similarly to the way nonlinear stability can be proven by constructing Lyapunov functions. Nonnegativity of these polynomials is enforced by requiring them to be representable as sums of squares, a condition that can be checked computationally using the convex optimization technique of semidefinite programming. The methods will be illustrated using the Lorenz equations, for which they produce novel bounds on various average quantities. I will then discuss work in progress on extending these computerassisted methods to dissipative PDEs, using the KuramotoSivashinky equation as an example.

2 March 2017
A multiscale model for vectorborne diseases
Hayriye Gulbudak
Arizona State UniversityThere is recent interest in mathematical models which connect the epidemiological aspects of infectious diseases to the withinhost dynamics of the pathogen and immune response. Multiscale modeling of infections allows for assessing how immunepathogen dynamics affect spread of the disease in the population. Here, I consider a withinhost model for immunepathogen dynamics nested in an agesinceinfection structured PDE system for vectorborne epidemics. First, we study pathogenhost coevolution by analytically establishing evolutionary stable strategies for parasite and host, and by utilizing computational methods to simulate the evolution in various settings. We find that vector inoculum size can contribute to virulence of vectorborne diseases in distinct ways. Next, we develop a robust methodology for identifiability and estimation of parameters with multiscale data, along with sensitivity analysis. The nested multiscale model is fit to combined withinhost and epidemiological data for Rift Valley Fever. An ultimate goal is to accurately model how control measures, such as vaccination and drug treatment, affect both scales of infection.

6 March 2017 (MONDAY 2:45)
From Picard groups to Picard categories
Michael Gurski
University of SheffieldThe Picard group of a commutative ring is a classical invariant that appears in a number of guises in algebraic geometry and number theory. This group can be enhanced to a more sophisticated invariant called a Picard category, and these appear naturally in the algebraic contexts previously mentioned as well as category theory and algebraic topology. While not a new structure, little work has been done on the algebra of Picard categories. I will introduce the basic notions, give examples and applications of Picard categories, and discuss how the algebra of Picard categories relates to that of both abelian groups and spectra via homological algebra.

8 March 2017 (WEDNESDAY 3:00)
Toric topology, polyhedral products and applications
Mentor Stafa
IndianaPurdue UniversityPolyhedral products are the central objects in the emerging field of toric topology, which stands at the crossroads of equivariant topology, algebraic and symplectic geometry, combinatorics and commutative algebra. In this talk I will give an introduction to these combinatorial constructions in topology, and give a few applications, including calculations of monodromy representations.

10 March 2017 (FRIDAY 11:15)
Configuration space integrals and integervalued cohomology classes in spaces of knots and links
Robin Koytcheff
University of MassachusettsConfiguration space integrals are a generalization of the Gauss linking integral which produce invariants of both knots and links. They can be used to construct all Vassiliev invariants, as well as nontrivial, realvalued “Vassiliev classes” in the cohomology of spaces of knots and links. I will review these ideas and then explain how configuration space integrals can be reinterpreted topologically to recover an integer lattice among the realvalued Vassiliev classes. This work also provides constructions of modp classes which need not be modp reductions of classes in this integer lattice.

28 March 2017 (TUESDAY 3:30)
Asymptotically Wellposed Boundary Conditions for Partitioned FluidStructure Algorithms
Longfei Li
Rensselaer Polytechnic InstituteA new partitioned algorithm is described for solving fluidstructure interaction (FSI) problems coupling incompressible flows with elastic structures undergoing finite deformations. The new algorithm, referred to as the AddedMass Partitioned (AMP) scheme, overcomes the addedmass instability that has for decades plagued partitioned FSI simulations of incompressible flows coupled to light structures. Within a FiniteDifference framework, the AMP scheme achieves fully secondorder accuracy and remains stable, without subtimestep iterations, even for very light structures when addedmass effects are strong. The stability and accuracy of the AMP scheme is validated through mode analysis and numerical experiments. Aiming to extend the AMP scheme to a FiniteElement framework, we also develop an accurate and efficient FiniteElement Method for solving the Incompressible NavierStokes Equations with highorder accuracy upto the boundary.

24 April 2017 (MONDAY 3:30)
Global Hopf bifurcation for differentialalgebraic equations with state dependent delay
Qingwen Hu
University of Texas at DallasWe discuss the type of differential equations with statedependent delays and the associated global Hopf bifurcation problems. In particular, we develop a global Hopf bifurcation theory for differential equations with a statedependent delay governed by an algebraic equation, using the $S^1$equivariant degree. We apply the global Hopf bifurcation theory to a model of genetic regulatory dynamics with degenerate threshold type statedependent delays, for a description of the global continuation of the periodic oscillations.