Refine
Has Fulltext
- no (29) (remove)
Year of publication
- 2011 (29) (remove)
Document Type
- Article (29) (remove)
Language
- English (29) (remove)
Is part of the Bibliography
- yes (29)
Keywords
Institute
- Institut für Mathematik (29) (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.