Nonadiabatic forces in ion-solid interactions:The initial stages of radiation damage

The Born-Oppenheimer approximation is the keystone for molecular dynamics simulations of radiation damage processes; however, actual materials response involves nonadiabatic energy exchange between nuclei and electrons. In this work, time dependent density functional theory is used to calculate the electronic excitations produced by energetic protons in Al. We study the influence of these electronic excitations on the interatomic forces and find that they differ substantially from the adiabatic case, revealing a nontrivial connection between electronic and nuclear stopping that is absent in the adiabatic case. These results unveil new effects in the early stages of radiation damage cascades.

The UK Renal Registry Advanced CKD Study: frequency of incorrect reporting of date of start of RRT

Background: A preliminary review of the UK Renal Registry (UKRR) pre-RRT study data revealed results suggesting that, for some patients, the date of start of renal replacement therapy (RRT), as reported to the UKRR, was incorrect and often significantly later than the true date of start. A more detailed study then aimed to validate a set of criteria to identify patients with an incorrect start date. Methods: Pre-RRT laboratory data were electronically extracted from 8,810 incident RRT patients from 9 UK renal centres. Any patient with a low urea (<15 mmol/L) at the start of RRT or with a substantial improvement in kidney function (either a fall in urea >10 mmol/L or rise in eGFR >2 ml/min/1.73 m) within the two months prior to RRT were considered to potentially have an incorrect date of start. In 4
selected centres, the electronic patient records of all patients flagged were reviewed to validate these criteria.
Results: Of 8,810 patients, 1,616 (18.3%) were flagged by the identification criteria as having a potentially incorrect date of start of RRT, although a single centre accounted for 41% of the total flagged cohort. Of these flagged patients, 61.7% had been assigned an incorrect date of start of haemodialysis (HD), 5.7% had evidence of acute RRT being given before the reported date of start of HD
and 9.2% had evidence of starting peritoneal dialysis exchanges prior to the reported date of start. Of
those flagged, 10.7% had a correct date of start of RRT.
Conclusions: Accurate reporting of RRT episodes is vital for the analysis of time dependent studies such as survival or time to transplantation. A proportion of patients starting RRT were assigned an incorrect start date. In order to improve the accuracy of this reporting the UK Renal Registry
must work with renal centres and clinical staff on improving data input for the start of RRT.

Uptake, translocation and transformation of arsenate and arsenite in sunflower (Helianthus annuus):formation of arsenic-phytochelatin complexes during exposure to high arsenic concentrations

The aim of the study was to determine the time-dependent formation of arsenic-phytochelatin (As-PC) complexes in the roots, stems and leaves of an arsenic-nontolerant plant (Helianthus annuus) during exposure to 66 mol l(-1) arsenite (As(III)) or arsenate (As(V)). We used our previously developed method of simultaneous element-specific (inductively coupled plasma mass spectrometry, ICP-MS) and molecular-specific (electrospray-ionization mass spectrometry, ES-MS) detection systems interfaced with a suitable chromatographic column and eluent conditions, which enabled us to identify and quantify As-PC complexes directly. Roots of As-exposed H. annuus contained up to 14 different arsenic species, including the complex of arsenite with two (gamma-Glu-Cys)(2)-Gly molecules [As((III))-(PC(2))(2)], the newly identified monomethylarsonic phytochelatin-2 or (gamma-Glu-Cys)(2)-Gly CH(3)As (MA((III))-PC(2)) and at least eight not yet identified species. The complex of arsenite with (gamma-Glu-Cys)(3)-Gly (As((III))-PC(3)) and the complex of arsenite with glutathione (GSH) and (gamma-Glu-Cys)(2)-Gly (GS-As((III))-PC(2)) were present in all samples (roots, stems and leaves) taken from plants exposed to As. The GS-As((III))-PC(2) complex was the dominant complex after 1 h of exposure. As((III))-PC(3) became the predominant As-PC complex after 3 h, binding up to 40% of the As present in the exposed plants. No As-PC complexes were found in sap (mainly xylem sap from the root system), in contrast to roots, stems and leaves, which is unequivocal evidence that As-PC complexes are not involved in the translocation of As from root to leaves of H. annuus.

A complete analytic theory for structure and dynamics of populations and communities spanning wide ranges in body size

The prediction and management of ecosystem responses to global environmental change would profit from a clearer understanding of the mechanisms determining the structure and dynamics of ecological communities. The analytic theory presented here develops a causally closed picture for the mechanisms controlling community and population size structure, in particular community size spectra, and their dynamic responses to perturbations, with emphasis on marine ecosystems. Important implications are summarised in non-technical form. These include the identification of three different responses of community size spectra to size-specific pressures (of which one is the classical trophic cascade), an explanation for the observed slow recovery of fish communities from exploitation, and clarification of the mechanism controlling predation mortality rates. The theory builds on a community model that describes trophic interactions among size-structured populations and explicitly represents the full life cycles of species. An approximate time-dependent analytic solution of the model is obtained by coarse graining over maturation body sizes to obtain a simple description of the model steady state, linearising near the steady state, and then eliminating intraspecific size structure by means of the quasi-neutral approximation. The result is a convolution equation for trophic interactions among species of different maturation body sizes, which is solved analytically using a novel technique based on a multiscale expansion.

