Evaluating kernels on Xeon Phi to accelerate Gysela application

This work describes the challenges presented by porting parts ofthe Gysela code to the Intel Xeon Phi coprocessor, as well as techniques used for optimization, vectorization and tuning that can be applied to other applications. We evaluate the performance of somegeneric micro-benchmark on Phi versus Intel Sandy Bridge. Several interpolation kernels useful for the Gysela application are analyzed and the performance are shown. Some memory-bound and compute-bound kernels are accelerated by a factor 2 on the Phi device compared to Sandy architecture. Nevertheless, it is hard, if not impossible, to reach a large fraction of the peek performance on the Phi device,especially for real-life applications as Gysela. A collateral benefit of this optimization and tuning work is that the execution time of Gysela (using 4D advections) has decreased on a standard architecture such as Intel Sandy Bridge.Comment: submitted to ESAIM proceedings for CEMRACS 2014 summer school version reviewe

Contribution of the screened self-energy to the Lamb shift of quasidegenerate states

Expressions for the effective Quantum Electrodynamics (QED) Hamiltonian due to self-energy screening (self-energy correction to the electron-electron interaction) are presented. We use the method of the two-time Green's function, which handles quasidegenerate atomic states. From these expression one can evaluate energy corrections to, e.g., 1s2p 3P1 and 1s2p 1P1 in helium and two-electron ions, to all orders in Z\alph

Relativistic and Radiative Energy Shifts for Rydberg States

We investigate relativistic and quantum electrodynamic effects for highly-excited bound states in hydrogenlike systems (Rydberg states). In particular, hydrogenic one-loop Bethe logarithms are calculated for all circular states (l = n-1) in the range 20 <= n <= 60 and successfully compared to an existing asymptotic expansion for large principal quantum number n. We provide accurate expansions of the Bethe logarithm for large values of n, for S, P and circular Rydberg states. These three expansions are expected to give any Bethe logarithms for principal quantum number n > 20 to an accuracy of five to seven decimal digits, within the specified manifolds of atomic states. Within the numerical accuracy, the results constitute unified, general formulas for quantum electrodynamic corrections whose validity is not restricted to a single atomic state. The results are relevant for accurate predictions of radiative shifts of Rydberg states and for the description of the recently investigated laser-dressed Lamb shift, which is observable in a strong coherent-wave light field.Comment: 8 pages; RevTeX

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

Line shape analysis of the Kβ\beta transition in muonic hydrogen

The K$\beta$ transition in muonic hydrogen was measured with a high-resolution crystal spectrometer. The spectrum is shown to be sensitive to the ground-state hyperfine splitting, the corresponding triplet-to-singlet ratio, and the kinetic energy distribution in the $3p$ state. The hyperfine splitting and triplet-to-singlet ratio are found to be consistent with the values expected from theoretical and experimental investigations and, therefore, were fixed accordingly in order to reduce the uncertainties in the further reconstruction of the kinetic energy distribution. The presence of high-energetic components was established and quantified in both a phenomenological, i.e. cascade-model-free fit, and in a direct deconvolution of the Doppler broadening based on the Bayesian approach.Comment: 22 pages, 21 figure

The first evidence for multiple pulsation axes: a new roAp star in the Kepler field, KIC 10195926

We have discovered a new rapidly oscillating Ap star among the Kepler Mission target stars, KIC 10195926. This star shows two pulsation modes with periods that are amongst the longest known for roAp stars at 17.1 min and 18.1 min, indicating that the star is near the terminal age main sequence. The principal pulsation mode is an oblique dipole mode that shows a rotationally split frequency septuplet that provides information on the geometry of the mode. The secondary mode also appears to be a dipole mode with a rotationally split triplet, but we are able to show within the improved oblique pulsator model that these two modes cannot have the same axis of pulsation. This is the first time for any pulsating star that evidence has been found for separate pulsation axes for different modes. The two modes are separated in frequency by 55 microHz, which we model as the large separation. The star is an alpha^2 CVn spotted magnetic variable that shows a complex rotational light variation with a period of Prot = 5.68459 d. For the first time for any spotted magnetic star of the upper main sequence, we find clear evidence of light variation with a period of twice the rotation period; i.e. a subharmonic frequency of $\nu_{\rm rot}/2$. We propose that this and other subharmonics are the first observed manifestation of torsional modes in an roAp star. From high resolution spectra we determine Teff = 7400 K, log g = 3.6 and v sin i = 21 km/s. We have found a magnetic pulsation model with fundamental parameters close to these values that reproduces the rotational variations of the two obliquely pulsating modes with different pulsation axes. The star shows overabundances of the rare earth elements, but these are not as extreme as most other roAp stars. The spectrum is variable with rotation, indicating surface abundance patches.Comment: 17 pages; 16 figures; MNRA

Theory of exciton-exciton correlation in nonlinear optical response

We present a systematic theory of Coulomb interaction effects in the nonlinear optical processes in semiconductors using a perturbation series in the exciting laser field. The third-order dynamical response consists of phase-space filling correction, mean-field exciton-exciton interaction, and two-exciton correlation effects expressed as a force-force correlation function. The theory provides a unified description of effects of bound and unbound biexcitons, including memory-effects beyond the Markovian approximation. Approximations for the correlation function are presented.Comment: RevTex, 35 pages, 10 PostScript figs, shorter version submitted to Physical Review

Fundamental properties of the Population II fiducial stars HD 122563 and Gmb 1830 from CHARA interferometric observations

We have determined the angular diameters of two metal-poor stars, HD 122563 and Gmb 1830, using CHARA and Palomar Testbed Interferometer observations. For the giant star HD 122563, we derive an angular diameter theta_3D = 0.940 +- 0.011 milliarcseconds (mas) using limb-darkening from 3D convection simulations and for the dwarf star Gmb 1830 (HD 103095) we obtain a 1D limb-darkened angular diameter theta_1D = 0.679 +- 0.007 mas. Coupling the angular diameters with photometry yields effective temperatures with precisions better than 55 K (Teff = 4598 +- 41 K and 4818 +- 54 K --- for the giant and the dwarf star, respectively). Including their distances results in very well-determined luminosities and radii (L = 230 +- 6 L_sun, R = 23.9 +- 1.9 R_sun and L = 0.213 +- 0.002 L_sun, R = 0.664 +- 0.015 R_sun, respectively). We used the CESAM2k stellar structure and evolution code in order to produce models that fit the observational data. We found values of the mixing-length parameter alpha (which describes 1D convection) that depend on the mass of the star. The masses were determined from the models with precisions of <3% and with the well-measured radii excellent constraints on the surface gravity are obtained (log g = 1.60 +- 0.04, 4.59 +- 0.02, respectively). The very small errors on both log g and Teff provide stringent constraints for spectroscopic analyses given the sensitivity of abundances to both of these values. The precise determination of Teff for the two stars brings into question the photometric scales for metal-poor stars.Comment: accepted A&A, 8 dbl-column pages, incl. 7 tables and 4 figure