Adaptive Eigenvalue Computation - Complexity Estimates

This paper is concerned with the design and analysis of a fully adaptive eigenvalue solver for linear symmetric operators. After transforming the original problem into an equivalent one formulated on $\ell_2$, the space of square summable sequences, the problem becomes sufficiently well conditioned so that a gradient type iteration can be shown to reduce the error by some fixed factor per step. It then remains to realize these (ideal) iterations within suitable dynamically updated error tolerances. It is shown under which circumstances the adaptive scheme exhibits in some sense asymptotically optimal complexity.Comment: submitted to Math. Com

Domain decomposition and locality optimization for large-scale lattice Boltzmann simulations

We present a simple, parallel and distributed algorithm for setting up and partitioning a sparse representation of a regular discretized simulation domain. This method is scalable for a large number of processes even for complex geometries and ensures load balance between the domains, reasonable communication interfaces, and good data locality within the domain. Applying this scheme to a list-based lattice Boltzmann flow solver can achieve similar or even higher flow solver performance than widely used standard graph partition based tools such as METIS and PT-SCOTCH

Gamma-widths, lifetimes and fluctuations in the nuclear quasi-continuum

Statistical $\gamma$-decay from highly excited states is determined by the nuclear level density (NLD) and the $\gamma$-ray strength function ($\gamma$SF). These average quantities have been measured for several nuclei using the Oslo method. For the first time, we exploit the NLD and $\gamma$SF to evaluate the $\gamma$-width in the energy region below the neutron binding energy, often called the quasi-continuum region. The lifetimes of states in the quasi-continuum are important benchmarks for a theoretical description of nuclear structure and dynamics at high temperature. The lifetimes may also have impact on reaction rates for the rapid neutron-capture process, now demonstrated to take place in neutron star mergers.Comment: CGS16, Shanghai 2017, Proceedings, 5 pages, 3 figure

Statistical properties of 243^{243}Pu, and 242^{242}Pu(n,γ\gamma) cross section calculation

The level density and gamma-ray strength function (gammaSF) of 243Pu have been measured in the quasi-continuum using the Oslo method. Excited states in 243Pu were populated using the 242Pu(d,p) reaction. The level density closely follows the constant-temperature level density formula for excitation energies above the pairing gap. The gammaSF displays a double-humped resonance at low energy as also seen in previous investigations of actinide isotopes. The structure is interpreted as the scissors resonance and has a centroid of omega_{SR}=2.42(5)MeV and a total strength of B_{SR}=10.1(15)mu_N^2, which is in excellent agreement with sum-rule estimates. The measured level density and gammaSF were used to calculate the 242Pu(n,gamma) cross section in a neutron energy range for which there were previously no measured data.Comment: 9 pages, 8 figure

Statistical properties of the well deformed 153,155^{153,155}Sm nuclei and the scissors resonance

The Nuclear Level Densities (NLDs) and the $\gamma$-ray Strength Functions ($\gamma$SFs) of $^{153,155}$Sm have been extracted from (d,p$\gamma$) coincidences using the Oslo method. The experimental NLD of $^{153}$Sm is higher than the NLD of $^{155}$Sm, in accordance with microscopic calculations. The $\gamma$SFs of $^{153,155}$Sm are in fair agreement with QRPA calculations based on the D1M Gogny interaction. An enhancement is observed in the $\gamma$SF for both $^{153,155}$Sm nuclei around 3 MeV in excitation energy and is attributed to the M1 Scissors Resonance (SR). Their integrated strengths were found to be in the range 1.3 - 2.1 and 4.4 - 6.4 $\mu^{2}_{N}$ for $^{153}$Sm and $^{155}$Sm, respectively. The strength of the SR for $^{155}$Sm is comparable to those for deformed even-even Sm isotopes from nuclear resonance fluorescence measurements, while that of $^{153}$Sm is lower than expected

Nuclear Level Density and γ\gamma-ray Strength Function of 67Ni^{67}\mathrm{Ni} and the impact on the i-process

Proton-$\gamma$ coincidences from $(\mathrm{d},\mathrm{p})$ reactions between a $^{66}\mathrm{Ni}$ beam and a deuterated polyethylene target have been analyzed with the inverse Oslo method to find the nuclear level density (NLD) and $\gamma$-ray strength function ($\gamma$SF) of $^{67}\mathrm{Ni}$. The $^{66}\mathrm{Ni}(\mathrm{n},\gamma)$ capture cross section has been calculated using the Hauser-Feshbach model in TALYS using the measured NLD and $\gamma$SF as constraints. We confirm that $^{66}\mathrm{Ni}(\mathrm{n},\gamma)$ acts as a bottleneck when relying on one-zone nucleosynthesis calculations. However, we find that the impact of this reaction is strongly damped in multi-zone low-metallicity AGB stellar models experiencing i-process nucleosynthesis.Comment: Submitted to Phys. Rev.