Electronic excitations in solids studied using inelastic x-ray scattering

Inelastic x-ray scattering can be used to study the electronic structure of matter. The x rays scattered from the target both induce and carry information on the electronic excitations taking place in the system. These excitations are the manifestations of the electronic structure and the physics governing the many-body system. This work presents results of non-resonant inelastic x-ray scattering experiments on a range of materials including metallic, insulating and semiconducting compounds as well as an organic polymer. The experiments were carried out at the National Synchrotron Light Source, USA and at the European Synchrotron Radiation Facility, France. The momentum transfer dependence of the experimental valence- and core-electron excitation spectra is compared with the results of theoretical first principles computations that incorporate the electron-hole interaction. A recently developed method for analyzing the momentum transfer dependence of core-electron excitation spectra is studied in detail. This method is based on real space multiple scattering calculations and is used to extract the angular symmetry components of the local unoccupied density of final states.

Fast calculation technique for scattering in T-matrix method

In scattering calculations using the T-matrix method, the calculation of the T-matrix involves multiplication and inversion of matrices. These two types of matrix operations are time-consuming, especially for the matrices with large size. Petrov et al. [D. Petrov, Y. Shkuratov, G. Videen, Opt. Lett. 32 (2007) 1168] proposed an optimized matrix inversion technique, which suggests the inversion of two matrices, each of which contains half the number of rows. This technique reduces time-consumption significantly. On the basis of this approach, we propose another fast calculation technique for scattering in the T-matrix method, which obtains the scattered fields through carrying out only the operations between matrices and the incident field coefficient. Numerical results show that this technique can decrease time-consumption by more than half that of the optimized matrix inversion technique by Petrov et al. (c) 2008 Elsevier B.V. All rights reserved.

Low energy elastic and electronically inelastic electron scattering from biomolecules.

Reactions initiated by collisions with low-energy secondary electrons has been found to be the prominent mechanism toward the radiation damage on living tissues through DNA strand breaks. Now it is widely accepted that during the interaction with these secondary species the selective breaking of chemical bonds is triggered by dissociative electron attachment (DEA), that is, the capture of the incident electron and the formation of temporary negative ion states [1,2,3]. One of the approaches largely used toward a deeper understanding of the radiation damage to DNA is through modeling of DEA with its basic constituents (nucleotide bases, sugar and other subunits). We have tried to simplify this approach and attempt to make it comprehensible at a more fundamental level by looking at even simple molecules. Studies involving organic systems such as carboxylic acids, alcohols and simple ¯ve-membered heterocyclic compounds are taken as starting points for these understanding. In the present study we investigate the role played by elastic scattering and electronic excitation of molecules on electron-driven chemical processes. Special attention is focused on the analysis of the in°uence of polarization and multichannel coupling e®ects on the magnitude of elastic and electronically inelastic cross-sections. Our aim is also to investigate the existence of resonances in the elastic and electronically inelastic channels as well as to characterize them with respect to its type (shape, core-excited or Feshbach), symmetry and position. The relevance of these issues is evaluated within the context of possible applications for the modeling of discharge environments and implications in the understanding of mutagenic rupture of DNA chains. The scattering calculations were carried out with the Schwinger multichannel method (SMC) [4] and its implementation with pseudopotentials (SMCPP) [5] at di®erent levels of approximation for impact energies ranging from 0.5 eV to 30 eV. References [1] B. Boudai®a, P. Cloutier, D. Hunting, M. A. Huels and L. Sanche, Science 287, 1658 (2000). [2] X. Pan, P. Cloutier, D. Hunting and L. Sanche, Phys. Rev. Lett. 90, 208102 (2003). [3] F. Martin, P. D. Burrow, Z. Cai, P. Cloutier, D. Hunting and L. Sanche, Phys. Rev. Lett. 93, 068101 (2004). [4] K. Takatsuka and V. McKoy, Phys. Rev. A 24, 2437 (1981); ibid. Phys. Rev. A 30, 1734 (1984). [5] M. H. F. Bettega, L. G. Ferreira and M. A. P. Lima, Phys. Rev. A 47, 1111 (1993).

First-principles study of ferroelectricity and isotope effects in H-bonded KH2PO4 crystals

