Spin gaps and spin-flip energies in density-functional theory

Energy gaps are crucial aspects of the electronic structure of finite and extended systems. Whereas much is known about how to define and calculate charge gaps in density-functional theory (DFT), and about the relation between these gaps and derivative discontinuities of the exchange-correlation functional, much less is known about spin gaps. In this paper we give density-functional definitions of spin-conserving gaps, spin-flip gaps and the spin stiffness in terms of many-body energies and in terms of single-particle (Kohn-Sham) energies. Our definitions are as analogous as possible to those commonly made in the charge case, but important differences between spin and charge gaps emerge already on the single-particle level because unlike the fundamental charge gap spin gaps involve excited-state energies. Kohn-Sham and many-body spin gaps are predicted to differ, and the difference is related to derivative discontinuities that are similar to, but distinct from, those usually considered in the case of charge gaps. Both ensemble DFT and time-dependent DFT (TDDFT) can be used to calculate these spin discontinuities from a suitable functional. We illustrate our findings by evaluating our definitions for the Lithium atom, for which we calculate spin gaps and spin discontinuities by making use of near-exact Kohn-Sham eigenvalues and, independently, from the single-pole approximation to TDDFT. The many-body corrections to the Kohn-Sham spin gaps are found to be negative, i.e., single-particle calculations tend to overestimate spin gaps while they underestimate charge gaps.

Density-functional description of superconducting and magnetic proximity effects across a tunneling barrier

A density-functional formalism for superconductivity and magnetism is presented. The resulting relations unify previously derived Kohn-Sham equations for superconductors and for noncollinear magnetism. The formalism, which discriminates Cooper-pair singlets from triplets, is applied to two quantum liquids coupled by tunneling through a barrier. An exact expression is derived, relating the eigenstates and eigenvalues of the Kohn-Sham equations, unperturbed by tunneling, on one side of the barrier to the proximity-induced ordering potential on the other.

Density functional theory and ab initio molecular dynamics study of NO adsorption on Pd(111) and Pt(111) surfaces

The origin of the unique geometry for nitric oxide (NO) adsorption on Pd(111) and Pt(111) surfaces as well as the effect of temperature were studied by density functional theory calculations and ab initio molecular dynamics at finite temperature. We found that at low coverage, the adsorption geometry is determined by electronic interactions, depending sensitively on the adsorption sites and coverages, and the effect of temperature on geometries is significant. At coverage of 0.25 monolayer (ML), adsorbed NO at hollow sites prefer an upright configuration, while NO adsorbed at top sites prefer a tilting configuration. With increase in the coverage up to 0.50 ML, the enhanced steric repulsion lead to the tilting of hollow NO. We found that the tilting was enhanced by the thermal effects. At coverage of 0.75 ML with p(2 x 2)-3NO(fcc+hcp+top) structure, we found that there was no preferential orientation for tilted top NO. The interplay of the orbital hybridization, thermal effects, steric repulsion, and their effects on the adsorption geometries were highlighted at the end.

Theory of nitride oxide adsorption on transition metal (111) surfaces: a first-principles investigation

In this work, we report a density functional theory study of nitric oxide (NO) adsorption on close-packed transition metal (TM) Rh(111), Ir(111), Pd(111) and Pt(111) surfaces in terms of adsorption sites, binding mechanism and charge transfer at a coverage of Theta(NO) = 0.25, 0.50, 0.75 monolayer (ML). Based on our study, an unified picture for the interaction between NO and TM(111) and site preference is established, and valuable insights are obtained. At low coverage (0.25 ML), we find that the interaction of NO/TM(111) is determined by an electron donation and back-donation process via the interplay between NO 5 sigma/2 pi* and TM d-bands. The extent of the donation and back-donation depends critically on the coordination number (adsorption sites) and TM d-band filling, and plays an essential role for NO adsorption on TM surfaces. DFT calculations shows that for TMs with high d-band filling such as Pd and Pt, hollow-site NO is energetically the most favorable, and top-site NO prefers to tilt away from the normal direction. While for TMs with low d-band filling (Rh and Ir), top-site NO perpendicular to the surfaces is energetically most favorable. Electronic structure analysis show that irrespective of the TM and adsorption site, there is a net charge transfer from the substrate to the adsorbate due to overwhelming back-donation from the TM substrate to the adsorbed NO molecules. The adsorption-induced change of the work function with respect to bare surfaces and dipole moment is however site dependent, and the work function increases for hollow-site NO, but decreases for top-site NO, because of differences in the charge redistribution. The interplay between the energetics, lateral interaction and charge transfer, which is element dependent, rationalizes the structural evolution of NO adsorption on TM(111) surfaces in the submonolayer regime.