Gamma Irradiation and Targeted Radionuclides Enhance the Expression of the Noradrenaline Transporter Transgene Controlled by the Radio-Inducible p21(WAF1/CIP1) Promoter

The use of radiation-inducible promoters to drive transgene expression offers the possibility of temporal and spatial regulation of gene activation. This study assessed the potential of one such promoter element, p21(WAF1/CIP1) (WAF1), to drive expression of the noradrenaline transporter (NAT) gene, which conveys sensitivity to radioiodinated meta-iodobenzylguanidine (MIBG). An expression vector containing NAT under the control of the radiation-inducible WAF1 promoter (pWAF/NAT) was produced. The non-NAT expressing cell lines UVW (glioma) and HCT116 (colorectal cancer) were transfected with this construct to assess radiation-controlled WAF1 activation of the NAT gene. Transfection of UVW and HCT cells with pWAF/NAT conferred upon them the ability to accumulate [(131)I]MIBG, which led to increased sensitivity to the radiopharmaceutical. Pretreatment of transfected cells with ? radiation or the radiopharmaceuticals [(123)I]MIBG or [(131)I]MIBG induced dose- and time-dependent increases in subsequent [(131)I]MIBG uptake and led to enhanced efficacy of [(131)I]MIBG-mediated cell kill. Gene therapy using WAF1-driven expression of NAT has the potential to expand the use of this therapeutic modality to tumors that lack a radio-targetable feature.

Influence of multiple ionization thresholds on harmonic generation from Ar+

We apply time-dependent R-matrix theory to investigate harmonic generation from ground-state Ar+ with M = 0 at a wavelength of 390 nm. Contributions associated with the different 3s(2)3p(4) ionization thresholds are assessed, including the interference between these. The dominant contribution originates from the second ionization threshold, 3s(2)3p(4 1)D. Changes to the harmonic yields arising from the higher 3s3p(5) thresholds are also assessed. We further confirm that Ar+ has a higher harmonic yield than He for the same laser pulse, despite having a higher ionization threshold. 

Electronic stopping power of H and He in Al and LiF from first principles

Non-linearities in the electronic stopping power of light projectiles in bulk Al and LiF are addressed from first principles using time-evolving time-dependent density functional theory. In the case of Al, the agreement of the calculations with experiments for H and He projectiles is fair, but a recently observed transition for He from one value of the electronic friction coefficient to a higher value at v ~ 0.3 a.u. is not reproduced by the calculations. For LiF, better accuracy is obtained as compared with previously published simulations, albeit the threshold remains overestimated. © 2013 Elsevier B.V.

yambo: An ab initio tool for excited state calculations

yambo is an ab initio code for calculating quasiparticle energies and optical properties of electronic systems within the framework of many-body perturbation theory and time-dependent density functional theory. Quasiparticle energies are calculated within the GW approximation for the self-energy. Optical properties are evaluated either by solving the Bethe-Salpeter equation or by using the adiabatic local density approximation. yambo is a plane-wave code that, although particularly suited for calculations of periodic bulk systems, has been applied to a large variety of physical systems. yambo relies on efficient numerical techniques devised to treat systems with reduced dimensionality, or with a large number of degrees of freedom. The code has a user-friendly command-line based interface, flexible 110 procedures and is interfaced to several publicly available density functional ground-state codes.

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.

On the required shape corrections to the local density and generalized gradient approximations to the Kohn-Sham potentials for molecular response calculations of (hyper)polarizabilities and excitation energies

It is well known that shape corrections have to be applied to the local-density (LDA) and generalized gradient (GGA) approximations to the Kohn-Sham exchange-correlation potential in order to obtain reliable response properties in time dependent density functional theory calculations. Here we demonstrate that it is an oversimplified view that these shape corrections concern primarily the asymptotic part of the potential, and that they affect only Rydberg type transitions. The performance is assessed of two shape-corrected Kohn-Sham potentials, the gradient-regulated asymptotic connection procedure applied to the Becke-Perdew potential (BP-GRAC) and the statistical averaging of (model) orbital potentials (SAOP), versus LDA and GGA potentials, in molecular response calculations of the static average polarizability alpha, the Cauchy coefficient S-4, and the static average hyperpolarizability beta. The nature of the distortions of the LDA/GGA potentials is highlighted and it is shown that they introduce many spurious excited states at too low energy which may mix with valence excited states, resulting in wrong excited state compositions. They also lead to wrong oscillator strengths and thus to a wrong spectral structure of properties like the polarizability. LDA, Becke-Lee-Yang-Parr (BLYP), and Becke-Perdew (BP) characteristically underestimate contributions to alpha and S-4 from bound Rydberg-type states and overestimate those from the continuum. Cancellation of the errors in these contributions occasionally produces fortuitously good results. The distortions of the LDA, BLYP, and BP spectra are related to the deficiencies of the LDA/GGA potentials in both the bulk and outer molecular regions. In contrast, both SAOP and BP-GRAC potentials produce high quality polarizabilities for 21 molecules and also reliable Cauchy moments and hyperpolarizabilities for the selected molecules. The analysis for the N-2 molecule shows, that both SAOP and BP-GRAC yield reliable energies omega(i) and oscillator strengths f(i) of individual excitations, so that they reproduce well the spectral structure of alpha and S-4.(C) 2002 American Institute of Physics.

