Coordinate-Space Hartree-Fock-Bogoliubov Solvers for Superfluid Fermi Systems in Large Boxes

The self-consistent Hartree-Fock-Bogoliubov problem in large boxes can be solved accurately in the coordinate space with the recently developed solvers HFB-AX (2D) and MADNESS-HFB (3D). This is essential for the description of superfluid Fermi systems with complicated topologies and significant spatial extend, such as fissioning nuclei, weakly-bound nuclei, nuclear matter in the neutron star rust, and ultracold Fermi atoms in elongated traps. The HFB-AX solver based on B-spline techniques uses a hybrid MPI and OpenMP programming model for parallel computation for distributed parallel computation, within a node multi-threaded LAPACK and BLAS libraries are used to further enable parallel calculations of large eigensystems. The MADNESS-HFB solver uses a novel multi-resolution analysis based adaptive pseudo-spectral techniques to enable fully parallel 3D calculations of very large systems. In this work we present benchmark results for HFB-AX and MADNESS-HFB on ultracold trapped fermions.Comment: Conference on Computational Physics (CCP 2011) Proceedin

Nonequilibrium Electron Interactions in Metal Films

Ultrafast relaxation dynamics of an athermal electron distribution is investigated in silver films using a femtosecond pump-probe technique with 18 fs pulses in off-resonant conditions. The results yield evidence for an increase with time of the electron-gas energy loss rate to the lattice and of the free electron damping during the early stages of the electron-gas thermalization. These effects are attributed to transient alterations of the electron average scattering processes due to the athermal nature of the electron gas, in agreement with numerical simulations

Size-Dependent Surface Plasmon Dynamics in Metal Nanoparticles

We study the effect of Coulomb correlations on the ultrafast optical dynamics of small metal particles. We demonstrate that a surface-induced dynamical screening of the electron-electron interactions leads to quasiparticle scattering with collective surface excitations. In noble-metal nanoparticles, it results in an interband resonant scattering of d-holes with surface plasmons. We show that this size-dependent many-body effect manifests itself in the differential absorption dynamics for frequencies close to the surface plasmon resonance. In particular, our self-consistent calculations reveal a strong frequency dependence of the relaxation, in agreement with recent femtosecond pump-probe experiments.Comment: 8 pages + 4 figures, final version accepted to PR

Hole dynamics in noble metals

We present a detailed analysis of hole dynamics in noble metals (Cu and Au), by means of first-principles many-body calculations. While holes in a free-electron gas are known to live shorter than electrons with the same excitation energy, our results indicate that d-holes in noble metals exhibit longer inelastic lifetimes than excited sp-electrons, in agreement with experiment. The density of states available for d-hole decay is larger than that for the decay of excited electrons; however, the small overlap between d- and sp-states below the Fermi level increases the d-hole lifetime. The impact of d-hole dynamics on electron-hole correlation effects, which are of relevance in the analysis of time-resolved two-photon photoemission experiments, is also addressed.Comment: 4 pages, 2 figures, to appear in Phys. Rev. Let

Deformed Coordinate-Space Hartree-Fock-Bogoliubov Approach to Weakly Bound Nuclei and Large Deformations

The coordinate space formulation of the Hartree-Fock-Bogoliubov (HFB) method enables self-consistent treatment of mean-field and pairing in weakly bound systems whose properties are affected by the particle continuum space. Of particular interest are neutron-rich, deformed drip-line nuclei which can exhibit novel properties associated with neutron skin. To describe such systems theoretically, we developed an accurate 2D lattice Skyrme-HFB solver {\hfbax} based on B-splines. Compared to previous implementations, we made a number of improvements aimed at boosting the solver's performance. These include: explicit imposition of axiality and space inversion, use of the modified Broyden's method to solve self-consistent equations, and a partial parallelization of the code. {\hfbax} has been benchmarked against other HFB solvers, both spherical and deformed, and the accuracy of the B-spline expansion was tested by employing the multiresolution wavelet method. Illustrative calculations are carried out for stable and weakly bound nuclei at spherical and very deformed shapes, including constrained fission pathways. In addition to providing new physics insights, {\hfbax} can serve as a useful tool to assess the reliability and applicability of coordinate-space and configuration-space HFB solvers, both existing and in development.Comment: 12 pages,7 figs, submitted to Phys. Rev.

Analytical solutions to the third-harmonic generation in trans-polyacetylene: Application of dipole-dipole correlation on the single electron models

The analytical solutions for the third-harmonic generation (THG) on infinite chains in both Su-Shrieffer-Heeger (SSH) and Takayama-Lin-Liu-Maki (TLM) models of trans-polyacetylene are obtained through the scheme of dipole-dipole ($DD$) correlation. They are not equivalent to the results obtained through static current-current ($J_0J_0$) correlation or under polarization operator $\hat{P}$. The van Hove singularity disappears exactly in the analytical forms, showing that the experimentally observed two-photon absorption peak (TPA) in THG may not be directly explained by the single electron models.Comment: 10 pages, 4 figures, submitted to Phys. Rev.

MADNESS: A Multiresolution, Adaptive Numerical Environment for Scientific Simulation

MADNESS (multiresolution adaptive numerical environment for scientific simulation) is a high-level software environment for solving integral and differential equations in many dimensions that uses adaptive and fast harmonic analysis methods with guaranteed precision based on multiresolution analysis and separated representations. Underpinning the numerical capabilities is a powerful petascale parallel programming environment that aims to increase both programmer productivity and code scalability. This paper describes the features and capabilities of MADNESS and briefly discusses some current applications in chemistry and several areas of physics

Learning and Memory Alterations Are Associated with Hippocampal N-acetylaspartate in a Rat Model of Depression as Measured by 1H-MRS

It is generally accepted that cognitive processes, such as learning and memory, are affected in depression. The present study used a rat model of depression, chronic unpredictable mild stress (CUMS), to determine whether hippocampal volume and neurochemical changes were involved in learning and memory alterations. A further aim was to determine whether these effects could be ameliorated by escitalopram treatment, as assessed with the non-invasive techniques of structural magnetic resonance imaging (MRI) and magnetic resonance spectroscopy (MRS). Our results demonstrated that CUMS had a dramatic influence on spatial cognitive performance in the Morris water maze task, and CUMS reduced the concentration of neuronal marker N-acetylaspartate (NAA) in the hippocampus. These effects could be significantly reversed by repeated administration of escitalopram. However, neither chronic stress nor escitalopram treatment influenced hippocampal volume. Of note, the learning and memory alterations of the rats were associated with right hippocampal NAA concentration. Our results indicate that in depression, NAA may be a more sensitive measure of cognitive function than hippocampal volume