Nuclear Hartree-Fock calculations on parallel computers

Hartree-Fock (HF) calculations have had remarkable success in describing large nuclei at high spin, temperature and deformation. To allow full range of possible deformations, the Skyrme HF equations can be discretized on a three-dimensional mesh. However, such calculations are currently limited by the computational resources provided by traditional supercomputers. To take advantage of recent developments in massively parallel computing technology, we have implemented the LLNL Skyrme-force static and rotational HF codes on Intel's DELTA and GAMMA systems at Caltech.

We decomposed the HF code by assigning a portion of the mesh to each node, with nearest neighbor meshes assigned to nodes connected by communication· channels. This kind of decomposition is well-suited for the DELTA and the GAMMA architecture because the only non-local operations are wave function orthogonalization and the boundary conditions of the Poisson equation for the Coulomb field.

Our first application of the HF code on parallel computers has been the study of identical superdeformed (SD) rotational bands in the Hg region. In the last ten years, many SD rotational bands have been found experimentally. One very surprising feature found in these SD rotational bands is that many pairs of bands in nuclei that differ by one or two mass units have nearly identical deexcitation gamma-ray energies. Our calculations of the five rotational bands in ^(192)Hg and ^(194)Pb show that the filling of specific orbitals can lead to bands with deexcitation gamma-ray energies differing by at most 2 keV in nuclei differing by two mass units and over a range of angular momenta comparable to that observed experimentally. Our calculations of SD rotational bands in the Dy region also show that twinning can be achieved by filling or emptying some specific orbitals.

The interpretation of future precise experiments on atomic parity nonconservation (PNC) in terms of parameters of the Standard Model could be hampered by uncertainties in the atomic and nuclear structure. As a further application of the massively parallel HF calculations, we calculated the proton and neutron densities of the Cesium isotopes from A = 125 to A = 139. Based on our good agreement with experimental charge radii, binding energies, and ground state spins, we conclude that the uncertainties in the ratios of weak charges are less than 10^(-3), comfortably smaller than the anticipated experimental error.

Hyperfine interactions and nuclear structure

An automatic experimental apparatus for perturbed angular correlation measurements, capable of incorporating Ge(Li) detectors as well as scintillation counters, has been constructed.

The gamma-gamma perturbed angular correlation technique has been used to measure magnetic dipole moments of several nuclear excited states in the osmium transition region. In addition, the hyperfine magnetic fields, experienced by nuclei of 'impurity' atoms embedded in ferromagnetic host lattices, have been determined for several '4d' and '5d' impurity atoms.

The following magnetic dipole moments were obtained in the osmium transition region μ_2⁺(¹⁹⁰Os) = 0.54 ± 0.06 nm μ_4⁺(¹⁹⁰Os) = 0.88 ± 0.48 nm μ_2⁺(¹⁹²Os) = 0.56 ± 0.08 nm μ_2⁺(¹⁹²Pt) = 0.56 ± 0.06 nm μ_2^+’(¹⁹²Pt) = 0.62 ± 0.14 nm.

These results are discussed in terms of three collective nuclear models; the cranking model, the rotation-vibration model and the pairing-plus-quadrupole model. The measurements are found to be in satisfactory agreement with collective descriptions of low lying nuclear states in this region.

The following hyperfine magnetic fields of 'impurities' in ferromagnetic hosts were determined; H_int(Cd Ni) = - (64.0 ± 0.8)kG H_int(Hg Fe) = - (440 ± 105)kG H_int(Hg Co) = - (370 ± 78)kG H_int(Hg Ni) = - (86 ± 22)kG H_int(Tl Fe) = - (185 ± 70)kG H_int(Tl Co) = - (90 ± 35)kG H_int(Ra Fe) = - (105 ± 20)kG H_int(Ra Co) = - (80 ± 16)kG H_int(Ra Ni) = - (30 ± 10)kG, where in H_int(AB); A is the impurity atom embedded in the host lattice B. No quantitative theory is available for comparison. However, these results are found to obey the general systematics displayed by these fields. Several mechanisms which may be responsible for the appearance of these fields are mentioned.

Finally, a theoretical expression for time-differential perturbed angular correlation measurement, which duplicates experimental conditions is developed and its importance in data analysis is discussed.

Chemical and electrochemical investigations of cobalt cyanide and ruthenium ammine complexes

Part I.

The stoichiometry and kinetics of the reaction between Co(CN₅H^3- and HgX₂ (X = CN, OH) have been investigated. The products of the reaction are two new complexes, [(NC)₅Co-HgX]^3- and [(NC)₅Co-Hg-Co(CN)₅]^6-, whose spectra are reported. The kinetic measurements produced a value for the forward rate constant of the reaction Co(CN)₅H^3- + OH^- k₁/k_-1 Co(CN)₅^4- +H₂O, k1 = (9.7 ± 0.8) x 10^-2 M^-1 sec^-1 at 24°C, and an equilibrium constant for the reaction K = 10^-6 M^-1.

Part II.

Unusually large and sharp "adsorption waves" appear in cyclic voltammograms of Co(CN)₅^3- and several cobalt(III) pentacyano complexes at stationary mercury electrodes. The nature of the adsorbed species and the reasons for the absence of the adsorption waves in polarograms taken with a d.m.e. have been examined. The data are compatible with the adsorption, in all cases, of a coordinatively unsaturated cobalt(II) complex, Co(CN)₄^2-, by means of a cobalt-mercury bond. When the resulting adsorbed complex is reduced, a series of subsequent chemical and electrode reactions is initiated in which three faradays of charge are consumed for each mole of adsorbed complex. The adsorption of the anionic complex strongly retards the reduction of other negatively charged complexes.

Part III.

A number of formal redox potentials for Ru^III (NH₃)₅L + e = Ru^II (NH₃)₅L and Ru^III(NH₃)₄L₂ + e = Ru^II (NH₃)₄L₂ (where L is various ligands) has been measured by cyclic voltammetry, potentiometry, and polarography and are discussed in terms of the properties of the ligands, such as π-accepting capability. Reduction of coordinated pyrazine in the complexes, Ru(NH₃)₅ Pz²⁺, cis- and trans-Ru(NH₃)₄Pz₂²⁺, on a mercury electrode has been observed. The behavior of this reduction in various acidity of the solution as well as the reoxidation of the reduction products are discussed.