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.

Implementation of Fixed-Nuclei Polyatomic MCTDHF Capability and the Future with Nuclear Motion

Multiconfigurational time-dependent density functional theory for atomic nuclei: Technical and numerical aspects

Nuclear-Electronic Orbital Time-Dependent Configuration Interaction Method

A multi-layer multi-configurational time-dependent Hartree approach to lattice models beyond one dimension

Adaptive Multi-Resolution 3D Hartree-Fock-Bogoliubov Solver for Nuclear Structure

Multi-layer Gaussian-based multi-configuration time-dependent Hartree (ML-GMCTDH) simulations of ultrafast charge separation in a donor-acceptor complex

Multiconfiguration Dirac-Hartree-Fock and configuration-interaction study of $4d-3p$ x-ray transitions in Cu- and Ni-like tungsten ions

Lecture Notes: many-body quantum dynamics with MCTDH-X

Efficient time-dependent vibrational coupled cluster computations with time-dependent basis sets at the two-mode coupling level: full and hybrid TDMVCC[2]

Molecular hyperdynamics coupled with the nonorthogonal tight-binding approach: Implementation and validation

Tunable-fidelity wave functions for the ab initio description of scattering and reactions

An implementation of nuclear many-body wave functions by the superposition of localized Gaussians

Density fitting in periodic systems: application to TDHF in diamond and oxides

Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) Approach to the Correlated Exciton-Vibrational Dynamics in the FMO Complex

Multicomponent Density-Functional Theory for Electrons and Nuclei

Review of the finite difference Hartree-Fock method for atoms and diatomic molecules, and its implementation in the x2dhf program

The TDHF code Sky3D version 1.2

A Coupled Cluster Approach with a Hybrid Treatment of Connected Triple Excitations: Implementation and Applications for Open-Shell Systems.

Covariant density functional theory for nuclear fission based on a two-center harmonic oscillator basis

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm-Dancoff approximation