Refine
Year of publication
- 2009 (40) (remove)
Document Type
- Article (27)
- Preprint (8)
- Doctoral Thesis (4)
- Monograph/Edited Volume (1)
Language
- English (40) (remove)
Keywords
Institute
- Institut für Mathematik (40) (remove)
The generalized hybrid Monte Carlo (GHMC) method combines Metropolis corrected constant energy simulations with a partial random refreshment step in the particle momenta. The standard detailed balance condition requires that momenta are negated upon rejection of a molecular dynamics proposal step. The implication is a trajectory reversal upon rejection, which is undesirable when interpreting GHMC as thermostated molecular dynamics. We show that a modified detailed balance condition can be used to implement GHMC without momentum flips. The same modification can be applied to the generalized shadow hybrid Monte Carlo (GSHMC) method. Numerical results indicate that GHMC/GSHMC implementations with momentum flip display a favorable behavior in terms of sampling efficiency, i.e., the traditional GHMC/GSHMC implementations with momentum flip got the advantage of a higher acceptance rate and faster decorrelation of Monte Carlo samples. The difference is more pronounced for GHMC. We also numerically investigate the behavior of the GHMC method as a Langevin-type thermostat. We find that the GHMC method without momentum flip interferes less with the underlying stochastic molecular dynamics in terms of autocorrelation functions and it to be preferred over the GHMC method with momentum flip. The same finding applies to GSHMC.
The generalized hybrid Monte Carlo (GHMC) method combines Metropolis corrected constant energy simulations with a partial random refreshment step in the particle momenta. The standard detailed balance condition requires that momenta are negated upon rejection of a molecular dynamics proposal step. The implication is a trajectory reversal upon rejection, which is undesirable when interpreting GHMC as thermostated molecular dynamics. We show that a modified detailed balance condition can be used to implement GHMC without momentum flips. The same modification can be applied to the generalized shadow hybrid Monte Carlo (GSHMC) method. Numerical results indicate that GHMC/GSHMC implementations with momentum flip display a favorable behavior in terms of sampling efficiency, i.e., the traditional GHMC/GSHMC implementations with momentum flip got the advantage of a higher acceptance rate and faster decorrelation of Monte Carlo samples. The difference is more pronounced for GHMC. We also numerically investigate the behavior of the GHMC method as a Langevin-type thermostat. We find that the GHMC method without momentum flip interferes less with the underlying stochastic molecular dynamics in terms of autocorrelation functions and it to be preferred over the GHMC method with momentum flip. The same finding applies to GSHMC.
The spectral theory of the Friedrichs model on the positive half line with Hilbert-Schmidt perturbations, equipped with distinguished analytic properties, is presented. In general, the (separable) multiplicity Hilbert space is assumed to be infinite-dimensional. The results include a spectral characterization of its resonances and the association of so-called Gamov vectors. Sufficient conditions are presented such that all resonances are simple poles of the scattering matrix. The connection between their residual terms and the associated Gamov vectors is pointed out.
The resonances (poles of the scattering matrix) of quantum mechanical scattering by central-symmetric potentials with compact support and zero angular momentum are spectrally characterized directly in terms of the Hamiltonian by a (generalized) eigenvalue problem distinguished by an additional condition (called boundary condition). The connection between the (generalized) eigenspace of a resonance and corresponding Gamov vectors is pointed out. A condition is presented such that a relation between special transition probabilities and infinite sums of residual terms for all complex-conjugated pairs of resonances can be proved. In the case of the square well potential the condition is satisfied.
We consider the problem of propagating an ensemble of solutions and its characterization in terms of its mean and covariance matrix. We propose differential equations that lead to a continuous matrix factorization of the ensemble into a generalized singular value decomposition (SVD). The continuous factorization is applied to ensemble propagation under periodic rescaling (ensemble breeding) and under periodic Kalman analysis steps (ensemble Kalman filter). We also use the continuous matrix factorization to perform a re-orthogonalization of the ensemble after each time-step and apply the resulting modified ensemble propagation algorithm to the ensemble Kalman filter. Results from the Lorenz-96 model indicate that the re-orthogonalization of the ensembles leads to improved filter performance.
In this paper we present a family of supersymmetric Wilson loops of N=4 supersymmetric Yang-Mills theory in Minkowski space. Our examples focus on curves restricted to hyperbolic submanifolds, H-3 and H-2, of space-time. Generically they preserve two supercharges, but in special cases more, including a case which has not been discussed before, of the hyperbolic line, conformal to the straight line and circle, which is 1/2 BPS. We discuss some general properties of these Wilson loops and their string duals and study special examples in more detail. Generically the string duals propagate on a complexification of AdS(5)xS(5) and in some specific examples the compact sphere is effectively replaced by a de Sitter space.
We introduce a new mixed finite element for solving the 2- and 3-dimensional wave equations and equations of incompressible flow. The element, which we refer to as P1(D)-P2, uses discontinuous piecewise linear functions for velocity and continuous piecewise quadratic functions for pressure. The aim of introducing the mixed formulation is to produce a new flexible element choice for triangular and tetrahedral meshes which satisfies the LBB stability condition and hence has no spurious zero-energy modes. The advantage of this particular element choice is that the mass matrix for velocity is block diagonal so it can be trivially inverted; it also allows the order of the pressure to be increased to quadratic whilst maintaining LBB stability which has benefits in geophysical applications with Coriolis forces. We give a normal mode analysis of the semi-discrete wave equation in one dimension which shows that the element pair is stable, and demonstrate that the element is stable with numerical integrations of the wave equation in two dimensions, an analysis of the resultant discrete Laplace operator in two and three dimensions on various meshes which shows that the element pair does not have any spurious modes. We provide convergence tests for the element pair which confirm that the element is stable since the convergence rate of the numerical solution is quadratic.
Using complex interpolation we prove new inclusion and coincidence theorems for multiple (fully) summing multilinear and holomorphic mappings. Among several other results we show that continuous n- linear forms on cotype 2 spaces are multiple (2; q(k),..., q(k))-summing, where 2(k-1) < n <= 2(k), q(0) = 2 and q(k+1) = 2q(k)/1+q(k) for k >= 0.
An n-ary cooperation is a mapping from a nonempty set A to the nth copower of A. A clone of cooperations is a set of cooperations which is closed under superposition and contains all injections. Coalgebras are pairs consisting of a set and a set of cooperations defined on this set. We define terms for coalgebras, coidentities and cohyperidentities. These concepts will be applied to give a new solution of the completeness problem for clones of cooperations defined on a two-element set and to separate clones of cooperations by coidentities.
The dimension of a variety V of algebras of a given type was introduced by E. Graczynska and D. Schweigert in [7] as the cardinality of the set of all derived varieties of V which are properly contained in V. In this paper, we characterize all solid varieties of dimensions 0, 1, and 2; prove that the dimension of a variety of finite type is at most N-0; give an example of a variety which has infinite dimension; and show that for every n is an element of N there is a variety with dimension n. Finally, we show that the dimension of a variety is related to the concept of the semantical kernel of a hypersubstitution and apply this connection to calculate the dimension of the class of all algebras of type tau = (n).
Ternutator identities
(2009)
The ternary commutator or ternutator, defined as the alternating sum of the product of three operators, has recently drawn much attention as an interesting structure generalizing the commutator. The ternutator satisfies cubic identities analogous to the quadratic Jacobi identity for the commutator. We present various forms of these identities and discuss the possibility of using them to define ternary algebras.
We consider an infinite system of non overlaping globules undergoing Brownian motions in R3. The term globules means that the objects we are dealing with are spherical, but with a radius which is random and time-dependent. The dynamics is modelized by an infinitedimensional Stochastic Differential Equation with local time. Existence and uniqueness of a strong solution is proven for such an equation with fixed deterministic initial condition. We also find a class of reversible measures.
We give an example of a commutative Prufer domain R with field of fractions F and a quaternion division algebra D with centre F such that R cannot be extended to a Prufer order in D in the sense of [AD]. This shows, that a general extension theorem for Prufer orders in central simple algebras does not exist and finally answers a question given in [MMU]. Moreover, in our example R is a Bezout domain which is the intersection of a countable number of (non-discrete) real valuation rings.
We construct a family of admissible analysis reconstruction pairs of wavelet families on the sphere. The construction is an extension of the isotropic Poisson wavelets. Similar to those, the directional wavelets allow a finite expansion in terms of off-center multipoles. Unlike the isotropic case, the directional wavelets are not a tight frame. However, at small scales, they almost behave like a tight frame. We give an explicit formula for the pseudodifferential operator given by the combination analysis-synthesis with respect to these wavelets. The Euclidean limit is shown to exist and an explicit formula is given. This allows us to quantify the asymptotic angular resolution of the wavelets.
Predictive population models designed to assist managers and policy makers require an explicit treatment of inherent uncertainty and variability. These are particular concerns when modelling non-native and reintroduced species, when data have been collected within one geographical or ecological context but predictions are required for another, or when extending models to predict the consequences of environmental change (e.g., climate or land-use). We present an aspatial, probabilistic framework of hierarchical process models for predicting population growth even when data are sparse or of poor quality. Insight into the factors affecting population dynamics in real landscapes can be provided and Kullback-Leibier distances are used to compare the relative output of models. This flexible yet robust framework gives easily interpretable results, allowing managers as well as modellers to invalidate anomalous models and apply others to real-life scenarios. We illustrate the framework's power with a meta-analysis of European wild boar (Sus scrofa) data. We test hypotheses about the effect of geographic region, hunting and mast years on wild boar population growth, to build models of wild boar dynamics for the UK. The framework quantifies the importance of hunting pressure as a driver of population growth, and confirms that reproductive success is greatly decreased in poor mast years, suggesting that the key to predicting wild boar dynamics is to ascertain local hunting pressure and to better understand changing food availability. Geography had no significant effect, indicating that it is not a good proxy for modelling the impact of change in climate or land-use on wild boar populations at the European scale. We use the framework to predict population abundance 9 years after an isolated population of wild boar established in the UK; in a comparison with the only field data and two independent modelling exercises, our framework provides the most robust and informative results.
We extend a classification of irreducible almost-commutative geometries, whose spectral action is dynamically nondegenerate, to internal algebras that have six simple summands. We find essentially four particle models: an extension of the standard model by a new species of fermions with vectorlike coupling to the gauge group and gauge invariant masses, two versions of the electrostrong model, and a variety of the electrostrong model with Higgs mechanism.
In this paper we present a method to recover symmetric and non-symmetric potential functions of inverse Sturm- Liouville problems from the knowledge of eigenvalues. The linear multistep method coupled with suitable boundary conditions known as boundary value method (BVM) is the main tool to approximate the eigenvalues in each iteration step of the used Newton method. The BVM was extended to work for Neumann-Neumann boundary conditions. Moreover, a suitable approximation for the asymptotic correction of the eigenvalues is given. Numerical results demonstrate that the method is able to give good results for both symmetric and non-symmetric potentials.
The Cauchy problem of the vacuum Einstein's equations aims to find a semimetric g(αβ) of a spacetime with vanishing Ricci curvature Rα,β and prescribed initial data. Under the harmonic gauge condition, the equations Rα,β = 0 are transferred into a system of quasi-linear wave equations which are called the reduced Einstein equations. The initial data for Einstein's equations are a proper Riemannian metric h(αβ) and a second fundamental form K(αβ). A necessary condition for the reduced Einstein equation to satisfy the vacuum equations is that the initial data satisfy Einstein constraint equations. Hence the data (h(αβ),K(αβ)) cannot serve as initial data for the reduced Einstein equations. Previous results in the case of asymptotically flat spacetimes provide a solution to the constraint equations in one type of Sobolev spaces, while initial data for the evolution equations belong to a different type of Sobolev spaces. The goal of our work is to resolve this incompatibility and to show that under the harmonic gauge the vacuum Einstein equations are well-posed in one type of Sobolev spaces.
For a general class of difference operators H-epsilon = T-epsilon + V-epsilon on l(2) ((epsilon Z)(d)), where V- epsilon is a multi-well potential and a is a small parameter. we analyze the asymptotic behavior as epsilon -> 0 of the (low-lying) eigenvalues and eigenfunctions. We show that the first it eigenvalues of H converge to the first it eigenvalues of the direct suns of harmonic oscillators oil R-d located at the several wells. Our proof is microlocal.
We present a Monte Carlo technique for sampling from the canonical distribution in molecular dynamics. The method is built upon the Nose-Hoover constant temperature formulation and the generalized hybrid Monte Carlo method. In contrast to standard hybrid Monte Carlo methods only the thermostat degree of freedom is stochastically resampled during a Monte Carlo step.
Studying the influence of the updating scheme for MCMC algorithm on spatially extended models is a well known problem. For discrete-time interacting particle systems we study through simulations the effectiveness of a synchronous updating scheme versus the usual sequential one. We compare the speed of convergence of the associated Markov chains from the point of view of the time-to-coalescence arising in the coupling-from-the-past algorithm. Unlike the intuition, the synchronous updating scheme is not always the best one. The distribution of the time-to-coalescence for these spatially extended models is studied too.
We construct elliptic elements in the algebra of (classical pseudo-differential) operators on a manifold M with conical singularities. The ellipticity of any such operator A refers to a pair of principal symbols (σ0, σ1) where σ0 is the standard (degenerate) homogeneous principal symbol, and σ1 is the so-called conormal symbol, depending on the complex Mellin covariable z. The conormal symbol, responsible for the conical singularity, is operator-valued and acts in Sobolev spaces on the base X of the cone. The σ1-ellipticity is a bijectivity condition for all z of real part (n + 1)/2 − γ, n = dimX, for some weight γ. In general, we have to rule out a discrete set of exceptional weights that depends on A. We show that for every operator A which is elliptic with respect to σ0, and for any real weight γ there is a smoothing Mellin operator F in the cone algebra such that A + F is elliptic including σ1. Moreover, we apply the results to ellipticity and index of (operator-valued) edge symbols from the calculus on manifolds with edges.
Multitype branching processes and Feller diffusion processes are conditioned on very late extinction. The conditioned laws are expressed as Doob h-transforms of the unconditioned laws, and an interpretation of the conditioned paths for the branching process is given, via the immortal particle. We study different limits for the conditioned process (increasing delay of extinction, long-time behavior, scaling limit) and provide an exhaustive list of exchangeability results.
This thesis considers on the one hand the construction of point processes via conditional intensities, motivated by the partial Integration of the Campbell measure of a point process. Under certain assumptions on the intensity the existence of such a point process is shown. A fundamental example turns out to be the Pólya sum process, whose conditional intensity is a generalisation of the Pólya urn dynamics. A Cox process representation for that point process is shown. A further process considered is a Poisson process of Gaussian loops, which represents a noninteracting particle system derived from the discussion of indistinguishable particles. Both processes are used to define particle systems locally, for which thermodynamic limits are determined.
Zufällige Punktprozesse beschreiben eine (zufällige) zeitliche Abfolge von Ereignissen oder eine (zufällige) räumliche Anordnung von Objekten. Deren wichtigster Vertreter ist der Poissonprozess. Der Poissonprozess zum Intensitätsmaß, das Lebesgue-Maß ordnet jedem Gebiet sein Volumen zu, erzeugt lokal, d.h in einem beschränkten Gebiet B, gerade eine mit dem Volumen von B poissonverteilte Anzahl von Punkten, die identisch und unabhängig voneinander in B plaziert werden; im Mittel ist diese Anzahl (B). Ersetzt man durch ein Vielfaches a, so wird diese Anzahl mit dem a-fachen Mittelwert erzeugt. Poissonprozesse, die im gesamten Raum unendlich viele Punkte realisieren, enthalten bereits in einer einzigen Stichprobe genügend Informationen, um Statistik betreiben zu können: Bedingt man lokal bzgl. der Anzahl der Teilchen einer Stichprobe, so fragt man nach allen Punktprozessen, die eine solche Beobachtung hätten liefern können. Diese sind Limespunktprozesse zu dieser Beobachtung. Kommt mehr als einer in Frage, spricht man von einem Phasenübergang. Da die Menge dieser Limespunktprozesse konvex ist, fragt man nach deren Extremalpunkten, dem Rand. Im ersten Teil wird ein Poissonprozess für ein physikalisches Teilchenmodell für Bosonen konstruiert. Dieses erzeugt sogenannte Loops, das sind geschlossene Polygonzüge, die dadurch charakterisiert sind, dass man an einem Ort mit einem Punkt startet, den mit einem normalverteilten Schritt läuft und dabei nach einer gegebenen, aber zufälligen Anzahl von Schritten zum Ausgangspunkt zurückkehrt. Für verschiedene Beobachtungen von Stichproben werden zugehörige Limespunktprozesse diskutiert. Diese Beobachtungen umfassen etwa das Zählen der Loops gemäaß ihrer Länge, das Zählen der Loops insgesamt, oder das Zählen der von den Loops gemachten Schritte. Jede Wahl zieht eine charakteristische Struktur der invarianten Punktprozesse nach sich. In allen hiesigen Fällen wird ein charakteristischer Phasenübergang gezeigt und Extremalpunkte werden als spezielle Poissonprozesse identifiziert. Insbesondere wird gezeigt, wie die Wahl der Beobachtung die Länge der Loops beeinflusst. Geometrische Eigenschaften dieser Poissonprozesse sind der Gegenstand des zweiten Teils der Arbeit. Die Technik der Palmschen Verteilungen eines Punktprozesses erlaubt es, unter den unendlich vielen Loops einer Realisierung den typischen Loop herauszupicken, dessen Geometrie dann untersucht wird. Eigenschaften sind unter anderem die euklidische Länge eines Schrittes oder, nimmt man mehrere aufeinander folgende Schritte, das Volumen des von ihnen definierten Simplex. Weiterhin wird gezeigt, dass der Schwerpunkt eines typischen Loops normalverteilt ist mit einer festen Varianz. Der dritte und letzte Teil befasst sich mit der Konstruktion, den Eigenschaften und der Statistik eines neuartigen Punktprozesses, der Polyascher Summenprozess genannt wird. Seine Konstruktion verallgemeinert das Prinzip der Polyaschen Urne: Im Gegensatz zum Poissonprozess, der alle Punkte unabhängig und vor allem identisch verteilt, werden hier die Punkte nacheinander derart verteilt, dass der Ort, an dem ein Punkt plaziert wird, eine Belohnung auf die Wahrscheinlichkeit bekommt, nach der nachfolgende Punkte verteilt werden. Auf diese Weise baut der Polyasche Summenprozess "Türmchen", indem sich verschiedene Punkte am selben Ort stapeln. Es wird gezeigt, dass dennoch grundlegende Eigenschaften mit denjenigen des Poissonprozesses übereinstimmen, dazu gehören unendliche Teilbarkeit sowie Unabhängigkeit der Zuwächse. Zudem werden sein Laplace-Funktional sowie seine Palmsche Verteilung bestimmt. Letztere zeigt, dass die Höhe der Türmchen gerade geometrisch verteilt ist. Abschließend werden wiederum Statistiken, nun für den Summenprozess, diskutiert. Je nach Art der Beobachtung von der Stichprobe, etwa Anzahl, Gesamthöhe der Türmchen oder beides, gibt es in jedem der drei Fälle charakteristische Limespunktprozesse und es stellt sich heraus, dass die zugehörigen Extremalverteilungen wiederum Polyasche Summenprozesse sind.
We give a survey on the calculus of (pseudo-differential) boundary value problems with the transmision property at the boundary, and ellipticity in the Shapiro-Lopatinskij sense. Apart from the original results of the work of Boutet de Monvel we present an approach based on the ideas of the edge calculus. In a final section we introduce symbols with the anti-transmission property.
The ellipticity of boundary value problems on a smooth manifold with boundary relies on a two-component principal symbolic structure (sigma(psi), sigma(partial derivative)), consisting of interior and boundary symbols. In the case of a smooth edge on manifolds with boundary, we have a third symbolic component, namely, the edge symbol sigma(boolean AND), referring to extra conditions on the edge, analogously as boundary conditions. Apart from such conditions 'in integral form' there may exist singular trace conditions, investigated in Kapanadze et al., Internal Equations and Operator Theory, 61, 241-279, 2008 on 'closed' manifolds with edge. Here, we concentrate on the phenomena in combination with boundary conditions and edge problem.
In this paper we will implement the inverse seesaw mechanism into the noncommutative framework on the basis of the AC extension of the standard model. The main difference from the classical AC model is the chiral nature of the AC fermions with respect to a U(1)(X) extension of the standard model gauge group. It is this extension which allows us to couple the right-handed neutrinos via a gauge invariant mass term to left-handed A particles. The natural scale of these gauge invariant masses is of the order of 10(17) GeV while the Dirac masses of the neutrino and the AC particles are generated dynamically and are therefore much smaller (similar to 1 to similar to 10(6) GeV). From this configuration, a working inverse seesaw mechanism for the neutrinos is obtained.
This paper provides a complete list of Krajewski diagrams representing the standard model of particle physics. We will give the possible representations of the algebra and the anomaly free lifts which provide the representation of the standard model gauge group on the fermionic Hilbert space. The algebra representations following from the Krajewski diagrams are not complete in the sense that the corresponding spectral triples do not necessarily obey to the axiom of Poincare duality. This defect may be repaired by adding new particles to the model, i.e., by building models beyond the standard model. The aim of this list of finite spectral triples (up to Poincare duality) is therefore to provide a basis for model building beyond the standard model.
In this publication we present an extension of the standard model within the framework of Connes' noncommutative geometry. The model presented here is based on a minimal spectral triple which contains the standard model particles, new vectorlike fermions, and a new U(1) gauge subgroup. Additionally a new complex scalar field appears that couples to the right-handed neutrino, the new fermions, and the standard Higgs particle. The bosonic part of the action is given by the spectral action which also determines relations among the gauge couplings, the quartic scalar couplings, and the Yukawa couplings at a cutoff energy of similar to 10(17) GeV. We investigate the renormalization group flow of these relations. The low energy behavior allows to constrain the Higgs mass, the mass of the new scalar, and the mixing between these two scalar fields.