Real-time approach to the optical properties of solids and nanostructures: Time-dependent Bethe-Salpeter equation

Many-body effects are known to play a crucial role in the electronic and optical properties of solids and nanostructures. Nevertheless, the majority of theoretical and numerical approaches able to capture the influence of Coulomb correlations are restricted to the linear response regime. In this work, we introduce an approach based on a real-time solution of the electronic dynamics. The proposed approach reduces to the well-known Bethe-Salpeter equation in the linear limit regime and it makes it possible, at the same time, to investigate correlation effects in nonlinear phenomena. We show the flexibility and numerical stability of the proposed approach by calculating the dielectric constants and the effect of a strong pulse excitation in bulk h-BN.

Improved description of chemical barriers with generalized gradient approximations (GGAs) and meta-GGAs

The performance of exchange and correlation (xc) functionals of the generalized gradient approximation (GGA) type and of the meta-GGA type in the calculation of chemical reactions is related to topological features of the electron density which, in turn, are connected to the orbital structure of chemical bonds within the Kohn-Sham (KS) theory. Seventeen GGA and meta-GGA xc functionals are assessed for 15 hydrogen abstraction reactions and 3 symmetrical S(N)2 reactions. Systems that are problematic for standard GGAs characteristically have enhanced values of the dimensionless gradient argument s(sigma)(2) with local maxima in the bonding region. The origin of this topological feature is the occupation of valence KS orbitals with an antibonding or essentially nonbonding character. The local enhancement of s(sigma)(2) yields too negative exchange-correlation energies with standard GGAs for the transition state of the S(N)2 reaction, which leads to the reduced calculated reaction barriers. The unwarranted localization of the effective xc hole of the standard GGAs, i.e., the nondynamical correlation that is built into them but is spurious in this case, wields its effect by their s(sigma)(2) dependence. Barriers are improved for xc functionals with the exchange functional OPTX as x component, which has a modified dependence on s(sigma)(2). Standard GGAs also underestimate the barriers for the hydrogen abstraction reactions. In this case the barriers are improved by correlation functionals, such as the Laplacian-dependent (LAP3) functional, which has a modified dependence on the Coulomb correlation of the opposite- and like-spin electrons. The best overall performance is established for the combination OLAP3 of OPTX and LAP3.

Exchange-correlation energy and potential as approximate functionals of occupied and virtual Kohn-Sham orbitals: Application to dissociating H-2

The standard local density approximation and generalized gradient approximations fail to properly describe the dissociation of an electron pair bond, yielding large errors (on the order of 50 kcal/mol) at long bond distances. To remedy this failure, a self-consistent Kohn-Sham (KS) method is proposed with the exchange-correlation (xc) energy and potential depending on both occupied and virtual KS orbitals. The xc energy functional of Buijse and Baerends [Mol. Phys. 100, 401 (2002); Phys. Rev. Lett. 87, 133004 (2001)] is employed, which, based on an ansatz for the xc-hole amplitude, is able to reproduce the important dynamical and nondynamical effects of Coulomb correlation through the efficient use of virtual orbitals. Self-consistent calculations require the corresponding xc potential to be obtained, to which end the optimized effective potential (OEP) method is used within the common energy denominator approximation for the static orbital Green's function. The problem of the asymptotic divergence of the xc potential of the OEP when a finite number of virtual orbitals is used is addressed. The self-consistent calculations reproduce very well the entire H-2 potential curve, describing correctly the gradual buildup of strong left-right correlation in stretched H-2. (C) 2003 American Institute of Physics.

Shape corrections to exchange-correlation potentials by gradient-regulated seamless connection of model potentials for inner and outer region

