On the Existence of Focal Points near Closed Geodesics on Surfaces

We establish some criteria for the existence or nonexistence of focal points near closed geodesics on surfaces. These criteria are in terms of the curvature of the manifold along the closed geodesic and the average values of the partial derivatives of the curvature in the direction perpendicular to the geodesic. Our criteria lead to a new family of examples of surfaces with no focal points. We also show that if S is a compact surface with no focal points and an inequality relating the curvature of the surface to the curvature of the horocycles holds, then the horocycles (considered as curves in S) are uniformly C ^2+Lipschitz.     

Rational approximation to trigonometric operators

We consider the approximation of trigonometric operator functions that arise in the numerical solution of wave equations by trigonometric integrators. It is well known that Krylov subspace methods for matrix functions without exponential decay show superlinear convergence behavior if the number of steps is larger than the norm of the operator. Thus, Krylov approximations may fail to converge for unbounded operators. In this paper, we propose and analyze a rational Krylov subspace method which converges not only for finite element or finite difference approximations to differential operators but even for abstract, unbounded operators. In contrast to standard Krylov methods, the convergence will be independent of the norm of the operator and thus of its spatial discretization. We will discuss efficient implementations for finite element discretizations and illustrate our analysis with numerical experiments. AMS subject classification (2000)  65F10, 65L60, 65M60, 65N22     

"Compression-only" behavior: a second-order nonlinear response of ultrasound contrast agent microbubbles

Oscillating phospholipid-coated ultrasound contrast agent microbubbles display a so-called "compression-only" behavior, where it is observed that the bubbles compress efficiently while their expansion is suppressed. Here, a theoretical understanding of the source of this nonlinear behavior is provided through a weakly nonlinear analysis of the shell buckling model proposed by Marmottant et al. [J. Acoust. Soc. Am. 118, 3499-3505 (2005)]. It is shown that the radial dynamics of the bubble can be considered as a superposition of a linear response at the fundamental driving frequency and a second-order nonlinear low-frequency response that describes the negative offset of the mean bubble radius. The analytical solution deduced from the weakly nonlinear analysis shows that the compression-only behavior results from a rapid change of the shell elasticity with bubble radius. In addition, the radial dynamics of single phospholipid-coated microbubbles was recorded as a function of both the amplitude and the frequency of the driving pressure pulse. The comparison between the experimental data and the theory shows that the magnitude of compression-only behavior is mainly determined by the initial phospholipids concentration on the bubble surface, which slightly varies from bubble to bubble.     

Exponential multistep methods of Adams-type

The paper is concerned with the construction, implementation and numerical analysis of exponential multistep methods. These methods are related to explicit Adams methods but, in contrast to the latter, make direct use of the exponential and related matrix functions of a (possibly rough) linearization of the vector field. This feature enables them to integrate stiff problems explicitly in time.     

On the use of two QMR algorithms for solving singular systems and applications in Markov chain modeling

Recently, Freund and Nachtigal proposed the quasi-minimal residual algorithm (QMR) for solving general nonsingular non-Hermitian linear systems. The method is based on the Lanczos process, and thus it involves matrix—vector products with both the coefficient matrix of the linear system and its transpose. Freund developed a variant of QMR, the transpose-free QMR algorithm (TFQMR), that only requires products with the coefficient matrix. In this paper, the use of QMR and TFQMR for solving singular systems is explored. First, a convergence result for the general class of Krylov-subspace methods applied to singular systems is presented. Then, it is shown that QMR and TFQMR both converge for consistent singular linear systems with coefficient matrices of index 1. Singular systems of this type arise in Markov chain modeling. For this particular application, numerical experiments are reported.     

Beantwortung einer Frage von J. Arason und A. Pfister als Korollar eines Ergebnisses von L. Bröcker

Ohne Zusammenfassung     

Quantum calculation of ro-vibrational states: methodology and DOCl application results

An improved Lanczos eigenvalue analysis method has been developed to compute the bound ro-vibrational states for the DOCl system at a total angular momentum of J = 0 and J = 30. In this method, the error norm is used to identify all the true eigenvalues, using the Lanczos algorithm without re-orthogonalization. For ro-vibrational spectroscopy calculations, the comparisons among experimental results, the exact quantum mechanical calculations, and the widely used approximate adiabatic rotation method have been made for J = 30. For J = 0, the density of states (DOS) in both the bound and unimolecular dissociation regime have been computed, whereas for the J = 30 case, only the DOS in the lower portion of the bound spectrum has been reported, because of substantial computational tasks.     

State-to-state reaction probabilities for the H+O2(v,j)-->O+OH(v',j') reaction on three potential energy surfaces

We report state-to-state and total reaction probabilities for J=0 and total reaction probabilities for J=2 and 4 for the title reaction, both for ground-state and initially rovibrationally excited reactants. The results for three different potential energy surfaces are compared and contrasted. The potential energy surfaces employed are the DMBE IV surface by Pastrana et al. [J. Phys. Chem. 94, 8073 (1990)], the surface by Troe and Ushakov (TU) [J. Chem. Phys. 115, 3621 (2001)], and the new XXZLG ab initio surface by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. Our results show that the total reaction probabilities from both the TU and XXZLG surfaces are much smaller in magnitude for collision energies above 1.2 eV compared to the DMBE IV surface. The three surfaces also show different behavior with regards to the effect of initial state excitation. The reactivity is increased on the XXZLG and the TU surfaces and decreased on the DMBE IV surface. Vibrational and rotational product state distributions for the XXZLG and the DMBE IV surface show different behaviors for both types of distributions. Our results show that for energies above 1.25 eV the dynamics on the DMBE IV surface are not statistical. However, there is also evidence that the dynamics on the XXZLG surface are not purely statistical for energies above the onset of the first excited product vibrational state v'=1. The magnitude of the total reaction probability is decreased for J>0 for the DMBE IV and the XXZLG surfaces for ground-state reactants. However, for initially rovibrationally excited reactants, the total reaction probability does not decrease as expected for both surfaces. As a result the total cross section averaged over all Boltzmann accessible rotational states may well be larger than the cross section reported in the literature for j=1.     

Isolation and structure elucidation of two new cytotoxic metabolites from red yeast rice

In this study, 10 already described secondary metabolites and 2 unknown metabolites were identified in an extract of Monascus purpureus by high-performance liquid chromatography-diode array detection. The unknown metabolites were isolated and their chemical structures were elucidated. The new metabolites possess the molecular formulas C(21)H(27)NO(4) and C(23)H(31)NO(4). They were named monascopyridines E and F due to their pyridine backbone. The cytotoxicity of the new compounds was studied using immortalised human kidney epithelial cells displaying IC(50) values in the micromolar range.