By means of extensive first-principles calculations we studied the ferroelectric phase transition and the associated isotope effect in KH2PO4 (KDP). Our calculations revealed that the spontaneous polarization of the ferroelectric phase is due to electronic charge redistributions and ionic displacements which are a consequence of proton ordering, and not vice versa. The experimentally observed double-peaked proton distribution in the paraelectric phase cannot be explained by a dynamics of only protons. This requires, instead, collective displacements within clusters that include also the heavier ions. These tunneling clusters can explain the recent evidence of tunneling obtained from Compton scattering measurements. The sole effect of mass change upon deuteration is not sufficient to explain the huge isotope effect. Instead, we find that structural modifications deeply connected with the chemistry of the H bonds produce a feedback effect on tunneling that strongly enhances the phenomenon. The resulting influence of the geometric changes on the isotope effect agrees with experimental data from neutron scattering. Calculations under pressure allowed us to analyze the issue of universality in the disappearance of ferroelectricity upon compression. Compressing DKDP so that the distance between the two peaks in the deuteron distribution is the same as for protons in KDP, corresponds to a modification of the underlying double-well potential, which becomes 23 meV shallower. This energy difference is what is required to modify the O-O distance in such a way as to have the same distribution for protons and deuterons. At the high pressures required experimentally, the above feedback mechanism is crucial to explain the magnitude of the geometrical effect.

Direct Observation of strong coupling in a laser-shock compressed plasma

In this Letter we report on a near collective x-ray scattering experiment on shock-compressed targets. A highly coupled Al plasma was generated and probed by spectrally resolving an x-ray source forward scattered by the sample. A significant reduction in the intensity of the elastic scatter was observed, which we attribute to the formation of an incipient long-range order. This speculation is confirmed by x-ray scattering calculations accounting for both electron degeneracy and strong coupling effects. Measurements from rear side visible diagnostics are consistent with the plasma parameters inferred from x-ray scattering data. These results give the experimental evidence of the strongly coupled ionic dynamics in dense plasmas.

Dissociative electron attachment to the highly reactive difluoromethylene molecule–importance of CF2 for negative ion formation in fluorocarbon plasmas

Dissociative electron attachment to the highly reactive difluoromethylene molecule, CF2, produced in a C3F6/He microwave plasma and stepwise via the fast atom reaction CF3I+H?CF3+HI and CF3+H?CF2+HF, has been investigated. The upper limit for the cross section of formation of F- via dissociative electron attachment to CF2 is estimated to be 5×10-4?Å2. This value is four orders of magnitude smaller than the cross section previously predicted from scattering calculations. It is concluded that difluoromethylene plays a negligible role in negative ion formation in fluorocarbon plasmas.

High-resolution electron energy-loss spectroscopy of Ba Ti O 3∕ Sr Ti O 3 multilayers

The dielectric properties of BaTiO3 thin films and multilayers are different from bulk materials because of nanoscale dimensions, interfaces, and stress-strain conditions. In this study, BaTiO3/SrTiO3 multilayers deposited on SrTiO3 substrates by pulsed laser deposition have been investigated by high-energy-resolution electron energy-loss spectroscopy. The fine structures in the spectra are discussed in terms of crystal-field splitting and the internal strain. The crystal-field splitting of the BaTiO3 thin layer is found to be a little larger than that of bulk BaTiO3, which has been interpreted by the presence of the internal strain induced by the misfit at the interface. This finding is consistent with the lattice parameters of the BaTiO3 thin layer determined by the selected area diffraction pattern. The near-edge structure of the oxygen K edge in BaTiO3 thin layers and in bulk BaTiO3 are simulated by first-principle self-consistent full multiple-scattering calculations. The results of the simulations are in a good agreement with the experimental results. Moreover, the aggregation of oxygen vacancies at the rough BaTiO3/SrTiO3 interface is indicated by the increased [Ti]/[O] element ratio, which dominates the difference of dielectric properties between BaTiO3 layer and bulk materials.

The 2p interloper resonances in collisions

Ab initio electron scattering calculations using the R -matrix approach have been performed for within a three-state valence configuration-interaction model (VCI). The lowest three electronic target states ( , and the ) of this molecular nitrogen cation are included in the close-coupling method, with each state being represented by a valence CI approximation. From a detailed analysis of the resonance structure found in our work for the symmetries we find four prominent Rydberg series of the type , , , and a interloper resonance. This interloper molecular resonance associated with the B state of is seen to cause distortions of the resulting resonance spectra. A comparison of our total cross sections for the X - B transition shows excellent agreement with the available experimental data.

The Structure of the chiral Pt{531} surface: a combined LEED and DFT study

