BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Gero Friesecke

Gero Friesecke appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math-phmath.MPquant-phmath.APphysics.chem-phcond-mat.str-elmath.OCmath.FAmath.NAnlin.SINumerical Analysisphysics.atom-phphysics.comp-phphysics.optics

Catalog footprint

What is connected

18works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2026arXiv

Mass splitting in the time-discrete generalized Euler equations and non-Monge solutions in multi-marginal optimal transport

The time-discretized, spatially continuous generalized Euler equations are a prototype example of multi-marginal optimal transport, yet the question whether they exhibit mass-splitting (or equivalently, whether they have solutions that are not of Monge form) has remained open. Here we resolve this question by giving a mass-splitting example in one spatial dimension. Moreover we present a related and very simple fully discrete example of mass-splitting which reveals a transparent underlying mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Density Functionals based on the mathematical structure of the strong-interaction limit of DFT

While in principle exact, Kohn-Sham density functional theory -- the workhorse of computational chemistry -- must rely on approximations for the exchange-correlation functional. Despite staggering successes, present-day approximations still struggle when the effects of electron-electron correlation play a prominent role. The limit in which the electronic Coulomb repulsion completely dominates the exchange-correlation functional offers a well-defined mathematical framework that provides insight for new approximations able to deal with strong correlation. In particular, the mathematical structure of this limit, which is now well-established thanks to its reformulation as an optimal transport problem, points to the use of very different ingredients (or features) with respect to the traditional ones used in present approximations. We focus on strategies to use these new ingredients to build approximations for computational chemistry and highlight future promising directions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Exact matrix product state representation and convergence of a fully correlated electronic wavefunction in the infinite basis limit

In this article we present the exact representation of a fully correlated electronic wavefunction as the single-particle basis approaches completeness. It consists of a half-infinite chain of matrices of exponentially increasing size. The complete basis limit is illustrated numerically using the density matrix renormalization group method by computing the core-valence entanglement in the C$_2$ ground state in increasing subsets of cc-pVTZ and pVQZ bases until convergence is reached.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the closedness and geometry of tensor network state sets

Tensor network states (TNS) are a powerful approach for the study of strongly correlated quantum matter. The curse of dimensionality is addressed by parametrizing the many-body state in terms of a network of partially contracted tensors. These tensors form a substantially reduced set of effective degrees of freedom. In practical algorithms, functionals like energy expectation values or overlaps are optimized over certain sets of TNS. Concerning algorithmic stability, it is important whether the considered sets are closed because, otherwise, the algorithms may approach a boundary point that is outside the TNS set and tensor elements diverge. We discuss the closedness and geometries of TNS sets, and we propose regularizations for optimization problems on non-closed TNS sets. We show that sets of matrix product states (MPS) with open boundary conditions, tree tensor network states (TTNS), and the multiscale entanglement renormalization ansatz (MERA) are always closed, whereas sets of translation-invariant MPS with periodic boundary conditions (PBC), heterogeneous MPS with PBC, and projected entangled-pair states (PEPS) are generally not closed. The latter is done using explicit examples like the W state, states that we call two-domain states, and fine-grained versions thereof.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Barycenters for the Hellinger--Kantorovich distance over $\mathbb{R}^d$

