Spectral estimates for saddle point matrices arising in weak constraint four‐dimensional variational data assimilation

We consider the large sparse symmetric linear systems of equations that arise in the solution of weak constraint four‐dimensional variational data assimilation, a method of high interest for numerical weather prediction. These systems can be written as saddle point systems with a 3 × 3 block structure but block eliminations can be performed to reduce them to saddle point systems with a 2 × 2 block structure, or further to symmetric positive definite systems. In this article, we analyse how sensitive the spectra of these matrices are to the number of observations of the underlying dynamical system. We also obtain bounds on the eigenvalues of the matrices. Numerical experiments are used to confirm the theoretical analysis and bounds.     

Smart Geodesic Tracing in GRworkbench

We have developed a new tool for numerical work in General Relativity: GRworkbench. We discuss how GRworkbench's implementation of a numerically-amenable analogue to Differential Geometry facilitates the development of robust and chart-independent numerical algorithms. We consider, as an example, geodesic tracing on two charts covering the exterior Schwarzschild space-time.     

Bringing the balance function into focus

Charge balance functions provide important insight concerning hadronization and charge transport in heavy ion collisions at RHIC. Unfortunately, this observable is clouded by several effects. One of these effects is resonance production. Here, we present a simple thermal model where resonance production is included via a Monte Carlo technique based on the canonical ensemble. Resonance production is not found to have a large effect on the width.     

Ultrafast creation and annihilation of space-charge domains in a semiconductor superlattice observed by use of Terahertz fields

We report an experimental study indicating ultrafast creation and annihilation of space-charge domains in a semiconductor superlattice under the action of a THz field. Our experiment was performed for an InGaAs/InAlAs superlattice with the conduction electrons undergoing miniband transport. We applied to a superlattice a dc bias that was slightly smaller than a critical bias necessary for the formation of space-charge domains caused by a static negative differential conductivity. Additionally subjecting the superlattice to a strong THz field, resulted in a dc transport governed by the formation of domains if the frequency of the field was smaller than an upper frequency limit (~3 THz). From this frequency limit for the creation and annihilation of domains we determined the characteristic time of the domain buildup. Our analysis shows that the buildup time of domains in a wide miniband and heavily doped superlattice is limited by the relaxation time due to scattering of the miniband electrons at polar optic phonons. Our results are of importance for both an understanding of ultrafast dynamics of pattern formation in nanostructures and the development of THz electronic devices.Received: 25 March 2004, Published online: 23 July 2004PACS: 72.20.Ht High-field and nonlinear effects - 72.30. + q High-frequency effects; plasma effects - 73.21.Cd SuperlatticesK.N. Alekseev: Permanent address: Department of Physical Sciences, P.O. Box 3000, University of Oulu FIN-90014, Finland.