Lower Bounds on the Exchange-Correlation Energy in Reduced Dimensions

Bounds on the exchange-correlation energy of many-electron systems are derived and tested. By using universal scaling properties of the electron-electron interaction, we obtain the exponent of the bounds in three, two, one, and quasione dimensions. From the properties of the electron gas in the dilute regime, the tightest estimate to date is given for the numerical prefactor of the bound, which is crucial in practical applications. Numerical tests on various low-dimensional systems are in line with the bounds obtained and give evidence of an interesting dimensional crossover between two and one dimensions.

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.

Contribution of the second Landau level to the exchange energy of the three-dimensional electron gas in a high magnetic field

We derive a closed analytical expression for the exchange energy of the three-dimensional interacting electron gas in strong magnetic fields, which goes beyond the quantum limit (L=0) by explicitly including the effect of the second, L=1, Landau level and arbitrary spin polarization. The inclusion of the L=1 level brings the fields to which the formula applies closer to the laboratory range, as compared to previous expressions, valid only for L=0 and complete spin polarization. We identify and explain two distinct regimes separated by a critical density n(c). Below n(c), the per particle exchange energy is lowered by the contribution of L=1, whereas above n(c) it is increased. As special cases of our general equation we recover various known more limited results for higher fields, and we identify and correct a few inconsistencies in some of these earlier expressions.

Antiferromagnetism with spin polarization of GaN-based diluted magnetic semiconductors

The diluted magnetic semiconductors are promising materials for spintronic applications. Usually one intents to find the ferromagnetic state but recently the antiferromagnetism (AFM) was proposed to have some advantages. In this work, we verify the possibility to obtain spin polarization with an AFM state. In particular, we studied GaN 5% double doped with two different transition metals atoms (Mn and Co or Cr and Ni), forming the Mn(x)Co(0.056-x)Ga(0.944)N and Cr(x)Ni(0.056-x)Ga(0.944)N quaternary alloys. In order to simulate these systems in a more realistic way, and take into account composition fluctuations, we adapted the generalized quasichemical approach to diluted alloys, which is used in combination with spin density-functional theory. We find that is possible to obtain an AFM ground state up to 70% spin polarization.

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]

Computer-aided drug design and ADMET predictions for identification and evaluation of novel potential farnesyltransferase inhibitors in cancer therapy

We have used various computational methodologies including molecular dynamics, density functional theory, virtual screening, ADMET predictions and molecular interaction field studies to design and analyze four novel potential inhibitors of farnesyltransferase (FTase). Evaluation of two proposals regarding their drug potential as well as lead compounds have indicated them as novel promising FTase inhibitors, with theoretically interesting pharmacotherapeutic profiles, when Compared to the very active and most cited FTase inhibitors that have activity data reported, which are launched drugs or compounds in clinical tests. One of our two proposals appears to be a more promising drug candidate and FTase inhibitor, but both derivative molecules indicate potentially very good pharmacotherapeutic profiles in comparison with Tipifarnib and Lonafarnib, two reference pharmaceuticals. Two other proposals have been selected with virtual screening approaches and investigated by LIS, which suggest novel and alternatives scaffolds to design future potential FTase inhibitors. Such compounds can be explored as promising molecules to initiate a research protocol in order to discover novel anticancer drug candidates targeting farnesyltransferase, in the fight against cancer. (C) 2009 Elsevier Inc. All rights reserved.

Molecular dynamics, flexible docking, virtual screening, ADMET predictions, and molecular interaction field studies to design novel potential MAO-B inhibitors