Shape corrections to the standard approximate Kohn-Sham exchange-correlation (xc) potentials are considered with the aim to improve the excitation energies (especially for higher excitations) calculated with time-dependent density functional perturbation theory. A scheme of gradient-regulated connection (GRAC) of inner to outer parts of a model potential is developed. Asymptotic corrections based either on the potential of Fermi and Amaldi or van Leeuwen and Baerends (LB) are seamlessly connected to the (shifted) xc potential of Becke and Perdew (BP) with the GRAC procedure, and are employed to calculate the vertical excitation energies of the prototype molecules N-2, CO, CH2O, C2H4, C5NH5, C6H6, Li-2, Na-2, K-2. The results are compared with those of the alternative interpolation scheme of Tozer and Handy as well as with the results of the potential obtained with the statistical averaging of (model) orbital potentials. Various asymptotically corrected potentials produce high quality excitation energies, which in quite a few cases approach the benchmark accuracy of 0.1 eV for the electronic spectra. Based on these results, the potential BP-GRAC-LB is proposed for molecular response calculations, which is a smooth potential and a genuine "local" density functional with an analytical representation. (C) 2001 American Institute of Physics.

Cycloidal magnetic order in the compound IrMnSi

A new compound, IrMnSi, has been synthesized, and its crystal structure and magnetic properties have been investigated by means of neutron powder diffraction, magnetization measurements, and first-principles theory. The crystal structure is found to be of the TiNiSi type (ordered Co2P, space group Pnma). The Mn-projected electronic states are situated at the Fermi level, giving rise to metallic binding, whereas a certain degree of covalent character is observed for the chemical bond between the It and Si atoms. A cycloidal, i.e., noncollinear, magnetic structure was observed below 460 K, with the propagation vector q=[0,0,0.4530(5)] at 10 K. The magnetism is dominated by large moments on the Mn sites, 3.8 mu(B)/atom from neutron diffraction. First-principles theory reproduces the propagation vector of the experimental magnetic structure as well as the angles between the Mn moments. The calculations further result in a magnetic moment of 3.21 mu(B) for the Mn atoms, whereas the Ir and Si moments are negligible, in agreement with observations. A calculation that more directly incorporates electron-electron interactions improves the agreement between the theoretical and experimental magnetic moments. A band mechanism is suggested to explain the observed magnetic order.

Deployment on GPUs of an application in computational atomic physics

This paper describes the deployment on GPUs of PROP, a program of the 2DRMP suite which models electron collisions with H-like atoms and ions. Because performance on GPUs is better in single precision than in double precision, the numerical stability of the PROP program in single precision has been studied. The numerical quality of PROP results computed in single precision and their impact on the next program of the 2DRMP suite has been analyzed. Successive versions of the PROP program on GPUs have been developed in order to improve its performance. Particular attention has been paid to the optimization of data transfers and of linear algebra operations. Performance obtained on several architectures (including NVIDIA Fermi) are presented.

GIANT LOBES OF CENTAURUS A RADIO GALAXY OBSERVED WITH THE SUZAKU X-RAY SATELLITE

We report on Suzaku observations of selected regions within the southern giant lobe of the radio galaxy Centaurus A. In our analysis we focus on distinct X-ray features detected with the X-ray Imaging Spectrometer within the range 0.5-10 keV, some of which are likely associated with fine structure of the lobe revealed by recent high-quality radio intensity and polarization maps. With the available photon statistics, we find that the spectral properties of the detected X-ray features are equally consistent with thermal emission from hot gas with temperatures kT > 1 keV, or with a power-law radiation continuum characterized by photon indices Gamma similar to 2.0 +/- 0.5. However, the plasma parameters implied by these different models favor a synchrotron origin for the analyzed X-ray spots, indicating that a very efficient acceleration of electrons up to greater than or similar to 10 TeV energies is taking place within the giant structure of Centaurus A, albeit only in isolated and compact regions associated with extended and highly polarized radio filaments. We also present a detailed analysis of the diffuse X-ray emission filling the whole field of view of the instrument, resulting in a tentative detection of a soft excess component best fitted by a thermal model with a temperature of kT similar to 0.5 keV. The exact origin of the observed excess remains uncertain, although energetic considerations point to thermal gas filling the bulk of the volume of the lobe and mixed with the non-thermal plasma, rather than to the alternative scenario involving a condensation of the hot intergalactic medium around the edges of the expanding radio structure. If correct, this would be the first detection of the thermal content of the extended lobes of a radio galaxy in X-rays. The corresponding number density of the thermal gas in such a case is n(g) similar to 10(-4) cm(-3), while its pressure appears to be in almost exact equipartition with the volume-averaged non-thermal pressure provided by the radio-emitting electrons and the lobes' magnetic field. A prominent large-scale fluctuation of the Galactic foreground emission, resulting in excess foreground X-ray emission aligned with the lobe, cannot be ruled out. Although tentative, our findings potentially imply that the structure of the extended lobes in active galaxies is likely to be highly inhomogeneous and non-uniform, with magnetic reconnection and turbulent acceleration processes continuously converting magnetic energy to internal energy of the plasma particles, leading to possibly significant spatial and temporal variations in the plasma beta parameter around the volume-averaged equilibrium condition beta similar to 1.

