Low-lying states in near-magic odd-odd nuclei and the effective interaction

The iterative quasi-particle-random-phase approximation (QRPA) method we previously developed to accurately calculate properties of individual nuclear states is extended so that it can be applied for nuclei with odd numbers of neutrons and protons. The approach is based on the proton-neutron-QRPA (pnQRPA) and uses an iterative non-hermitian Arnoldi diagonalization method where the QRPA matrix does not have to be explicitly calculated and stored. The method is used to calculate excitation energies of proton-neutron multiplets for several nuclei. The influence of a pairing interaction in the $T=0$ channel is studied

Effective pseudopotential for energy density functionals with higher order derivatives

We derive a zero-range pseudopotential that includes all possible terms up to sixth order in derivatives. Within the Hartree-Fock approximation, it gives the average energy that corresponds to a quasi-local nuclear Energy Density Functional (EDF) built of derivatives of the one-body density matrix up to sixth order. The direct reference of the EDF to the pseudopotential acts as a constraint that divides the number of independent coupling constants of the EDF by two. This allows, e.g., for expressing the isovector part of the functional in terms of the isoscalar part, or vice versa. We also derive the analogous set of constraints for the coupling constants of the EDF that is restricted by spherical, space-inversion, and time-reversal symmetries.Comment: 18 LaTeX pages, 2 EPS Figures, 27 Tables, and 18 files of the supplemental material (LaTeX, Mathematica, and Fortran), introduction rewritten, table XXVII and figure 2 corrected, in press in Physical Review

Convergence of density-matrix expansions for nuclear interactions

We extend density-matrix expansions in nuclei to higher orders in derivatives of densities and test their convergence properties. The expansions allow for converting the interaction energies characteristic to finite- and short-range nuclear effective forces into quasi-local density functionals. We also propose a new type of expansion that has excellent convergence properties when benchmarked against the binding energies obtained for the Gogny interaction.Comment: 4 pages, 3 figure

Fluctuating parts of nuclear ground state correlation energies

Background: Heavy atomic nuclei are often described using the Hartree-Fock-Bogoliubov (HFB) method. In principle, this approach takes into account Pauli effects and pairing correlations while other correlation effects are mimicked through the use of effective density-dependent interactions. Purpose: Investigate the influence of higher order correlation effects on nuclear binding energies using Skyrme's effective interaction. Methods: A cut-off in relative momenta is introduced in order to remove ultraviolet divergences caused by the zero-range character of the interaction. Corrections to binding energies are then calculated using the quasiparticle-random-phase approximation (QRPA) and second order many-body perturbation theory (MBPT2). Result: Contributions to the correlation energies are evaluated for several isotopic chains and an attempt is made to disentangle which parts give rise to fluctuations that may be difficult to incorporate on the HFB level. The dependence of the results on the cut-off is also investigated. Conclusions: The improved interaction allows explicit summations of perturbation series which is useful for the description of some nuclear observables. However, refits of the interaction parameters are needed to obtain more quantitative results

Parallel Mapper

The construction of Mapper has emerged in the last decade as a powerful and effective topological data analysis tool that approximates and generalizes other topological summaries, such as the Reeb graph, the contour tree, split, and joint trees. In this paper, we study the parallel analysis of the construction of Mapper. We give a provably correct parallel algorithm to execute Mapper on multiple processors and discuss the performance results that compare our approach to a reference sequential Mapper implementation. We report the performance experiments that demonstrate the efficiency of our method

Hard X-ray and UV Observations of the 2005 January 15 Two-ribbon Flare

In this paper, we present comprehensive analysis of a two-ribbon flare observed in UV 1600{\AA} by Transition Region and Coronal Explorer and in HXRs by Reuven Ramaty High Energy Solar Spectroscopic Imager. HXR (25-100 keV) imaging observations show two kernels of size (FWHM) 15?? moving along the two UV ribbons. We find the following results. (1) UV brightening is substantially enhanced wherever and whenever the compact HXR kernel is passing, and during the HXR transit across a certain region, the UV count light curve in that region is temporally correlated with the HXR total flux light curve. After the passage of the HXR kernel, the UV light curve exhibits smooth monotonical decay. (2)We measure the apparent motion speed of the HXR sources and UV ribbon fronts, and decompose the motion into parallel and perpendicular motions with respect to the magnetic polarity inversion line (PIL). It is found that HXR kernels and UV fronts exhibit similar apparent motion patterns and speeds. The parallel motion dominates during the rise of the HXR emission, and the perpendicular motion starts and dominates at the HXR peak, the apparent motion speed being 10-40 km s-1. (3) We also find that UV emission is characterized by a rapid rise correlated with HXRs, followed by a long decay on timescales of 15-30 minutes. The above analysis provides evidence that UV brightening is primarily caused by beam heating, which also produces thick-target HXR emission. The thermal origin of UV emission cannot be excluded, but would produce weaker heating by one order of magnitude. The extended UV ribbons in this event are most likely a result of sequential reconnection along the PIL, which produces individual flux tubes (post-flare loops), subsequent non-thermal energy release and heating in these flux tubes, and then the very long cooling time of the transition region at the feet of these flux tubes.Comment: 8 figure

Modeling of Covalent Bonding in Solids by Inversion of Cohesive Energy Curves

We provide a systematic test of empirical theories of covalent bonding in solids using an exact procedure to invert ab initio cohesive energy curves. By considering multiple structures of the same material, it is possible for the first time to test competing angular functions, expose inconsistencies in the basic assumption of a cluster expansion, and extract general features of covalent bonding. We test our methods on silicon, and provide the direct evidence that the Tersoff-type bond order formalism correctly describes coordination dependence. For bond-bending forces, we obtain skewed angular functions that favor small angles, unlike existing models. As a proof-of-principle demonstration, we derive a Si interatomic potential which exhibits comparable accuracy to existing models.Comment: 4 pages revtex (twocolumn, psfig), 3 figures. Title and some wording (but no content) changed since original submission on 24 April 199