Source author record

Jürgen Gauss

Jürgen Gauss appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

10 published item(s)

preprint2026arXiv

The Direct-Product Decomposition Approach for Symmetry Exploitation in Many-Body Methods in Case of Non-Abelian Point Groups

We demonstrate for the specific case of $C_{3v}$ how the direct-product decomposition scheme for the treatment of symmetry in coupled-cluster (CC) calculations can be extended to non-Abelian point groups. We show that for the two-electron integrals and CC amplitudes a block structure can be obtained by resolving the reducible products of two irreducible representations into their irreducible representations. To deal with the necessary resorts of the ordering of the two-electron integrals and amplitudes, spin-adaptation, and the O(M$^5$) contractions (with M as the number of basis functions) of a CC calculation, we suggest a strategy that uses both the reduced and non-reduced representation of the corresponding quantities and switches back and forth between them. While the reduced representations are the ones used in the O(M$^6$) contractions, the other steps are better carried out in the non-reduced representation. Our pilot implementation of the CC singles and doubles method confirms in test calculations for NH$_3$ and PH$_3$ using different basis sets that significant savings (of more than 20 compared to treatments without symmetry and about 5 compared to treatments using $C_s$ symmetry) are possible and suggest that the exploitation of non-Abelian symmetry would render CC computations on large highly symmetric molecules possible

preprint2022arXiv

Computation of NMR shieldings at the CASSCF level using gauge-including atomic orbitals and Cholesky decomposition

We present an implementation of coupled-perturbed complete active space self-consistent field (CP-CASSCF) theory for the computation of nuclear magnetic resonance chemical shifts using gauge-including atomic orbitals and Cholesky decomposed two-electron integrals. The CP-CASSCF equations are solved using a direct algorithm where the magnetic Hessian matrix-vector product is expressed in terms of one-index transformed quantities. Numerical tests on systems with up to about 1300 basis functions provide information regarding both the computational efficiency and limitations of our implementation.

preprint2021arXiv

A black-box, general purpose quadratic self-consistent field code with and without Cholesky Decomposition of the two-electron integrals

