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

Christian Schilling

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

ResearcherUnclaimed source record
quant-phmath-phmath.MPSystems and Controlphysics.chem-phArtificial Intelligencecond-mat.quant-gasLogic in Computer Sciencecond-mat.str-eleess.SYMachine Learningmath.DSmath.NAmath.OCNumerical Analysis

Catalog footprint

What is connected

33works
15topics
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

33 published item(s)

preprint2026arXiv

A Geometric Approach to Strongly Correlated Bosons: From $N$-Representability to the Generalized BEC Force

Building on recent advances in reduced density matrix theory, we develop a geometric framework for describing strongly correlated lattice bosons. We first establish that translational symmetry, together with a fixed pair interaction, enables an exact functional formulation expressed solely in terms of momentum occupation numbers. Employing the constrained-search formalism and exploiting a geometric correspondence between $N$-boson configuration states and their one-particle reduced density matrices, we derive the general form of the ground-state functional. Its structure highlights the omnipresent significance of one-body $N$-representability: (i) the domain is exactly determined by the $N$-representability conditions; (ii) at its boundary, the gradient of the functional diverges repulsively, thereby generalizing the recently discovered Bose-Einstein condensate (BEC) force; and (iii) an explicit expression for this boundary force follows directly from geometric arguments. These key results are demonstrated analytically for few-site lattice systems, and we illustrate the broader significance of our functional form in defining a systematic hierarchy of functional approximations.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Implementations of two Algorithms for the Threshold Synthesis Problem

A linear pseudo-Boolean constraint (LPB) is an expression of the form $a_1 \cdot \ell_1 + \dots + a_m \cdot \ell_m \geq d$, where each $\ell_i$ is a literal (it assumes the value 1 or 0 depending on whether a propositional variable $x_i$ is true or false) and $a_1, \dots, a_m, d$ are natural numbers. An LPB represents a Boolean function, and those Boolean functions that can be represented by exactly one LPB are called threshold functions. The problem of finding an LPB representation of a Boolean function if possible is called threshold recognition problem or threshold synthesis problem. The problem has an $O(m^7 t^5)$ algorithm using linear programming, where $m$ is the dimension and $t$ the number of terms in the DNF input. It has been an open question whether one can recognise threshold functions through an entirely combinatorial procedure. Smaus has developed such a procedure for doing this, which works by decomposing the DNF and "counting" the variable occurrences in it. We have implemented both algorithms as a thesis project. We report here on this experience. The most important insight was that the algorithm by Smaus is, unfortunately, incomplete.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Comment on "Self-Consistent-Field Method for Correlated Many-Electron Systems with an Entropic Cumulant Energy"

In [Phys. Rev. Lett. 128, 013001 (2022)] a novel ground state method was proposed. It has been suggested that this $i$-DMFT would be a method within one-particle reduced density matrix functional theory (DMFT), capable of describing accurately molecules at various geometries with an information-theoretical nature. We reassess this work and its suggestions from a conceptual and practical point of view, leading to the following conclusions: i) A method which assigns to each molecule $\mathcal{M}$ its own functional $\mathcal{F}_{\!\mathcal{M}}$ is not a functional theory (striking violation of "universality") ii) even for the simplest systems $i$-DMFT yields incorrect one-particle reduced density matrices and iii) the use of an information-theoretical concept to describe molecular dissociation limits was not essential. The latter insight may help to fix the deficiency of $i$-DMFT to not reproduce correctly the smaller occupation numbers and thus to not recover the important dynamic correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Conservative Time Discretization: A Comparative Study

We present the first review of methods to overapproximate the set of reachable states of linear time-invariant systems subject to uncertain initial states and input signals for short time horizons. These methods are fundamental to state-of-the-art reachability algorithms for long time horizons, which proceed in two steps: First they use such a method to discretize the system for a short time horizon, and then they efficiently obtain a solution of the new discrete system for the long time horizon. Traditionally, both qualitative and quantitative comparison between different reachability algorithms has only considered the combination of both steps. In this paper we study the first step in isolation. We perform a variety of numerical experiments for six fundamental discretization methods from the literature. As we show, these methods have different trade-offs regarding accuracy and computational cost and, depending on the characteristics of the system, some methods may be preferred over others. We also discuss preprocessing steps to improve the results and efficient implementation strategies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Efficient reachability analysis of parametric linear hybrid systems with time-triggered transitions