Density Functional study of the photoactive yellow protein's chromophore

We have investigated the structural and electronic properties of p-coumaric acid, the chromophore of the photoactive yellow protein (PYP), by means of first-principles molecular dynamics based on density functional theory (DFT). We have studied the chromophore both in the vacuum and in an extended model which includes the nearest residues in the binding pocket of PYP, as derived from crystallographic data. We have characterized the ground state of the isolated chromophore in its protonated and deprotonated forms and computed the energy barrier involved in the trans to cis isomerization process around the carbon-carbon double bond. A comparison of the optimized structures of the chromophore in the vacuum and in the extended protein model, both in the trans (ground state of PYP in the dark) and cis (first light-activated intermediate) configuration, shows how the protein environment affects the chromophore in the first step of the photocycle. Our model gives an energy storage of 25 kcal/mol associated with the trans-to-cia photoisomerization. Finally, we have elucidated the nature of the electronic excitation relevant for the photochemistry of PYP by means of time-dependent DFT calculations.

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.

Double-detonation sub-Chandrasekhar supernovae: Synthetic observables for minimum helium shell mass models

In the double-detonation scenario for Type Ia supernovae, it is suggested that a detonation initiates in a shell of helium-rich material accreted from a companion star by a sub-Chandrasekhar-mass white dwarf. This shell detonation drives a shock front into the carbon-oxygen white dwarf that triggers a secondary detonation in the core. The core detonation results in a complete disruption of the white dwarf. Earlier studies concluded that this scenario has difficulties in accounting for the observed properties of Type Ia supernovae since the explosion ejecta are surrounded by the products of explosive helium burning in the shell. Recently, however, it was proposed that detonations might be possible for much less massive helium shells than previously assumed (Bildsten et al.). Moreover, it was shown that even detonations of these minimum helium shell masses robustly trigger detonations of the carbon-oxygen core (Fink et al.). Therefore, it is possible that the impact of the helium layer on observables is less than previously thought. Here, we present time-dependent multi-wavelength radiative transfer calculations for models with minimum helium shell mass and derive synthetic observables for both the optical and ? -ray spectral regions. These differ strongly from those found in earlier simulations of sub-Chandrasekhar-mass explosions in which more massive helium shells were considered. Our models predict light curves that cover both the range of brightnesses and the rise and decline times of observed Type Ia supernovae. However, their colors and spectra do not match the observations. In particular, their B - V colors are generally too red. We show that this discrepancy is mainly due to the composition of the burning products of the helium shell of the Fink et al. models which contain significant amounts of titanium and chromium. Using a toy model, we also show that the burning products of the helium shell depend crucially on its initial composition. This leads us to conclude that good agreement between sub-Chandrasekhar-mass explosions and observed Type Ia supernovae may still be feasible but further study of the shell properties is required.

Light curves for off-centre ignition models of Type Ia supernovae

Motivated by recent models involving off-centre ignition of Type Ia supernova explosions, we undertake three-dimensional time-dependent radiation transport simulations to investigate the range of bolometric light-curve properties that could be observed from supernovae in which there is a lop-sided distribution of the products from nuclear burning. We consider both a grid of artificial toy models which illustrate the conceivable range of effects and a recent three-dimensional hydrodynamical explosion model. We find that observationally significant viewing angle effects are likely to arise in such supernovae and that these may have important ramifications for the interpretation of the observed diversity of Type Ia supernova and the systematic uncertainties which relate to their use as standard candles in contemporary cosmology. © 2007 RAS.

On the γ-ray emission of Type Ia supernovae

A multidimension, time-dependent Monte Carlo code is used to compute sample ?-ray spectra to explore whether unambiguous constraints could be obtained from ?-ray observations of Type Ia supernovae. Both spherical and aspherical geometries are considered and it is shown that moderate departures from sphericity can produce viewing-angle effects that are at least as significant as those caused by the variation of key parameters in 1D models. Thus, ?-ray data could, in principle, carry some geometrical information, and caution should be applied when discussing the value of ?-ray data based only on 1D explosion models. In light of the limited sensitivity of current ?-ray observatories, the computed theoretical spectra are studied to revisit the issue of whether useful constraints could be obtained for moderately nearby objects. The most useful ?-ray measurements are likely to be of the light curve and time-dependent hardness ratios, but sensitivity higher than currently available, particularly at relatively hard energies (~2-3 MeV), is desirable. © 2008 The Authors. Journal compilation © 2008 RAS.