Michał Lesiuk

Reexamination of the calculation of two-center, two-electron integrals over Slater-type orbitals. III. Case study of the beryllium dimer

In this paper we present results of ab-initio calculations for the beryllium dimer with basis set of Slater-type orbitals (STOs). Nonrelativistic interaction energy of the system is determined using the frozen-core full configuration interaction calculations combined with high-level coupled cluster correction for inner-shell effects. Newly developed STOs basis sets, ranging in quality from double to sextuple zeta, are used in these computations. Principles of their construction are discussed and several atomic benchmarks are presented. Relativistic effects of order

Crossover between few and many fermions in a harmonic trap

{\alpha}^2

Reexamination of the calculation of two-center, two-electron integrals over Slater-type orbitals. I. Coulomb and hybrid integrals.

are calculated perturbatively by using the Breit-Pauli Hamiltonian and are found to be significant. We also estimate the leading-order QED effects. Influence of the adiabatic correction is found to be negligible. Finally, the interaction energy of the beryllium dimer is determined to be 929.0

Many interacting fermions in a one-dimensional harmonic trap: a quantum-chemical treatment

\,\pm\,

Reexamination of the calculation of two-center, two-electron integrals over Slater-type orbitals. II. Neumann expansion of the exchange integrals.

1.9

A first-principles-based correlation functional for harmonious connection of short-range correlation and long-range dispersion

cm^{-1}

Molecular electrostatic potential at the atomic sites in the effective core potential approximation.

, in a very good agreement with the recent experimental value. The results presented here appear to be the most accurate ab-initio calculations for the beryllium dimer available in the literature up to date and probably also one of the most accurate calculations for molecular systems containing more than four electrons.

Analytical two-center integrals over Slater geminal functions

The properties of a balanced two-component Fermi gas in a one-dimensional harmonic trap are studied by means of the coupled-cluster method. For few fermions we recover the results of exact diagonalization, yet with this method we are able to study much larger systems. We compute the energy, the chemical potential, the pairing gap, and the density profile of the trapped clouds, smoothly mapping the crossover between the few-body and many-body limits. The energy is found to converge surprisingly rapidly to the many-body result for every value of the interaction strength. Many more particles are instead needed to give rise to the nonanalytic behavior of the pairing gap, and to smoothen the pronounced even-odd oscillations of the chemical potential induced by the shell structure of the trap.

Combining Slater-type orbitals and effective core potentials

In this paper, which constitutes the first part of the series, we consider calculation of two-center Coulomb and hybrid integrals over Slater-type orbitals. General formulas for these integrals are derived with no restrictions on the values of the quantum numbers and nonlinear parameters. Direct integration over the coordinates of one of the electrons leaves us with the set of overlaplike integrals which are evaluated by using two distinct methods. The first one is based on the transformation to the ellipsoidal coordinates system and the second utilizes a recursive scheme for consecutive increase of the angular momenta in the integrand. In both methods simple one-dimensional numerical integrations are used in order to avoid severe digital erosion connected with the straightforward use of the alternative analytical formulas. It is discussed that the numerical integration does not introduce a large computational overhead since the integrands are well-behaved functions, calculated recursively with decent speed. Special attention is paid to the numerical stability of the algorithms. Applicability of the resulting scheme over a large range of the nonlinear parameters is tested on examples of the most difficult integrals appearing in the actual calculations including, at most, 7i-type functions (l=6).

Exploring Chemical Space with Alchemical Derivatives: BN-Simultaneous Substitution Patterns in C60

We employ \textit{ab initio} methods of quantum chemistry to investigate spin-1/2 fermions interacting via a two-body contact potential in a one-dimensional harmonic trap. The convergence of the total energy with the size of the one-particle basis set is analytically investigated for the two-body problem and the same form of the convergence formula is numerically confirmed to be valid for the many-body case. Benchmark calculations for two to six fermions with the full configuration interaction method equivalent to the exact diagonalization approach, and the coupled cluster method including single, double, triple, and quadruple excitations are presented. The convergence of the correlation energy with the level of excitations included in the coupled cluster model is analyzed. The range of the interaction strength for which single-reference coupled cluster methods work is examined. Next, the coupled cluster method restricted to single, double, and noniterative triple excitations, CCSD(T), is employed to study a two-component Fermi gas composed of 6 to 80 atoms in a one-dimensional harmonic trap. The density profiles of trapped atomic clouds are also reported. Finally, a comparison with experimental results for few-fermion systems is presented. Upcoming possible applications and extensions of the presented approach are discussed.