Seminarios Primavera 2017

Ejemplos de Matemáticas en Biomedicina

Viernes 19 de Mayo, 13:00 - 14:00, Salón B-3

Miguel Angel Moreles Vázquez

La charla en tres partes. Describimos de manera panorámica aplicaciones de bioimpedancia como método de diagnostico no invasivo.

Mostramos que el fenómeno subyacente se puede modelar con Cálculo Fraccionario, por lo cual hacemos un paréntesis técnico para su discusión.

Terminamos el Frankestein con una propuesta de diagnóstico de diabetes mellitus a partir de mediciones de la prueba oral de tolerancia a la glucosa.

Estructuras aritméticas de gráficas

Viernes 28 de Abril, 13:00 - 14:00, Salón B-3

Carlos Valencia
CINVESTAV

Una estructura aritmética sobre una gráfica G con n vértices es una dupla (d,r) de vectores de dimensión n con entradas enteras positivas tal que r es primitivo y r esta en el kernel de la matriz obtenida al restarle a la matriz diagonal de d la matriz de adyacencia de G. El concepto de estructura aritmética fue introducido en 1989 por Lorenzini como una matriz de intersección, las cuales aparecen en el estudio de curvas de degeneración en geometría algebraica.

Bajo ciertas condiciones se puede demostrar que el conjunto de estructuras aritméticas A(G) de una gráfica es finito. En esta platica presentaremos algunos de los conceptos básicos de las estructuras aritméticas en el contexto de M-matrices y algunos de los últimos resultados que hemos obtenido. Exploraremos como se comportan el conjunto de estructuras aritméticas de una gráfica bajo algunas operaciones de gráficas como la subdivision de aristas, la duplicación de vértices, etc. Estos resultados nos permitirán describir el conjunto de estructuras aritméticas para algunas familias de gráficas. Por ejemplo, se puede demostrar que el número de estructura aritméticas del camino es el número de Catalán.

Simetrías y coreografías en familias que bifurcan de equilibrios relativos poligonales del problema de los n cuerpos

Viernes 21 de Abril, 13:00 - 14:00, Salón B-3

Renato Calleja
IIMAS-UNAM

En esta plática describiré la continuación numérica y técnicas de bifurcación en un marco de valores de frontera utilizadas para seguir familias de Lyapunov de órbitas periódicas. Estas familias surgen de un sistema poligonal de n cuerpos en un marco de referencia rotatorio. Cuando la frecuencia de una órbita de Lyapunov y la frecuencia del marco rotatorio tiene una relación racional, la órbita es también periódica en el marco de referencia inercial. Mostraré que un conjunto denso de órbitas de Lyapunov con frecuencias que satisfacen una ecuación diofántica corresponden a coreografías. Se presentarán ejemplos de varias coreografías que encontramos numéricamente para n = 4,6,7, 8 y 9 cuerpos. Este es trabajo conjunto con Eusebius Doedel y Carlos García Azpeitia.

El triángulo hiperbólico

Viernes 5 de Mayo, 13:00 - 14:00, Salón B-3

Leobardo Fernández Roman
ITAM

A partir de la negación del quinto postulado de Euclides aparecen lo que se conocen como geometrías no Euclidianas. En esta plática presentaremos un modelo para la Geometría Hiperbólica, el modelo del semi-plano superior de Poincaré, y discutiremos algunas de las concurrencias, en este modelo, de líneas importantes en un triángulo.

La dinámica combinatoria (forcing) para funciones tipo ’skew products’ en el cilindro

Miércoles 19 de Abril, 13:00 - 14:00, Salón 104

Leopoldo Morales López
Universidad Autónoma de Barcelona

En el año 1964 Sharkovskiı̆ enunció y demostró un célebre teorema que supuso, entre otros aspectos, el inicio del estudio de lo que hoy conocemos como dinámica combinatoria en el intervalo. Este resultado afirma que la existencia de orbitas periódicas de un determinado periodo en una aplicación del intervalo “fuerza” la existencia de órbitas periódicas de otros periodos. Un refinamiento de este teorema es lo que conocemos como teorı́a del forcing de órbitas periódicas en el intervalo.

En el artı́culo Roberta Fabbri, Tobias Jäger, Russell Johnson y Gerhard Keller, A Sharkovskii-type Theorem for Minimally
Forced Interval Maps. Topol. Methods Nonlinear Anal.26 (2005),no. 1, 163-188, el Teorema de Sharkovskiı̆ fue extendido a una clase de funciones triangulares en el cilindro. Concretamente funciones continuas T : S^1 × I −→ S^1 × I donde T (θ, x) = (θ + ω, f (θ, x)) con ω ∈ R \ Q. A esta clase de funciones se las conoce en la literatura como skew-product en el cilindro.

