Reverse engineering in many-body quantum physics: Correspondence between many-body systems and effective single-particle equations

The mapping, exact or approximate, of a many-body problem onto an effective single-body problem is one of the most widely used conceptual and computational tools of physics. Here, we propose and investigate the inverse map of effective approximate single-particle equations onto the corresponding many-particle system. This approach allows us to understand which interacting system a given single-particle approximation is actually describing, and how far this is from the original physical many-body system. We illustrate the resulting reverse engineering process by means of the Kohn-Sham equations of density-functional theory. In this application, our procedure sheds light on the nonlocality of the density-potential mapping of density-functional theory, and on the self-interaction error inherent in approximate density functionals.

Physical signatures of discontinuities of the time-dependent exchange-correlation potential

The exact exchange-correlation (XC) potential in time-dependent density-functional theory (TDDFT) is known to develop steps and discontinuities upon change of the particle number in spatially confined regions or isolated subsystems. We demonstrate that the self-interaction corrected adiabatic local-density approximation for the XC potential has this property, using the example of electron loss of a model quantum well system. We then study the influence of the XC potential discontinuity in a real-time simulation of a dissociation process of an asymmetric double quantum well system, and show that it dramatically affects the population of the resulting isolated single quantum wells. This indicates the importance of a proper account of the discontinuities in TDDFT descriptions of ionization, dissociation or charge transfer processes.

Intermediate motions and dipolar couplings as studied by Lee-Goldburg cross-polarization NMR: Hartmann-Hahn matching profiles

In this article, we evaluate the use of simple Lee-Goldburg cross-polarization (LG-CP) NMR experiments for obtaining quantitative information of molecular motion in the intermediate regime. In particular, we introduce the measurement of Hartmann-Hahn matching profiles for the assessment of heteronuclear dipolar couplings as well as dynamics as a reliable and robust alternative to the more common analysis of build-up curves. We have carried out dynamic spin dynamics simulations in order to test the method's sensitivity to intermediate motion and address its limitations concerning possible experimental imperfections. We further demonstrate the successful use of simple theoretical concepts, most prominently Anderson-Weiss (AW) theory, to analyze the data. We further propose an alternative way to estimate activation energies of molecular motions, based upon the acquisition of only two LG-CP spectra per temperature at different temperatures. As experimental tests, molecular jumps in imidazole methyl sulfonate, trimethylsulfoxonium iodide, and bisphenol A polycarbonate were investigated with the new method.

Determination of the molecular weight of proteins in solution from a single small-angle X-ray scattering measurement on a relative scale

This paper describes a new and simple method to determine the molecular weight of proteins in dilute solution, with an error smaller than similar to 10%, by using the experimental data of a single small-angle X-ray scattering (SAXS) curve measured on a relative scale. This procedure does not require the measurement of SAXS intensity on an absolute scale and does not involve a comparison with another SAXS curve determined from a known standard protein. The proposed procedure can be applied to monodisperse systems of proteins in dilute solution, either in monomeric or multimeric state, and it has been successfully tested on SAXS data experimentally determined for proteins with known molecular weights. It is shown here that the molecular weights determined by this procedure deviate from the known values by less than 10% in each case and the average error for the test set of 21 proteins was 5.3%. Importantly, this method allows for an unambiguous determination of the multimeric state of proteins with known molecular weights.

Accurate band gaps of AlGaN, InGaN, and AlInN alloys calculations based on LDA-1/2 approach

We present parameter-free calculations of electronic properties of InGaN, InAlN, and AlGaN alloys. The calculations are based on a generalized quasichemical approach, to account for disorder and composition effects, and first-principles calculations within the density functional theory with the LDA-1/2 approach, to accurately determine the band gaps. We provide precise results for AlGaN, InGaN, and AlInN band gaps for the entire range of compositions, and their respective bowing parameters. (C) 2011 American Institute of Physics. [doi:10.1063/1.3576570]

Low-energy electron scattering from methanol and ethanol