Exact nonequilibrium work generating function for a small classical system

We obtain the exact nonequilibrium work generating function (NEWGF) for a small system consisting of a massive Brownian particle connected to internal and external springs. The external work is provided to the system for a finite-time interval. The Jarzynski equality, obtained in this case directly from the NEWGF, is shown to be valid for the present model, in an exact way regardless of the rate of external work.

TESTING STATISTICAL HYPOTHESIS ON RANDOM TREES AND APPLICATIONS TO THE PROTEIN CLASSIFICATION PROBLEM

Efficient automatic protein classification is of central importance in genomic annotation. As an independent way to check the reliability of the classification, we propose a statistical approach to test if two sets of protein domain sequences coming from two families of the Pfam database are significantly different. We model protein sequences as realizations of Variable Length Markov Chains (VLMC) and we use the context trees as a signature of each protein family. Our approach is based on a Kolmogorov-Smirnov-type goodness-of-fit test proposed by Balding et at. [Limit theorems for sequences of random trees (2008), DOI: 10.1007/s11749-008-0092-z]. The test statistic is a supremum over the space of trees of a function of the two samples; its computation grows, in principle, exponentially fast with the maximal number of nodes of the potential trees. We show how to transform this problem into a max-flow over a related graph which can be solved using a Ford-Fulkerson algorithm in polynomial time on that number. We apply the test to 10 randomly chosen protein domain families from the seed of Pfam-A database (high quality, manually curated families). The test shows that the distributions of context trees coming from different families are significantly different. We emphasize that this is a novel mathematical approach to validate the automatic clustering of sequences in any context. We also study the performance of the test via simulations on Galton-Watson related processes.

Multiclass Hammersley-Aldous-Diaconis process and multiclass-customer queues

