Hamiltonian-preserving schemes for the Liouville equation with discontinuous potentials

Speaker: 

Prof. Shi Jin

Institution: 

University of Wisconsin-Madison

Time: 

Monday, November 14, 2005 - 4:00pm

Location: 

MSTB 254

When numerically solving the Liouville equation with a discontinuous potential, one faces the problem of severe time step restriction,and the inconsistency to the constant Hamiltonian which is related to the problem of how the weak solution should be defined for such linear hyperbolic equations with singular coefficients. In this talk,we present a class of Hamiltonian-preserving schemes that are able to overcome these numerical deficiencies. The key idea is to build into the numerical flux the behavior of a classical particle at a potential barrier. We establish the stability theory of these new schemes, and analyze their numerical accuracy. Numerical experiments are carried out to verify the theoretical results.

This method can also be applied to the level
set methods for the computations of multivalued physical observables in the semiclassical limit of the linear Schrodinger equation with a discontinuous potential. For wave equations with
discontinuous local speeds, this leads to numerical schemes consistent
with Snell's Law of Refraction.

Transition Pathways in Complex Systems: Reaction Coordinates, Isocommittor Surfaces, and Transition

Speaker: 

Prof. Eric Vanden-Eijnden

Institution: 

Courant Institute

Time: 

Monday, November 21, 2005 - 4:00pm

Location: 

MSTB 254

I will present a statistical theory of reactive trajectories between two pre-specified sets A and B, i.e. the portions of the path of a
Markov process during which the path makes a transition from A to B. This
problem is relevant e.g. in the context of metastability, in which case
the two sets A and B are metastable sets, though the formalism we propose is independent of any such assumptions on A and B. Various probability distributions on the reactive trajectories can be expressed in terms of the equilibrium distribution of the process and the so-called committor functions which give the probability that the process reaches first B before reaching A, either backward or forward in time. Using these objects, the following can be obtained: (i) the distribution of reactive trajectories, which gives the proportion of time reactive trajectories spend in sets outside of A and B; (ii) the hitting point distribution of the reactive trajectories the integral of which on a surface gives the net average flux
of reactive trajectories across this surface; and (iv) the average frequency of reactive trajectories, which gives the average number of transitions between A and B per unit of time. Time permitting, I will also illustrated how this theoretical framework can be used to design a numerical algorithm, termed the string method, to identify the committor function in high
dimensional systems.

Mathematical modeling of gene structures for prokaryotic genomes

Speaker: 

Prof. Zhen-Su She

Institution: 

UCLA

Time: 

Monday, November 28, 2005 - 4:00pm

Location: 

MSTB 254

The computational prediction for prokaryotic genes have been carried out either with very simple assumption of long open reading frame (ORF) or with Markov models with many parameters determined through machine learning algorithms. The latter usually gives highly reliable prediction of genes, but it shed little light to the structure of genes under study because of an enormous amount of parameters. We attempt to develop a mathematical model that explicitly accounts for universal features related to genes and its translation mechanisms. Our model contains relatively fewer parameters that have clear biological meaning. We design also an algorithm MED 2.0 that enables an unsupervised learning process for genome-specific parameters before the prediction of genes. The MED 2.0 not only predicts a set of genes for any newly sequenced prokaryotic genome, but also yields related parameters characterizing gene starts. We report the performance of the algorithm and discuss how mathematical models may help to gain biological insight.

On the $L^p$ Spectrum of the Hodge Laplacian on Non-Compact Manifolds

Speaker: 

Professor Nelia Charalambous

Institution: 

UCI

Time: 

Tuesday, September 27, 2005 - 4:00pm

Location: 

MSTB 254

One of the central questions in Geometric Analysis is the
interplay between the curvature of the manifold and the spectrum
of an operator.

In this talk, we will be considering the Hodge Laplacian on
differential forms of any order $k$ in the Banach Space $L^p$. In
particular, under sufficient curvature conditions, it will be
demonstrated that the $L^p\,$ spectrum is independent of $p$ for
$1\!\leq\!p\!\leq\! \infty.$ The underlying space is a
$C^{\infty}$-smooth non-compact manifold $M^n$ with a lower bound
on its Ricci Curvature and the Weitzenb\"ock Tensor. The further
assumption on subexponential growth of the manifold is also
necessary. We will see that in the case of Hyperbolic space the
$L^p$ spectrum does in fact depend on $p.$

As an application, we will show that the spectrum of the Laplacian
on one-forms has no gaps on certain manifolds with a pole and on
manifolds that are in a warped product form. This will be done
under weaker curvature restrictions than what have been used
previously; it will be achieved by finding the $L^1$ spectrum of
the Laplacian.

Time permitting, we will take a short look at an alternative
method for finding the Gaussian Heat kernel bounds for the Hodge
Laplacian via Logarithmic Sobolev Inequalities. Such bounds are
necessary in the proof of the $L^p$ independence.

Essential spectrum of the linearized Euler equations of an ideal fluid

Speaker: 

Prof. Yuri Latushkin

Institution: 

University of Missouri-Columbia

Time: 

Monday, December 5, 2005 - 4:00pm

Location: 

MSTB 254

We will discuss recent results on the
description of the essential spectrum of the
linearized Euler operator using Lyapunov exponents
and the Multiplicative Ergodic Theorem. Connections
to the evolution semigroups and to bicharacteristic
amplitude equations will be also mentioned.

Global convergence of the Yamabe flow in dimension 6 and higher

Speaker: 

Professor Simon Brendle

Institution: 

Stanford University

Time: 

Tuesday, October 18, 2005 - 3:00pm

Location: 

AP&M 7421 (UCSD)

Let $M$ be a compact manifold of dimension $n \geq 3$. Along the > Yamabe flow, a Riemannian metric on $M$ is deformed according to the > equation $\frac{\partial g}{\partial t} = -(R_g - r_g) \, g$, where $R_g$ > is the scalar curvature associated with the metric $g$ and $r_g$ denotes > the mean value of $R_g$. > > It is known that the Yamabe flow exists for all time. Moreover, if $3 \leq > n \leq 5$ or $M$ is locally conformally flat, then the solution approaches > a metric of constant scalar curvature as $t \to \infty$. I will describe > how this result can be generalized to dimensions $6$ and higher under a > technical condition on the Weyl tensor. The proof requires the > construction of a suitable family of test functions.

Inverse Scattering Theory for Jacobi Operators with Quasi-Periodic Background

Speaker: 

gerald Teschl

Institution: 

University of Vienna

Time: 

Thursday, August 25, 2005 - 2:00pm

Location: 

MSTB 254

I will consider inverse scattering theory for Jacobi operators which
are
short range perturbation of a quasi-periodic finite-gap background
operator.
In particular I want to investigate the algebraic constraints on the
scattering
data in this situation.

The Resonance Counting Function for Schrodinger Operators

Speaker: 

Professor Peter Hislop

Institution: 

University of Kentucky

Time: 

Tuesday, August 16, 2005 - 2:00pm

Location: 

MSTB 254

I will discuss the resonance counting function for Schrodinger operators with compactly-supported, $L^\infty$, real-, or complex-valued potentials, in odd dimensions $d \geq 3$. In joint work with T. Christiansen, we prove that the set of such potentials for which the resonance counting function has maximal order of growth $d$ is generic.

Pages

Subscribe to UCI Mathematics RSS