Refine
Year of publication
- 2011 (37) (remove)
Document Type
- Article (29)
- Monograph/Edited Volume (3)
- Doctoral Thesis (3)
- Preprint (2)
Language
- English (37) (remove)
Keywords
Institute
- Institut für Mathematik (37) (remove)
We study a new approach to determine the asymptotic behaviour of quantum many-particle systems near coalescence points of particles which interact via singular Coulomb potentials. This problem is of fundamental interest in electronic structure theory in order to establish accurate and efficient models for numerical simulations. Within our approach, coalescence points of particles are treated as embedded geometric singularities in the configuration space of electrons. Based on a general singular pseudo-differential calculus, we provide a recursive scheme for the calculation of the parametrix and corresponding Green operator of a nonrelativistic Hamiltonian. In our singular calculus, the Green operator encodes all the asymptotic information of the eigenfunctions. Explicit calculations and an asymptotic representation for the Green operator of the hydrogen atom and isoelectronic ions are presented.
We deduce a new formula for the perihelion advance $Theta$ of a test particle in the Schwarzschild black hole by applying a newly developed non-linear transformation within the Schwarzschild space-time. By this transformation we are able to apply the well-known formula valid in the weak-field approximation near infinity also to trajectories in the strong-field regime near the horizon of the black hole. The resulting formula has the structure $Theta = c_1 - c_2 ln(c^2_3 - e^2) $ with positive constants $c_{1,2,3}$ depending on the angular momentum of the test particle. It is especially useful for orbits with large eccentricities $e < c_3 < 1$ showing that $Theta o infty$ as $e o c_3$.
For the Lagrangian L = G ln G where G is the Gauss-Bonnet curvature scalar we deduce the field equation and solve it in closed form for 3-flat Friedman models using a statefinder parametrization. Further we show, that among all lagrangians F(G) this L is the only one not having the form G^r with a real constant r but possessing a scale-invariant field equation. This turns out to be one of its analogies to f(R)-theories in 2-dimensional space-time. In the appendix, we systematically list several formulas for the decomposition of the Riemann tensor in arbitrary dimensions n, which are applied in the main deduction for n=4.
In this article, the fractional variational iteration method is employed for computing the approximate analytical solutions of degenerate parabolic equations with fractional time derivative. The time-fractional derivatives are described by the use of a new approach, the so-called Jumarie modified Riemann-Liouville derivative, instead in the sense of Caputo. The approximate solutions of our model problem are calculated in the form of convergent series with easily computable components. Moreover, the numerical solution is compared with the exact solution and the quantitative estimate of accuracy is obtained. The results of the study reveal that the proposed method with modified fractional Riemann-Liouville derivatives is efficient, accurate, and convenient for solving the fractional partial differential equations in multi-dimensional spaces without using any linearization, perturbation or restrictive assumptions.
The morphological features in the deviations of the total electron content (TEC) of the ionosphere from the background undisturbed state as possible precursors of the earthquake of January 12, 2010 (21:53 UT (16:53 LT), 18.46A degrees N, 72.5A degrees W, 7.0 M) in Haiti are analyzed. To identify these features, global and regional differential TEC maps based on global 2-h TEC maps provided by NASA in the IONEX format were plotted. For the considered earthquake, long-lived disturbances, presumably of seismic origin, were localized in the near-epicenter area and were accompanied by similar effects in the magnetoconjugate region. Both decreases and increases in the local TEC over the period from 22 UT of January 10 to 08 UT of January 12, 2010 were observed. The horizontal dimensions of the anomalies were similar to 40A degrees in longitude and similar to 20A degrees in latitude, with the magnitude of TEC disturbances reaching similar to 40% relative to the background near the epicenter and more than 50% in the magnetoconjugate area. No significant geomagnetic disturbances within January 1-12, 2010 were observed, i.e., the detected TEC anomalies were manifestations of interplay between processes in the lithosphere-atmosphere-ionosphere system.
The thermospheric crosswind velocities at an altitude of 400 km measured by an accelerometer on board of the CHAMP satellite are compared with the results of model calculations performed using the Upper Atmosphere Model (UAM). The results of measurements averaged over the year in 2003 reveal a two-vortex structure of high-latitude winds corresponding to magnetospheric-ionospheric convection of ions in the F2 ionosphere region. A similar picture with similar speed values was obtained in model calculations. A comparison of the crosswind speed obtained in individual measurements on October 28, 2003 with the corresponding model values revealed close agreement between them in some flights and differences in others. Taking into account the dependence of convection electric field on the B (y) component of interplanetary magnetic field sometimes improved agreement between thermospheric crosswind speeds obtained in model calculations and measured using the satellite.
We consider a resonantly perturbed system of coupled nonlinear oscillators with small dissipation and outer periodic perturbation. We show that for the large time t similar to s(-2) one component of the system is described for the most part by the inhomogeneous Mathieu equation while the other component represents pulsation of large amplitude. A Hamiltonian system is obtained which describes for the most part the behavior of the envelope in a special case. The analytic results agree with numerical simulations.
We discuss to what extent a given earthquake catalog and the assumption of a doubly truncated Gutenberg-Richter distribution for the earthquake magnitudes allow for the calculation of confidence intervals for the maximum possible magnitude M. We show that, without further assumptions such as the existence of an upper bound of M, only very limited information may be obtained. In a frequentist formulation, for each confidence level alpha the confidence interval diverges with finite probability. In a Bayesian formulation, the posterior distribution of the upper magnitude is not normalizable. We conclude that the common approach to derive confidence intervals from the variance of a point estimator fails. Technically, this problem can be overcome by introducing an upper bound (M) over tilde for the maximum magnitude. Then the Bayesian posterior distribution can be normalized, and its variance decreases with the number of observed events. However, because the posterior depends significantly on the choice of the unknown value of (M) over tilde, the resulting confidence intervals are essentially meaningless. The use of an informative prior distribution accounting for pre-knowledge of M is also of little use, because the prior is only modified in the case of the occurrence of an extreme event. Our results suggest that the maximum possible magnitude M should be better replaced by M(T), the maximum expected magnitude in a given time interval T, for which the calculation of exact confidence intervals becomes straightforward. From a physical point of view, numerical models of the earthquake process adjusted to specific fault regions may be a powerful alternative to overcome the shortcomings of purely statistical inference.
We develop a multigrid, multiple time stepping scheme to reduce computational efforts for calculating complex stress interactions in a strike-slip 2D planar fault for the simulation of seismicity. The key elements of the multilevel solver are separation of length scale, grid-coarsening, and hierarchy. In this study the complex stress interactions are split into two parts: the first with a small contribution is computed on a coarse level, and the rest for strong interactions is on a fine level. This partition leads to a significant reduction of the number of computations. The reduction of complexity is even enhanced by combining the multigrid with multiple time stepping. Computational efficiency is enhanced by a factor of 10 while retaining a reasonable accuracy, compared to the original full matrix-vortex multiplication. The accuracy of solution and computational efficiency depend on a given cut-off radius that splits multiplications into the two parts. The multigrid scheme is constructed in such a way that it conserves stress in the entire half-space.
We propose a novel strategy for global sensitivity analysis of ordinary differential equations. It is based on an error-controlled solution of the partial differential equation (PDE) that describes the evolution of the probability density function associated with the input uncertainty/variability. The density yields a more accurate estimate of the output uncertainty/variability, where not only some observables (such as mean and variance) but also structural properties (e.g., skewness, heavy tails, bi-modality) can be resolved up to a selected accuracy. For the adaptive solution of the PDE Cauchy problem we use the Rothe method with multiplicative error correction, which was originally developed for the solution of parabolic PDEs. We show that, unlike in parabolic problems, conservation properties necessitate a coupling of temporal and spatial accuracy to avoid accumulation of spatial approximation errors over time. We provide convergence conditions for the numerical scheme and suggest an implementation using approximate approximations for spatial discretization to efficiently resolve the coupling of temporal and spatial accuracy. The performance of the method is studied by means of low-dimensional case studies. The favorable properties of the spatial discretization technique suggest that this may be the starting point for an error-controlled sensitivity analysis in higher dimensions.
This note is a revised and enlarged version of the german article [16] in a slightly different framework. We here correct a serious mistake in the first version and generalize the class of Polya sum processes considered there. (A corrected version of the same results can be found already in the thesis of Mathias Rafler [12].) Moreover, the class of Polya difference processes is constructed here for the first time. In analogy to classical statistical mechanics we propose a theory of interacting Bosons and Fermions. We consider Papangelou processes. These are point processes specified by some kernel which represents the conditional intensity of the process. The main result is a general construction of a large class of such processes which contains Cox, Gibbs processes of classical statistical mechanics, but also interacting Bose and Fermi processes.
We prove a theorem on separation of boundary null points for generators of continuous semigroups of holomorphic self-mappings of the unit disk in the complex plane. Our construction demonstrates rather strikingly the particular role of the binary operation au broken vertical bar given by 1/ f au broken vertical bar g = 1/f + 1/g on generators.
We define the Dirichlet to Neumann operator for an elliptic complex of first order differential operators on a compact Riemannian manifold with boundary. Under reasonable conditions the Betti numbers of the complex prove to be completely determined by the Dirichlet to Neumann operator on the boundary.
We study the averaged macroscopic strain tensor for a sand pile consisting of soft convex polygonal particles numerically, using the discrete-element method (DEM). First, we construct two types of "sand piles" by two different pouring protocols. Afterwards, we deform the sand piles, relaxing them under a 10% reduction of gravity. Four different types of methods, three best-fit strains and a derivative strain, are adopted for determining the strain distribution under a sand pile. The results of four different versions of strains obtained from DEM simulation are compared with each other. Moreover, we compare the vertical normal strain tensor between two types of sand piles qualitatively and show how the construction history of the piles affects their strain distribution.