Strong-field approximation for Coulomb explosion of H-2(+) and D-2(+) by short intense laser pulses

We present a simple quantum mechanical model to describe Coulomb explosion of H-2(+) and D-2(+) by short, intense infrared laser pulses. The model is based on the length gauge version of the molecular strong-field approximation and is valid when the process of dissociation prior to ionization is negligible. The results are compared with recent experimental data for the proton kinetic energy spectrum [Th. Ergler , Phys. Rev. Lett. 95, 093001 (2005); D. S. Murphy , J. Phys. B 40, S359 (2007)]. Using a Franck-Condon distribution over initial vibrational states, the theory reproduces the overall shape of the spectrum with only a small overestimation of slow protons. The agreement between theory and experiment can be made perfect by using a non-Frank-Condon initial distribution characteristic for H-2(+) (D-2(+)) targets produced by strong-field ionization of H-2 (D-2). For comparison, we also present results obtained by two different tunneling models for this process.

Testing numerical implementations of strong-field electrodynamics

We test current numerical implementations of laser-matter interactions by comparison with exact analytical results. Focusing on photon emission processes, it is found that the numerics accurately reproduce analytical emission spectra in all considered regimes, except for the harmonic structures often singled out as the most significant high-intensity (multiphoton) effects. We find that this discrepancy originates in the use of the locally constant field approximation.

Interactions of 1-Methylimidazole with UO2(CH3CO2)2 and UO2(NO3)2: Structural, Spectroscopic and Theoretical Evidence for Imidazole Binding to the Uranyl Ion

The first definitive high-resolution single-crystal X-ray structure for the coordination of the 1-methylimidazole (Meimid) ligand to UO2(Ac)2 (Ac = CH3CO2) is reported. The crystal structure evidence is confirmed by IR, Raman, and UV-vis spectroscopic data. Direct participation of the nitrogen atom of the Meimid ligand in binding to the uranium center is confirmed. Structural analysis at the DFT (B3LYP) level of theory showed a conformational difference of the Meimid ligand in the free gas-phase complex versus the solid state due to small energetic differences and crystal packing effects. Energetic analysis at the MP2 level in the gas phase supported stronger Meimid binding over H2O binding to both UO2(Ac)2 and UO2(NO3)2. In addition, self-consistent reaction field COSMO calculations were used to assess the aqueous phase energetics of combination and displacement reactions involving H2O and Meimid ligands to UO2R2 (R = Ac, NO3). For both UO2(NO3)2 and UO2(Ac)2, the displacement of H2O by Meimid was predicted to be energetically favorable, consistent with experimental results that suggest Meimid may bind uranyl at physiological pH. Also, log(Knitrate/KAc) calculations supported experimental evidence that the binding stoichiometry of the Meimid ligand is dependent upon the nature of the reactant uranyl complex. These results clearly demonstrate that imidazole binds to uranyl and suggest that binding of histidine residues to uranyl could occur under normal biological conditions.

Asymptotic modeling of short plane-parallel high-pressure glow discharge

Here a self-consistent continuum model is presented for a narrow gap plane-parallel dc glow discharge. The set of governing equations consisting of continuity and momentum equations for positive ions, fast (emitted by the cathode) and slow electrons (generated by fast electron impact ionization) coupled with Poisson's equation is treated by the technique of matched asymptotic expansions. Explicit results are obtained in the asymptotic limit: (chi delta) much less than 1, where chi = e Phi(a)/kT, delta = (r(D)/L)(2) (Phi(a) is the applied voltage, r(D) is the Debye radius) and pL much greater than 1(Hg mm cm), where p is the gas pressure and L is the gap length. In the case of high pressure, the electron energy relaxation length is much smaller than the gap length, and so the local field approximation is valid. The discharge space divides naturally into a cathode fall sheath, a quasineutral plasma region, and an anode fall sheath. The electric potential distribution obtained for each region in a (semi)analytical form is asymptotically matched to the adjoining regions in the region of overlap. The effects of the gas pressure, gap length, and applied voltage on the length of each region are investigated. (C) 2000 American Institute of Physics. [S1070-664X(00)01302-1].

Ionization energies of beryllium in strong magnetic fields: a frozen core approximation

An effective frozen core approximation has been developed and applied to the calculation of energy levels and ionization energies of the beryllium atom in magnetic field strengths up to 2.35 x 10(5) T. Systematic improvement over the existing results for the beryllium ground and low-lying states has been accomplished by taking into account most of the correlation effects in the four-electron system. To our knowledge, this is the first calculation of the electronic properties of the beryllium atom in a strong magnetic field carried out using a configuration interaction approximation and thus allowing a treatment beyond that of Hartree-Fock. Differing roles played by strong magnetic fields in intrashell correlation within different states are observed. In addition, possible ways to gain further improvement in the energies of the states of interest are proposed and discussed briefly.