Efficiently handling time-triggered and possibly nondeterministic switches for hybrid systems reachability is a challenging task. In this paper we present an approach based on conservative set-based enclosure of the dynamics that can handle systems with uncertain parameters and inputs, where the uncertainties are bound to given intervals. The method is evaluated on the plant model of an experimental electro-mechanical braking system with periodic controller. In this model, the fast-switching controller dynamics requires simulation time scales of the order of nanoseconds. Accurate set-based computations for relatively large time horizons are known to be expensive. However, by appropriately decoupling the time variable with respect to the spatial variables, and enclosing the uncertain parameters using interval matrix maps acting on zonotopes, we show that the computation time can be lowered to 5000 times faster with respect to previous works. This is a step forward in formal verification of hybrid systems because reduced run-times allow engineers to introduce more expressiveness in their models with a relatively inexpensive computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Excitations of Quantum Many-Body Systems via Purified Ensembles: A Unitary-Coupled-Cluster-based Approach

State-average calculations based on mixture of states are increasingly being exploited across chemistry and physics as versatile procedures for addressing excitations of quantum many-body systems. If not too many states should need to be addressed, calculations performed on individual states is also a common option. Here we show how the two approaches can be merged into one method, dealing with a generalized yet single pure state. Implications in electronic structure calculations are discussed and for quantum computations are pointed out.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Foundation of one-particle reduced density matrix functional theory for excited states

In [Phys. Rev. Lett. 127, 023001 (2021)] a reduced density matrix functional theory (RDMFT) has been proposed for calculating energies of selected eigenstates of interacting many-fermion systems. Here, we develop a solid foundation for this so-called $\boldsymbol{w}$-RDMFT and present the details of various derivations. First, we explain how a generalization of the Ritz variational principle to ensemble states with fixed weights $\boldsymbol{w}$ in combination with the constrained search would lead to a universal functional of the one-particle reduced density matrix. To turn this into a viable functional theory, however, we also need to implement an exact convex relaxation. This general procedure includes Valone's pioneering work on ground state RDMFT as the special case $\boldsymbol{w}=(1,0,\ldots)$. Then, we work out in a comprehensive manner a methodology for deriving a compact description of the functional's domain. This leads to a hierarchy of generalized exclusion principle constraints which we illustrate in great detail. By anticipating their future pivotal role in functional theories and to keep our work self-contained, several required concepts from convex analysis are introduced and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantifying Electron Entanglement Faithfully

Entanglement is one of the most fascinating concepts of modern physics. In striking contrast to its abstract, mathematical foundation, its practical side is, however, remarkably underdeveloped. Even for systems of just two orbitals or sites no faithful entanglement measure is known yet. By exploiting the spin symmetries of realistic many-electron systems, we succeed in deriving a closed formula for the relative entropy of entanglement between electron orbitals. Its broad applicability in the quantum sciences is demonstrated: (i) in light of the second quantum revolution, it quantifies the true physical entanglement by incorporating the crucial fermionic superselection rule (ii) an analytic description of the long-distance entanglement in free electron chains is found, refining Kohn's locality principle (iii) the bond-order wave phase in the extended Hubbard model can be confirmed, and (iv) the quantum complexity of common molecular bonding structures could be marginalized through orbital transformations, thus rationalizing zero-seniority wave function ansatzes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

SpecRepair: Counter-Example Guided Safety Repair of Deep Neural Networks

Deep neural networks (DNNs) are increasingly applied in safety-critical domains, such as self-driving cars, unmanned aircraft, and medical diagnosis. It is of fundamental importance to certify the safety of these DNNs, i.e. that they comply with a formal safety specification. While safety certification tools exactly answer this question, they are of no help in debugging unsafe DNNs, requiring the developer to iteratively verify and modify the DNN until safety is eventually achieved. Hence, a repair technique needs to be developed that can produce a safe DNN automatically. To address this need, we present SpecRepair, a tool that efficiently eliminates counter-examples from a DNN and produces a provably safe DNN without harming its classification accuracy. SpecRepair combines specification-based counter-example search and resumes training of the DNN, penalizing counter-examples and certifying the resulting DNN. We evaluate SpecRepair's effectiveness on the ACAS Xu benchmark, a DNN-based controller for unmanned aircraft, and two image classification benchmarks. The results show that SpecRepair is more successful in producing safe DNNs than comparable methods, has a shorter runtime, and produces safe DNNs while preserving their classification accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Verification of Neural-Network Control Systems by Integrating Taylor Models and Zonotopes