The structure of the chiral kinked Pt{531} surface has been determined by low-energy electron diffraction intensity-versus-energy (LEED-IV) analysis and density functional theory (DFT). Large contractions and expansions of the vertical interlayer distances with respect to the bulk-terminated surface geometry were found for the first six layers (LEED: d(12) = 0.44 angstrom, d(23) = 0.69 angstrom, d(34) = 0.49 angstrom, d(45) = 0.95 angstrom, d(56) = 0.56 angstrom; DFT: d(12) = 0.51 angstrom, d(23) = 0.55 angstrom, d(34) = 0.74 angstrom, d(45) = 0.78 angstrom, d(56) = 0.63 angstrom; d(bulk) = 0.66 angstrom). Energy-dependent cancellations of LEED spots over unusually large energy ranges, up to 100 eV, can be explained by surface roughness and reproduced by applying a model involving 0.25 ML of vacancies and adatoms in the scattering calculations. The agreement between the results from LEED and DFT is not as good as in other cases, which could be due to this roughness of the real surface.

A FORTRAN-90 low-energy electron diffraction program (LEED90 v1.1)

We describe a FORTRAN-90 program to compute low-energy electron diffraction I(V) curves. Plane-waves and layer doubling are used to compute the inter-layer multiple-scattering, while the intra-layer multiple-scattering is computed in the standard way expanding the wavefield on a basis of spherical waves. The program is kept as general as possible, in order to allow testing different parts of multiple-scattering calculations. In particular, it can handle non-diagonal t-matrices describing the scattering of non-spherical potentials, anisotropic vibrations, anharmonicity, etc. The program does not use old FORTRAN flavours, and has been written keeping in mind the advantage for parallelism brought forward by FORTRAN-90.

RELATIVISTIC 3-PARTICLE DYNAMICAL EQUATIONS .1. THEORETICAL DEVELOPMENT

Starting from the two-particle Bethe-Salpeter equation in the ladder approximation and integrating over the time component of momentum, we rederive three-dimensional scattering integral equations satisfying constraints of relativistic unitarity and convariance, first derived by Weinberg and by Blankenbecler and Sugar. These two-particle equations are shown to be related by a transformation of variables. Hence we show how to perform and relate identical dynamical calculation using these two equations. Similarly, starting from the Bethe-Salpeter-Faddeev equation for the three-particle system and integrating over the time component of momentum, we derive several three-dimensional three-particle scattering equations satisfying constraints of relativistic unitarity and convariance. We relate two of these three-particle equations by a transformation of variables as in the two-particle case. The three-particle equations we derive are very practical and suitable for performing relativistic scattering calculations. (C) 1994 Academic Press, Inc.

Dissociative electron attachment to chloromethane.

A low energy electron may attach to a molecule, forming a metastable resonance, which may dissociate into a stable anion and a neutral radical. Chloromethane has been a good target for dissociative electron attachment studies, since it is a small molecule with a clear dissociative ‘sigma*’ shape resonance. We present potential energy curves for CH3Cl and its anion, as a function of the C-Cl distance. Due to the resonant nature of the anion, a correct description requires a treatment based on scattering calculations. In order to compute elastic cross sections and phase shifts we employed the Schwinger multichannel method, implemented with pseudopotentials of Bachelet, Hamann and Schlüter, at the static-exchange plus polarization approximation. At the equilibrium geometry, the resonance was found arround 3.3 eV, in accordance to experience. The incoming electron is captured by a ‘sigma*’ orbital located at the C-Cl bond, which will relax in the presence of this extra electron. We took this bond as the reaction coordinate, and performed several scattering calculations for a series of nuclear conformations. The phase shift obtained in each calculation was fitted by a two component function, consisting in the usual Breit-Wigner profile, which captures the resonant character, and a second order polynomial in the wave number, which accounts for the background contribution. That way, we obtained position and width of the resonance, which allowed us to build the potential energy curve. For larger distances, the anion becomes stable and usual electronic structure calculations suffice. Furthermore, the existence of a dipole-bound anion state is revealed when we employed a set of very diffuse functions. The knowledge on the behaviour of the neutral and anionic electronic states helps us in elucidating how the dissociation takes place.

Shape resonance spectra of uracil, 5-fluorouracil, and 5-chlorouracil.