In the Hammersley-Aldous-Diaconis process, infinitely many particles sit in R and at most one particle is allowed at each position. A particle at x, whose nearest neighbor to the right is at y, jumps at rate y - x to a position uniformly distributed in the interval (x, y). The basic coupling between trajectories with different initial configuration induces a process with different classes of particles. We show that the invariant measures for the two-class process can be obtained as follows. First, a stationary M/M/1 queue is constructed as a function of two homogeneous Poisson processes, the arrivals with rate, and the (attempted) services with rate rho > lambda Then put first class particles at the instants of departures (effective services) and second class particles at the instants of unused services. The procedure is generalized for the n-class case by using n - 1 queues in tandem with n - 1 priority types of customers. A multi-line process is introduced; it consists of a coupling (different from Liggett's basic coupling), having as invariant measure the product of Poisson processes. The definition of the multi-line process involves the dual points of the space-time Poisson process used in the graphical construction of the reversed process. The coupled process is a transformation of the multi-line process and its invariant measure is the transformation described above of the product measure.

Can we make a Bohmian electron reach the speed of light, at least for one instant?

In Bohmian mechanics, a version of quantum mechanics that ascribes world lines to electrons, we can meaningfully ask about an electron's instantaneous speed relative to a given inertial frame. Interestingly, according to the relativistic version of Bohmian mechanics using the Dirac equation, a massive particle's speed is less than or equal to the speed of light, but not necessarily less. That is, there are situations in which the particle actually reaches the speed of light-a very nonclassical behavior. That leads us to the question of whether such situations can be arranged experimentally. We prove a theorem, Theorem 5, implying that for generic initial wave functions the probability that the particle ever reaches the speed of light, even if at only one point in time, is zero. We conclude that the answer to the question is no. Since a trajectory reaches the speed of light whenever the quantum probability current (psi) over bar gamma(mu)psi is a lightlike 4-vector, our analysis concerns the current vector field of a generic wave function and may thus be of interest also independently of Bohmian mechanics. The fact that the current is never spacelike has been used to argue against the possibility of faster-than-light tunneling through a barrier, a somewhat similar question. Theorem 5, as well as a more general version provided by Theorem 6, are also interesting in their own right. They concern a certain property of a function psi : R(4) -> C(4) that is crucial to the question of reaching the speed of light, namely being transverse to a certain submanifold of C(4) along a given compact subset of space-time. While it follows from the known transversality theorem of differential topology that this property is generic among smooth functions psi : R(4) -> C(4), Theorem 5 asserts that it is also generic among smooth solutions of the Dirac equation. (C) 2010 American Institute of Physics. [doi:10.1063/1.3520529]

Identification of archaeal proteins that affect the exosome function in vitro

Background: The archaeal exosome is formed by a hexameric RNase PH ring and three RNA binding subunits and has been shown to bind and degrade RNA in vitro. Despite extensive studies on the eukaryotic exosome and on the proteins interacting with this complex, little information is yet available on the identification and function of archaeal exosome regulatory factors. Results: Here, we show that the proteins PaSBDS and PaNip7, which bind preferentially to poly-A and AU-rich RNAs, respectively, affect the Pyrococcus abyssi exosome activity in vitro. PaSBDS inhibits slightly degradation of a poly-rA substrate, while PaNip7 strongly inhibits the degradation of poly-A and poly-AU by the exosome. The exosome inhibition by PaNip7 appears to depend at least partially on its interaction with RNA, since mutants of PaNip7 that no longer bind RNA, inhibit the exosome less strongly. We also show that FITC-labeled PaNip7 associates with the exosome in the absence of substrate RNA. Conclusions: Given the high structural homology between the archaeal and eukaryotic proteins, the effect of archaeal Nip7 and SBDS on the exosome provides a model for an evolutionarily conserved exosome control mechanism.

The Essential Nucleolar Yeast Protein Nop8p Controls the Exosome Function during 60S Ribosomal Subunit Maturation

The yeast nucleolar protein Nop8p has previously been shown to interact with Nip7p and to be required for 60S ribosomal subunit formation. Although depletion of Nop8p in yeast cells leads to premature degradation of rRNAs, the biochemical mechanism responsible for this phenotype is still not known. In this work, we show that the Nop8p amino-terminal region mediates interaction with the 5.8S rRNA, while its carboxyl-terminal portion interacts with Nip7p and can partially complement the growth defect of the conditional mutant strain Dnop8/GAL:NOP8. Interestingly, Nop8p mediates association of Nip7p to pre-ribosomal particles. Nop8p also interacts with the exosome subunit Rrp6p and inhibits the complex activity in vitro, suggesting that the decrease in 60S ribosomal subunit levels detected upon depletion of Nop8p may result from degradation of pre-rRNAs by the exosome. These results strongly indicate that Nop8p may control the exosome function during pre-rRNA processing.

Structure-Function Analysis of the HrpB2-HrcU Interaction in the Xanthomonas citri Type III Secretion System

Bacterial type III secretion systems deliver protein virulence factors to host cells. Here we characterize the interaction between HrpB2, a small protein secreted by the Xanthomonas citri subsp. citri type III secretion system, and the cytosolic domain of the inner membrane protein HrcU, a paralog of the flagellar protein FlhB. We show that a recombinant fragment corresponding to the C-terminal cytosolic domain of HrcU produced in E. coli suffers cleavage within a conserved Asn264-Pro265-Thr266-His267 (NPTH) sequence. A recombinant HrcU cytosolic domain with N264A, P265A, T266A mutations at the cleavage site (HrcU(AAAH)) was not cleaved and interacted with HrpB2. Furthermore, a polypeptide corresponding to the sequence following the NPTH cleavage site also interacted with HrpB2 indicating that the site for interaction is located after the NPTH site. Non-polar deletion mutants of the hrcU and hrpB2 genes resulted in a total loss of pathogenicity in susceptible citrus plants and disease symptoms could be recovered by expression of HrpB2 and HrcU from extrachromossomal plasmids. Complementation of the Delta hrcU mutant with HrcU(AAAH) produced canker lesions similar to those observed when complemented with wild-type HrcU. HrpB2 secretion however, was significantly reduced in the Delta hrcU mutant complemented with HrcU(AAAH), suggesting that an intact and cleavable NPTH site in HrcU is necessary for total functionally of T3SS in X. citri subsp. citri. Complementation of the Delta hrpB2 X. citri subsp. citri strain with a series of hrpB2 gene mutants revealed that the highly conserved HrpB2 C-terminus is essential for T3SS-dependent development of citrus canker symptoms in planta.

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]