We study the verification problem for closed-loop dynamical systems with neural-network controllers (NNCS). This problem is commonly reduced to computing the set of reachable states. When considering dynamical systems and neural networks in isolation, there exist precise approaches for that task based on set representations respectively called Taylor models and zonotopes. However, the combination of these approaches to NNCS is non-trivial because, when converting between the set representations, dependency information gets lost in each control cycle and the accumulated approximation error quickly renders the result useless. We present an algorithm to chain approaches based on Taylor models and zonotopes, yielding a precise reachability algorithm for NNCS. Because the algorithm only acts at the interface of the isolated approaches, it is applicable to general dynamical systems and neural networks and can benefit from future advances in these areas. Our implementation delivers state-of-the-art performance and is the first to successfully analyze all benchmark problems of an annual reachability competition for NNCS.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Concept of orbital entanglement and correlation in quantum chemistry

A recent development in quantum chemistry has established the quantum mutual information between orbitals as a major descriptor of electronic structure. This has already facilitated remarkable improvements of numerical methods and may lead to a more comprehensive foundation for chemical bonding theory. Building on this promising development, our work provides a refined discussion of quantum information theoretical concepts by introducing the physical correlation and its separation into classical and quantum parts as distinctive quantifiers of electronic structure. In particular, we succeed in quantifying the entanglement. Intriguingly, our results for different molecules reveal that the total correlation between orbitals is mainly classical, raising questions about the general significance of entanglement in chemical bonding. Our work also shows that implementing the fundamental particle number superselection rule, so far not accounted for in quantum chemistry, removes a major part of correlation and entanglement previously seen. In that respect, realizing quantum information processing tasks with molecular systems might be more challenging than anticipated.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

How creating one additional well can generate Bose-Einstein condensation

