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.

Exchange potential from the common energy denominator approximation for the Kohn-Sham Green's function: Application to (hyper)polarizabilities of molecular chains

An approximate Kohn-Sham (KS) exchange potential v(xsigma)(CEDA) is developed, based on the common energy denominator approximation (CEDA) for the static orbital Green's function, which preserves the essential structure of the density response function. v(xsigma)(CEDA) is an explicit functional of the occupied KS orbitals, which has the Slater v(Ssigma) and response v(respsigma)(CEDA) potentials as its components. The latter exhibits the characteristic step structure with "diagonal" contributions from the orbital densities \psi(isigma)\(2), as well as "off-diagonal" ones from the occupied-occupied orbital products psi(isigma)psi(j(not equal1)sigma). Comparison of the results of atomic and molecular ground-state CEDA calculations with those of the Krieger-Li-Iafrate (KLI), exact exchange (EXX), and Hartree-Fock (HF) methods show, that both KLI and CEDA potentials can be considered as very good analytical "closure approximations" to the exact KS exchange potential. The total CEDA and KLI energies nearly coincide with the EXX ones and the corresponding orbital energies epsilon(isigma) are rather close to each other for the light atoms and small molecules considered. The CEDA, KLI, EXX-epsilon(isigma) values provide the qualitatively correct order of ionizations and they give an estimate of VIPs comparable to that of the HF Koopmans' theorem. However, the additional off-diagonal orbital structure of v(xsigma)(CEDA) appears to be essential for the calculated response properties of molecular chains. KLI already considerably improves the calculated (hyper)polarizabilities of the prototype hydrogen chains H-n over local density approximation (LDA) and standard generalized gradient approximations (GGAs), while the CEDA results are definitely an improvement over the KLI ones. The reasons of this success are the specific orbital structures of the CEDA and KLI response potentials, which produce in an external field an ultranonlocal field-counteracting exchange potential. (C) 2002 American Institute of Physics.

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.

Free energy and molecular dynamics calculations for the cubic-tetragonal phase transition in zirconia

The high-temperature cubic-tetragonal phase transition of pure stoichiometric zirconia is studied by molecular dynamics (MD) simulations and within the framework of the Landau theory of phase transformations. The interatomic forces are calculated using an empirical, self-consistent, orthogonal tight-binding model, which includes atomic polarizabilities up to the quadrupolar level. A first set of standard MD calculations shows that, on increasing temperature, one particular vibrational frequency softens. The temperature evolution of the free-energy surfaces around the phase transition is then studied with a second set of calculations. These combine the thermodynamic integration technique with constrained MD simulations. The results seem to support the thesis of a second-order phase transition but with unusual, very anharmonic behavior above the transition temperature.

Relative energetics and structural properties of zirconia using a self-consistent tight-binding model

We describe an empirical, self-consistent, orthogonal tight-binding model for zirconia, which allows for the polarizability of the anions at dipole and quadrupole levels and for crystal field splitting of the cation d orbitals, This is achieved by mixing the orbitals of different symmetry on a site with coupling coefficients driven by the Coulomb potentials up to octapole level. The additional forces on atoms due to the self-consistency and polarizabilities are exactly obtained by straightforward electrostatics, by analogy with the Hellmann-Feynman theorem as applied in first-principles calculations. The model correctly orders the zero temperature energies of all zirconia polymorphs. The Zr-O matrix elements of the Hamiltonian, which measure covalency, make a greater contribution than the polarizability to the energy differences between phases. Results for elastic constants of the cubic and tetragonal phases and phonon frequencies of the cubic phase are also presented and compared with some experimental data and first-principles calculations. We suggest that the model will be useful for studying finite temperature effects by means of molecular dynamics.

Reduced-size polarized basis sets for calculations of molecular electric properties. III. Second-row atoms

Reduced-size polarized (ZmPolX) basis sets are developed for the second-row atoms X = Si, P, S, and Cl. The generation of these basis sets follows from a simple physical model of the polarization effect of the external electric field which leads to highly compact polarization functions to be added to the chosen initial basis set. The performance of the ZmPolX sets has been investigated in calculations of molecular dipole moments and polarizabilities. Only a small deterioration of the quality of the calculated molecular electric properties has been found. Simultaneously the size of the present reduced-size ZmPolX basis sets is about one-third smaller than that of the usual polarized (PolX) sets. This reduction considerably widens the range of applications of the ZmPolX sets in calculations of molecular dipole moments, dipole polarizabilities, and related properties.

Quantum circuits for spin and flavor degrees of freedom of quarks forming nucleons

We discuss the quantum-circuit realization of the state of a nucleon in the scope of simple simmetry groups. Explicit algorithms are presented for the preparation of the state of a neutron or a proton as resulting from the composition of their quark constituents. We estimate the computational resources required for such a simulation and design a photonic network for its implementation. Moreover, we highlight that current work on three-body interactions in lattices of interacting qubits, combined with the measurement-based paradigm for quantum information processing, may also be suitable for the implementation of these nucleonic spin states.

