Orthogonal Polynomials on the Unit Circle. Spectral transformations and their applications to integrable systems
Series: Analysis Seminar
In this talk we will present some recent results about the matrix representation of the multiplication operator in terms of a basis of either orthogonal polynomials (OPUC) or orthogonal Laurent polynomials (OLPUC) with respect to a nontrivial probability measure supported on the unit circle. These are the so called GGT and CMV matrices.When spectral linear transformations of the measure are introduced, we will find the GGT and CMV matrices associated with the new sequences of OPUC and OLPUC, respectively. A connection with the QR factorization of such matrices will be stated. A conjecture about the generator system of such spectral transformations will be discussed.Finally, the Lax pair for the GGT and CMV matrices associated with some special time-depending deformations of the measure will be analyzed. In particular, we will study the Schur flow, which is characterized by a complex semidiscrete modified KdV equation and where a discrete analogue of the Miura transformation appears. Some open problems for time-depending deformations related to spectral linear transformations will be stated.This is a joint work with K. Castillo (Universidad Carlos III de Madrid) and L. Garza (Universidad Autonoma de Tamaulipas, Mexico).
Wednesday, February 3, 2010 - 11:00 , Location: Skiles 255 , Ilya Nemenman , Emory University , Organizer:
Even the simplest biochemical networks often have more degrees of freedoms than one can (or should!) analyze. Can we ever hope to do the physicists' favorite trick of coarse-graining, simplifying the networks to a much smaller set of effective dynamical variables that still capture the relevant aspects of the kinetics? I will argue then that methods of statistical physics provide hints at the existence of rigorous coarse-grained methodologies in modeling biological information processing systems, allowing to identify features of the systems that are relevant to their functions. While a general solution is still far away, I will focus on a specific example illustrating the approach. Namely, for a a general stochastic network exhibiting the kinetic proofreading behavior, I will show that the microscopic parameters of the system are largely important only to the extent that they contribute to a single aggregate parameter, the mean first passage time through the network, and the higher cumulants of the escape time distribution are related to this parameter uniquely. Thus a phenomenological model with a single parameter does a good job explaining all of the observable data generated by this complex system.
Series: PDE Seminar
Couette flows are shear flows with a linear velocity profile. Known by Orr in 1907, the vertical velocity of the linearized Euler equations at Couette flows is known to decay in time, for L^2 vorticity. It is interesting to know if the perturbed Euler flow near Couette tends to a nearby shear flow. Such problems of nonlinear inviscid damping also appear for other stable flows and are important to understand the appearance of coherent structures in 2D turbulence. With Chongchun Zeng, we constructed non-parallel steady flows arbitrarily near Couette flows in H^s (s<3/2) norm of vorticity. Therefore, the nonlinear inviscid damping is not true in (vorticity) H^s (s<3/2) norm. We also showed that in (vorticity) H^s (s>3/2) neighborhood of Couette flows, the only steady structures (including travelling waves) are stable shear flows. This suggests that the long time dynamics near Couette flows in (vorticity) H^s (s>3/2) space might be simpler. Similar results will also be discussed for the problem of nonlinear Landau damping in 1D electrostatic plasmas.
Series: Research Horizons Seminar
Hosted by: Huy Huynh and Yao Li
I will discuss some theorems and conjectures in the relatively new field of arithmetic dynamics, focusing in particular on some methods from number theory which can be used to study the orbits of points in algebraic dynamical systems.
Series: Analysis Working Seminar
Monday, February 1, 2010 - 13:00 , Location: Skiles 255 , Manu O. Platt , Biomedical Engineering (BME), Georgia Tech , Organizer:
Tissue remodeling involves the activation of proteases, enzymes capable of degrading the structural proteins of tissue and organs. The implications of the activation of these enzymes span all organ systems and therefore, many different disease pathologies, including cancer metastasis. This occurs when local proteolysis of the structural extracellular matrix allows for malignant cells to break free from the primary tumor and spread to other tissues. Mathematical models add value to this experimental system by explaining phenomena difficult to test at the wet lab bench and to make sense of complex interactions among the proteases or the intracellular signaling changes leading to their expression. The papain family of cysteine proteases, the cathepsins, is an understudied class of powerful collagenases and elastases implicated in extracellular matrix degradation that are secreted by macrophages and cancer cells and shown to be active in the slightly acidic tumor microenvironment. Due to the tight regulatory mechanisms of cathepsin activity and their instability outside of those defined spaces, detection of the active enzyme is difficult to precisely quantify, and therefore challenging to target therapeutically. Using valid assumptions that consider these complex interactions we are developing and validating a system of ordinary differential equations to calculate the concentrations of mature, active cathepsins in biological spaces. The system of reactions considers four enzymes (cathepsins B, K, L, and S, the most studied cathepsins with reaction rates available), three substrates (collagen IV, collagen I, and elastin) and one inhibitor (cystatin C) and comprise more than 30 differential equations with over 50 specified rate constants. Along with the mathematical model development, we have been developing new ways to quantify proteolytic activity to provide further inputs. This predictive model will be a useful tool in identifying the time scale and culprits of proteolytic breakdown leading to cancer metastasis and angiogenesis in malignant tumors.
Series: CDSNS Colloquium
We improved a computational model of leukemia development from stem cells to terminally differentiated cells by replacing the probabilistic, agent-based model of Roeder et al. (2006) with a system of deterministic, difference equations. The model is based on the relatively recent theory that cancer originates from cancer stem cells that reside in a microenvironment, called the stem cell niche. Depending on a stem cell’s location within the stem cell niche, the stem cell may remain quiescent or begin proliferating. This emerging theory states that leukemia (and potentially other cancers) is caused by the misregulation of the cycle ofproliferation and quiescence within the stem cell niche.Unlike the original agent-based model, which required seven hours per simulation, our model could be numerically evaluated in less than five minutes. The results of our numerical simulations showed that our model closely replicated the average behavior of the original agent-based model. We then extended our difference equation model to a system of age-structured partial differential equations (PDEs), which also reproduced the behavior of the Roeder model. Furthermore, the PDE model was amenable to mathematical stability analysis, which revealed three modes of behavior: stability at 0 (cancer dies out), stability at a nonzero equilibrium (a scenario akin to chronic myelogenous leukemia), and periodic oscillations (a scenario akin to accelerated myelogenous leukemia).The PDE formulation not only makes the model suitable for analysis, but also provides an effective mathematical framework for extending the model to include other aspects, such as the spatial distribution of stem cells within the niche.
Series: Combinatorics Seminar
One of the biggest hurdles in high performance computing today is the analysis of massive quantities of data. As the size of the datasets grows to petascale (and beyond), new techniques are needed to efficiently compute meaningful information from the raw data. Graph-based data (which is ubiquitous in social networks, biological interaction networks, etc) poses additional challenges due to the difficulty of parallelizing many common graph algorithms. A key component in success is the generation of "realistic" random data sets for testing and benchmarking new algorithms. The R-MAT graph generator introduced by Chakrabarti, Faloutsos, and Zhan (2004) offers a simple, fast method for generating very large directed graphs. One commonly held belief regarding graphs produced by R-MAT is that they are "scale free"; in other words, their degree distribution follows a power law as is observed in many real world networks. These properties have made R-MAT a popular choice for generating graphs for use in a variety of research disciplines including graph theoretic benchmarks, social network analysis, computational biology, and network monitoring. However, despite its wide usage and elegant, parsimonius design, our recent work provides the first rigorous mathematical analysis of the degree distributions of the generated graphs. Applying results from occupancy problems in probability theory, we derive exact expressions for the degree distributions and other parameters. We also prove that in the limit (as the number of vertices tends to infinity), graphs generated with R-MAT have degree distributions that can be expressed as a mixture of normal distributions. This talk will focus on the techniques used in solving this applied problem in terms of classical "ball and urn" results, including a minor extension of Chistyakov's theorem.
Friday, January 29, 2010 - 14:00 , Location: Skiles 269 , Mohammad Ghomi , School of Mathematics, Georgia Tech , Organizer:
We prove that convex hypersurfaces M in R^n which are level sets of functions f: R^n --> R are C^1-regular if f has a nonzero partial derivative of some order at each point of M. Furthermore, applying this result, we show that if f is algebraic and M is homeomorphic to R^(n-1), then M is an entire graph, i.e., there exists a line L in R^n such that M intersects every line parallel L at precisely one point. Finally we will give a number of examples to show that these results are sharp.