The realization of Bose-Einstein condensation in ultracold trapped gases has led to a revival of interest in that fascinating quantum phenomenon. This experimental achievement necessitated both extremely low temperatures and sufficiently weak interactions. Particularly in reduced spatial dimensionality even an infinitesimal interaction immediately leads to a departure to quasi-condensation. We propose a system of strongly interacting bosons which overcomes those obstacles by exhibiting a number of intriguing related features: (i) The tuning of just a single control parameter drives a transition from quasi-condensation to complete condensation, (ii) the destructive influence of strong interactions is compensated by the respective increased mobility, (iii) topology plays a crucial role since a crossover from one- to `infinite'-dimensionality is simulated, (iv) a ground state gap opens which makes the condensation robust to thermal noise. Remarkably, all these features can be derived by analytical and exact numerical means despite the non-perturbative character of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Into the Unknown: Active Monitoring of Neural Networks

Neural-network classifiers achieve high accuracy when predicting the class of an input that they were trained to identify. Maintaining this accuracy in dynamic environments, where inputs frequently fall outside the fixed set of initially known classes, remains a challenge. The typical approach is to detect inputs from novel classes and retrain the classifier on an augmented dataset. However, not only the classifier but also the detection mechanism needs to adapt in order to distinguish between newly learned and yet unknown input classes. To address this challenge, we introduce an algorithmic framework for active monitoring of a neural network. A monitor wrapped in our framework operates in parallel with the neural network and interacts with a human user via a series of interpretable labeling queries for incremental adaptation. In addition, we propose an adaptive quantitative monitor to improve precision. An experimental evaluation on a diverse set of benchmarks with varying numbers of classes confirms the benefits of our active monitoring framework in dynamic scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Synthesis of Hybrid Automata with Affine Dynamics from Time-Series Data

Formal design of embedded and cyber-physical systems relies on mathematical modeling. In this paper, we consider the model class of hybrid automata whose dynamics are defined by affine differential equations. Given a set of time-series data, we present an algorithmic approach to synthesize a hybrid automaton exhibiting behavior that is close to the data, up to a specified precision, and changes in synchrony with the data. A fundamental problem in our synthesis algorithm is to check membership of a time series in a hybrid automaton. Our solution integrates reachability and optimization techniques for affine dynamical systems to obtain both a sufficient and a necessary condition for membership, combined in a refinement framework. The algorithm processes one time series at a time and hence can be interrupted, provide an intermediate result, and be resumed. We report experimental results demonstrating the applicability of our synthesis approach.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Correlation paradox of the dissociation limit: A quantum information perspective

The interplay between electron interaction and geometry in a molecular system can lead to rather paradoxical situations. The prime example is the dissociation limit of the hydrogen molecule: While a significant increase of the distance $r$ between the two nuclei marginalizes the electron-electron interaction, the exact ground state does, however, not take the form of a single Slater determinant. By first reviewing and then employing concepts from quantum information theory, we resolve this paradox and its generalizations to more complex systems in a quantitative way. To be more specific, we illustrate and prove that thermal noise due to finite, possibly even just infinitesimally low, temperature $T$ will destroy the entanglement beyond a critical separation distance $r_{\mathrm{crit}}$($T$) entirely. Our analysis is comprehensive in the sense that we simultaneously discuss both total correlation and entanglement in the particle picture as well as in the orbital/mode picture. Our results reveal a conceptually new characterization of static and dynamical correlation in ground states by relating them to the (non)robustness of correlation with respect to thermal noise.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Outside the Box: Abstraction-Based Monitoring of Neural Networks

Neural networks have demonstrated unmatched performance in a range of classification tasks. Despite numerous efforts of the research community, novelty detection remains one of the significant limitations of neural networks. The ability to identify previously unseen inputs as novel is crucial for our understanding of the decisions made by neural networks. At runtime, inputs not falling into any of the categories learned during training cannot be classified correctly by the neural network. Existing approaches treat the neural network as a black box and try to detect novel inputs based on the confidence of the output predictions. However, neural networks are not trained to reduce their confidence for novel inputs, which limits the effectiveness of these approaches. We propose a framework to monitor a neural network by observing the hidden layers. We employ a common abstraction from program analysis - boxes - to identify novel behaviors in the monitored layers, i.e., inputs that cause behaviors outside the box. For each neuron, the boxes range over the values seen in training. The framework is efficient and flexible to achieve a desired trade-off between raising false warnings and detecting novel inputs. We illustrate the performance and the robustness to variability in the unknown classes on popular image-classification benchmarks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reachability analysis of linear hybrid systems via block decomposition

Reachability analysis aims at identifying states reachable by a system within a given time horizon. This task is known to be computationally expensive for linear hybrid systems. Reachability analysis works by iteratively applying continuous and discrete post operators to compute states reachable according to continuous and discrete dynamics, respectively. In this paper, we enhance both of these operators and make sure that most of the involved computations are performed in low-dimensional state space. In particular, we improve the continuous-post operator by performing computations in high-dimensional state space only for time intervals relevant for the subsequent application of the discrete-post operator. Furthermore, the new discrete-post operator performs low-dimensional computations by leveraging the structure of the guard and assignment of a considered transition. We illustrate the potential of our approach on a number of challenging benchmarks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reduced Density Matrix Functional Theory for Bosons

Based on a generalization of Hohenberg-Kohn's theorem, we propose a ground state theory for bosonic quantum systems. Since it involves the one-particle reduced density matrix $γ$ as a natural variable but still recovers quantum correlations in an exact way it is particularly well-suited for the accurate description of Bose-Einstein condensates. As a proof of principle we study the building block of optical lattices. The solution of the underlying $v$-representability problem is found and its peculiar form identifies the constrained search formalism as the ideal starting point for constructing accurate functional approximations: The exact functionals for this $N$-boson Hubbard dimer and general Bogoliubov-approximated systems are determined. The respective gradient forces are found to diverge in the regime of Bose-Einstein condensation, $\nabla_γ \mathcal{F} \propto 1/\sqrt{1-N_{\mathrm{BEC}}/N}$, providing a natural explanation for the absence of complete BEC in nature.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Implications of pinned occupation numbers for natural orbital expansions. I: Generalizing the concept of active spaces

The concept of active spaces simplifies the description of interacting quantum many-body systems by restricting to a neighbourhood of active orbitals around the Fermi level. The respective wavefunction ansatzes which involve all possible electron configurations of active orbitals can be characterized by the saturation of a certain number of Pauli constraints $0 \leq n_i \leq 1$, identifying the occupied core orbitals ($n_i=1$) and the inactive virtual orbitals ($n_j=0$). In Part I, we generalize this crucial concept of active spaces by referring to the generalized Pauli constraints. To be more specific, we explain and illustrate that the saturation of any such constraint on fermionic occupation numbers characterizes a distinctive set of active electron configurations. A converse form of this selection rule establishes the basis for corresponding multiconfigurational wavefunction ansatzes. In Part II, we provide rigorous derivations of those findings. Moroever, we extend our results to non-fermionic multipartite quantum systems, revealing that extremal single-body information has always strong implications for the multipartite quantum state. In that sense, our work also confirms that pinned quantum systems define new physical entities and the presence of pinnings reflect the existence of (possibly hidden) ground state symmetries.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Implications of pinned occupation numbers for natural orbital expansions. II: Rigorous derivation and extension to non-fermionic systems

We have explained and comprehensively illustrated in Part I that the generalized Pauli constraints suggest a natural extension of the concept of active spaces. In the present Part II, we provide rigorous derivations of the theorems involved therein. This will offer in particular deeper insights into the underlying mathematical structure and will explain why the saturation of generalized Pauli constraints implies a specific simplified structure of the corresponding many-fermion quantum state. Moreover, we extend the results of Part I to non-fermionic multipartite quantum systems, revealing that extremal single-body information has always strong implications for the multipartite quantum state. In that sense, our work also confirms that pinned quantum systems define new physical entities and the presence of pinnings reflect the existence of (possibly hidden) ground state symmetries.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Reach Set Approximation through Decomposition with Low-dimensional Sets and High-dimensional Matrices

Approximating the set of reachable states of a dynamical system is an algorithmic yet mathematically rigorous way to reason about its safety. Although progress has been made in the development of efficient algorithms for affine dynamical systems, available algorithms still lack scalability to ensure their wide adoption in the industrial setting. While modern linear algebra packages are efficient for matrices with tens of thousands of dimensions, set-based image computations are limited to a few hundred. We propose to decompose reach set computations such that set operations are performed in low dimensions, while matrix operations like exponentiation are carried out in the full dimension. Our method is applicable both in dense- and discrete-time settings. For a set of standard benchmarks, it shows a speed-up of up to two orders of magnitude compared to the respective state-of-the art tools, with only modest losses in accuracy. For the dense-time case, we show an experiment with more than 10.000 variables, roughly two orders of magnitude higher than possible with previous approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Instrumenting an SMT Solver to Solve Hybrid Network Reachability Problems

PDDL+ planning has its semantics rooted in hybrid automata (HA) and recent work has shown that it can be modeled as a network of HAs. Addressing the complexity of nonlinear PDDL+ planning as HAs requires both space and time efficient reasoning. Unfortunately, existing solvers either do not address nonlinear dynamics or do not natively support networks of automata. We present a new algorithm, called HNSolve, which guides the variable selection of the dReal Satisfiability Modulo Theories (SMT) solver while reasoning about network encodings of nonlinear PDDL+ planning as HAs. HNSolve tightly integrates with dReal by solving a discrete abstraction of the HA network. HNSolve finds composite runs on the HA network that ignore continuous variables, but respect mode jumps and synchronization labels. HNSolve admissibly detects dead-ends in the discrete abstraction, and posts conflict clauses that prune the SMT solver's search. We evaluate the benefits of our HNSolve algorithm on PDDL+ benchmark problems and demonstrate its performance with respect to prior work.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Invariant Clusters for Hybrid Systems

In this paper, we propose an approach to automatically compute invariant clusters for semialgebraic hybrid systems. An invariant cluster for an ordinary differential equation (ODE) is a multivariate polynomial invariant g(u,x)=0, parametric in u, which can yield an infinite number of concrete invariants by assigning different values to u so that every trajectory of the system can be overapproximated precisely by a union of concrete invariants. For semialgebraic systems, which involve ODEs with multivariate polynomial vector flow, invariant clusters can be obtained by first computing the remainder of the Lie derivative of a template multivariate polynomial w.r.t. its Groebner basis and then solving the system of polynomial equations obtained from the coefficients of the remainder. Based on invariant clusters and sum-of-squares (SOS) programming, we present a new method for the safety verification of hybrid systems. Experiments on nonlinear benchmark systems from biology and control theory show that our approach is effective and efficient.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Number-parity effect for confined fermions in one dimension

For $N$ spin-polarized fermions with harmonic pair interactions in a $1$-dimensional trap an odd-even effect is found. The spectrum of the $1$-particle reduced density matrix of the system's ground state differs qualitatively for $N$ odd and $N$ even. This effect does only occur for strong attractive and repulsive interactions. Since it does not exists for bosons, it must originate from the repulsive nature implied by the fermionic exchange statistics. In contrast to the spectrum, the $1$-particle density and correlation function for strong attractive interactions do not show any sensitivity on the number parity. This also suggests that reduced-density-matrix-functional theory has a more subtle $N$-dependency than density functional theory.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pinning of Fermionic Occupation Numbers: General Concepts and One Dimension

Analytical evidence for the physical relevance of generalized Pauli constraints (GPCs) has recently been provided in [PRL 110, 040404]: Natural occupation numbers $\vecλ\equiv (λ_i)$ of the ground state of a model system in the regime of weak couplings $κ$ of three spinless fermions in one spatial dimension were found extremely close, in a distance $D_{min}\sim κ^8$ to the boundary of the allowed region. We provide a self-contained and complete study of this quasipinning phenomenon. In particular, we develop tools for its systematic exploration and quantification. We confirm that quasipinning in one dimension occurs also for larger particle numbers and extends to intermediate coupling strengths, but vanishes for very strong couplings. We further explore the non-triviality of our findings by comparing quasipinning by GPCs to potential quasipinning by the less restrictive Pauli exclusion principle constraints. This allows us to eventually confirm the significance of GPCs beyond Pauli's exclusion principle.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pinning of fermionic occupation numbers: Higher spatial dimensions and spin

The role of the generalized Pauli constraints (GPCs) in higher spatial dimensions and by incorporating spin degrees of freedom is systematically explored for a system of interacting fermions confined by a harmonic trap. Physical relevance of the GPCs is confirmed by analytical means for the ground state in the regime of weak couplings by finding its vector of natural occupation numbers close to the boundary of the allowed region. Such quasipinning is found to become weaker in the intermediate and strong coupling regime. The study of crossovers between different spatial dimensions by detuning the harmonic trap frequencies suggests that quasipinning is essentially an effect for systems with reduced spatial dimensionality. In addition, we find that quasipinning becomes stronger by increasing the degree of spin-polarization. Consequently, the number of states available around the Fermi level plays a key role for the occurrence of quasipinning. This suggests that quasipinning emerges from the conflict between energy minimization and fermionic exchange symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hubbard model: Pinning of occupation numbers and role of symmetries

Fermionic natural occupation numbers do not only obey Pauli's exclusion principle, but are even further restricted by so-called generalized Pauli constraints. Such restrictions are particularly relevant whenever they are saturated by given natural occupation numbers $\vecλ=(λ_i)$. For few-site Hubbard models we explore the occurrence of this pinning effect. By varying the on-site interaction $U$ for the fermions we find sharp transitions from pinning of $\vecλ$ to the boundary of the allowed region to nonpinning. We analyze the origin of this phenomenon which turns out be either a crossing of natural occupation numbers $λ_{i}(U), λ_{i+1}(U)$ or a crossing of $N$-particle energies. Furthermore, we emphasize the relevance of symmetries for the occurrence of pinning. Based on recent progress in the field of ultracold atoms our findings suggest an experimental set-up for the realization of the pinning effect.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum Marginal Problem and its Physical Relevance

The Pauli exclusion principle as constraint on fermionic occupation numbers is a consequence of the much deeper fermionic exchange statistics. Just recently, it was shown by Klyachko that this antisymmetry of fermionic wave functions leads to further restrictions on natural occupation numbers. These so-called generalized Pauli constraints (GPC) significantly strengthen Pauli's exclusion principle. Our first goal is to develop an understanding of the mathematical concepts behind Klyachko's work, in the context of quantum marginal problems. Afterwards, we explore the physical relevance of GPC and study concrete physical systems from that new viewpoint. In the first part of this thesis we review Klyachko's solution of the univariate quantum marginal problem. In particular we break his abstract derivation based on algebraic topology down to a more elementary level and reveal the geometrical picture behind it. The second part explores the possible physical relevance of GPC. We review the effect of pinning, i.e. the saturation of some GPC by given natural occupation numbers and explain its consequences. Although this effect would be quite spectacular we argue that pinning is unnatural. Instead, we conjecture the effect of quasipinning, defined by occupation numbers close to (but not exactly on) the boundary of the allowed region. In the third part we study concrete fermionic quantum systems from the new viewpoint of GPC. In particular, we compute the natural occupation numbers for the ground state of a family of interacting fermions in a harmonic potential. Intriguingly, we find that the occupation numbers are strongly quasipinned, even up to medium interaction strengths. We identify this as an effect of the lowest few energy eigenstates, which provides first insights into the mechanism behind quasipinning.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Duality of reduced density matrices and their eigenvalues

For states of quantum systems of $N$ particles with harmonic interactions we prove that each reduced density matrix $ρ$ obeys a duality condition. This condition implies duality relations for the eigenvalues $λ_k$ of $ρ$ and relates a harmonic model with length scales $l_1,l_2, \ldots, l_N $ with another one with inverse lengths $1/l_1, 1/l_2,\ldots, 1/l_N$. Entanglement entropies and correlation functions inherit duality from $ρ$. Self-duality can only occur for noninteracting particles in an isotropic harmonic trap.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quasipinning and its relevance for $N$-Fermion quantum states

Fermionic natural occupation numbers (NON) do not only obey Pauli's famous exclusion principle but are even further restricted to a polytope by the generalized Pauli constraints, conditions which follow from the fermionic exchange statistics. Whenever given NON are pinned to the polytope's boundary the corresponding $N$-fermion quantum state $|Ψ_N\rangle$ simplifies due to a selection rule. We show analytically and numerically for the most relevant settings that this rule is stable for NON close to the boundary, if the NON are non-degenerate. In case of degeneracy a modified selection rule is conjectured and its validity is supported. As a consequence the recently found effect of quasipinning is physically relevant in the sense that its occurrence allows to approximately reconstruct $|Ψ_N\rangle$, its entanglement properties and correlations from 1-particle information. Our finding also provides the basis for a generalized Hartree-Fock method by a variational ansatz determined by the selection rule.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Quantum Marginal Problem

The question of whether given density operators for subsystems of a multipartite quantum system are compatible to one common total density operator is known as the quantum marginal problem. We briefly review the solution of a subclass of such problems found just recently. In particular, this provides the solution of the $1$-body $N$-representability problem. Its solution, the so-called generalized Pauli constraints, restrict the set of mathematically possible fermionic occupation numbers significantly, and strengthens Pauli's exclusion principle. Moreover, we review the study of a concrete physical model of interacting fermions confined to a harmonic trap. There, we found occupation numbers close, but not exactly on the boundary of the allowed region. This new effect of quasipinning is physically relevant since it corresponds to a simplified structure of the corresponding $N$-fermion quantum state.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Natural Orbitals and Occupation Numbers for Harmonium: Fermions vs. Bosons

For a quantum system of N identical, harmonically interacting particles in a one-dimensional harmonic trap we calculate for the bosonic and fermionic ground state the corresponding 1-particle reduced density operator $ρ_1$ analytically. In case of bosons $ρ_1$ is a Gibbs state for an effective harmonic oscillator. Hence the natural orbitals are Hermite functions and their occupation numbers obey a Boltzmann distribution. Intriguingly, for fermions with not too large couplings the natural orbitals coincide up to just a very small error with the bosonic ones. In case of strong coupling this still holds qualitatively. Moreover, the decay of the decreasingly ordered fermionic natural occupation numbers is given by the bosonic one, but modified by an algebraic prefactor. Significant differences to bosons occur only for the largest occupation numbers. After all the "discontinuity" at the "Fermi level" decreases with increasing coupling strength but remains well pronounced even for strong interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pinning of Fermionic Occupation Numbers

The Pauli exclusion principle is a constraint on the natural occupation numbers of fermionic states. It has been suspected since at least the 1970's, and only proved very recently, that there is a multitude of further constraints on these numbers, generalizing the Pauli principle. Here, we provide the first analytic analysis of the physical relevance of these constraints. We compute the natural occupation numbers for the ground states of a family of interacting fermions in a harmonic potential. Intriguingly, we find that the occupation numbers are almost, but not exactly, pinned to the boundary of the allowed region (quasi-pinned). The result suggests that the physics behind the phenomenon is richer than previously appreciated. In particular, it shows that for some models, the generalized Pauli constraints play a role for the ground state, even though they do not limit the ground-state energy. Our findings suggest a generalization of the Hartree-Fock approximation.

0 citations0 reviewsNo code linkedOpen access