CO oxidation on Pd(100) and Pd(111): A comparative study of reaction pathways and reactivity at low and medium coverages

We have performed density functional theory calculations with the generalized gradient approximation to investigate CO oxidation on a close-packed transition metal surface, Pd(lll), and a more open surface, Pd(100), aiming to shed light on surface structure effects on reaction pathways and reactivity, an important issue in catalysis. Reaction pathways on both surfaces at two different coverages have been studied. It is found that the reaction pathways on both surfaces possess crucial common features despite the fact that they have different surface symmetries. Having determined reaction barriers in these systems, we find that the reaction on Pd(lll) is strongly coverage dependent. Surface coverages, however, have little effect on the reaction on Pd(100). Calculations also reveal that the low coverage reactions are structure sensitive while the medium coverage reactions are not. Detailed discussions on these results are given.

One-dimensional particle simulation of the filamentation instability: Electrostatic field driven by the magnetic pressure gradient force

Two counterpropagating cool and equally dense electron beams are modeled with particle-in-cell simulations. The electron beam filamentation instability is examined in one spatial dimension, which is an approximation for a quasiplanar filament boundary. It is confirmed that the force on the electrons imposed by the electrostatic field, which develops during the nonlinear stage of the instability, oscillates around a mean value that equals the magnetic pressure gradient force. The forces acting on the electrons due to the electrostatic and the magnetic field have a similar strength. The electrostatic field reduces the confining force close to the stable equilibrium of each filament and increases it farther away, limiting the peak density. The confining time-averaged total potential permits an overlap of current filaments with an opposite flow direction.

Statistical-mechanical description of classical test-particle dynamics in the presence of an external force field: modelling noise and damping from first principles

Aiming to establish a rigorous link between macroscopic random motion (described e.g. by Langevin-type theories) and microscopic dynamics, we have undertaken a kinetic-theoretical study of the dynamics of a classical test-particle weakly coupled to a large heat-bath in thermal equilibrium. Both subsystems are subject to an external force field. From the (time-non-local) generalized master equation a Fokker-Planck-type equation follows as a "quasi-Markovian" approximation. The kinetic operator thus defined is shown to be ill-defined; in specific, it does not preserve the positivity of the test-particle distribution function f(x, v; t). Adopting an alternative approach, previously introduced for quantum open systems, is proposed to lead to a correct kinetic operator, which yields all the expected properties. A set of explicit expressions for the diffusion and drift coefficients are obtained, allowing for modelling macroscopic diffusion and dynamical friction phenomena, in terms of an external field and intrinsic physical parameters.

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.

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.

Dielectrics in a time-dependent electric field: a real-time approach based on density-polarization functional theory

In the presence of a (time-dependent) macroscopic electric field the electron dynamics of dielectrics cannot be described by the time-dependent density only. We present a real-time formalism that has the density and the macroscopic polarization P as key quantities. We show that a simple local function of P already captures long-range correlation in linear and nonlinear optical response functions. Specifically, after detailing the numerical implementation, we examine the optical absorption, the second- and third-harmonic generation of bulk Si, GaAs, AlAs and CdTe at different level of approximation. We highlight links with ultranonlocal exchange-correlation functional approximations proposed within linear response time-dependent density functional theory framework.

Tracing the Chromospheric and Coronal Magnetic Field with AIA, IRIS, IBIS, and ROSA Data

The aim of this study is to explore the suitability of chromospheric images for magnetic modeling of active regions. We use high-resolutionimages (≈0.2"-0.3"), from the Interferometric Bidimensional Spectrometer in the Ca II 8542 Å line, the Rapid Oscillations in the Solar Atmosphere instrument in the Hα 6563Å line, the Interface Region Imaging Spectrograph in the 2796Å line, and compare non-potential magnetic field models obtainedfrom those chromospheric images with those obtained from images of the Atmospheric Imaging Assembly in coronal (171 Å, etc.) and inchromospheric (304 Å) wavelengths. Curvi-linear structures are automatically traced in those images with the OCCULT-2 code, to which we forward-fitted magnetic field lines computed with the Vertical-current Approximation Nonlinear Force Free Field code. We find that the chromospheric images: (1) reveal crisp curvi-linear structures (fibrils, loop segments, spicules) that are extremely well-suited for constraining magnetic modeling; (2) that these curvi-linear structures arefield-aligned with the best-fit solution by a median misalignment angle of μ2 ≈ 4°–7° (3) the free energy computed from coronal data may underestimate that obtained from chromospheric data by a factor of ≈2–4, (4) the height range of chromospheric features is confined to h≲4000 km, while coronal features are detected up to h = 35,000 km; and (5) the plasma-β parameter is β ≈ 10^-5 - 10^-1 for all traced features. We conclude that chromospheric images reveal important magnetic structures that are complementary to coronal images and need to be included in comprehensive magnetic field models, something that is currently not accomodated in standard NLFFF codes.