532 research outputs found
Tensor Numerical Methods in Quantum Chemistry: from Hartree-Fock Energy to Excited States
We resume the recent successes of the grid-based tensor numerical methods and
discuss their prospects in real-space electronic structure calculations. These
methods, based on the low-rank representation of the multidimensional functions
and integral operators, led to entirely grid-based tensor-structured 3D
Hartree-Fock eigenvalue solver. It benefits from tensor calculation of the core
Hamiltonian and two-electron integrals (TEI) in complexity using
the rank-structured approximation of basis functions, electron densities and
convolution integral operators all represented on 3D
Cartesian grids. The algorithm for calculating TEI tensor in a form of the
Cholesky decomposition is based on multiple factorizations using algebraic 1D
``density fitting`` scheme. The basis functions are not restricted to separable
Gaussians, since the analytical integration is substituted by high-precision
tensor-structured numerical quadratures. The tensor approaches to
post-Hartree-Fock calculations for the MP2 energy correction and for the
Bethe-Salpeter excited states, based on using low-rank factorizations and the
reduced basis method, were recently introduced. Another direction is related to
the recent attempts to develop a tensor-based Hartree-Fock numerical scheme for
finite lattice-structured systems, where one of the numerical challenges is the
summation of electrostatic potentials of a large number of nuclei. The 3D
grid-based tensor method for calculation of a potential sum on a lattice manifests the linear in computational work, ,
instead of the usual scaling by the Ewald-type approaches
Tensor-Structured Coupled Cluster Theory
We derive and implement a new way of solving coupled cluster equations with
lower computational scaling. Our method is based on decomposition of both
amplitudes and two electron integrals, using a combination of tensor
hypercontraction and canonical polyadic decomposition. While the original
theory scales as with respect to the number of basis functions, we
demonstrate numerically that we achieve sub-millihartree difference from the
original theory with scaling. This is accomplished by solving directly
for the factors that decompose the cluster operator. The proposed scheme is
quite general and can be easily extended to other many-body methods
Block Circulant and Toeplitz Structures in the Linearized Hartree–Fock Equation on Finite Lattices: Tensor Approach
This paper introduces and analyses the new grid-based tensor approach to
approximate solution of the elliptic eigenvalue problem for the 3D
lattice-structured systems. We consider the linearized Hartree-Fock equation
over a spatial lattice for both periodic and
non-periodic problem setting, discretized in the localized Gaussian-type
orbitals basis. In the periodic case, the Galerkin system matrix obeys a
three-level block-circulant structure that allows the FFT-based
diagonalization, while for the finite extended systems in a box (Dirichlet
boundary conditions) we arrive at the perturbed block-Toeplitz representation
providing fast matrix-vector multiplication and low storage size. The proposed
grid-based tensor techniques manifest the twofold benefits: (a) the entries of
the Fock matrix are computed by 1D operations using low-rank tensors
represented on a 3D grid, (b) in the periodic case the low-rank tensor
structure in the diagonal blocks of the Fock matrix in the Fourier space
reduces the conventional 3D FFT to the product of 1D FFTs. Lattice type systems
in a box with Dirichlet boundary conditions are treated numerically by our
previous tensor solver for single molecules, which makes possible calculations
on rather large lattices due to reduced numerical
cost for 3D problems. The numerical simulations for both box-type and periodic
lattice chain in a 3D rectangular "tube" with up to
several hundred confirm the theoretical complexity bounds for the
block-structured eigenvalue solvers in the limit of large .Comment: 30 pages, 12 figures. arXiv admin note: substantial text overlap with
arXiv:1408.383
A literature survey of low-rank tensor approximation techniques
During the last years, low-rank tensor approximation has been established as
a new tool in scientific computing to address large-scale linear and
multilinear algebra problems, which would be intractable by classical
techniques. This survey attempts to give a literature overview of current
developments in this area, with an emphasis on function-related tensors
Magnus and Dyson Series for Master Integrals
We elaborate on the method of differential equations for evaluating Feynman
integrals. We focus on systems of equations for master integrals having a
linear dependence on the dimensional parameter. For these systems we identify
the criteria to bring them in a canonical form, recently identified by Henn,
where the dependence of the dimensional parameter is disentangled from the
kinematics. The determination of the transformation and the computation of the
solution are obtained by using Magnus and Dyson series expansion. We apply the
method to planar and non-planar two-loop QED vertex diagrams for massive
fermions, and to non-planar two-loop integrals contributing to 2 -> 2
scattering of massless particles. The extension to systems which are polynomial
in the dimensional parameter is discussed as well.Comment: 32 pages, 6 figures, 2 ancillary files. v2: references added, typos
corrected in the text and in the ancillary file
In-plane magnetoelectric response in bilayer graphene
A graphene bilayer shows an unusual magnetoelectric response whose magnitude
is controlled by the valley-isospin density, making it possible to link
magnetoelectric behavior to valleytronics. Complementary to previous studies,
we consider the effect of static homogeneous electric and magnetic fields that
are oriented parallel to the bilayer's plane. Starting from a tight-binding
description and using quasi-degenerate perturbation theory, the low-energy
Hamiltonian is derived including all relevant magnetoelectric terms whose
prefactors are expressed in terms of tight-binding parameters. We confirm the
existence of an expected axion-type pseudoscalar term, which turns out to have
the same sign and about twice the magnitude of the previously obtained
out-of-plane counterpart. Additionally, small anisotropic corrections to the
magnetoelectric tensor are found that are fundamentally related to the skew
interlayer hopping parameter . We discuss possible ways to identify
magnetoelectric effects by distinctive features in the optical conductivity.Comment: 14 pages, 7 figure
- …