As a musical instrument construction material, wood is both musically and aesthetically pleasing. Easy to work and abundant, it has traditionally been the material of choice. Unfortunately, wood is also a very inconsistent material. Due to great environmental and climatic variations, wooden specimens present large variations in their mechanical properties, even within species of a similar region. Surprisingly, an industry based entirely on acoustics has done very little to account for these variations. For this reason, manufactured wooden guitars are acoustically inconsistent. Previous work has shown that varying the dimensions of a guitar soundboard brace is a good method for taking into account variations in the mechanical properties of the wooden soundboard plate. In this thesis, the effects of a scalloped-shaped brace on the natural frequencies of a brace-plate system have been studied and tools have been developed in order to calculate the dimensions of the brace required to account for variations in the mechanical properties of the plate. It has been shown that scalloped braces can be used to modify two natural frequencies of a brace-plate system simultaneously. Furthermore, the most important criteria in modifying any given frequency of a brace-plate system is the mass and stiffness properties of the brace at the antinode of the given frequency’s associated modeshape. Subsequently, designing a brace for desired system natural frequencies, by taking into account the mechanical properties of the wooden plate, is an inverse eigenvalue problem. Since few methods exist for solving the inverse eigenvalue problem of general matrices, a new method based on the generalized Cayley-Hamilton theorem was proposed in the thesis. A further method, based on the determinant of the generalized eigenvalue problem was also presented. Both methods work well, although the determinant method is shown to be more efficient for partially described systems. Finally, experimental results were obtained for the natural frequencies of simply supported wooden plates, with and without a brace, as well as the inverse eigenvalue determinant method. Good correlation was found between theoretical and experimental results.

Cette thèse est consacrée à l'étude de la méthode PAW (projector augmented-wave) et d'une de ses modifications, baptisée méthode PAW variationnelle (VPAW), pour le calcul de l'état fondamental d'Hamiltoniens en géométrie périodique. Ces méthodes visent à améliorer la vitesse de convergence des méthodes d'ondes planes (ou méthodes de Fourier) en appliquant une transformation inversible au problème aux valeurs propres initial agissant au voisinage de chaque site atomique. Cette transformation permet de capter une partie des difficultés dues aux singularités coulombiennes. La méthode VPAW est analysée pour un opérateur de Schr\"odinger unidimensionnel avec des potentiels de Dirac. Les fonctions propres de ce modèle comprennent des sauts de dérivées similaires aux cusps électroniques. Le saut de dérivée des fonctions propres du problème aux valeurs propres issu de la méthode VPAW est réduit de façon importante. Cela entraîne une accélération de convergence en ondes planes du calcul des valeurs propres corroborée par une étude numérique. Une étude de la méthode VPAW est conduite pour des Hamiltoniens 3D périodiques avec des singularités coulombiennes, parvenant à des conclusions similaires. Pour la méthode PAW, la transformation inversible comporte des sommes infinies qui sont tronquées en pratique. Ceci introduit une erreur, qui est rarement quantifiée en pratique. Elle est analysée dans le cas de l'opérateur de Schrödinger unidimensionnel avec des potentiels de Dirac. Des bornes sur la plus basse valeur propre en fonction des paramètres PAW sont prouvées conformes aux tests numériques. / This thesis is devoted to the study of the PAW method (projector augmented-wave) and of a variant called the variational PAW method (VPAW). These methods aim to accelerate the convergence of plane-wave methods in electronic structure calculations. They rely on an invertible transformation applied to the eigenvalue problem, which acts in a neighborhood of each atomic site. The transformation captures some difficulties caused by the Coulomb singularities. The VPAW method is applied to a periodic one-dimensional Schr\"odinger operator with Dirac potentials and analyzed in this setting. Eigenfunctions of this model have derivative jumps similar to the electronic cusps. The derivative jumps of eigenfunctions of the VPAW eigenvalue problem are significantly reduced. Hence, a smaller plane-wave cut-off is required for a given accuracy level. The study of the VPAW method is also carried out for 3D periodic Hamiltonians with Coulomb singularities yielding similar results. In the PAW method, the invertible transformation has infinite sums that are truncated in practice. The induced error is analyzed in the case of the periodic one-dimensional Schrödinger operator with Dirac potentials. Error bounds on the lowest eigenvalue are proved depending on the PAW parameters.