Thermo-solvatochromism in binary mixtures of water and ionic liquids: on the relative importance of solvophobic interactions

The thermo-solvatochromism of 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr(2), has been studied in mixtures of water, W, with ionic liquids, ILs, in the temperature range of 10 to 60 degrees C, where feasible. The objectives of the study were to test the applicability of a recently introduced solvation model, and to assess the relative importance of solute-solvent solvophobic interactions. The ILs were 1-allyl-3-alkylimidazolium chlorides, where the alkyl groups are methyl, 1-butyl, and 1-hexyl, respectively. The equilibrium constants for the interaction of W and the ILs were calculated from density data; they were found to be linearly dependent on N(C), the number of carbon atoms of the alkyl group; van't Hoff equation (log K versus 1/T) applied satisfactorily. Plots of the empirical solvent polarities, E(T) (MePMBr(2)) in kcal mol(-1), versus the mole fraction of water in the binary mixture, chi(w), showed non-linear, i.e., non-ideal behavior. The dependence of E(T) (MePMBr(2)) on chi(w), has been conveniently quantified in terms of solvation by W, IL, and the ""complex"" solvent IL-W. The non-ideal behavior is due to preferential solvation by the IL and, more efficiently, by IL-W. The deviation from linearity increases as a function of increasing N(C) of the IL, and is stronger than that observed for solvation of MePMBr(2) by aqueous 1-propanol, a solvent whose lipophilicity is 12.8 to 52.1 times larger than those of the ILs investigated. The dependence on N(C) is attributed to solute-solvent solvophobic interactions, whose relative contribution to solvation are presumably greater than that in mixtures of water and 1-propanol.

Structure, stability, depolarized light scattering, and vibrational spectra of fullerenols from all-electron density-functional-theory calculations

We have investigated the stability, electronic properties, Rayleigh (elastic), and Raman (inelastic) depolarization ratios, infrared and Raman absorption vibrational spectra of fullerenols [C(60)(OH)(n)] with different degrees of hydroxylation by using all-electron density-functional-theory (DFT) methods. Stable arrangements of these molecules were found by means of full geometry optimizations using Becke's three-parameter exchange functional with the Lee, Yang, and Parr correlation functional. This DFT level has been combined with the 6-31G(d,p) Gaussian-type basis set, as a compromise between accuracy and capability to treat highly hydroxylated fullerenes, e.g., C(60)(OH)(36). Thus, the molecular properties of fullerenols were systematically analyzed for structures with n=1, 2, 3, 4, 8, 10, 16, 18, 24, 32, and 36. From the electronic structure analysis of these molecules, we have evidenced an important effect related to the weak chemical reactivity of a possible C(60)(OH)(24) isomer. To investigate Raman scattering and the vibrational spectra of the different fullerenols, frequency calculations are carried out within the harmonic approximation. In this case a systematic study is only performed for n=1-4, 8, 10, 16, 18, and 24. Our results give good agreements with the expected changes in the spectral absorptions due to the hydroxylation of fullerenes.

Excitation energies from ground-state density-functionals by means of generator coordinates

The generator-coordinate method is a flexible and powerful reformulation of the variational principle. Here we show that by introducing a generator coordinate in the Kohn-Sham equation of density-functional theory, excitation energies can be obtained from ground-state density functionals. As a viability test, the method is applied to ground-state energies and various types of excited-state energies of atoms and ions from the He and the Li isoelectronic series. Results are compared to a variety of alternative DFT-based approaches to excited states, in particular time-dependent density-functional theory with exact and approximate potentials.