Abstract:The Sky3D code has been widely used to describe nuclear ground states, collective vibrational excitations, and heavy-ion collisions. The approach is based on Skyrme forces or related energy density functionals. The static and dynamic equations are solved on a three-dimensional grid, and pairing is been implemented in the BCS approximation. This updated version of the code aims to facilitate the calculation of nuclear strength functions in the regime of linear response theory, while retaining all existing functionality and use cases. The strength functions are benchmarked against available RPA codes, and the user has the freedom of choice when selecting the nature of external excitation (from monopole to hexadecapole and more). Some utility programs are also provided that calculate the strength function from the time-dependent output of the dynamic calculations of the Sky3D code. New version program summary Program Title: Sky3D CPC Library link to program files: https://doi.org/10.17632/vzbrzvyrn4.2 Developer's repository link: https://github.com/manybody/sky3d Licensing provisions: GPLv3 Programming language: Fortran, with one post-processing utility in Python. Journal reference of previous version:: Schuetrumpf, B., Reinhard, P.G., Stevenson, P.D., Umar, A.S., and Maruhn, J.A. (2018). The TDHF code Sky3D version 1.1. Comput. Phys. Commun. 229 (2018) 211–213. Does the new version supersede the previous version?: Yes. Reasons for the new version: The capability of reproducing the nuclear strength function for a variety of newly-coded external boosts has been added. Nature of problem: Calculating nuclear multipole strength functions is a crucial probe that can help model the nuclear system and its structure properties. A variety of models exist for this task, such as QRPA (Quasiparticle Random Phase Approximation) and its variants, but such approaches are often limited due to symmetry constraints. Time-dependent Hartree Fock (TDHF) has been used to simulate nuclear vibrations and collisions between nuclei for low energies without assuming any symmetry in the system. This code extends the TDHF to calculate the multi-pole strength functions of atomic nuclei. We showcase its reliability by comparing it with the established RPA codes for the calculation of such strength functions. Solution method: We extended previous versions of the Sky3D code [1] , [2] to include an external boost of multipole type where the user can provide custom input that decides the nature of the multipole (monopole, quadrupole, octupole, and so on) boost. The principal aim is to calculate the multipole strength function, which is the Fourier transform of the time-dependent expectation value of the multipole operator, which has the same form as the external boost. The proper unit conversion is done so that one can extract the exact unit of the thus calculated strength function, which is comparable to the available literature in the field. The boundary conditions are chosen such that a Woods-Saxon-like function cuts off the external field, driving it to zero at the boundary. Summary of revisions: The ability to give a more general external field for excitation has been implemented, which can be of a chosen multipolarity and isospin nature. Then, as a function of time, moments for both isospins and for multipolarities up to L=5 are followed. A new analysis program is included to calculate strength functions and energy weighted sum rules. The ability to specify a mixing parameter in the pairing force is included. This allows a mixed surface-volume pairing to be selected in a way more compatable with QRPA codes for comparison. An updated Makefile is included which makes compilation on Apple silicon computers easier.

Conditional Aharonov-Bohm Phases with Double Quantum Dots

Secondly quantized multi-configurational approach for atomic databases

A program for simplifying summation of Wigner 3 j -symbols

Molecular structure calculations: a unified quantum mechanical description of electrons and nuclei using explicitly correlated Gaussian functions and the global vector representation

LibSC: Library for Scaling Correction Methods in Density Functional Theory

Atomic basis functions for molecular electronic structure calculations

High-order finite element method for atomic structure calculations

JRAF: A Julia Package for Computation of the Relativistic Molecular Auxiliary Functions

Relativistic quantum theory and algorithms: a toolbox for modeling many-fermion systems in different scenarios

High Accurate Computational Method and Program for Approximate Relativistic and Non-relativistic Density Functional Calculations

Libcint: an Efficient General Integral Library for Gaussian Basis Functions

Developing Orbital-Dependent Corrections for the Non-Additive Kinetic Energy in Subsystem Density Functional Theory

Ab initio Calculations in Atoms, Molecules, and Solids, Treating Spin-Orbit Coupling and Electron Interaction on Equal Footing

The LISE Package: Solvers for Static and Time-Dependent Superfluid Local Density Approximation Equations in Three Dimensions.

A Sub-Linear Scaling Algorithm for Computing the Electronic Structure of Materials

eQE 2.0: Subsystem DFT beyond GGA functionals

The TDHF code Sky3D version 1.2

Instruction of my personal computing library

On matrix elements of the vector physical quantities and methods of angular momentum in quantum mechanics

The Alexandria library, a quantum-chemical database of molecular properties for force field development

ZKCM: A C++ library for multiprecision matrix computation with applications in quantum information