We study the barycenter of the Hellinger--Kantorovich metric over non-negative measures on compact, convex subsets of $\mathbb{R}^d$. The article establishes existence, uniqueness (under suitable assumptions) and equivalence between a coupled-two-marginal and a multi-marginal formulation. We analyze the HK barycenter between Dirac measures in detail, and find that it differs substantially from the Wasserstein barycenter by exhibiting a local `clustering' behaviour, depending on the length scale of the input measures. In applications it makes sense to simultaneously consider all choices of this scale, leading to a 1-parameter family of barycenters. We demonstrate the usefulness of this family by analyzing point clouds sampled from a mixture of Gaussians and inferring the number and location of the underlying Gaussians.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inversion symmetry of singular values and a new orbital ordering method in tensor train approximations for quantum chemistry

The tensor train approximation of electronic wave functions lies at the core of the QC-DMRG (Quantum Chemistry Density Matrix Renormalization Group) method, a recent state-of-the-art method for numerically solving the $N$-electron Schrödinger equation. It is well known that the accuracy of TT approximations is governed by the tail of the associated singular values, which in turn strongly depends on the ordering of the one-body basis. Here we find that the singular values $s_1\ge s_2\ge ... \ge s_d$ of tensors representing ground states of noninteracting Hamiltonians possess a surprising inversion symmetry, $s_1s_d=s_2s_{d-1}=s_3s_{d-2}=...$, thus reducing the tail behaviour to a single hidden invariant, which moreover depends explicitly on the ordering of the basis. For correlated wavefunctions, we find that the tail is upper bounded by a suitable superposition of the invariants. Optimizing the invariants or their superposition thus provides a new ordering scheme for QC-DMRG. Numerical tests on simple examples, i.e. linear combinations of a few Slater determinants, show that the new scheme reduces the tail of the singular values by several orders of magnitudes over existing methods, including the widely used Fiedler order.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sparse Control of Quantum Systems

A new class of cost functionals for optimal control of quantum systems which produces controls which are sparse in frequency and smooth in time is proposed. This is achieved by penalizing a suitable time-frequency representation of the control field, rather than the control field itself, and by employing norms which are of $L^1$ or measure form with respect to frequency but smooth with respect to time. We prove existence of optimal controls for the resulting nonsmooth optimization problem, derive necessary optimality conditions, and rigorously establish the frequency-sparsity of the optimizers. More precisely, we show that the time-frequency representation of the control field, which a priori admits a continuum of frequencies, is supported on only \textit{ finitely many} frequencies. These results cover important systems of physical interest, including (infinite-dimensional) Schrödinger dynamics on multiple potential energy surfaces as arising in laser control of chemical reactions. Numerical simulations confirm that the optimal controls, unlike those obtained with the usual $L^2$ costs, concentrate on just a few frequencies, even in the infinite-dimensional case of laser-controlled chemical reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pair densities in density functional theory

The exact interaction energy of a many-electron system is determined by the electron pair density, which is not well-approximated in standard Kohn-Sham density functional models. Here we study the (complicated but well-defined) exact universal map from density to pair density. We survey how many common functionals, including the most basic version of the LDA (Dirac exchange with no correlation contribution), arise from particular approximations of this map. We develop an algorithm to compute the map numerically, and apply it to one-parameter families {a*rho(a*x)} of one-dimensional homogeneous and inhomogeneous single-particle densities. We observe that the pair density develops remarkable multiscale patterns which strongly depend on both the particle number and the "width" 1/a of the single-particle density. The simulation results are confirmed by rigorous asymptotic results in the limiting regimes a>>1 and a<<1. For one-dimensional homogeneous systems, we show that the whole spectrum of patterns is reproduced surprisingly well by a simple asymptotics-based ansatz which slowly smoothens out the "strictly correlated" a=0 pair density while slowly turning on the a=infty "exchange" terms as a increases. Our findings lend theoretical support to the celebrated semi-empirical idea [Becke93] to mix in a fractional amount of exchange, albeit not to assuming the mixing to be additive and taking the fraction to be a system independent constant.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Twisted X-rays: incoming waveforms yielding discrete diffraction patterns for helical structures

Conventional X-ray methods use incoming plane waves and result in discrete diffraction patterns when scattered at crystals. Here we find, by a systematic method, incoming waveforms which exhibit discrete diffraction patterns when scattered at helical structures. As examples we present simulated diffraction patterns of carbon nanotubes and tobacco mosaic virus. The new incoming waveforms, which we call twisted waves due to their geometric shape, are found theoretically as closed-form solutions to Maxwell's equations. The theory of the ensuing diffraction patterns is developed in detail. A twisted analogue of the Von Laue condition is seen to hold, with the peak locations encoding the symmetry and the helix parameters, and the peak intensities indicating the electronic structure in the unit cell. If suitable twisted X-ray sources can in the future be realized experimentally, it appears from our mathematical results that they will provide a powerful tool for directly determining the detailed atomic structure of numerous biomolecules and nanostructures with helical symmetries. This would eliminate the need to crystallize those structures or their subunits.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical methods for a Kohn-Sham density functional model based on optimal transport

In this paper, we study numerical discretizations to solve density functional models in the "strictly correlated electrons" (SCE) framework. Unlike previous studies our work is not restricted to radially symmetric densities. In the SCE framework, the exchange-correlation functional encodes the effects of the strong correlation regime by minimizing the pairwise Coulomb repulsion, resulting in an optimal transport problem. We give a mathematical derivation of the self-consistent Kohn-Sham-SCE equations, construct an efficient numerical discretization for this type of problem for N = 2 electrons, and apply it to the H2 molecule in its dissociating limit. Moreover, we prove that the SCE density functional model is correct for the H2 molecule in its dissociating limit.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Infinite-body optimal transport with Coulomb Cost

We introduce and analyze symmetric infinite-body optimal transport (OT) problems with cost function of pair potential form. We show that for a natural class of such costs, the optimizer is given by the independent product measure all of whose factors are given by the one-body marginal. This is in striking contrast to standard finite-body OT problems, in which the optimizers are typically highly correlated, as well as to infinite-body OT problems with Gangbo-Swiech cost. Moreover, by adapting a construction from the study of exchangeable processes in probability theory, we prove that the corresponding $N$-body OT problem is well approximated by the infinite-body problem. To our class belongs the Coulomb cost which arises in many-electron quantum mechanics. The optimal cost of the Coulombic N-body OT problem as a function of the one-body marginal density is known in the physics and quantum chemistry literature under the name SCE functional, and arises naturally as the semiclassical limit of the celebrated Hohenberg-Kohn functional. Our results imply that in the inhomogeneous high-density limit (i.e. $N\to\infty$ with arbitrary fixed inhomogeneity profile $ρ/N$), the SCE functional converges to the mean field functional. We also present reformulations of the infinite-body and N-body OT problems as two-body OT problems with representability constraints and give a dual characterization of representable two-body measures which parallels an analogous result by Kummer on quantum representability of two-body density matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

N-density representability and the optimal transport limit of the Hohenberg-Kohn functional

We derive and analyze a hierarchy of approximations to the strongly correlated limit of the Hohenberg-Kohn functional. These "density representability approximations" are obtained by first noting that in the strongly correlated limit, N-representability of the pair density reduces to the requirement that the pair density must come from a symmetric N-point density. One then relaxes this requirement to the existence of a representing symmetric k-point density with k < N. The approximate energy can be computed by simulating a fictitious k-electron system. We investigate the approximations by deriving analytically exact results for a 2-site model problem, and by incorporating them into a self-consistent Kohn-Sham calculation for small atoms. We find that the low order representability conditions already capture the main part of the correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cnoidal Waves on Fermi-Pasta-Ulam Lattices

We study a chain of infinitely many particles coupled by nonlinear springs, obeying the equations of motion [\ddot{q}_n = V'(q_{n+1}-q_n) - V'(q_n-q_{n-1})] with generic nearest-neighbour potential $V$. We show that this chain carries exact spatially periodic travelling waves whose profile is asymptotic, in a small-amlitude long-wave regime, to the KdV cnoidal waves. The discrete waves have three interesting features: (1) being exact travelling waves they keep their shape for infinite time, rather than just up to a timescale of order wavelength$^{-3}$ suggested by formal asymptotic analysis, (2) unlike solitary waves they carry a nonzero amount of energy per particle, (3) analogous behaviour of their KdV continuum counterparts suggests long-time stability properties under nonlinear interaction with each other. Connections with the Fermi-Pasta-Ulam recurrence phenomena are indicated. Proofs involve an adaptation of the renormalization approach of Friesecke and Pego (1999) to a periodic setting and the spectral theory of the periodic Schrödinger operator with KdV cnoidal wave potential.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Density functional theory and optimal transportation with Coulomb cost

We present here novel insight into exchange-correlation functionals in density functional theory, based on the viewpoint of optimal transport. We show that in the case of two electrons and in the semiclassical limit, the exact exchange-correlation functional reduces to a very interesting functional of novel form, which depends on an optimal transport map $T$ associated with a given density $ρ$. Since the above limit is strongly correlated, the limit functional yields insight into electron correlations. We prove the existence and uniqueness of such an optimal map for any number of electrons and each $ρ$, and determine the map explicitly in the case when $ρ$ is radially symmetric.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Efficient Algorithm for Asymptotics-Based Configuration-Interaction Methods and Electronic Structure of Transition Metal Atoms

Asymptotics-based configuration-interaction (CI) methods [G. Friesecke and B. D. Goddard, Multiscale Model. Simul. 7, 1876 (2009)] are a class of CI methods for atoms which reproduce, at fixed finite subspace dimension, the exact Schrödinger eigenstates in the limit of fixed electron number and large nuclear charge. Here we develop, implement, and apply to 3d transition metal atoms an efficient and accurate algorithm for asymptotics-based CI. Efficiency gains come from exact (symbolic) decomposition of the CI space into irreducible symmetry subspaces at essentially linear computational cost in the number of radial subshells with fixed angular momentum, use of reduced density matrices in order to avoid having to store wavefunctions, and use of Slater-type orbitals (STO's). The required Coulomb integrals for STO's are evaluated in closed form, with the help of Hankel matrices, Fourier analysis, and residue calculus. Applications to 3d transition metal atoms are in good agreement with experimental data. In particular we reproduce the anomalous magnetic moment and orbital filling of Chromium in the otherwise regular series Ca, Sc, Ti, V, Cr.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Semiclassical limit of quantum dynamics with rough potentials and well posedness of transport equations with measure initial data

In this paper we study the semiclassical limit of the Schrödinger equation. Under mild regularity assumptions on the potential $U$ which include Born-Oppenheimer potential energy surfaces in molecular dynamics, we establish asymptotic validity of classical dynamics globally in space and time for "almost all" initial data, with respect to an appropriate reference measure on the space of initial data. In order to achieve this goal we prove existence, uniqueness and stability results for the flow in the space of measures induced by the continuity equation.

0 citations0 reviewsNo code linkedOpen access