This thesis describes improvements of the periodic QZ algorithm and several variants of the Schur algorithm for block Toeplitz matrices. Documentation of the available software is included.

Cette thèse considère la solution numérique d'un problème aux valeurs propres de grandes dimensions, dans lequel l'opérateur est dérivé d'un problème de transfert radiatif. Ainsi, cette thèse étudie l'utilisation de matrices hiérarchiques, une représentation efficace de tableaux, très intéressante pour une utilisation avec des problèmes de grandes dimensions. Les matrices sont des représentations hiérarchiques de structures de données efficaces pour les matrices denses, l'idée de base étant la division d'une matrice en une hiérarchie de blocs et l´approximation de certains blocs par une matrice de petite caractéristique. Son utilisation permet de diminuer la mémoire nécessaire tout en réduisant les coûts informatiques. L'application de l'utilisation de matrices hiérarchique est analysée dans le contexte de la solution numérique d'un problème aux valeurs propres de grandes dimensions résultant de la discrétisation d'un opérateur intégral. L'opérateur est de convolution et est défini par la première fonction exponentielle intégrale, donc faiblement singulière. Pour le calcul informatique, nous avons accès à HLIB (Hierarchical matrices LIBrary) qui fournit des routines pour la construction de la structure hiérarchique des matrices et des algorithmes pour les opérations approximative avec ces matrices. Nous incorporons certaines routines comme la multiplication matrice-vecteur ou la decomposition LU, en SLEPc (Hierarchical matrices LIBrary) pour explorer les algorithmes existants afin de résoudre les problèmes de valeur propre. Nous développons aussi des expressions analytiques pour l'approximation des noyaux dégénérés utilisés dans la thèse et déduire ainsi les limites supérieures d'erreur pour ces approximations. Les résultats numériques obtenus avec d'autres techniques pour résoudre le problème en question sont utilisés pour la comparaison avec ceux obtenus avec la nouvelle technique, illustrant l'efficacité de ce dernier / In this thesis, we consider the numerical solution of a large eigenvalue problem in which the integral operator comes from a radiative transfer problem. It is considered the use of hierarchical matrices, an efficient data-sparse representation of matrices, especially useful for large dimensional problems. It consists on low-rank subblocks leading to low memory requirements as well as cheap computational costs. We discuss the use of the hierarchical matrix technique in the numerical solution of a large scale eigenvalue problem arising from a finite rank discretization of an integral operator. The operator is of convolution type, it is defined through the first exponential-integral function and hence it is weakly singular. We access HLIB (Hierarchical matrices LIBrary) that provides, among others, routines for the construction of hierarchical matrix structures and arithmetic algorithms to perform approximative matrix operations. Moreover, it is incorporated the matrix-vector multiply routines from HLIB, as well as LU factorization for preconditioning, into SLEPc (Scalable Library for Eigenvalue Problem Computations) in order to exploit the available algorithms to solve eigenvalue problems. It is also developed analytical expressions for the approximate degenerate kernels and deducted error upper bounds for these approximations. The numerical results obtained with other approaches to solve the problem are used to compare with the ones obtained with this technique, illustrating the efficiency of the techniques developed and implemented in this work

In this thesis, a general pseudospectral formulation for a class of Sturm-Liouville eigenvalue problems is consructed. It is shown that almost all, regular or singular, Sturm-Liouville eigenvalue problems in the Schr&ouml / dinger form may be transformed into a more tractable form. This tractable form will be called here a weighted equation of hypergeometric type with a perturbation (WEHTP) since the non-weighted and unperturbed part of it is known as the equation of hypergeometric type (EHT). It is well known that the EHT has polynomial solutions which form a basis for the Hilbert space of square integrable functions. Pseudospectral methods based on this natural expansion basis are constructed to approximate the eigenvalues of WEHTP, and hence the energy eigenvalues of the Schr&ouml / dinger equation. Exemplary computations are performed to support the convergence numerically.

This thesis is concerned with the finite element analysis and the a posteriori error estimation for eigenvalue problems for general operator pencils on two-dimensional manifolds. A specific application of the presented theory is the computation of corner singularities. Engineers use the knowledge of the so-called singularity exponents to predict the onset and the propagation of cracks. All results of this thesis are explained for two model problems, the Laplace and the linear elasticity problem, and verified by numerous numerical results.