Measured and calculated differential cross sections for elastic (rotationally unresolved) electron scattering from two primary alcohols, methanol (CH(3)OH) and ethanol (C(2)H(5)OH), are reported. The measurements are obtained using the relative flow method with helium as the standard gas and a thin aperture as the collimating target gas source. The relative flow method is applied without the restriction imposed by the relative flow pressure conditions on helium and the unknown gas. The experimental data were taken at incident electron energies of 1, 2, 5, 10, 15, 20, 30, 50, and 100 eV and for scattering angles of 5 degrees-130 degrees. There are no previous reports of experimental electron scattering differential cross sections for CH(3)OH and C(2)H(5)OH in the literature. The calculated differential cross sections are obtained using two different implementations of the Schwinger multichannel method, one that takes all electrons into account and is adapted for parallel computers, and another that uses pseudopotentials and considers only the valence electrons. Comparison between theory and experiment shows that theory is able to describe low-energy electron scattering from these polyatomic targets quite well.

Transition-metal 13-atom clusters assessed with solid and surface-biased functionals

First-principles density-functional theory studies have reported open structures based on the formation of double simple-cubic (DSC) arrangements for Ru(13), Rh(13), Os(13), and Ir(13), which can be considered an unexpected result as those elements crystallize in compact bulk structures such as the face-centered cubic and hexagonal close-packed lattices. In this work, we investigated with the projected augmented wave method the dependence of the lowest-energy structure on the local and semilocal exchange-correlation (xc) energy functionals employed in density-functional theory. We found that the local-density approximation (LDA) and generalized-gradient formulations with different treatment of the electronic inhomogeneities (PBE, PBEsol, and AM05) confirm the DSC configuration as the lowest-energy structure for the studied TM(13) clusters. A good agreement in the relative total energies are obtained even for structures with small energy differences, e. g., 0.10 eV. The employed xc functionals yield the same total magnetic moment for a given structure, i.e., the differences in the bond lengths do not affect the moments, which can be attributed to the atomic character of those clusters. Thus, at least for those systems, the differences among the LDA, PBE, PBEsol, and AM05 functionals are not large enough to yield qualitatively different results. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3577999]

The role of electron localization in the atomic structure of transition-metal 13-atom clusters: the example of Co(13), Rh(13), and Hf(13)

The crystalline structure of transition-metals (TM) has been widely known for several decades, however, our knowledge on the atomic structure of TM clusters is still far from satisfactory, which compromises an atomistic understanding of the reactivity of TM clusters. For example, almost all density functional theory (DFT) calculations for TM clusters have been based on local (local density approximation-LDA) and semilocal (generalized gradient approximation-GGA) exchange-correlation functionals, however, it is well known that plain DFT fails to correct the self-interaction error, which affects the properties of several systems. To improve our basic understanding of the atomic and electronic properties of TM clusters, we report a DFT study within two nonlocal functionals, namely, the hybrid HSE (Heyd, Scuseria, and Ernzerhof) and GGA + U functionals, of the structural and electronic properties of the Co(13), Rh(13), and Hf(13) clusters. For Co(13) and Rh(13), we found that improved exchange-correlation functionals decrease the stability of open structures such as the hexagonal bilayer (HBL) and double simple-cubic (DSC) compared with the compact icosahedron (ICO) structure, however, DFT-GGA, DFT-GGA + U, and DFT-HSE yield very similar results for Hf(13). Thus, our results suggest that the DSC structure obtained by several plain DFT calculations for Rh(13) can be improved by the use of improved functionals. Using the sd hybridization analysis, we found that a strong hybridization favors compact structures, and hence, a correct description of the sd hybridization is crucial for the relative energy stability. For example, the sd hybridization decreases for HBL and DSC and increases for ICO in the case of Co(13) and Rh(13), while for Hf(13), the sd hybridization decreases for all configurations, and hence, it does not affect the relative stability among open and compact configurations.

Direct evidence of singlet molecular oxygen generation from peroxynitrate, a decomposition product of peroxynitrite