Monoamine oxidase is a flavoenzyme bound to the mitochondrial outer membranes of the cells, which is responsible for the oxidative deamination of neurotransmitter and dietary amines. It has two distinct isozymic forms, designated MAO-A and MAO-B, each displaying different substrate and inhibitor specificities. They are the well-known targets for antidepressant, Parkinson`s disease, and neuroprotective drugs. Elucidation of the x-ray crystallographic structure of MAO-B has opened the way for the molecular modeling studies. In this work we have used molecular modeling, density functional theory with correlation, virtual screening, flexible docking, molecular dynamics, ADMET predictions, and molecular interaction field studies in order to design new molecules with potential higher selectivity and enzymatic inhibitory activity over MAO-B.

Evaluation of the enthalpy of formation, proton affinity, and gas-phase basicity of gamma-butyrolactone and 2-pyrrolidinone by isodesmic reactions

The knowledge of thermochemical parameters such as the enthalpy of formation, gas-phase basicity, and proton affinity may be the key to understanding molecular reactivity. The obtention of these thermochemical parameters by theoretical chemical models may be advantageous when experimental measurements are difficult to accomplish. The development of ab initio composite models represents a major advance in the obtention of these thermochemical parameters,. but these methods do not always lead to accurate values. Aiming at achieving a comparison between the ab initio models and the hybrid models based on the density functional theory (DFT), we have studied gamma-butyrolactone and 2-pyrrolidinone with a goal of obtaining high-quality thermochemical parameters using the composite chemical models G2, G2MP2, MP2, G3, CBS-Q, CBS-4, and CBS-QB3; the DFT methods B3LYP, B3P86, PW91PW91, mPW1PW, and B98; and the basis sets 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d), 6-311+G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ. Values obtained for the enthalpies of formation, proton affinity, and gas-phase basicity of the two target molecules were compared to the experimental data reported in the literature. The best results were achieved with the use of DFT models, and the B3LYP method led to the most accurate data.

Continuum, Discrete, and Explicit Solvation Models for Describing the Low-Lying Absorption Spectrum of the Pterin Acid in Aqueous Environment

The absorption spectrum of the acid form of pterin in water was investigated theoretically. Different procedures using continuum, discrete, and explicit models were used to include the solvation effect on the absorption spectrum, characterized by two bands. The discrete and explicit models used Monte Carlo simulation to generate the liquid structure and time-dependent density functional theory (B3LYP/6-31G+(d)) to obtain the excitation energies. The discrete model failed to give the correct qualitative effect on the second absorption band. The continuum model, in turn, has given a correct qualitative picture and a semiquantitative description. The explicit use of 29 solvent molecules, forming a hydration shell of 6 angstrom, embedded in the electrostatic field of the remaining solvent molecules, gives absorption transitions at 3.67 and 4.59 eV in excellent agreement with the S(0)-S(1) and S(0)-S(2) absorption bands at of 3.66 and 4.59 eV, respectively, that characterize the experimental spectrum of pterin in water environment. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110: 2371-2377, 2010

Many-body energy decomposition of hydrogen-bonded glycine clusters in gas-phase

A detailed analysis of the many-body contribution to the interaction energies of the gas-phase hydrogen-bonded glycine clusters, (Gly)(N), N = 1-4 is presented. The energetics of the hydrogen-bonded dimer, trimer and tetramer complexes have been analyzed using density-functional theory. The magnitude of the two-through four-body energy terms have been calculated and compared. The relaxation energy and the two-body energy terms are the principal contributors to the total binding energy. Four-body contribution is negligible. However, the three-body contribution is found to be sizable and the formation of the cyclic glycine trimer presents geometric strains that make it less favorable. (C) 2010 Elsevier B.V. All rights reserved.

Ab initio study of the EFG at the N sites in imidazole

We study the nuclear quadrupole interaction at the nitrogen sites in the molecular and crystalline phases of the imidazole compound. We use PAW which is a state-of-the-art method to calculate the electronic structure and electric field gradient at the nucleus in the framework of the density functional theory. The quadrupole frequencies at both imino and amino N sites are in excellent agreement with measurements. This is the first time that the electric field gradient at crystalline imidazole is correctly treated by an ab initio theoretical approach.

Electronic transport properties of ascorbic acid and nicotinamide adsorbed on single-walled carbon nanotubes

Ab initio simulations of carbon nanotubes interacting with ascorbic acid and nicotinamide are reported. The electronic transport properties of these systems are studied using a combination of density functional theory and non-equilibrium Green`s functions methods. The adsorptions of both molecules are observed to depend strongly on their functionalization. The interaction through the appropriate functionalized species modifies the structural and electronic properties of the original system, resulting in a chemisorption regime. Changes in the electronic transport properties are also observed, with reductions on the total electronic transmission probabilities. Nevertheless, when the molecules interact through the pristine form, a physisorption interaction is observed with insignificant structural and electronic transport changes. (c) 2011 Elsevier B.V. All rights reserved.

Driving forces for the adsorption of cyclopentene on InP(001)

In this work the interaction of cyclopentene with a set of InP(001) surfaces is investigated by means of the density functional theory. We propose a simple approach for evaluating the surface strain and based on it we have found a linear relation between bond and strain energies and the adsorption energy. Our results also indicate that the higher the bond energy, the more disperse the charge distribution is around the adsorption site associated to the high occupied state, a key feature that characterizes the adsorption process. Different adsorption coverages are used to evaluate the proposed equation. Our results suggest that the proposed approach might be extended to other systems where the interaction of the semiconductor surface and the molecule is restricted to first neighbor sites. (C) 2011 Elsevier B.V. All rights reserved.

IxV Curves of Boron and Nitrogen Doping Zigzag Graphene Nanoribbons

We investigate the transport properties (IxV curves and zero bias transmittance) of pristine graphene nanoribbons (GNRs) as well as doped with boron and nitrogen using an approach that combines nonequilibrium Green`s functions and density functional theory (DFT) [NEGF-DFT]. Even for a pristine nanoribbon we verify a spin-filter effect under finite bias voltage when the leads have an antiparallel magnetization. The presence of the impurities at the edges of monohydrogenated zigzag GNRs changes dramatically the charge transport properties inducing a spin-polarized conductance. The IxV curves for these systems show that depending on the bias voltage the spin polarization can be inverted. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 111: 1379-1386, 2011