L'objectif de cette thèse est de fournir des bornes d'erreur pour des problèmes aux valeurs propres linéaires et non linéaires issus du calcul de structure électronique, en particulier celui de l'état fondamental avec la théorie de la fonctionnelle de la densité. Ces bornes d'erreur reposent principalement sur des estimations a posteriori. D'abord, nous étudions un phénomène de compensation d'erreur de discrétisation pour un problème linéaire aux valeurs propres, grâce à une analyse a priori de l'erreur sur l'énergie. Ensuite, nous présentons une analyse a posteriori pour le problème du laplacien aux valeurs propres discrétisé par une large classe d'éléments finis. Les bornes d'erreur proposées pour les valeurs propres simples et leurs vecteurs propres associés sont garanties, calculables et efficaces. Nous nous concentrons alors sur des problèmes aux valeurs propres non linéaires. Nous proposons des bornes d'erreur pour l'équation de Gross-Pitaevskii, valables sous des hypothèses vérifiables numériquement, et pouvant être séparées en deux composantes venant respectivement de la discrétisation et de l'algorithme itératif utilisé pour résoudre le problème non linéaire aux valeurs propres. L'équilibrage de ces composantes d'erreur permet d'optimiser les ressources numériques. Enfin, nous présentons une méthode de post-traitement pour le problème de Kohn-Sham discrétisé en ondes planes, améliorant la précision des résultats à un faible coût de calcul. Les solutions post-traitées peuvent être utilisées soit comme solutions plus précises du problème, soit pour calculer une estimation de l'erreur de discrétisation, qui n'est plus garantie, mais néanmoins proche de l'erreur. / The objective of this thesis is to provide error bounds for linear and nonlinear eigenvalue problems arising from electronic structure calculation. We focus on ground-state calculations based on Density Functional Theory, including Kohn-Sham models. Our bounds mostly rely on a posteriori error analysis. More precisely, we start by studying a phenomenon of discretization error cancellation for a simple linear eigenvalue problem, for which analytical solutions are available. The mathematical study is based on an a priori analysis for the energy error. Then, we present an a posteriori analysis for the Laplace eigenvalue problem discretized with finite elements. For simple eigenvalues of the Laplace operator and their corresponding eigenvectors , we provide guaranteed, fully computable and efficient error bounds. Thereafter, we focus on nonlinear eigenvalue problems. First, we provide an a posteriori analysis for the Gross-Pitaevskii equation. The error bounds are valid under assumptions that can be numerically checked, and can be separated in two components coming respectively from the discretization and the iterative algorithm used to solve the nonlinear eigenvalue problem. Balancing these error components allows to optimize the computational resources. Second, we present a post-processing method for the Kohn-Sham problem, which improves the accuracy of planewave computations of ground state orbitals at a low computational cost. The post-processed solutions can be used either as a more precise solution of the problem, or used for computing an estimation of the discretization error. This estimation is not guaranteed, but in practice close to the real error.