Ion Acceleration in Multispecies Targets Driven by Intense Laser Radiation Pressure

The acceleration of ions from ultrathin foils has been investigated by using 250 TW, subpicosecond laser pulses, focused to intensities of up to 3 X 10(20) W cm(-2). The ion spectra show the appearance of narrow-band features for protons and carbon ions peaked at higher energies (in the 5-10 MeV/nucleon range) and with significantly higher flux than previously reported. The spectral features and their scaling with laser and target parameters provide evidence of a multispecies scenario of radiation pressure acceleration in the light sail mode, as confirmed by analytical estimates and 2D particle-in-cell simulations. The scaling indicates that monoenergetic peaks with more than 100 MeV/nucleon are obtainable with moderate improvements of the target and laser characteristics, which are within reach of ongoing technical developments.

STO and GTO field-induced polarization functions for H to Kr

Field-induced polarization (FIP) functions were proposed over two decades ago to improve the accuracy of calculated response properties, and the FIP functions in GTO form for H and C to F were tested on small molecules, with encouraging results. The concept of FIP,is now extended to all atoms up to Kr. New simplifying approximations for the description of asymptotic highest occupied atomic orbitals. (HOAOs) are introduced in this study. They provide the basis for STO and GTO exponents of a complete set of FIP functions from H to Kr, which are both listed for the convenience of the users. Tests on the polarizabilities of a series of atoms and molecules demonstrate that addition of the FIP basis functions to a series' of standard basis sets drastically improves the performance of all these basis sets compared to converged results. Moreover, the byproduct of this study (approximate asymptotic HOAOs) provides information for the construction of accurate basis sets for long-range ground state properties. (C) 2003 Wiley Periodicals, Inc.

Charge equilibrium of a laser-generated carbon-ion beam in warm dense matter

Using ion carbon beams generated by high intensity short pulse lasers we perform measurements of single shot mean charge equilibration in cold or isochorically heated solid density aluminum matter. We demonstrate that plasma effects in such matter heated up to 1 eV do not significantly impact the equilibration of carbon ions with energies 0.045-0.5 MeV/nucleon. Furthermore, these measurements allow for a first evaluation of semiempirical formulas or ab initio models that are being used to predict the mean of the equilibrium charge state distribution for light ions passing through warm dense matter.

Antiproton induced DNA damage: proton like in flight, carbon-ion like near rest

Biological validation of new radiotherapy modalities is essential to understand their therapeutic potential. Antiprotons have been proposed for cancer therapy due to enhanced dose deposition provided by antiproton-nucleon annihilation. We assessed cellular DNA damage and relative biological effectiveness (RBE) of a clinically relevant antiproton beam. Despite a modest LET (,19 keV/mm), antiproton spread out Bragg peak (SOBP) irradiation caused significant residual c-H2AX foci compared to X-ray, proton and antiproton plateau irradiation. RBE of ,1.48 in the SOBP and ,1 in the plateau were measured and used for a qualitative effective dose curve comparison with proton and carbon-ions. Foci in the antiproton SOBP were larger and more structured compared to X-rays, protons and carbon-ions. This is likely due to overlapping particle tracks near the annihilation vertex, creating spatially correlated DNA lesions. No biological effects were observed at 28–42 mm away from the primary beam suggesting minimal risk from long-range secondary particles.

Experimental investigation of hole boring and light sail regimes of RPA by varying laser and target parameters

Temporal evolution of plasma jets from micrometre-scale thick foils following the interaction of intense (3 × 10 W cm ) laser pulses is studied systematically by time resolved optical interferometry. The fluid velocity in the plasma jets is determined by comparing the data with 2D hydrodynamic simulation, which agrees with the expected hole-boring (HB) velocity due to the laser radiation pressure. The homogeneity of the plasma density across the jets has been found to be improved substantially when irradiating the laser at circular polarization compared to linear polarization. While overdense plasma jets were formed efficiently for micrometre thick targets, decreasing the target areal density and/or increasing the irradiance on the target have provided indication of transition from the 'HB' to the 'light sail (LS)' regime of RPA, characterized by the appearance of narrow-band spectral features at several MeV/nucleon in proton and carbon spectra.

Modified Thomson spectrometer design for high energy, multi-species ion sources

A modification to the standard Thomson parabola spectrometer is discussed, which is designed to measure high energy (tens of MeV/nucleon), broad bandwidth spectra of multi-species ions accelerated by intense laser plasma interactions. It is proposed to implement a pair of extended, trapezoidal shaped electric plates, which will not only resolve ion traces at high energies, but will also retain the lower energy part of the spectrum. While a longer (along the axis of the undeflected ion beam direction) electric plate design provides effective charge state separation at the high energy end of the spectrum, the proposed new trapezoidal shape will enable the low energy ions to reach the detector, which would have been clipped or blocked by simply extending the rectangular plates to enhance the electrostatic deflection.