We report on the shape resonance spectra of uracil, 5-fluorouracil, and 5-chlorouracil, as obtained from fixed-nuclei elastic scattering calculations performed with the Schwinger multichannel method with pseudopotentials. Our results are in good agreement with the available electron transmission spectroscopy data, and support the existence of three π* resonances in uracil and 5-fluorouracil. As expected, the anion states are more stable in the substituted molecules than in uracil. Since the stabilization is stronger in 5-chlorouracil, the lowest π* resonance in this system becomes a bound anion state. The present results also support the existence of a low-lying σ ∗ CCl shape resonance in 5- chlorouracil. Exploratory calculations performed at selected C–Cl bond lengths suggest that the σ ∗ CCl resonance could couple to the two lowest π* states, giving rise to a very rich dissociation dynamics. These facts would be compatible with the complex branching of the dissociative electron attachment cross sections, even though we cannot discuss any details of the vibration dynamics based only on the present fixed-nuclei results.

Electron Interactions with Disulfide Bridges

Because of its electronic properties, sulfur plays a major role in a variety of metabolic processes and, more in general, in the chemistry of life. In particular, S-S bridges between cysteines are present in the amino acid backbone of proteins. Protein disulfur radical anions may decay following different paths through competing intra and intermolecular routes, including bond cleavage, disproportionation, protein-protein cross linking, and electron transfer. Indeed, mass spectrometry ECD (electron capture dissociation massspectroscopy) studies have shown that capture of low-energy (<0.2 eV) electrons by multiply protonated proteins is followed by dissociation of S-S bonds holding two peptide chains together. In view of the importance of organic sulfur chemistry, we report on electron interactions with disulphide bridges. To study these interactions we used as prototypes the molecules dimethyl sulfide [(CH3)2S] and dimethyl disulfide [(H3C)S2(CH3)]. We seek to better understand the electron-induced cleavage of the disulfide bond. To explore dissociative processes we performed electron scattering calculations with the Schwinger Multichannel Method with pseudopotentials (SMCPP), recently parallelized with OpenMP directives and optimized with subroutines for linear algebra (BLAS) and LAPACK routines. Elastic cross sections obtained for different S-S bond lengths indicate stabilization of the anion formed by electron attachment to a σ*SS antibonding orbital, such that dissociation would be expected.

The calculation of Afρ and mass loss rate for comets

Ab initio calculations of Afρ are presented using Mie scattering theory and a Direct Simulation Monte Carlo (DSMC) dust outflow model in support of the Rosetta mission and its target 67P/Churyumov-Gerasimenko (CG). These calculations are performed for particle sizes ranging from 0.010 μm to 1.0 cm. The present status of our knowledge of various differential particle size distributions is reviewed and a variety of particle size distributions is used to explore their effect on Afρ , and the dust mass production View the MathML sourcem˙. A new simple two parameter particle size distribution that curtails the effect of particles below 1 μm is developed. The contributions of all particle sizes are summed to get a resulting overall Afρ. The resultant Afρ could not easily be predicted a priori and turned out to be considerably more constraining regarding the mass loss rate than expected. It is found that a proper calculation of Afρ combined with a good Afρ measurement can constrain the dust/gas ratio in the coma of comets as well as other methods presently available. Phase curves of Afρ versus scattering angle are calculated and produce good agreement with observational data. The major conclusions of our calculations are: – The original definition of A in Afρ is problematical and Afρ should be: qsca(n,λ)×p(g)×f×ρqsca(n,λ)×p(g)×f×ρ. Nevertheless, we keep the present nomenclature of Afρ as a measured quantity for an ensemble of coma particles.– The ratio between Afρ and the dust mass loss rate View the MathML sourcem˙ is dominated by the particle size distribution. – For most particle size distributions presently in use, small particles in the range from 0.10 to 1.0 μm contribute a large fraction to Afρ. – Simplifying the calculation of Afρ by considering only large particles and approximating qsca does not represent a realistic model. Mie scattering theory or if necessary, more complex scattering calculations must be used. – For the commonly used particle size distribution, dn/da ∼ a−3.5 to a−4, there is a natural cut off in Afρ contribution for both small and large particles. – The scattering phase function must be taken into account for each particle size; otherwise the contribution of large particles can be over-estimated by a factor of 10. – Using an imaginary index of refraction of i = 0.10 does not produce sufficient backscattering to match observational data. – A mixture of dark particles with i ⩾ 0.10 and brighter silicate particles with i ⩽ 0.04 matches the observed phase curves quite well. – Using current observational constraints, we find the dust/gas mass-production ratio of CG at 1.3 AU is confined to a range of 0.03–0.5 with a reasonably likely value around 0.1.