947 research outputs found
Calculation of electron density of periodic systems using non-orthogonal localised orbitals
Methods for calculating an electron density of a periodic crystal constructed
using non-orthogonal localised orbitals are discussed. We demonstrate that an
existing method based on the matrix expansion of the inverse of the overlap
matrix into a power series can only be used when the orbitals are highly
localised (e.g. ionic systems). In other cases including covalent crystals or
those with an intermediate type of chemical bonding this method may be either
numerically inefficient or fail altogether. Instead, we suggest an exact and
numerically efficient method which can be used for orbitals of practically
arbitrary localisation. Theory is illustrated by numerical calculations on a
model system.Comment: 12 pages, 4 figure
Non-equilibrium statistical mechanics of classical nuclei interacting with the quantum electron gas
Kinetic equations governing time evolution of positions and momenta of atoms
in extended systems are derived using quantum-classical ensembles within the
Non-Equilibrium Statistical Operator Method (NESOM). Ions are treated
classically, while their electrons quantum mechanically; however, the
statistical operator is not factorised in any way and no simplifying
assumptions are made concerning the electronic subsystem. Using this method, we
derive kinetic equations of motion for the classical degrees of freedom (atoms)
which account fully for the interaction and energy exchange with the quantum
variables (electrons). Our equations, alongside the usual Newtonian-like terms
normally associated with the Ehrenfest dynamics, contain additional terms,
proportional to the atoms velocities, which can be associated with the
electronic friction. Possible ways of calculating the friction forces which are
shown to be given via complicated non-equilibrium correlation functions, are
discussed. In particular, we demonstrate that the correlation functions are
directly related to the thermodynamic Matsubara Green's functions, and this
relationship allows for the diagrammatic methods to be used in treating
electron-electron interaction perturbatively when calculating the correlation
functions. This work also generalises previous attempts, mostly based on model
systems, of introducing the electronic friction into Molecular Dynamics
equations of atoms.Comment: 18 page
Fr\'echet frames, general definition and expansions
We define an {\it -frame} with Banach spaces , , and a -space (\Theta, \snorm[\cdot]).
Then by the use of decreasing sequences of Banach spaces
and of sequence spaces , we define a general Fr\'
echet frame on the Fr\' echet space . We give
frame expansions of elements of and its dual , as well of some of
the generating spaces of with convergence in appropriate norms. Moreover,
we give necessary and sufficient conditions for a general pre-Fr\' echet frame
to be a general Fr\' echet frame, as well as for the complementedness of the
range of the analysis operator .Comment: A new section is added and a minor revision is don
Alternating Minimal Energy Methods for Linear Systems in Higher Dimensions
We propose algorithms for the solution of high-dimensional symmetrical positive definite (SPD) linear systems with the matrix and the right-hand side given and the solution sought in a low-rank format. Similarly to density matrix renormalization group (DMRG) algorithms, our methods optimize the components of the tensor product format subsequently. To improve the convergence, we expand the search space by an inexact gradient direction. We prove the geometrical convergence and estimate the convergence rate of the proposed methods utilizing the analysis of the steepest descent algorithm. The complexity of the presented algorithms is linear in the mode size and dimension, and the demonstrated convergence is comparable to or even better than the one of the DMRG algorithm. In the numerical experiment we show that the proposed methods are also efficient for non-SPD systems, for example, those arising from the chemical master equation describing the gene regulatory model at the mesoscopic scale
Approximate analytical description of the nonaffine response of amorphous solids
An approximation scheme for model disordered solids is proposed that leads to
the fully analytical evaluation of the elastic constants under explicit account
of the inhomogeneity (nonaffinity) of the atomic displacements. The theory is
in quantitative agreement with simulations for central-force systems and
predicts the vanishing of the shear modulus at the isostatic point with the
linear law {\mu} ~ (z - 2d), where z is the coordination number. The vanishing
of rigidity at the isostatic point is shown to be a consequence of the
canceling out of positive affine and negative nonaffine terms
Fredholm factorization of Wiener-Hopf scalar and matrix kernels
A general theory to factorize the Wiener-Hopf (W-H) kernel using Fredholm Integral Equations (FIE) of the second kind is presented. This technique, hereafter called Fredholm factorization, factorizes the W-H kernel using simple numerical quadrature. W-H kernels can be either of scalar form or of matrix form with arbitrary dimensions. The kernel spectrum can be continuous (with branch points), discrete (with poles), or mixed (with branch points and poles). In order to validate the proposed method, rational matrix kernels in particular are studied since they admit exact closed form factorization. In the appendix a new analytical method to factorize rational matrix kernels is also described. The Fredholm factorization is discussed in detail, supplying several numerical tests. Physical aspects are also illustrated in the framework of scattering problems: in particular, diffraction problems. Mathematical proofs are reported in the pape
Numerical solution and spectrum of boundary-domain integral equation for the Neumann BVP with variable coefficient
This is the post-print version of the Article. The official published version can be accessed from the link below - Copyright @ 2012 Taylor & Francis.In this paper, a numerical implementation of a direct united boundary-domain integral equation (BDIE) related to the Neumann boundary value problem for a scalar elliptic partial differential equation with a variable coefficient is discussed. The BDIE is reduced to a uniquely solvable one by adding an appropriate perturbation operator. The mesh-based discretization of the BDIEs with quadrilateral domain elements leads to a system of linear algebraic equations (discretized BDIE). Then, the system is solved by LU decomposition and Neumann iterations. Convergence of the iterative method is discussed in relation to the distribution of eigenvalues of the corresponding discrete operators calculated numerically.The work was supported by the grant EP/H020497/1 "Mathematical analysis of localised boundary-domain integral equations for BVPs with variable coefficients" of the EPSRC, UK
Curved Noncommutative Tori as Leibniz Quantum Compact Metric Spaces
We prove that curved noncommutative tori, introduced by Dabrowski and Sitarz,
are Leibniz quantum compact metric spaces and that they form a continuous
family over the group of invertible matrices with entries in the commutant of
the quantum tori in the regular representation, when this group is endowed with
a natural length function.Comment: 16 Pages, v3: accepted in Journal of Math. Physic
Suspensions of supracolloidal magnetic polymers: self-assembly properties from computer simulations
We study self-assembly in suspensions of supracolloidal polymer-like
structures made of crosslinked magnetic particles. Inspired by self-assembly
motifs observed for dipolar hard spheres, we focus on four different topologies
of the polymer-like structures: linear chains, rings, Y-shaped and X-shaped
polymers. We show how the presence of the crosslinkers, the number of beads in
the polymer and the magnetic interparticle interaction affect the structure of
the suspension. It turns out that for the same set of parameters, the rings are
the least active in assembling larger structures, whereas the system of Y- and
especially X-like magnetic polymers tend to form very large loose aggregates
Local Volatility Calibration by Optimal Transport
The calibration of volatility models from observable option prices is a
fundamental problem in quantitative finance. The most common approach among
industry practitioners is based on the celebrated Dupire's formula [6], which
requires the knowledge of vanilla option prices for a continuum of strikes and
maturities that can only be obtained via some form of price interpolation. In
this paper, we propose a new local volatility calibration technique using the
theory of optimal transport. We formulate a time continuous martingale optimal
transport problem, which seeks a martingale diffusion process that matches the
known densities of an asset price at two different dates, while minimizing a
chosen cost function. Inspired by the seminal work of Benamou and Brenier [1],
we formulate the problem as a convex optimization problem, derive its dual
formulation, and solve it numerically via an augmented Lagrangian method and
the alternative direction method of multipliers (ADMM) algorithm. The solution
effectively reconstructs the dynamic of the asset price between the two dates
by recovering the optimal local volatility function, without requiring any time
interpolation of the option prices
- …