We present the implementation of a quadratically convergent Self-consistent field (QCSCF) algorithm based on an adaptive trust-radius optimization scheme for restricted open-shell Hartree-Fock (ROHF), restricted Hartree-Fock (RHF), and unrestricted Hartree-Fock (UHF) references. The algorithm can exploit Cholesky decomposition (CD) of the two-electron integrals to allow calculations on larger systems. The most important feature of the QCSCF code lies in its black-box nature -- probably the most important quality desired by a generic user. As shown for pilot applications, it does not require one to tune the self-consistent field (SCF) parameters (damping, Pulay's DIIS, and other similar techniques) in difficult-to-converge molecules. Also, it can be used to obtain a very thigh convergence with extended basis set - a situation often needed when computing high-order molecular properties - where the standard SCF algorithm starts to oscillate. Nevertheless, trouble may appear even with a QCSCF solver. In this respect, we discuss what can go wrong, focusing on the multiple UHF solutions of ortho-benzyne

preprint2021arXiv

Complex Ground-State and Excitation Energies in Coupled-Cluster Theory

Since in coupled-cluster (CC) theory ground-state and excitation energies are eigenvalues of a non-Hermitian matrix, these energies can in principle take on complex values. In this paper we discuss the appearance of complex energy values in CC calculations from a mathematical perspective. We analyze the behaviour of the eigenvalues of Hermitian matrices that are perturbed (in a non-Hermitian manner) by a real parameter. Based on these results we show that for CC calculations with real-valued Hamiltonian matrices the ground-state energy generally takes a real value. Furthermore, we show that in the case of real-valued Hamiltonian matrices complex excitation energies only occur in the context of conical intersections. In such a case, unphysical consequences are encountered such as a wrong dimension of the intersection seam, large numerical deviations from full configuration-interaction (FCI) results, and the square-root-like behaviour of the potential surfaces near the conical intersection. In the case of CC calculations with complex-valued Hamiltonian matrix elements, it turns out that complex energy values are to be expected for ground and excited states when no symmetry is present. We confirm the occurrence of complex energies by sample calculations using a six-state model and by CC calculations for the H2O molecule in a strong magnetic field. We furthermore show that symmetry can prevent the occurrence of complex energy values. Lastly, we demonstrate that in most cases the real part of the complex energy values provides a very good approximation to the FCI energy.

preprint2021arXiv

Incremental Treatments of the Full Configuration Interaction Problem

The recent many-body expanded full configuration interaction (MBE-FCI) method is reviewed by critically assessing its advantages and drawbacks in the context of contemporary near-exact electronic structure theory. Besides providing a succinct summary of the history of MBE-FCI to date within a generalized and unified theoretical setting, its finer algorithmic details are discussed alongside our optimized computational implementation of the theory. A selected few of the most recent applications of MBE-FCI are revisited, before we close by outlining its future research directions as well as its place among modern near-exact wave function-based methods.

preprint2020arXiv

Determination of accurate rest frequencies and hyperfine structure parameters of cyanobutadiyne, HC$_5$N

Very accurate transition frequencies of HC$_5$N were determined between 5.3 and 21.4 GHz with a Fourier transform microwave spectrometer. The molecules were generated by passing a mixture of HC$_3$N and C$_2$H$_2$ highly diluted in neon through a discharge valve followed by supersonic expansion into the Fabry-Perot cavity of the spectrometer. The accuracies of the data permitted us to improve the experimental $^{14}$N nuclear quadrupole coupling parameter considerably and the first experimental determination of the $^{14}$N nuclear spin-rotation parameter. The transition frequencies are also well suited to determine in astronomical observations the local speed of rest velocities in molecular clouds with high fidelity. The same setup was used to study HC$_7$N, albeit with modest improvement of the experimental $^{14}$N nuclear quadrupole coupling parameter. Quantum chemical calculations were carried out to determine $^{14}$N nuclear quadrupole and spin-rotation coupling parameters of HC$_5$N, HC$_7$N, and related molecules. These calculations included evaluation of vibrational and relativistic corrections to the non-relativistic equilibrium quadrupole coupling parameters; their considerations improved the agreement between calculated and experimental values substantially.

preprint2016arXiv

Assessment of the accuracy of coupled cluster perturbation theory for open-shell systems. I. Triples expansions

The accuracy at which total energies of open-shell atoms and organic radicals may be calculated is assessed for selected coupled cluster perturbative triples expansions, all of which augment the coupled cluster singles and doubles (CCSD) energy by a non-iterative correction for the effect of triple excitations. Namely, the second- through sixth-order models of the recently proposed CCSD(T-n) triples series [J. Chem. Phys. 140, 064108 (2014)] are compared to the acclaimed CCSD(T) model for both unrestricted as well as restricted open-shell Hartree-Fock (UHF/ROHF) reference determinants. By comparing UHF- and ROHF-based statistical results for a test set of 18 modest-sized open-shell species with comparable RHF-based results, no behavioral differences are observed for the higher-order models of the CCSD(T-n) series in their correlated descriptions of closed- and open-shell species. In particular, we find that the convergence rate throughout the series towards the coupled cluster singles, doubles, and triples (CCSDT) solution is identical for the two cases. For the CCSD(T) model, on the other hand, not only its numerical consistency, but also its established, yet fortuitous cancellation of errors breaks down in the transition from closed- to open-shell systems. The higher-order CCSD(T-n) models (orders n>3) thus offer a consistent and significant improvement in accuracy relative to CCSDT over the CCSD(T) model, equally for RHF, UHF, and ROHF reference determinants, albeit at an increased computational cost.

preprint2016arXiv

Assessment of the accuracy of coupled cluster perturbation theory for open-shell systems. II. Quadruples expansions

We extend our assessment of the potential of perturbative coupled cluster (CC) expansions for a test set of open-shell atoms and organic radicals to the description of quadruple excitations. Namely, the second- through sixth-order models of the recently proposed CCSDT(Q-n) quadruples series [J. Chem. Phys. 140, 064108 (2014)] are compared to the prominent CCSDT(Q) and lambda-CCSDT(Q) models. From a comparison of the models in terms of their recovery of total CC singles, doubles, triples, and quadruples (CCSDTQ) energies, we find that the performance of the CCSDT(Q-n) models is independent of the reference used (unrestricted or restricted (open-shell) Hartree-Fock), in contrast to the CCSDT(Q) and lambda-CCSDT(Q) models, for which the accuracy is strongly dependent on the spin of the molecular ground state. By further comparing the ability of the models to recover relative CCSDTQ atomization energies, the discrepancy between them is found to be even more pronounced, stressing how a balanced description of both closed- and open-shell species - as found in the CCSDT(Q-n) models - is indeed of paramount importance if any perturbative CC model is to be of chemical relevance for high-accuracy applications. In particular, the third-order CCSDT(Q-3) model is found to offer an encouraging alternative to the existing choices of quadruples models used in modern computational thermochemistry, since the model is still only of moderate cost, albeit markedly more costly than, e.g., the CCSDT(Q) and lambda-CCSDT(Q) models.

preprint2015arXiv

Coupled-cluster theory for atoms and molecules in strong magnetic fields

An implementation of coupled-cluster (CC) theory to treat atoms and molecules in finite magnetic fields is presented. The main challenges stem from the magnetic-field dependence in the Hamiltonian, or, more precisely, the appearance of the angular momentum operator, due to which the wave function becomes complex and which introduces a gauge-origin dependence. For this reason, an implementation of a complex CC code is required together with the use of gauge-including atomic orbitals to ensure gauge-origin independence. Results of coupled-cluster singles--doubles--perturbative-triples (CCSD(T)) calculations are presented for atoms and molecules with a focus on the dependence of correlation and binding energies on the magnetic field.

preprint2009arXiv

Hyperfine structure in the J = 1-0 transitions of DCO+, DNC, and HN13C: astronomical observations and quantum-chemical calculations

We have observed the rotational ground-state (J = 1-0) transitions of DCO+, HN13C and DNC with the IRAM 30m telescope toward the dark cloud LDN 1512 which has exceptionally narrow lines permitting hyperfine splitting to be resolved in part. The measured splittings of 50-300 kHz are used to derive nuclear quadrupole and spin-rotation parameters for these species. The measurements are supplemented by high-level quantum-chemical calculations using coupled-cluster techniques and large atomic-orbital basis sets. We find eQq = +151.12 (400) kHz and C_I = -1.12 (43) kHz for DCO+, eQq = 272.5 (51) kHz for HN13C, and eQq(D) = 265.9 (83) kHz and eQq(N) = 288.2 (71) kHz for DNC. The numbers for DNC are consistent with previous laboratory data, while our constants for DCO+ are somewhat smaller than previous results based on astronomical data. For both DCO+ and DNC, our results are more accurate than previous determinations. Our results are in good agreement with the corresponding best theoretical estimates. We also derive updated rotational constants for HN13C: B = 43545.6000 (47) MHz and D = 93.7 (20) kHz. The hyperfine splittings of the DCO+, DNC and HN13C J = 1-0 lines range over 0.47-1.28 km/s, which is comparable to typical line widths in pre-stellar cores and to systematic gas motions on ~1000 AU scales in protostellar cores. We present tabular information to allow inclusion of the hyperfine splitting in astronomical data interpretation. The large differences in the 14N quadrupole parameters of DNC and HN13C have been traced to differences in the vibrational corrections caused by significant non-rigidity of these molecules, particularly along the bending coordinate.