Los objetos invariantes considerados en este caso, no son ya órbitas periódicas (ni tan solo objetos minimales) sino una generalización de curvas invariantes, que los autores llaman bandas periódicas. Intuitivamente una banda es un subconjunto compacto del cilindro tal que sus fibras en un conjunto residual de S^1. Una banda n-periódica es un conjunto de n bandas disjuntas que se aplican por la función de manera periódica entre ellas.

El trabajo que presentamos en esta charla, refina el resultado obtenido en Fabbri, Jäger et al, para obtener una teoría del forcing entre patterns de bandas periódicas. En particular demostraremos que la relación de forcing en el intervalo y en nuestra clase coinciden.

El tratamiento algunos sistemas de tipo Cauchy-Riemann inhomogéneos desde la perspectiva del análisis cuaterniónico.

Viernes 31 de Marzo, 13:00 - 14:00, Saló B-3

Marco Antonio Pérez de la Rosa
ITAM

 Se sabe que para una función de valores complejos g, definida en un dominio acotado con frontera suave, la ecuación de Cauchy-Riemann inhomogénea con dato g tiene solución para cualquier g continua y la solución general es de la forma:
f=T[g]+h, donde h es una función holomorfa en el dominio y T es el llamado operador de Theodoresco. En esta plática se mostrará una condición necesaria y suficiente para la solubilidad de algunos sistemas de tipo Cauchy-Riemann inhomogéneos y describiremos la solución general encajando el sistema en un marco cuaterniónico apropiado.

¿Qué tanto rotan tres vórtices puntuales?

Viernes 3 de Marzo, 13:00 - 14:00, Salón B-3

Antonio Hernández Garduño
Departamento de Matemáticas UAM-Iztapalapa

En su celebrado artículo “How much does the rigid body rotate?”, R. Montgomery deriva una fórmula  para la rotación de un cuerpo rígido libre durante un periodo de su dinámica reducida, interpretándola en términos de una fase de Berry.  En esta charla abordaremos el problema análogo para el problema de tres vórtices puntuales en el plano.  La rotación total tiene dos contribuciones:  la holonomía asociada a una conexión principal (fase geométrica) y un término de corrección que se calcula a partir de una versión modificada de las ecuaciones de movimiento (fase dinámica).  Discutiremos cómo una elección adecuada de la conexión permite obtener una fórmula notablemente simple para la fase dinámica, análoga a la solución de Montgomery para el cuerpo rígido.

Soluciones periódicas tipo herradura en el problema de 2n+1 cuerpos

Viernes 17 de Marzo, 13:00 - 14:00, Salón B-3

Abimael Bengochea
ITAM

Jano y Epimeteo, satélites de Saturno, siguen lo que se denomina órbita de herradura. El estudio de este tipo de órbitas se puede llevar a cabo con el modelo de 3 cuerpos. Mediante la identificación del mecanismo que da origen a dicha órbita, "construimos" las órbitas correspondientes al caso de 2n satélites con masas iguales. El estudio de estas órbitas lo realizamos con ayuda de un problema de contorno, asociado a las ecuaciones de movimiento (diferenciales ordinarias no lineales); determinamos numéricamente soluciones periódicas simétricas para n=1,2,3,4. También, mediante la continuación de las familias periódicas, mostramos que las soluciones de herradura bifurcan de soluciones homográficas (órbitas elípticas por parte de los satélites). En esta plática también se mencionará la relación entre las órbitas de herradura y la conjetura de Saari.

Quasi-periodic solutions for state dependent delay equations

Viernes 24 de Marzo, 13:00 - 14:00, Salón B-3

Rafael de la Llave
GeorgiaTech

We consider delay differential equations in which the delay depends on  the state of the system. These equations appear naturally in  electrodynamics for particles interacting with retarded potentials (the delay is proportional to the distance) as well as in several biological models.  We note that for these equations the phase space is in?finite dimensional and not easy to describe.  Questions such  as existence, uniqueness, dependence on parameters are still puzzling. We develop a theory of quasi-periodic solutions that bypasses the questions of existence for general initial data.  We develop a functional equation for the quasi-periodic equations and study them by functional analysis methods.  The main results are stated in an a-posteriori format that states that given approximate solutions that satisfy some explicit non-degeneracy solutions, there are true solutions nearby.  This can be used to justify some numerical solutions that have been produced.  We show that in a one-parameter family, there are smooth solutions.  Furthermore, we can find a large
measure set where the quasi-periodic solution are analytic.  We conjecture that this regularity is optimal.  We also develop a theory of stable/unstable manifolds.  This is joint work with Xiaolong  He.