[EN] The neutron diffusion equation is an approximation of the neutron transport equation that describes the neutron population in a nuclear reactor core. In particular, we will consider here VVER-type reactors which use the neutron diffusion equation discretized on hexagonal meshes. Most of the simulation codes of a nuclear power reactor use the multigroup neutron diffusion equation to describe the neutron distribution inside the reactor core.To study the stationary state of a reactor, the reactor criticality is forced in artificial way leading to a generalized differential eigenvalue problem, known as the Lambda modes equation, which is solved to obtain the dominant eigenvalues of the reactor and their corresponding eigenfunctions. To discretize this model a finite element method with h-p adaptivity is used. This method allows to use heterogeneous meshes, and allows different refinements such as the use of h-adaptive meshes, reducing the size of specific cells, and p-refinement, increasing the polynomial degree of the basic functions used in the expansions of the solution in the different cells. Once the solution for the steady state neutron distribution is obtained, it is used as initial condition for the time integration of the neutron diffusion equation. To simulate the behaviour of a nuclear power reactor it is necessary to be able to integrate the time-dependent neutron diffusion equation inside the reactor core. The spatial discretization of this equation is done using a finite element method that permits h-p refinements for different geometries. Transients involving the movement of the control rod banks have the problem known as the rod-cusping effect. Previous studies have usually approached the problem using a fixed mesh scheme defining averaged material properties and many techniques exist for the treatment of the rod cusping problem. The present work proposes the use of a moving mesh scheme that uses spatial meshes that change with the movement of the control rods avoiding the necessity of using equivalent material cross sections for the partially inserted cells. The performance of the moving mesh scheme is tested studying different benchmark problems. For reactor calculations, the accuracy of a diffusion theory solution is limited for for complex fuel assemblies or fine mesh calculations. To improve these results a method that incorporates higher-order approximations for the angular dependence, as the simplified spherical harmonics (SPN ) method must be employed. In this work an h-p Finite Element Method (FEM) is used to obtain the dominant Lambda mode associated with a configuration of a reactor core using the SPN approximation. The performance of the SPN (N= 1, 3, 5) approximations has been tested for different reactor benchmarks. / [ES] La ecuación de la difusión neutrónica es una aproximación de la ecuación del transporte de neutrones que describe la población de neutrones en el núcleo de un reactor nuclear. En particular, consideraremos reactores de tipo VVER y para simular su comportamiento se utilizará la ecuación de la difusión neutrónica para cuya discretización se hace uso de mallas hexagonales. La mayoría de los códigos de simulación de reactores nucleares utilizan aproximación multigrupo de energía de la ecuación de la difusión neutrónica para describir la distribución de neutrones en el interior del núcleo del reactor. Para estudiar el estado estacionario del reactor, es posible forzar la criticidad del reactor de forma artificial modificando las secciones eficaces de forma que se obtiene un problema de valores propios diferencial, conocido como el problema de los Modos Lambda, que se resuelve para obtener los valores propios dominantes del reactor y sus correspondientes funciones propias. Para discretizar este modelo se ha hecho uso de un método de elementos finitos con adaptabilidad h-p. Este método permite el uso de mallas heterogéneas, y de diferentes refinamientos como el uso mallas h-adaptativas, reduciendo el tamaño de los distintos nodos, y el p-refinado, aumentando el grado del polinomio de las funciones básicas utilizado en los desarrollos de la solución en los diferentes nodos. Se ha desarrollado un código basado en un método de elementos finitos de alto orden para resolver el problema de los Modos Lambda en un reactor con geometría hexagonal y se han obtenido los Modos dominantes para distintos problemas de referencia. Una vez que se ha obtenido la solución para la distribución de neutrones en estado estacionario, ésta se utiliza como condición inicial para la integración de la ecuación de difusión neutrónica dependiente del tiempo. Para simular el comportamiento de un reactor nuclear para un determinado transitorio, es necesario ser capaz de integrar la ecuación de la difusión neutrónica dependiente del tiempo en el interior del núcleo del reactor. La discretización espacial de esta ecuación se hace usando un método de elementos finitos de alto orden que permite refinados de tipo h-p para distintas geometrías. Los transitorios que implican el movimiento de los bancos de las barras de control tienen el problema conocido como el efecto 'rod-cusping'. Estudios anteriores, por lo general, han abordado este problema utilizando una malla fija y definiendo propiedades promedio para los materiales correspondientes a las celdas donde se tiene la barra de control parcialmente insertada. En el presente trabajo se propone el uso de un esquema de malla móvil, de forma que en mallado espacial va cambiando con el movimiento de la barra de control, evitando la necesidad de utilizar secciones eficaces equivalentes para las celdas parcialmente insertadas. El funcionamiento de este esquema de malla móvil propuesto se estudia resolviendo distintos problemas tipo. La precisión obtenida mediante de la teoría de la difusión en los cálculos de reactores es limitada cuando se tienen elementos de combustible complejos o se pretenden realizar cálculos en malla fina. Para mejorar estos resultados, es necesario disponer de un método que incorpore aproximaciones de orden superior de la ecuación del transporte de neutrones. Una posibilidad es hacer uso de las ecuaciones PN simplificadas (SPN ). En este trabajo se utiliza un método de elementos finitos h-p para obtener los modos dominantes asociados con una configuración dada del núcleo de un reactor nuclear con geometría hexagonal usando la aproximación SPN . El funcionamiento de las aproximaciones SPN (N = 1, 3, 5) se ha estudiado para distintos problemas de referencia. / [CAT] L'equació de la difusió neutrònica és una aproximació de l'equació del transport de neutrons que descriu la població de neutrons en el nucli de un reactor nuclear. En particular, considerarem reactors de tipus VVER i per a simular el seu comportament s'utilitzarà l'equació de la difusió neutrónica que es discretitza fent ús de malles hexagonals. La majoria dels codis de simulació de reactors nuclears utilitzen l'aproximació multigrup d'energia de l'equació de la difusió neutrónica per a descriure la distribució de neutrons a l'interior del nucli del reactor. Per a estudiar l'estat estacionari del reactor, és possible forçar la seua criticitat de forma artificial modificant les seccions eficaces de manera que s'obté un problema de valors propis diferencial, conegut com el problema dels Modes Lambda, que es resol per a obtenir els valors propis dominants del reactor i les seues corresponents funcions pròpies. Per a discretitzar aquest model s'ha fet ús d'un mètode d'elements finits amb adaptabilitat h-p. Aquest mètode permet l'ús de malles heterogènies, i de diferents refinaments com l'ús malles h-adaptatives, reduint la grandària dels diferents nodes, i el p-refinat, augmentant el grau del polinomi de les funcions bàsiques utilitzat en els desenvolupaments de la solució en els diferents nodes. S'ha desenvolupat un codi basat en un mètode d'elements finits d'alt ordre per a resoldre el problema dels Modes Lambda en un reactor amb geometria hexagonal i s'han obtingut els Modes dominants per a diferents problemes de referència. Una vegada que s'ha obtingut la solució per a la distribució de neutrons en estat estacionari, aquesta s'utilitza com a condició inicial per a la integració de l'equació de difusió neutrònica depenent del temps. Per a simular el comportament d'un reactor nuclear per a un determinat transitori, és necessari ser capaç d'integrar l'equació de la difusió neutrónica depenent del temps a l'interior del nucli del reactor. La discretitzación espacial d'aquesta equació es fa usant un mètode d'elements finits d'alt ordre que permet refinats de tipus h-p per a diferents geometries. Els transitoris que impliquen el moviment dels bancs de les barres de control tenen el problema conegut com l'efecte 'rod-cusping'. Estudis anteriors, en general, han abordat aquest problema utilitzant una malla fixa i definint propietats equivalents per als materials corresponents a les cel·les on es té la barra de control parcialment inserida. En el present treball es proposa l'ús d'un esquema de malla mòbil, de manera que en mallat espacial va canviant amb el moviment de la barra de control, evitant la necessitat d'utilitzar seccions eficaces equivalents per a les cel·les parcialment inserides. El funcionament de aquest esquema de malla mòbil s'estudia resolent diferents problemes tipus. La precisió obtinguda mitjançant de la teoria de la difusió en els càlculs de reactors és limitada quan es tenen elements de combustible complexos o es pretenen realitzar càlculs en malla fina. Per a millorar aquests resultats, és necessari disposar d'un mètode que incorpore aproximacions d'ordre superior de l'equació del transport de neutrons. Una possibilitat és fer ús de les equacions PN simplificades (SPN ). En aquest treball s'utilitza un mètode d'elements finits h- p per a obtenir els modes dominants associats amb una configuració donada del nucli de un reactor amb geometria hexagonal usant l'aproximació SPN . El funcionament de les aproximacions SPN (N = 1, 3, 5) s'ha estudiat per a diferents problemes de referència. / Fayez Moustafa Moawad, R. (2016). Approximation of The Neutron Diffusion Equation on Hexagonal Geometries Using a h-p finite element method [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/65353 / TESIS