Diffusive shock acceleration at laser-driven shocks: studying cosmic-ray accelerators in the laboratory

The non-thermal particle spectra responsible for the emission from many astrophysical systems are thought to originate from shocks via a first order Fermi process otherwise known as diffusive shock acceleration. The same mechanism is also widely believed to be responsible for the production of high energy cosmic rays. With the growing interest in collisionless shock physics in laser produced plasmas, the possibility of reproducing and detecting shock acceleration in controlled laboratory experiments should be considered. The various experimental constraints that must be satisfied are reviewed. It is demonstrated that several currently operating laser facilities may fulfil the necessary criteria to confirm the occurrence of diffusive shock acceleration of electrons at laser produced shocks. Successful reproduction of Fermi acceleration in the laboratory could open a range of possibilities, providing insight into the complex plasma processes that occur near astrophysical sources of cosmic rays.

RADIATIVE SIGNATURES OF RELATIVISTIC SHOCKS

Particle-in-cell simulations of relativistic, weakly magnetized collisionless shocks show that particles can gain energy by repeatedly crossing the shock front. This requires scattering off self-generated small length-scale magnetic fluctuations. The radiative signature of this first-order Fermi acceleration mechanism is important for models of both the prompt and afterglow emission in gamma-ray bursts and depends on the strength parameter a = lambda e/delta B/mc(2) of the fluctuations (lambda is the length scale and vertical bar delta B vertical bar is the magnitude of the fluctuations). For electrons (and positrons), acceleration saturates when the radiative losses produced by the scattering cannot be compensated by the energy gained on crossing the shock. We show that this sets an upper limit on both the electron Lorentz factor gamma <10(6) (n/1 cm(-3))(-1/6)(-1/6) and on the energy of the photons radiated during the scattering process h omega(max) <40Max(a, 1)(n/1 cm(-3))(1/6)(-1/6) eV, where n is the number density of the plasma and (gamma) over bar is the thermal Lorentz factor of the downstream plasma, provided a <a(crit) similar to 10(6). This rules out "jitter" radiation on self-excited fluctuations with a <I as a source of gamma rays, although high-energy photons might still be produced when the jitter photons are upscattered in an analog of the synchrotron self-Compton process. In fluctuations with a > 1, radiation is generated by the standard synchrotron mechanism, and the maximum photon energy rises linearly with a, until saturating at 70 MeV, when a = a(crit).

Stochastic particle acceleration in the lobes of giant radio galaxies

We investigate the acceleration of particles by Alfven waves via the second-order Fermi process in the lobes of giant radio galaxies. Such sites are candidates for the accelerators of ultra-high-energy cosmic rays (UHECR). We focus on the nearby Fanaroff-Riley type I radio galaxy Centaurus A. This is motivated by the coincidence of its position with the arrival direction of several of the highest energy Auger events. The conditions necessary for consistency with the acceleration time-scales predicted by quasi-linear theory are reviewed. Test particle calculations are performed in fields which guarantee electric fields with no component parallel to the local magnetic field. The results of quasi-linear theory are, to an order of magnitude, found to be accurate at low turbulence levels for non-relativistic Alfven waves and at both low and high turbulence levels in the mildly relativistic case. We conclude that for pure stochastic acceleration via Alfven waves to be plausible as the generator of UHECR in Cen A, the baryon number density would need to be several orders of magnitude below currently held upper limits.

Quantum chemistry and in situ FTir spectroscopy studies on potential-dependent properties of CO adsorbed on Pt electrodes