Estabilidad y bifurcación en la configuración elíptica de Maxwell

Viernes 10 de Marzo, 13:00 - 14:00, Salón B-3

Dr. Carlos García Azpeitia
Facultad de Ciencias, UNAM.

La configuración de Maxwell consiste de n cuerpos con masas iguales (anillo de Saturno) girando alrededor de una masa central (Saturno). También se puede considerar la misma configuración de Maxwell, pero cuando los n cuerpos siguen órbitas elípticas en lugar de órbitas circulares. Usando una simetría del tipo Hip-Hop, probamos existencia de lenguas de Arnold en el conjunto de parámetros para la masa central y la excentricidad de las órbitas elípticas. Dentro de las lenguas la configuración de Maxwell es inestable. Además, la configuración elíptica tiene una bifurcación de soluciones periódicas cada vez que la excentricidad cruza una de las lenguas. Este es un trabajo en progreso, en colaboración con Abimael Bengochea y Renato Calleja.

Formación de parches bioquímicos en pelos radiculares; dinámica de estructuras localizadas.

Viernes 24 de Febrero, 13:00 - 14:00, Salón B3

Víctor Breña Medina
Centro de Ciencias Matemáticas, UNAM Campus Morelia

Una amplia diversidad de fenómenos biológicos son susceptibles de ser estudiados desde el punto de las matemáticas. Desde esta manera, la biomatemática y la biología matemática, transdisciplina la primera e interdisciplina la segunda, son áreas del conocimiento que han cobrado particular interés en la comunidad científica. En particular, esto se debe a que la biología ofrece de manera natural un amplio espectro de problemas abiertos a las matemáticas; a su vez, las teorías y métodos matemáticos proporcionan una capacidad rigurosa de análisis de hipótesis fuera del alcance experimental.

En esta charla se expondrán los ingredientes clave del análisis numérico y analítico de las propiedades dinámicas de un modelo de iniciación en una célula localizada en la epidermis de la raíz de la planta Arabidopsis thaliana. Este modelo consiste en un sistema no homogéneo de reacción-difusión. Se mostrarán algunos resultados obtenidos por medio del análisis de bifurcación numérico en una y dos dimensiones espaciales. De igual manera, se presentará un bosquejo del análisis realizado a través de la teoría de interacciones semi-fuertes y de problemas de valores propios no locales y no lineales. Los resultados obtenidos arrojan luz en el papel que juegan algunas características físicas en interacciones que se consideran de origen puramente genético.

Numerical computation of coherent structures in stochastic neural networks: from interface dynamics to coarse-grained bifurcation analysis

Viernes 10 de Febrero, 13:00 - 14:00, Salón B3

Daniele Avitabile
University of Nottingham

I will discuss the formation of coherent structures in spatially-extended neural networks posed on a ring or a torus. Phenomenological models of neural networks have been intensively studied in the past and are known to support a variety of coherent structures observed experimentally (localised bumps of activity, travelling fronts, travelling bumps, lurching waves, rotating waves). These models are typically written as integro-differential equations, where the integral term is a convolution between a synaptic kernel, specifying the neural anatomical connectivity, and a sigmoidal firing rate. Successful strategies for the analysis of these systems includespecial choices of the synaptic kernels (leading to equivalent PDE formulations) and interface methods. The latter are obtained via approximating the firing rates with Heaviside distributions. In this case, the evolution of the system is described entirely by the loci of points where the neural activity attains the firing rate threshold value.

In this talk, I will initially consider a 1D model with heterogeneous synaptic kernel and Heaviside firing rate and show that interface methods allow for the explicit construction of a bifurcation equation for localised steady states, so that analytical expressions for a classical "snakes and ladders” bifurcation scenario can be derived.

I will then consider a 2D model which does not admit an equivalent PDE formulation.  In this (and other neural field models featuring a convolution structure) it is advantageous to combine FFT and Newton-Krylov solvers to perform numerical bifurcation analysis directly on the integral model. In particular, I will show that steep sigmoidal firing rates give rise to well-conditioned linear systems, thereby eliminating the need of a preconditioner in the Krylov step, even in large-dimensional discretisations of the network.