<p>Predmet istraživanja u doktorskoj disertaciji je metoda za konstrukciju<br />lokalizacionih skupova za spektar i pseudospektar nelinearnih problema<br />karakterističnih korena bazirana na Geršgorinovoj teoremi i njenim<br />generalizacijama koja koristi osobine poznatih podklasa H-matrica.<br />Navedena tvrđenja i primeri rasvetljavaju odnose između navedenih<br />lokalizacionih skupova, što je posebno značajno za primenu u praksi.<br />Sadržaj ovog rada time predstavlja polaznu tačku za dublja istraživanja na<br />temu konstrukcije lokalizacionih skupova za spektar i pseudospektar<br />nelinearnih problema karakterističnih korena Geršgorinovog tipa.</p> / <p>The subject of research in the doctoral dissertation is a method for constructing<br />spectra and pseudospectra localization sets for nonlinear eigenvalue problems<br />based on Ger&scaron;gorin theorem and its generalizations, that uses the properties of<br />well-known subclasses of H-matrices. Theorems and examples given in this<br />paper are showing relations between stated localization sets, which is very<br />important for practical applications. Therefore, the content of this paper represent<br />the starting point for deeper explorations on the subject of constructing spectra<br />and pseudospectra localization sets for Ger&scaron;gorin type nonlinear eigenvalue<br />problems.</p>