The electronic and vibrational properties of CO adsorbed on Pt electrodes at different potentials have been studied, by using methods of self-consistent-charge discrete variational Xa (SCC-DV-Xa) cluster calculations and in situ FTir spectroscopy. Two new models have been developed and verified to be successful: (1) using a "metallic state cluster" to imitate a metal (electrode) surface; and (2) charging the cluster and shifting its Fermi level (e{lunate}) to simulate, according to the relation of -d e{lunate}e dE, quantitatively the variation of the electrode potential (E). It is shown that the binding of PtCO is dominated by the electric charge transfer of dp ? 2p, while that of s ? Pt is less important in this binding. The electron occupancy of the 2p orbital of CO weakens the CO bond and decreases the v. Variation of E mainly influences the charge transfer process of dp ? 2p, but hardly influences that of s ? Pt. A linear potential-dependence of v has been shown and the calculated dv/dE = 35.0 cm V. All results of calculations coincide with the ir experimental data. © 1993.

Positron scattering and annihilation in hydrogenlike ions

Diagrammatic many-body theory is used to calculate the scattering phase shifts, normalized annihilation rates Zeff, and annihilation ? spectra for positron collisions with the hydrogenlike ions He+, Li2+, B4+, and F8+. Short-range electron-positron correlations and longer-range positron-ion correlations are accounted for by evaluating nonlocal corrections to the annihilation vertex and the exact positron self-energy. The numerical calculation of the many-body theory diagrams is performed using B-spline basis sets. To elucidate the role of the positron-ion repulsion, the annihilation rate is also estimated analytically in the Coulomb-Born approximation. It is found that the energy dependence and magnitude of Zeff are governed by the Gamow factor that characterizes the suppression of the positron wave function near the ion. For all of the H-like ions, the correlation enhancement of the annihilation rate is found to be predominantly due to corrections to the annihilation vertex, while the corrections to the positron wave function play only a minor role. Results of the calculations for s-, p-, and d-wave incident positrons of energies up to the positronium-formation threshold are presented. Where comparison is possible, our values are in excellent agreement with the results obtained using other, e.g., variational, methods. The annihilation-vertex enhancement factors obtained in the present calculations are found to scale approximately as 1+(1.6+0.46l)/Zi, where Zi is the net charge of the ion and l is the positron orbital angular momentum. Our results for positron annihilation in H-like ions provide insights into the problem of positron annihilation with core electrons in atoms and condensed matter systems, which have similar binding energies.

Acceleration of ions up to 20MeV/nucleon in the ultrashort, high-intensity regime

The measurements reported here provide scaling laws for the ion acceleration process in the regime of ultrashort (50 fs), ultrahigh contrast (10) and ultrahigh intensity (> 10W/cm ), never investigated previously. The scaling of the accelerated ion energies was studied by varying a number of parameters such as target thickness (down to 10nm), target material (C and Al) and laser light polar- ization (circular and linear) at 35° and normal laser incidence. A twofold increase in proton energy and an order of magnitude enhancement in ion flux have been observed over the investigated thickness range at 35° angle of incidence. Further- more, at normal laser incidence, measured peak proton energies of about 20 MeV are observed almost independently of the target thickness over a wide range (50nm- 10 µm). 1. © 2012 by Società Italiana di Fisica.

Laser-driven acceleration of ions:State of the art, perspectives and applications

The acceleration of ions with high-power lasers has been a very active field of research during the past 10 years. This paper summarizes the main results obtained in the field, detailing the mechanisms of the acceleration process and the main observed beam characteristics. Perspectives for future development of the field and current and future applications are also discussed. © 2012 by Società Italiana di Fisica.

High energy resolution fluorescence detection XANES - an in situ method to study the interaction of adsorbed molecules with metal catalysts in the liquid phase

The change in the Pt electronic structure following the adsorption of an a,ß-unsaturated aldehyde and ketone was followed by in situ HERFD-XANES in the liquid phase. The resulting shift in the Pt Fermi energy is in good agreement with the molecule adsorption energy trends calculated by DFT and provides insight into the reaction selectivity.