Implementation and testing of Lanczos-based algorithms for Random-Phase Approximation eigenproblems

The treatment of the Random-Phase Approximation Hamiltonians, encountered in different frameworks, like time-dependent density functional theory or Bethe-Salpeter equation, is complicated by their non-Hermicity. Compared to their Hermitian Hamiltonian counterparts, computational methods for the treatment of non-Hermitian Hamiltonians are often less efficient and less stable, sometimes leading to the breakdown of the method. Recently [Gruning et al. Nano Lett. 8 (2009) 28201, we have identified that such Hamiltonians are usually pseudo-Hermitian. Exploiting this property, we have implemented an algorithm of the Lanczos type for Random-Phase Approximation Hamiltonians that benefits from the same stability and computational load as its Hermitian counterpart, and applied it to the study of the optical response of carbon nanotubes. We present here the related theoretical grounds and technical details, and study the performance of the algorithm for the calculation of the optical absorption of a molecule within the Bethe-Salpeter equation framework. (C) 2011 Elsevier B.V. All rights reserved.

Use of a random phase plate as a KrF laser beam homogenizer for thin film deposition applications

Fabrication of devices based on thin film structures deposited using the pulsed laser deposition technique relies on reproducibility and control of deposition rates over substrate areas as large as possible. Here we present an application of the random phase plate technique to smooth and homogenize the intensity distribution of a KrF laser footprint on the surface of a target which is to be ablated. It is demonstrated that intensity distributions over millimeter-sized spots on the target can be made insensitive to the typical changes that occur in the near-field intensity distribution of the ultraviolet output from a KrF laser. (C) 1999 American Institute of Physics. [S0034-6748(99)02504-6].

Photo-double detachment from the F- ion

The correlated process of photodetaching two electrons from the F- ion following the absorption of a single photon has been investigated over an energy range 20-62 eV. In the experiment, a beam of photons from the Advanced Light Source was collinearly merged with a counter-propagating beam of F- ions from a sputter ion source. The F+ ions produced in the interaction region were detected, and the normalized signal was used to monitor the relative cross section for the double-detachment reaction. An absolute scale for the cross section was established by measuring the spatial overlap of the two beams and by determining the efficiency for collection and detection of the F+ ions. The measured cross section is compared with R-matrix and random phase approximation calculations. These calculations show that the Auger decay of the 2s2p(6) core-excited state of the F atom plays a minor role in the production of F+ ions and that double detachment is likely to be dominated by simultaneous correlated ejection of two valence electrons at energies well above threshold.

Local-field and excitonic effects in the optical response of alpha-alumina

This paper introduces key ingredients of the dielectric response of a-alumina that go beyond an independent-particle (IP) treatment of the valence-electron excitations. The optical-response functions were calculated from first-principles both at the Bethe-Salpeter and the random-phase approximation (RPA) levels. Excitonic effects obtained within the Bethe-Salpeter framework were found essential for reproducing the low-energy part of the experimental spectra (below 15 eV) and the bound exciton in particular. For higher energies, local-field effects introduced through the RPA modified considerably the IP results and provided a satisfactory account of the reflectivity spectra and of the position and shape of the dominant bulk plasmon resonance in the electron energy-loss spectra.

Density functionals from many-body perturbation theory: The band gap for semiconductors and insulators

Theoretically the Kohn-Sham band gap differs from the exact quasiparticle energy gap by the derivative discontinuity of the exchange-correlation functional. In practice for semiconductors and insulators the band gap calculated within any local or semilocal density approximations underestimates severely the experimental energy gap. On the other hand, calculations with an "exact" exchange potential derived from many-body perturbation theory via the optimized effective potential suggest that improving the exchange-correlation potential approximation can yield a reasonable agreement between the Kohn-Sham band gap and the experimental gap. The results in this work show that this is not the case. In fact, we add to the exact exchange the correlation that corresponds to the dynamical (random phase approximation) screening in the GW approximation. This accurate exchange-correlation potential provides band structures similar to the local density approximation with the corresponding derivative discontinuity that contributes 30%-50% to the energy gap. Our self-consistent results confirm substantially the results for Si and other semiconductors obtained perturbatively [R. W. Godby , Phys. Rev. B 36, 6497 (1987)] and extend the conclusion to LiF and Ar, a wide-gap insulator and a noble-gas solid. (c) 2006 American Institute of Physics.

Performance of polarisation functionals for linear and nonlinear optical properties of bulk zinc chalcogenides ZnX (X=S,Se,Te).

We calculated the frequency dependent macroscopic dielectric function and second-harmonic generation of cubic ZnS, ZnSe and ZnTe within time-dependent density-polarisation functional theory. The macroscopic dielectric function is calculated in a linear response framework, and second-harmonic generation in a real-time framework. The macroscopic exchange–correlation electric field that enters the time-dependent Kohn–Sham equations and accounts for long range correlation is approximated as a simple polarisation functional αP, where P is the macroscopic polarisation. Expressions for α are taken from the recent literature. The performance of the resulting approximations for the exchange–correlation electric field is analysed by comparing the theoretical spectra with experimental results and results obtained at the levels of the independent particle approximation and the random-phase approximation. For the dielectric function we also compare with state-of-the art calculations at the level of the Bethe–Salpeter equation.

van der Waals coefficients for positronium interactions with atoms

The random-phase approximation with exchange (RPAE) is used with a B-spline basis to compute dynamic dipole polarizabilities of noble-gas atoms and several other closed-shell atoms (Be, Mg, Ca, Zn, Sr, Cd, and Ba). From these, values of the van der Waals C₆ constants for positronium interactions with these atoms are determined and compared with existing data. After correcting the RPAE polarizabilities to fit the most accurate static polarizability data, our best predictions of C₆ for Ps–noble-gas pairs are expected to be accurate to within 1%, and to within a few percent for the alkaline-earth metals. We also used accurate dynamic dipole polarizabilities from the literature to compute the C₆ coefficients for the alkali-metal atoms. Implications of increased C₆ values for Ps scattering from more polarizable atoms are discussed.

Electron number density measurements in magnesium laser produced plumes

Interferometry has been used to investigate the spatio-temporal evolution of electron number density following 248 nm laser ablation of a magnesium target. Fringe shifts were measured as a function of laser power density for a circular spot obtained using a random phase plate. Line averaged electron number densities were obtained at delay times up to ∼100 ns after the laser pulse. Density profiles normal to the target surface were recorded for power densities on target in the range 125–300 MW cm−2.

Absolute photoionization cross sections for Xe4+, Xe5+, and Xe6+ near 13.5 nm: Experiment and theory

Absolute photoionization cross-section measurements for a mixture of ground and metastable states of Xe4+, Xe5+, and Xe6+ are reported in the photon energy range of 4d -> nf transitions, which occur within or adjacent to the 13.5 nm window for extreme ultraviolet lithography light source development. The reported values allow the quantification of opacity effects in xenon plasmas due to these 4d -> nf autoionizing states. The oscillator strengths for the 4d -> 4f and 4d -> 5f transitions in Xeq+ (q=1-6) ions are calculated using nonrelativistic Hartree-Fock and random phase approximations. These are compared with published experimental values for Xe+ to Xe3+ and with the values obtained from the present experimental cross-section measurements for Xe4+ to Xe6+. The calculations assisted in the determination of the metastable content in the ion beams for Xe5+ and Xe6+. The experiments were performed by merging a synchrotron photon beam generated by an undulator beamline of the Advanced Light Source with an ion beam produced by an electron cyclotron resonance ion source.

Macroscopic limit of time-dependent density-functional theory for adiabatic local approximations of the exchange-correlation kernel

Time-dependent density-functional theory is a rather accurate and efficient way to compute electronic excitations for finite systems. However, in the macroscopic limit (systems of increasing size), for the usual adiabatic random-phase, local-density, or generalized-gradient approximations, one recovers the Kohn-Sham independent-particle picture, and thus the incorrect band gap. To clarify this trend, we investigate the macroscopic limit of the exchange-correlation kernel in such approximations by means of an algebraical analysis complemented with numerical studies of a one-dimensional tight-binding model. We link the failure to shift the Kohn-Sham spectrum of these approximate kernels to the fact that the corresponding operators in the transition space act only on a finite subspace.

Simulation of quantum random walks using the interference of a classical field

We suggest a theoretical scheme for the simulation of quantum random walks on a line using beam splitters, phase shifters, and photodetectors. Our model enables us to simulate a quantum random walk using of the wave nature of classical light fields. Furthermore, the proposed setup allows the analysis of the effects of decoherence. The transition from a pure mean-photon-number distribution to a classical one is studied varying the decoherence parameters.

Phase and space resolved optical emission spectroscopic investigations of an inductively coupled RF plasma using an imaging acousto-optic spectrometer

A novel acousto-optic spectrometer (IfU Diagnostic Systems GmbH) for 2-dimensional (2D) optical emission spectroscopy with high spectral resolution has been developed. The spectrometer is based on acousto-optic tuneable filter technology with fast random wavelength access. Measurements for characterisation of the imaging quality, the spatial resolution, and the spectral resolution are presented. The applicability for 2D-space and phase resolved optical emission spectroscopy (2D-PROES) is shown. 2D-PROES has been applied to an inductively coupled plasma with radio frequency excitation at 13.56 MHz.

Systematic random sampling of the comet assay

The comet assay is a technique used to quantify DNA damage and repair at a cellular level. In the assay, cells are embedded in agarose and the cellular content is stripped away leaving only the DNA trapped in an agarose cavity which can then be electrophoresed. The damaged DNA can enter the agarose and migrate while the undamaged DNA cannot and is retained. DNA damage is measured as the proportion of the migratory ‘tail’ DNA compared to the total DNA in the cell. The fundamental basis of these arbitrary values is obtained in the comet acquisition phase using fluorescence microscopy with a stoichiometric stain in tandem with image analysis software. Current methods deployed in such an acquisition are expected to be both objectively and randomly obtained. In this paper we examine the ‘randomness’ of the acquisition phase and suggest an alternative method that offers both objective and unbiased comet selection. In order to achieve this, we have adopted a survey sampling approach widely used in stereology, which offers a method of systematic random sampling (SRS). This is desirable as it offers an impartial and reproducible method of comet analysis that can be used both manually or automated. By making use of an unbiased sampling frame and using microscope verniers, we are able to increase the precision of estimates of DNA damage. Results obtained from a multiple-user pooled variation experiment showed that the SRS technique attained a lower variability than that of the traditional approach. The analysis of a single user with repetition experiment showed greater individual variances while not being detrimental to overall averages. This would suggest that the SRS method offers a better reflection of DNA damage for a given slide and also offers better user reproducibility.

Random particle motion in magnetized plasma

A multivariate Fokker-Planck-type kinetic equation modeling a test - panicle weakly interacting with an electrostatic plasma. in the presence of a magnetic field B . is analytically solved in an Ornstein - Uhlenbeck - type approximation. A new set of analytic expressions are obtained for variable moments and panicle density as a function of time. The process is diffusive.

Current-induced atomic dynamics, instabilities, and Raman signals: Quasiclassical Langevin equation approach

We derive and employ a semiclassical Langevin equation obtained from path integrals to describe the ionic dynamics of a molecular junction in the presence of electrical current. The electronic environment serves as an effective nonequilibrium bath. The bath results in random forces describing Joule heating, current-induced forces including the nonconservative wind force, dissipative frictional forces, and an effective Lorentz-type force due to the Berry phase of the nonequilibrium electrons. Using a generic two-level molecular model, we highlight the importance of both current-induced forces and Joule heating for the stability of the system. We compare the impact of the different forces, and the wide-band approximation for the electronic structure on our result. We examine the current-induced instabilities (excitation of runaway "waterwheel" modes) and investigate the signature of these in the Raman signals.