10,342 research outputs found
Lecture Notes of Tensor Network Contractions
Tensor network (TN), a young mathematical tool of high vitality and great
potential, has been undergoing extremely rapid developments in the last two
decades, gaining tremendous success in condensed matter physics, atomic
physics, quantum information science, statistical physics, and so on. In this
lecture notes, we focus on the contraction algorithms of TN as well as some of
the applications to the simulations of quantum many-body systems. Starting from
basic concepts and definitions, we first explain the relations between TN and
physical problems, including the TN representations of classical partition
functions, quantum many-body states (by matrix product state, tree TN, and
projected entangled pair state), time evolution simulations, etc. These
problems, which are challenging to solve, can be transformed to TN contraction
problems. We present then several paradigm algorithms based on the ideas of the
numerical renormalization group and/or boundary states, including density
matrix renormalization group, time-evolving block decimation,
coarse-graining/corner tensor renormalization group, and several distinguished
variational algorithms. Finally, we revisit the TN approaches from the
perspective of multi-linear algebra (also known as tensor algebra or tensor
decompositions) and quantum simulation. Despite the apparent differences in the
ideas and strategies of different TN algorithms, we aim at revealing the
underlying relations and resemblances in order to present a systematic picture
to understand the TN contraction approaches.Comment: 134 pages, 68 figures. In this version, the manuscript has been
changed into the format of book; new sections about tensor network and
quantum circuits have been adde
Combining synchrosqueezed wave packet transform with optimization for crystal image analysis
We develop a variational optimization method for crystal analysis in atomic
resolution images, which uses information from a 2D synchrosqueezed transform
(SST) as input. The synchrosqueezed transform is applied to extract initial
information from atomic crystal images: crystal defects, rotations and the
gradient of elastic deformation. The deformation gradient estimate is then
improved outside the identified defect region via a variational approach, to
obtain more robust results agreeing better with the physical constraints. The
variational model is optimized by a nonlinear projected conjugate gradient
method. Both examples of images from computer simulations and imaging
experiments are analyzed, with results demonstrating the effectiveness of the
proposed method
Bridging Proper Orthogonal Decomposition methods and augmented Newton-Krylov algorithms: an adaptive model order reduction for highly nonlinear mechanical problems
This article describes a bridge between POD-based model order reduction
techniques and the classical Newton/Krylov solvers. This bridge is used to
derive an efficient algorithm to correct, "on-the-fly", the reduced order
modelling of highly nonlinear problems undergoing strong topological changes.
Damage initiation problems are addressed and tackle via a corrected
hyperreduction method. It is shown that the relevancy of reduced order model
can be significantly improved with reasonable additional costs when using this
algorithm, even when strong topological changes are involved
Optimizing the double description method for normal surface enumeration
Many key algorithms in 3-manifold topology involve the enumeration of normal
surfaces, which is based upon the double description method for finding the
vertices of a convex polytope. Typically we are only interested in a small
subset of these vertices, thus opening the way for substantial optimization.
Here we give an account of the vertex enumeration problem as it applies to
normal surfaces, and present new optimizations that yield strong improvements
in both running time and memory consumption. The resulting algorithms are
tested using the freely available software package Regina.Comment: 27 pages, 12 figures; v2: Removed the 3^n bound from Section 3.3,
fixed the projective equation in Lemma 4.4, clarified "most triangulations"
in the introduction to section 5; v3: replace -ise with -ize for Mathematics
of Computation (note that this changes the title of the paper
PyDEC: Software and Algorithms for Discretization of Exterior Calculus
This paper describes the algorithms, features and implementation of PyDEC, a
Python library for computations related to the discretization of exterior
calculus. PyDEC facilitates inquiry into both physical problems on manifolds as
well as purely topological problems on abstract complexes. We describe
efficient algorithms for constructing the operators and objects that arise in
discrete exterior calculus, lowest order finite element exterior calculus and
in related topological problems. Our algorithms are formulated in terms of
high-level matrix operations which extend to arbitrary dimension. As a result,
our implementations map well to the facilities of numerical libraries such as
NumPy and SciPy. The availability of such libraries makes Python suitable for
prototyping numerical methods. We demonstrate how PyDEC is used to solve
physical and topological problems through several concise examples.Comment: Revised as per referee reports. Added information on scalability,
removed redundant text, emphasized the role of matrix based algorithms,
shortened length of pape
Optimized Schwarz waveform relaxation for Primitive Equations of the ocean
In this article we are interested in the derivation of efficient domain
decomposition methods for the viscous primitive equations of the ocean. We
consider the rotating 3d incompressible hydrostatic Navier-Stokes equations
with free surface. Performing an asymptotic analysis of the system with respect
to the Rossby number, we compute an approximated Dirichlet to Neumann operator
and build an optimized Schwarz waveform relaxation algorithm. We establish the
well-posedness of this algorithm and present some numerical results to
illustrate the method
Drastic Reduction of Cutoff Effects in 2-d Lattice O(N) Models
We investigate the cutoff effects in 2-d lattice O(N) models for a variety of
lattice actions, and we identify a class of very simple actions for which the
lattice artifacts are extremely small. One action agrees with the standard
action, except that it constrains neighboring spins to a maximal relative angle
delta. We fix delta by demanding that a particular value of the step scaling
function agrees with its continuum result already on a rather coarse lattice.
Remarkably, the cutoff effects of the entire step scaling function are then
reduced to the per mille level. This also applies to the theta-vacuum effects
of the step scaling function in the 2-d O(3) model. The cutoff effects of other
physical observables including the renormalized coupling and the mass in the
isotensor channel are also reduced drastically. Another choice, the mixed
action, which combines the standard quadratic with an appropriately tuned large
quartic term, also has extremely small cutoff effects. The size of cutoff
effects is also investigated analytically in 1-d and at N = infinity in 2-d.Comment: 39 pages, 18 figure
- …