In dieser Arbeit werden verschiedene Ansätze untersucht, um Tensorzerlegungsmethoden auf die Full-Configuration-Interaction-Methode (FCI) anzuwenden. Das Ziel dieser Ansätze ist es, zuverlässig konvergierende Algorithmen zu erstellen, welche es erlauben, die Wellenfunktion effizient im Canonical-Product-Tensorformat (CP) zu approximieren. Hierzu werden drei Ansätze vorgestellt, um die FCI-Wellenfunktion zu repräsentieren und darauf basierend die benötigten Koeffizienten zu bestimmen. Der erste Ansatz beruht auf einer Entwicklung der Wellenfunktion als Linearkombination von Slaterdeterminanten, bei welcher in einer Hierarchie ausgehend von der Hartree-Fock-Slaterdeterminante sukzessive besetzte Orbitale durch virtuelle Orbitale ersetzt werden. Unter Nutzung von Tensorrepräsentationen im CP wird ein lineares Gleichungssystem gelöst, um die FCI-Koeffizienten zu bestimmen. Im darauf folgenden Ansatz, welcher an Direct-CI angelehnt ist, werden Tensorrepräsentationen der Hamiltonmatrix und des Koeffizientenvektors aufgestellt, welche zur Lösung des FCI-Eigenwertproblems erforderlich sind. Hier wird ein Algorithmus vorgestellt, mit welchem das Eigenwertproblem im CP gelöst wird. In einem weiteren Ansatz wird die Repräsentation der Hamiltonmatrix und des Koeffizientenvektors im Fockraum formuliert. Dieser Ansatz erlaubt die Lösung des FCI-Eigenwertproblems mit Hilfe verschiedener Algorithmen. Diese orientieren sich an den Rayleighquotienteniterationen oder dem Davidsonalgorithmus, wobei für den ersten Algorithmus eine zweite Version entwickelt wurde, wo die Rangreduktion teilweise durch Projektionen ersetzt wurde. Für den Davidsonalgorithmus ist ein breiteres Spektrum von Molekülen behandelbar und somit können erste Untersuchungen zur Skalierung und zu den zu erwartenden Fehlern vorgestellt werden. Schließlich wird ein Ausblick auf mögliche Weiterentwicklungen gegeben, welche eine effizientere Berechnung ermöglichen und somit FCI im CP auch für größere Moleküle zugänglich macht. / In this thesis, various approaches are investigated to apply tensor decomposition methods to the Full Configuration Interaction method (FCI). The aim of these approaches is the development of algorithms, which converge reliably and which permit to approximate the wave function efficiently in the Canonical Product format (CP). Three approaches are introduced to represent the FCI wave function and to obtain the corresponding coefficients. The first approach ist based on an expansion of the wave function as a linear combination of slater determinants. In this hierarchical expansion, starting from the Hartree Fock slater determinant, the occupied orbitals are substituted by virtual orbitals. Using tensor representations in the CP, a linear system of equations is solved to obtain the FCI coefficients. In a further approach, tensor representations of the Hamiltonian matrix and the coefficient vectors are set up, which are required to solve the FCI eigenvalue problem. The tensor contractions and an algorithm to solve the eigenvalue problem in the CP are explained her in detail. In the next approach, tensor representations of the Hamiltonian matrix and the coefficient vector are constructed in the Fock space. This approach allows the application of various algorithms. They are based on the Rayleight Quotient Algorithm and the Davidson algorithm and for the first one, there exists a second version, where the rank reduction algorithm is replaced by projections. The Davidson algorithm allows to treat a broader spectrum of molecules. First investigations regarding the scaling behaviour and the expectable errors can be shown for this approach. Finally, an outlook on the further development is given, that allows for more efficient calculations and makes FCI in the CP accessible for larger molecules.