Finally, I will discuss how the insight gained in the Heaviside limit may be used to perform coarse-grained bifurcation analysis on neural networks, even in cases where the network does not evolve according to an integro-differential equation. As an illustrative example, a heterogeneous neural network in the form of a discrete Markov chain with discrete ternary state space, posed on a lattice, will be considered. The model supports coarse bumps, multi-bumps and travelling waves, but the derivation of a coarse evolution equation is nontrivial. Hence, the emerging states are followed in parameter space, using numerical coarse bifurcation analysis. At the core of the coarse bifurcation analysis algorithm is a lifting step which provides an ensemble of
microscopic realisations compatible with a given macroscopic observable of the system. I will show that, by choosing the interface as the coarse variable, it is possible to construct efficient lifting steps for this system. Numerical results indicate that refractoriness in this network is responsible for the transition from stationary states to travelling states.

Efecto anarmónico en el solitón de Davydov en medios continuo y discreto.

7 de Abril 13:00 - 14:00, Salón B3

Luis Alberto Cisneros Ake
ESFM, IPN

Consideramos el modelo de Davydov para la transferencia de energía en un medio cristalino con interacciones no lineales del tipo
cúbico. En el límite continuo encontramos ecuaciones efectivas descritas por el acoplamiento entre las ecuaciones NLS y KdV. En dicho sistema acoplado encontramos una variedad de soluciones que van desde el solitón de Davydov de primera y segunda especies hasta los solitones obscuros (aquí llamados de primera y segunda especie, también). Finalmente, hacemos una búsqueda por medios variacionales y numéricos de éste tipo de soluciones en el medio discreto originalmente planteado.

Iso-Differential Calculus

Viernes 3 de Febrero, 13:00 - 14:00, Salón B-3

Svetlin G. Georgiev
Sorbonne University, Paris 6, France y Sofia University, Sofia, Bulgaria

In this lecture are introduced iso-numbers, iso-functions and isoderivatives. They are made comparisons with the real numbers, realvalued functions and their derivatives. They are given some applications of iso-functions.

Finite Elements on anisotropic meshes

Viernes 17 de Febrero, 13:00 - 14:00, Salón B-3

Andreas Wachtel
ITAM

Finite Element Methods (FEM) are used to approximate solutions of (partial) dierential equations. The field of applications is wide, for example, membrane equations, plate equations, fluid  flow problems, etc. However, a good quality of the approximation is in general only guaranteed for isotropic meshes. This is the most common limitation of such methods. We start this talk gently with a denition of such meshes and methods. Then, we dive into what happens when this condition is violated.

Metaheuristics for the Irregular Bin Packing Problem with free rotations

Viernes 20 de Enero, 13:00 - 14:00, Salón Por Confirmar

Antonio Martinez Sykora
University of Southampton

 We present a number of variants of a constructive algorithm able to solve a wide variety of variants of the Two-Dimensional Irregular Bin Packing Problem (2DIBPP). The aim of the 2DIBPP is to pack a set of irregular pieces,which may have concavities, into stock sheets (bins) with fixed dimensions in such a way that the utilization is maximized. This problem is inspired by a real application from a ceramic company in Spain. In addition, this problem arises in other industries such as the garment industry or ship building. The constructive procedure presented in this paper allows both free orientation for the pieces, as in the case of the ceramic industry, or a finite set of orientations as in the case of the garment industry. We explicitly model the assignment of pieces to bins and compare with the more common strategy of packing bins sequentially. There are very few papers in the literature that address the bin packing problem with irregular pieces and to our knowledge this is the first to additionally consider free rotation of pieces with bin packing. We propose several Integer Programming models to determine the association between pieces and bins and then we use a Mixed Integer Programming model for placing the pieces into the bins. The computational results show that the algorithm obtains high quality results in sets of instances with different properties. We have used both industry data and the available data in the literature of 2D irregular strip packing and bin packing problems.

Get packing! Key concepts and future directions in cutting and packing problems

Miércoles 25 de Enero, 13:00 - 14:00, Salón por Confirmar

Julia Bennell
University of Southampton

A huge number of products we use, wear and consume begin as raw material that requires cutting as part of the production process. Clothes, furniture, tools, pipes, shoes and windows are just a few examples. Moreover, the transportation of products efficiently, and safely, require an intelligent methodology for packing and loading taking into account many complex constraints. Cutting and packing problems cover a wide range of applications and with these comes many diverse and interesting challenges. Researchers have been tackling these problems since the 1960s and there is a strong legacy of methodological and application focused contributions. Despite this, there still remains many interesting open problems as well as rich opportunities for working at the interface of cutting and packing with other domains such as transportation and production planning. In this talk I will give a flavour of the diverse scope of cutting and packing problems and discuss some of the emerging application areas. The talk will review some of the key concepts and methodologies used for cutting and packing and highlight some of the current challenges. Moreover, I hope, the talk will enthuse researchers to engage with cutting and packing problems and its research community.