The decomposition of peroxynitrite to nitrite and dioxygen at neutral pH follows complex kinetics, compared to its isomerization to nitrate at low pH. Decomposition may involve radicals or proceed by way of the classical peracid decomposition mechanism. Peroxynitrite (ONOOH/ONOO(-)) decomposition has been proposed to involve formation of peroxynitrate (O(2)NOOH/O(2)NOO(-)) at neutral pH (D. Gupta, B. Harish, R. Kissner and W. H. Koppenol, Dalton Trans., 2009, DOI: 10.1039/b905535e, see accompanying paper in this issue). Peroxynitrate is unstable and decomposes to nitrite and dioxygen. This study aimed to investigate whether O(2)NOO(-) formed upon ONOOH/ONOO(-) decomposition generates singlet molecular oxygen [O(2) ((1)Delta(g))]. As unequivocally revealed by the measurement of monomol light emission in the near infrared region at 1270 nm and by chemical trapping experiments, the decomposition of ONOO(-) or O(2)NOOH at neutral to alkaline pH generates O(2) ((1)Delta(g)) at a yield of ca. 1% and 2-10%, respectively. Characteristic light emission, corresponding to O(2) ((1)Delta(g)) monomolecular decay was observed for ONOO(-) and for O(2)NOOH prepared by reaction of H(2)O(2) with NO(2)BF(4) and of H(2)O(2) with NO(2)(-) in HClO(4). The generation of O(2) ((1)Delta(g)) from ONOO(-) increased in a concentration-dependent manner in the range of 0.1-2.5 mM and was dependent on pH, giving a sigmoid pro. le with an apparent pK(a) around pD 8.1 (pH 7.7). Taken together, our results clearly identify the generation of O(2) ((1)Delta(g)) from peroxynitrate [O(2)NOO(-) -> NO(2)(-) + O(2) ((1)Delta(g))] generated from peroxynitrite and also from the reactions of H(2)O(2) with either NO(2)BF(4) or NO(2)(-) in acidic media.

Singlet molecular oxygen trapping by the fluorescent probe diethyl-3,3 '-(9,10-anthracenediyl)bisacrylate synthesized by the Heck reaction

Singlet molecular oxygen O(2)((1)Delta(g)) is a potent oxidant that can react with different biomolecules, including DNA, lipids and proteins. Many polycyclic aromatic hydrocarbons have been studied as O(2)((1)Delta(g)) chemical traps. Nevertheless, a suitable modification in the polycyclic aromatic ring must be made to increase the yield of O(2)((1)Delta(g)) chemical trapping. With this goal, an anthracene derivative, diethyl-3,3 '-(9,10-anthracenediyl)bisacrylate (DADB), was obtained from the reaction of 9,10-dibromoanthracene and ethyl acrylate through the Heck coupling reaction. The coupling of ethyl acrylate with the anthracene ring produced a new lipophilic, esterified, fluorescent probe reactive toward O(2)((1)Delta(g)). This compound reacts with O(2)((1)Delta(g)) at a rate of k(r) = 1.69 x 10(6) M(-1) s(-1) forming a stable endoperoxide (DADBO(2)), which was characterized by UV-Vis, fluorescence, HPLC/MS and (1)H and (13)C NMR techniques. The photophysical, photochemical and thermostability features of DADB were also evaluated. Furthermore, this compound has the potential for great application in biological systems because it is easily synthetized in large amount and generates specific endoperoxide (DADBO(2)), which can be easily detected by HPLC tandem mass spectrometry (HPLC/MS/MS).

Molecular dynamics simulation of liquid trimethylphosphine

