Refine
Has Fulltext
- no (40) (remove)
Year of publication
- 2015 (40) (remove)
Document Type
- Article (40) (remove)
Is part of the Bibliography
- yes (40)
Keywords
- Ellipticity of corner-degenerate operators (2)
- Fredholm property (2)
- Toeplitz operators (2)
- 31A25 (1)
- 65F18 (1)
- AERONET (1)
- AMS (1)
- Aerosols (1)
- Asymptotics of solutions (1)
- Averaging principle (1)
Institute
- Institut für Mathematik (40) (remove)
In this work we extract the microphysical properties of aerosols for a collection of measurement cases with low volume depolarization ratio originating from fire sources captured by the Raman lidar located at the National Institute of Optoelectronics (INOE) in Bucharest. Our algorithm was tested not only for pure smoke but also for mixed smoke and urban aerosols of variable age and growth. Applying a sensitivity analysis on initial parameter settings of our retrieval code was proved vital for producing semi-automatized retrievals with a hybrid regularization method developed at the Institute of Mathematics of Potsdam University. A direct quantitative comparison of the retrieved microphysical properties with measurements from a Compact Time of Flight Aerosol Mass Spectrometer (CToF-AMS) is used to validate our algorithm. Microphysical retrievals performed with sun photometer data are also used to explore our results. Focusing on the fine mode we observed remarkable similarities between the retrieved size distribution and the one measured by the AMS. More complicated atmospheric structures and the factor of absorption appear to depend more on particle radius being subject to variation. A good correlation was found between the aerosol effective radius and particle age, using the ratio of lidar ratios (LR: aerosol extinction to backscatter ratios) as an indicator for the latter. Finally, the dependence on relative humidity of aerosol effective radii measured on the ground and within the layers aloft show similar patterns. (C) 2015 Elsevier Inc. All rights reserved.
We present simulations of binary black-hole mergers in which, after the common outer horizon has formed, the marginally outer trapped surfaces (MOTSs) corresponding to the individual black holes continue to approach and eventually penetrate each other. This has very interesting consequences according to recent results in the theory of MOTSs. Uniqueness and stability theorems imply that two MOTSs which touch with a common outer normal must be identical. This suggests a possible dramatic consequence of the collision between a small and large black hole. If the penetration were to continue to completion, then the two MOTSs would have to coalesce, by some combination of the small one growing and the big one shrinking. Here we explore the relationship between theory and numerical simulations, in which a small black hole has halfway penetrated a large one.
We consider the volume- normalized Ricci flow close to compact shrinking Ricci solitons. We show that if a compact Ricci soliton (M, g) is a local maximum of Perelman's shrinker entropy, any normalized Ricci flowstarting close to it exists for all time and converges towards a Ricci soliton. If g is not a local maximum of the shrinker entropy, we showthat there exists a nontrivial normalized Ricci flow emerging from it. These theorems are analogues of results in the Ricci- flat and in the Einstein case (Haslhofer and Muller, arXiv:1301.3219, 2013; Kroncke, arXiv: 1312.2224, 2013).
We find necessary conditions for a second order ordinary differential equation to be equivalent to the Painleve III equation under a general point transformation. Their sufficiency is established by reduction to known results for the equations of the form y ' = f (x, y). We consider separately the generic case and the case of reducibility to an autonomous equation. The results are illustrated by the primary resonance equation.
The primary motivation for systematic bases in first principles electronic structure simulations is to derive physical and chemical properties of molecules and solids with predetermined accuracy. This requires a detailed understanding of the asymptotic behaviour of many-particle Coulomb systems near coalescence points of particles. Singular analysis provides a convenient framework to study the asymptotic behaviour of wavefunctions near these singularities. In the present work, we want to introduce the mathematical framework of singular analysis and discuss a novel asymptotic parametrix construction for Hamiltonians of many-particle Coulomb systems. This corresponds to the construction of an approximate inverse of a Hamiltonian operator with remainder given by a so-called Green operator. The Green operator encodes essential asymptotic information and we present as our main result an explicit asymptotic formula for this operator. First applications to many-particle models in quantum chemistry are presented in order to demonstrate the feasibility of our approach. The focus is on the asymptotic behaviour of ladder diagrams, which provide the dominant contribution to short-range correlation in coupled cluster theory. Furthermore, we discuss possible consequences of our asymptotic analysis with respect to adaptive wavelet approximation.
Processes having the same bridges as a given reference Markov process constitute its reciprocal class. In this paper we study the reciprocal class of compound Poisson processes whose jumps belong to a finite set . We propose a characterization of the reciprocal class as the unique set of probability measures on which a family of time and space transformations induces the same density, expressed in terms of the reciprocal invariants. The geometry of plays a crucial role in the design of the transformations, and we use tools from discrete geometry to obtain an optimal characterization. We deduce explicit conditions for two Markov jump processes to belong to the same class. Finally, we provide a natural interpretation of the invariants as short-time asymptotics for the probability that the reference process makes a cycle around its current state.
By edge algebra we understand a pseudo-differential calculus on a manifold with edge. The operators have a two-component principal symbolic hierarchy which determines operators up to lower order terms. Those belong to a filtration of the corresponding operator spaces. We give a new characterisation of this structure, based on an alternative representation of edge amplitude functions only containing holomorphic edge-degenerate Mellin symbols.
For an irreducible continuous time Markov chain, we derive the distribution of the first passage time from a given state i to another given state j and the reversed passage time from j to i, each under the condition of no return to the starting point. When these two distributions are identical, we say that i and j are in time duality. We introduce a new condition called permuted balance that generalizes the concept of reversibility and provides sufficient criteria, based on the structure of the transition graph of the Markov chain. Illustrative examples are provided.
Certain curvature conditions for the stability of Einstein manifolds with respect to the Einstein-Hilbert action are given. These conditions are given in terms of quantities involving the Weyl tensor and the Bochner tensor. In dimension six, a stability criterion involving the Euler characteristic is given.
We consider the signal detection problem in the Gaussian design trace regression model with low rank alternative hypotheses. We derive the precise (Ingster-type) detection boundary for the Frobenius and the nuclear norm. We then apply these results to show that honest confidence sets for the unknown matrix parameter that adapt to all low rank sub-models in nuclear norm do not exist. This shows that recently obtained positive results in [5] for confidence sets in low rank recovery problems are essentially optimal.
We study systematically the estimation of Earth's core angular momentum (CAM) variation between 1962.0 and 2008.0 by using core surface flow models derived from the recent geomagnetic field model C(3)FM2. Various flow models are derived by changing four parameters that control the least squares flow inversion. The parameters include the spherical harmonic (SH) truncation degree of the flow models and two Lagrange multipliers that control the weights of two additional constraints. The first constraint forces the energy spectrum of the flow solution to follow a power law l-p, where l is the SH degree and p is the fourth parameter. The second allows to modulate the solution continuously between the dynamical states of tangential geostrophy (TG) and tangential magnetostrophy (TM). The calculated CAM variations are examined in reference to two features of the observed length-of-day (LOD) variation, namely, its secular trend and 6year oscillation. We find flow models in either TG or TM state for which the estimated CAM trends agree with the LOD trend. It is necessary for TM models to have their flows dominate at planetary scales, whereas TG models should not be of this scale; otherwise, their CAM trends are too steep. These two distinct types of flow model appear to correspond to the separate regimes of previous numerical dynamos that are thought to be applicable to the Earth's core. The phase of the subdecadal CAM variation is coherently determined from flow models obtained with extensively varying inversion settings. Multiple sources of model ambiguity need to be allowed for in discussing whether these phase estimates properly represent that of Earth's CAM as an origin of the observed 6year LOD oscillation.
We study infinitesimal Einstein deformations on compact flat manifolds and on product manifolds. Moreover, we prove refinements of results by Koiso and Bourguignon which yield obstructions on the existence of infinitesimal Einstein deformations under certain curvature conditions. (C) 2014 Elsevier B.V. All rights reserved.
Performance of the generalized shadow hybrid Monte Carlo (GSHMC) method [1], which proved to be superior in sampling efficiency over its predecessors [2-4], molecular dynamics and hybrid Monte Carlo, can be further improved by combining it with multi-time-stepping (MTS) and mollification of slow forces. We demonstrate that the comparatively simple modifications of the method not only lead to better performance of GSHMC itself but also allow for beating the best performed methods, which use the similar force splitting schemes. In addition we show that the same ideas can be successfully applied to the conventional generalized hybrid Monte Carlo method (GHMC). The resulting methods, MTS-GHMC and MTS-GSHMC, provide accurate reproduction of thermodynamic and dynamical properties, exact temperature control during simulation and computational robustness and efficiency. MTS-GHMC uses a generalized momentum update to achieve weak stochastic stabilization to the molecular dynamics (MD) integrator. MTS-GSHMC adds the use of a shadow (modified) Hamiltonian to filter the MD trajectories in the HMC scheme. We introduce a new shadow Hamiltonian formulation adapted to force-splitting methods. The use of such Hamiltonians improves the acceptance rate of trajectories and has a strong impact on the sampling efficiency of the method. Both methods were implemented in the open-source MD package ProtoMol and were tested on a water and a protein systems. Results were compared to those obtained using a Langevin Molly (LM) method [5] on the same systems. The test results demonstrate the superiority of the new methods over LM in terms of stability, accuracy and sampling efficiency. This suggests that putting the MTS approach in the framework of hybrid Monte Carlo and using the natural stochasticity offered by the generalized hybrid Monte Carlo lead to improving stability of MTS and allow for achieving larger step sizes in the simulation of complex systems.
Metastability in a class of hyperbolic dynamical systems perturbed by heavy-tailed Levy type noise
(2015)
We consider a finite dimensional deterministic dynamical system with finitely many local attractors K-iota, each of which supports a unique ergodic probability measure P-iota, perturbed by a multiplicative non-Gaussian heavy-tailed Levy noise of small intensity epsilon > 0. We show that the random system exhibits a metastable behavior: there exists a unique epsilon-dependent time scale on which the system reminds of a continuous time Markov chain on the set of the invariant measures P-iota. In particular our approach covers the case of dynamical systems of Morse-Smale type, whose attractors consist of points and limit cycles, perturbed by multiplicative alpha-stable Levy noise in the Ito, Stratonovich and Marcus sense. As examples we consider alpha-stable Levy perturbations of the Duffing equation and Pareto perturbations of a biochemical birhythmic system with two nested limit cycles.
In this study, we analyze acoustic emission (AE) data recorded at the Morsleben salt mine, Germany, to assess the catalog completeness, which plays an important role in any seismicity analysis. We introduce the new concept of a magnitude completeness interval consisting of a maximum magnitude of completeness (M-c(max)) in addition to the well-known minimum magnitude of completeness. This is required to describe the completeness of the catalog, both for the smallest events (for which the detection performance may be low) and for the largest ones (which may be missed because of sensors saturation). We suggest a method to compute the maximum magnitude of completeness and calculate it for a spatial grid based on (1) the prior estimation of saturation magnitude at each sensor, (2) the correction of the detection probability function at each sensor, including a drop in the detection performance when it saturates, and (3) the combination of detection probabilities of all sensors to obtain the network detection performance. The method is tested using about 130,000 AE events recorded in a period of five weeks, with sources confined within a small depth interval, and an example of the spatial distribution of M-c(max) is derived. The comparison between the spatial distribution of M-c(max) and of the maximum possible magnitude (M-max), which is here derived using a recently introduced Bayesian approach, indicates that M-max exceeds M-c(max) in some parts of the mine. This suggests that some large and important events may be missed in the catalog, which could lead to a bias in the hazard evaluation.
We study the possibility of obtaining a computational turbulence model by means of non-dissipative regularisation of the compressible atmospheric equations for climate-type applications. We use an -regularisation (Lagrangian averaging) of the atmospheric equations. For the hydrostatic and compressible atmospheric equations discretised using a finite volume method on unstructured grids, deterministic and non-deterministic numerical experiments are conducted to compare the individual solutions and the statistics of the regularised equations to those of the original model. The impact of the regularisation parameter is investigated. Our results confirm the principal compatibility of -regularisation with atmospheric dynamics and encourage further investigations within atmospheric model including complex physical parametrisations.
We compute explicitly, and without any extra regularity assumptions, the large time limit of the fibrewise heat operator for Bismut-Lott type superconnections in the L-2-setting. This is motivated by index theory on certain non-compact spaces (families of manifolds with cocompact group action) where the convergence of the heat operator at large time implies refined L-2-index formulas.
As applications, we prove a local L-2-index theorem for families of signature operators and an L-2-Bismut-Lott theorem, expressing the Becker-Gottlieb transfer of flat bundles in terms of Kamber-Tondeur classes. With slightly stronger regularity we obtain the respective refined versions: we construct L-2-eta forms and L-2-torsion forms as transgression forms.