Structural and dynamical properties of liquid trimethylphosphine (TMP), (CH(3))(3)P, as a function of temperature is investigated by molecular dynamics (MD) simulations. The force field used in the MD simulations, which has been proposed from molecular mechanics and quantum chemistry calculations, is able to reproduce the experimental density of liquid TMP at room temperature. Equilibrium structure is investigated by the usual radial distribution function, g(r), and also in the reciprocal space by the static structure factor, S(k). On the basis of center of mass distances, liquid TMP behaves like a simple liquid of almost spherical particles, but orientational correlation due to dipole-dipole interactions is revealed at short-range distances. Single particle and collective dynamics are investigated by several time correlation functions. At high temperatures, diffusion and reorientation occur at the same time range as relaxation of the liquid structure. Decoupling of these dynamic properties starts below ca. 220 K, when rattling dynamics of a given TMP molecules due to the cage effect of neighbouring molecules becomes important. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3624408]

Polarization effects in molecular dynamics simulations of glass-formers Ca(NO(3))(2)center dot nH(2)O, n=4, 6, and 8

Thermodynamics, equilibrium structure, and dynamics of glass-forming liquids Ca(NO(3))(2)center dot nH(2)O, n=4, 6, and 8, have been investigated by molecular dynamics (MD) simulations. A polarizable model was considered for H(2)O and NO(3)- on the basis of previous fluctuating charge models for pure water and the molten salt 2Ca(NO(3))(2)center dot 3KNO(3). Similar thermodynamic properties have been obtained with nonpolarizable and polarizable models. The glass transition temperature, T(g), estimated from MD simulations was dependent on polarization, in particular the dependence of T(g) with electrolyte concentration. Significant polarization effects on equilibrium structure were observed in cation-cation, cation-anion, and water-water structures. Polarization increases the diffusion coefficient of H(2)O, but does not change significantly the diffusion coefficients of ions. Viscosity decreases upon inclusion of polarization, but the conductivity calculated with the polarizable model is smaller than the nonpolarizable model because polarization enhances anion-cation interactions.

Full-dimensional analytical ab initio potential energy surface of the ground state of HOI

Extensive ab initio calculations using a complete active space second-order perturbation theory wavefunction, including scalar and spin-orbit relativistic effects with a quadruple-zeta quality basis set were used to construct an analytical potential energy surface (PES) of the ground state of the [H, O, I] system. A total of 5344 points were fit to a three-dimensional function of the internuclear distances, with a global root-mean-square error of 1.26 kcal mol(-1). The resulting PES describes accurately the main features of this system: the HOI and HIO isomers, the transition state between them, and all dissociation asymptotes. After a small adjustment, using a scaling factor on the internal coordinates of HOI, the frequencies calculated in this work agree with the experimental data available within 10 cm(-1). (C) 2011 American Institute of Physics. [doi: 10.1063/1.3615545]

High performance gold nanorods and silver nanocubes in surface-enhanced Raman spectroscopy of pesticides

The behavior of Au nanorods and Ag nanocubes as analytical sensors was evaluated for three different classes of herbicides. The use of such anisotropic nanoparticles in surface-enhanced Raman scattering (SERS) experiments allows the one to obtain the spectrum of crystal violet dye in the single molecule regime, as well as the pesticides dichlorophenoxyacetic acid (2,4-D), trichlorfon and ametryn. Such metallic substrates show high SERS performance at low analyte concentrations making them adequate for use as analytical sensors. Density functional theory (DFT) calculations of the geometries and vibrational wavenumbers of the adsorbates in the presence of silver or gold atoms were used to elucidate the nature of adsorbate-nanostructure bonding in each case and support the enhancement patterns observed in each SERS spectrum.

Fragility and glassy dynamics of 2Ca(NO(3))(2)center dot 3KNO(3) under pressure: Molecular dynamics simulations

Molecular dynamics simulations of the glass-forming liquid 2Ca(NO(3))(2)center dot 3KNO(3) (CKN) were performed from high temperature liquid states down to low temperature glassy states at six different pressures from 10(-4) to 5.0 GPa. The temperature dependence of the structural relaxation time indicates that the fragility of liquid CKN changes with pressure. In line with recent proposal [Scopigno , Science 302, 849 (2003)], the change on liquid fragility is followed by a proportional change of the nonergodicity factor of the corresponding glass at low temperature. (c) 2008 American Institute of Physics.