Science.gov

Sample records for combinatorial density functional

  1. Combinatorial density functional theory-based screening of surface alloys for the oxygen reduction reaction.

    SciTech Connect

    Greeley, J.; Norskov, J.; Center for Nanoscale Materials; Technical Univ. of Denmark

    2009-03-26

    A density functional theory (DFT) -based, combinatorial search for improved oxygen reduction reaction (ORR) catalysts is presented. A descriptor-based approach to estimate the ORR activity of binary surface alloys, wherein alloying occurs only in the surface layer, is described, and rigorous, potential-dependent computational tests of the stability of these alloys in aqueous, acidic environments are presented. These activity and stability criteria are applied to a database of DFT calculations on nearly 750 binary transition metal surface alloys; of these, many are predicted to be active for the ORR but, with few exceptions, they are found to be thermodynamically unstable in the acidic environments typical of low-temperature fuel cells. The results suggest that, absent other thermodynamic or kinetic mechanisms to stabilize the alloys, surface alloys are unlikely to serve as useful ORR catalysts over extended periods of operation.

  2. Fuzzy electron density fragments in macromolecular quantum chemistry, combinatorial quantum chemistry, functional group analysis, and shape-activity relations.

    PubMed

    Mezey, Paul G

    2014-09-16

    Conspectus Just as complete molecules have no boundaries and have "fuzzy" electron density clouds approaching zero density exponentially at large distances from the nearest nucleus, a physically justified choice for electron density fragments exhibits similar behavior. Whereas fuzzy electron densities, just as any fuzzy object, such as a thicker cloud on a foggy day, do not lend themselves to easy visualization, one may partially overcome this by using isocontours. Whereas a faithful representation of the complete fuzzy density would need infinitely many such isocontours, nevertheless, by choosing a selected few, one can still obtain a limited pictorial representation. Clearly, such images are of limited value, and one better relies on more complete mathematical representations, using, for example, density matrices of fuzzy fragment densities. A fuzzy density fragmentation can be obtained in an exactly additive way, using the output from any of the common quantum chemical computational techniques, such as Hartree-Fock, MP2, and various density functional approaches. Such "fuzzy" electron density fragments properly represented have proven to be useful in a rather wide range of applications, for example, (a) using them as additive building blocks leading to efficient linear scaling macromolecular quantum chemistry computational techniques, (b) the study of quantum chemical functional groups, (c) using approximate fuzzy fragment information as allowed by the holographic electron density theorem, (d) the study of correlations between local shape and activity, including through-bond and through-space components of interactions between parts of molecules and relations between local molecular shape and substituent effects, (e) using them as tools of density matrix extrapolation in conformational changes, (f) physically valid averaging and statistical distribution of several local electron densities of common stoichiometry, useful in electron density databank mining, for

  3. Neurotransmitters drive combinatorial multistate postsynaptic density networks.

    PubMed

    Coba, Marcelo P; Pocklington, Andrew J; Collins, Mark O; Kopanitsa, Maksym V; Uren, Rachel T; Swamy, Sajani; Croning, Mike D R; Choudhary, Jyoti S; Grant, Seth G N

    2009-01-01

    The mammalian postsynaptic density (PSD) comprises a complex collection of approximately 1100 proteins. Despite extensive knowledge of individual proteins, the overall organization of the PSD is poorly understood. Here, we define maps of molecular circuitry within the PSD based on phosphorylation of postsynaptic proteins. Activation of a single neurotransmitter receptor, the N-methyl-D-aspartate receptor (NMDAR), changed the phosphorylation status of 127 proteins. Stimulation of ionotropic and metabotropic glutamate receptors and dopamine receptors activated overlapping networks with distinct combinatorial phosphorylation signatures. Using peptide array technology, we identified specific phosphorylation motifs and switching mechanisms responsible for the integration of neurotransmitter receptor pathways and their coordination of multiple substrates in these networks. These combinatorial networks confer high information-processing capacity and functional diversity on synapses, and their elucidation may provide new insights into disease mechanisms and new opportunities for drug discovery. PMID:19401593

  4. ωB97M-V: A combinatorially optimized, range-separated hybrid, meta-GGA density functional with VV10 nonlocal correlation.

    PubMed

    Mardirossian, Narbe; Head-Gordon, Martin

    2016-06-01

    A combinatorially optimized, range-separated hybrid, meta-GGA density functional with VV10 nonlocal correlation is presented. The final 12-parameter functional form is selected from approximately 10 × 10(9) candidate fits that are trained on a training set of 870 data points and tested on a primary test set of 2964 data points. The resulting density functional, ωB97M-V, is further tested for transferability on a secondary test set of 1152 data points. For comparison, ωB97M-V is benchmarked against 11 leading density functionals including M06-2X, ωB97X-D, M08-HX, M11, ωM05-D, ωB97X-V, and MN15. Encouragingly, the overall performance of ωB97M-V on nearly 5000 data points clearly surpasses that of all of the tested density functionals. In order to facilitate the use of ωB97M-V, its basis set dependence and integration grid sensitivity are thoroughly assessed, and recommendations that take into account both efficiency and accuracy are provided. PMID:27276948

  5. ωB97M-V: A combinatorially optimized, range-separated hybrid, meta-GGA density functional with VV10 nonlocal correlation

    NASA Astrophysics Data System (ADS)

    Mardirossian, Narbe; Head-Gordon, Martin

    2016-06-01

    A combinatorially optimized, range-separated hybrid, meta-GGA density functional with VV10 nonlocal correlation is presented. The final 12-parameter functional form is selected from approximately 10 × 109 candidate fits that are trained on a training set of 870 data points and tested on a primary test set of 2964 data points. The resulting density functional, ωB97M-V, is further tested for transferability on a secondary test set of 1152 data points. For comparison, ωB97M-V is benchmarked against 11 leading density functionals including M06-2X, ωB97X-D, M08-HX, M11, ωM05-D, ωB97X-V, and MN15. Encouragingly, the overall performance of ωB97M-V on nearly 5000 data points clearly surpasses that of all of the tested density functionals. In order to facilitate the use of ωB97M-V, its basis set dependence and integration grid sensitivity are thoroughly assessed, and recommendations that take into account both efficiency and accuracy are provided.

  6. Some functional metrics in algebraic and combinatorial coding

    NASA Astrophysics Data System (ADS)

    Choen, G.

    1980-06-01

    Three approaches to coding problems can be systematically distinguished: probabilistic (essentially existential), algebraic, and combinatorial. This last approach searches for optimal configurations and relegates to the second order, the problems of complexity related to decoding. Enumeration, graphs, designs, and the extreme theory of groups are used. The optimization of a functional metric was used with the combinatorial approach in order to define the space considered and the distance. The codes then become particular groups of the metric space, which is defined by parameters such as length, number of words, and capacity for correction. Some of these parameters are imposed.

  7. Combinatorial function of ETS transcription factors in the developing vasculature

    PubMed Central

    Pham, Van N.; Lawson, Nathan D.; Mugford, Joshua W.; Dye, Louis; Castranova, Daniel; Lo, Brigid; Weinstein, Brant M.

    2007-01-01

    Members of the ETS family of transcription factors are among the first genes expressed in the developing vasculature, but loss-of-function experiments for individual ETS factors in mice have not uncovered important early functional roles for these genes. However, multiple ETS factors are expressed in spatially and temporally overlapping patterns in the developing vasculature, suggesting possible functional overlap. We have taken a comprehensive approach to exploring the function of these factors during vascular development by employing the genetic and experimental tools available in the zebrafish to analyze four ETS family members expressed together in the zebrafish vasculature; fli1, fli1b, ets1, and etsrp. We isolated and characterized an ENU-induced mutant with defects in trunk angiogenesis and positionally cloned the defective gene from this mutant, etsrp. Using the etsrp morpholinos targeting each of the four genes, we show that the four ETS factors function combinatorially during vascular and hematopoietic development. Reduction of etsrp or any of the other genes alone results in either partial or no defects in endothelial differentiation, while combined reduction in the function of all four genes causes dramatic loss of endothelial cells. Our results demonstrate that combinatorial ETS factor function is essential for early endothelial specification and differentiation. PMID:17125762

  8. One and two quasiparticle state densities in the ESM: Combinatorial approach vs. exact results

    SciTech Connect

    Noy, R.C.; Silvera, M.P.; Martinez, R.P.

    1994-12-31

    An exact method proposed by Zhang and Yang to calculate the effect of Pauli exclusion principle in particle-hole state densities is used to derive closed formulas for one and two quasiparticle state densities. Analytical results are tested against exact combinatorial calculation in the framework of BCS theory using equidistant shell model spectra. Relationship between particle-hole and quasiparticle state densities is discussed.

  9. Combinatorial Gene Regulatory Functions Underlie Ultraconserved Elements in Drosophila.

    PubMed

    Warnefors, Maria; Hartmann, Britta; Thomsen, Stefan; Alonso, Claudio R

    2016-09-01

    Ultraconserved elements (UCEs) are discrete genomic elements conserved across large evolutionary distances. Although UCEs have been linked to multiple facets of mammalian gene regulation their extreme evolutionary conservation remains largely unexplained. Here, we apply a computational approach to investigate this question in Drosophila, exploring the molecular functions of more than 1,500 UCEs shared across the genomes of 12 Drosophila species. Our data indicate that Drosophila UCEs are hubs for gene regulatory functions and suggest that UCE sequence invariance originates from their combinatorial roles in gene control. We also note that the gene regulatory roles of intronic and intergenic UCEs (iUCEs) are distinct from those found in exonic UCEs (eUCEs). In iUCEs, transcription factor (TF) and epigenetic factor binding data strongly support iUCE roles in transcriptional and epigenetic regulation. In contrast, analyses of eUCEs indicate that they are two orders of magnitude more likely than the expected to simultaneously include protein-coding sequence, TF-binding sites, splice sites, and RNA editing sites but have reduced roles in transcriptional or epigenetic regulation. Furthermore, we use a Drosophila cell culture system and transgenic Drosophila embryos to validate the notion of UCE combinatorial regulatory roles using an eUCE within the Hox gene Ultrabithorax and show that its protein-coding region also contains alternative splicing regulatory information. Taken together our experiments indicate that UCEs emerge as a result of combinatorial gene regulatory roles and highlight common features in mammalian and insect UCEs implying that similar processes might underlie ultraconservation in diverse animal taxa. PMID:27247329

  10. Combinatorial Gene Regulatory Functions Underlie Ultraconserved Elements in Drosophila

    PubMed Central

    Warnefors, Maria; Hartmann, Britta; Thomsen, Stefan; Alonso, Claudio R.

    2016-01-01

    Ultraconserved elements (UCEs) are discrete genomic elements conserved across large evolutionary distances. Although UCEs have been linked to multiple facets of mammalian gene regulation their extreme evolutionary conservation remains largely unexplained. Here, we apply a computational approach to investigate this question in Drosophila, exploring the molecular functions of more than 1,500 UCEs shared across the genomes of 12 Drosophila species. Our data indicate that Drosophila UCEs are hubs for gene regulatory functions and suggest that UCE sequence invariance originates from their combinatorial roles in gene control. We also note that the gene regulatory roles of intronic and intergenic UCEs (iUCEs) are distinct from those found in exonic UCEs (eUCEs). In iUCEs, transcription factor (TF) and epigenetic factor binding data strongly support iUCE roles in transcriptional and epigenetic regulation. In contrast, analyses of eUCEs indicate that they are two orders of magnitude more likely than the expected to simultaneously include protein-coding sequence, TF-binding sites, splice sites, and RNA editing sites but have reduced roles in transcriptional or epigenetic regulation. Furthermore, we use a Drosophila cell culture system and transgenic Drosophila embryos to validate the notion of UCE combinatorial regulatory roles using an eUCE within the Hox gene Ultrabithorax and show that its protein-coding region also contains alternative splicing regulatory information. Taken together our experiments indicate that UCEs emerge as a result of combinatorial gene regulatory roles and highlight common features in mammalian and insect UCEs implying that similar processes might underlie ultraconservation in diverse animal taxa. PMID:27247329

  11. A Code For Combinatorial Calculation of Level Densities.

    Energy Science and Technology Software Center (ESTSC)

    1990-11-26

    Version 00 ICAR performs state and level density calculations for a specified type of nucleon. Calculations can be carried out in terms of the free gas model or in the frame of the BCS formalism to account for the pairing interaction.

  12. New optical field and its Wigner function obtained by partial tracing over one- and two-mode combinatorial squeezed state

    NASA Astrophysics Data System (ADS)

    Wang, Tong-Tong; Fan, Hong-Yi

    2016-12-01

    Based on the one- and two-mode combinatorial squeezed state (H.Y. Fan, Phys. Rev. A. 41(3), 1526 (1990))which can enhance squeezing effect, we derive a new optical field by using partial tracing method, we not only obtain its density operator but also deduce its Wigner function by virtue of operators' Weyl ordering property. This new photon field possesses more photon numbers than the corresponding chaotic field, and can be applied to quantum controlling and quantum information processing.

  13. Partition density functional theory

    NASA Astrophysics Data System (ADS)

    Nafziger, Jonathan

    Partition density functional theory (PDFT) is a method for dividing a molecular electronic structure calculation into fragment calculations. The molecular density and energy corresponding to Kohn Sham density-functional theory (KS-DFT) may be exactly recovered from these fragments. Each fragment acts as an isolated system except for the influence of a global one-body 'partition' potential which deforms the fragment densities. In this work, the developments of PDFT are put into the context of other fragment-based density functional methods. We developed three numerical implementations of PDFT: One within the NWChem computational chemistry package using basis sets, and the other two developed from scratch using real-space grids. It is shown that all three of these programs can exactly reproduce a KS-DFT calculation via fragment calculations. The first of our in-house codes handles non-interacting electrons in arbitrary one-dimensional potentials with any number of fragments. This code is used to explore how the exact partition potential changes for different partitionings of the same system and also to study features which determine which systems yield non-integer PDFT occupations and which systems are locked into integer PDFT occupations. The second in-house code, CADMium, performs real-space calculations of diatomic molecules. Features of the exact partition potential are studied for a variety of cases and an analytical formula determining singularities in the partition potential is derived. We introduce an approximation for the non-additive kinetic energy and show how this quantity can be computed exactly. Finally a PDFT functional is developed to address the issues of static correlation and delocalization errors in approximations within DFT. The functional is applied to the dissociation of H2 + and H2.

  14. Combinatorial atmospheric pressure chemical vapor deposition (cAPCVD): a route to functional property optimization.

    PubMed

    Kafizas, Andreas; Parkin, Ivan P

    2011-12-21

    We demonstrate how combinatorial atmospheric pressure chemical vapor deposition (cAPCVD) can be used as a synthetic tool for rapidly optimizing the functional properties of thin-films, by analyzing the self-cleaning properties of tungsten doped anatase as an example. By introducing reagents at separate points inside the reactor, a tungsten/titanium compositional gradient was formed and a diverse range of film growth conditions were obtained. By partially mixing the metal sources, a combinatorial film with a compositional profile that varied primarily in the lateral plane was synthesized. A combinatorial thin-film of anatase TiO(2) doped with an array of tungsten levels as a solid solution ranging from 0.38-13.8 W/Ti atom % was formed on a single glass substrate. The compositional-functional relationships were understood through comprehensively analyzing combinatorial phase space, with 200 positions investigated by high-throughput methods in this study. Physical and functional properties, and their compositional dependencies, were intercorrelated. It was found that increases in photocatalytic activity and conductivity were most highly dependent on film crystallinity within the 0.38-13.8 atom % W/Ti doping regime. However, enhancements in photoinduced surface wetting were primarily dependent on increases in preferred growth in the (211) crystal plane. PMID:22050427

  15. A drop-on-demand ink-jet printer for combinatorial libraries and functionally graded ceramics.

    PubMed

    Mohebi, Mohammad Masoud; Evans, Julian R G

    2002-01-01

    A printer has been designed and built for the preparation of combinatorial libraries of ceramics and for solid freeforming of functionally graded ceramics with three-dimensionally programmable spatial variation in composition. Several ceramic suspensions (as inks) can be subjected to micromixing behind the nozzle and printed at precise positions. Both mixing and positioning are computer-controlled. The machine consists of an XY table to control the geometry, a set of electromagnetic valves that manage the mixing, a combined electromagnetic valve and sapphire nozzle that form the print head, and a computer that controls the whole system. The mixing valves can eject as little as 1 mg/s ink into the mixing chamber. The printer has been controlled, run, calibrated and tested; the composition and geometry of printed mixtures can be controlled precisely. This method for the controlled mixing of powders facilitates the advance of combinatorial methods within the materials sciences. PMID:12099843

  16. A combinatorial optogenetic approach to medial habenula function

    NASA Astrophysics Data System (ADS)

    Turner, Eric E.; Hsu, Yun-Wei; Wang, Si; Morton, Glenn; Zeng, Hongkui

    2013-03-01

    The habenula is a brain region found in all vertebrate species. It consists of medial and lateral subnuclei which make complex descending connections to the brainstem. Although the medial habenula (MHb) and its projection, the fasciculus retroflexus (FR), have been recognized for decades, their function remains obscure. The small size of the MHb in rodents, and the cellular and molecular complexity of this region, have made it difficult to study the function of this region with high specificity. Here we describe a Cre-mediated genetic system for expressing the microbial opsin channelrhodopsin (ChR2) specifically in the dorsal (dMHb) and ventral (vMHb) medial habenula. Genetically targeted expression of ChR2 allows MHb neurons to be selectively activated with light in acute brain slices with electrophysiological readouts, and in vivo by means of custom-built fiber optic cannulas. These tools will allow highly specific studies of MHb circuitry and the role of the MHb in behaviors related to addiction and mood regulation.

  17. Density-dependent covariant energy density functionals

    SciTech Connect

    Lalazissis, G. A.

    2012-10-20

    Relativistic nuclear energy density functionals are applied to the description of a variety of nuclear structure phenomena at and away fromstability line. Isoscalar monopole, isovector dipole and isoscalar quadrupole giant resonances are calculated using fully self-consistent relativistic quasiparticle randomphase approximation, based on the relativistic Hartree-Bogoliubovmodel. The impact of pairing correlations on the fission barriers in heavy and superheavy nuclei is examined. The role of pion in constructing desnity functionals is also investigated.

  18. Alloy Design, Combinatorial Synthesis, and Microstructure-Property Relations for Low-Density Fe-Mn-Al-C Austenitic Steels

    NASA Astrophysics Data System (ADS)

    Raabe, D.; Springer, H.; Gutierrez-Urrutia, I.; Roters, F.; Bausch, M.; Seol, J.-B.; Koyama, M.; Choi, P.-P.; Tsuzaki, K.

    2014-09-01

    We present recent developments in the field of austenitic steels with up to 18% reduced mass density. The alloys are based on the Fe-Mn-Al-C system. Here, two steel types are addressed. The first one is a class of low-density twinning-induced plasticity or single phase austenitic TWIP (SIMPLEX) steels with 25-30 wt.% Mn and <4-5 wt.% Al or even <8 wt.% Al when naturally aged. The second one is a class of κ-carbide strengthened austenitic steels with even higher Al content. Here, κ-carbides form either at 500-600°C or even during quenching for >10 wt.% Al. Three topics are addressed in more detail, namely, the combinatorial bulk high-throughput design of a wide range of corresponding alloy variants, the development of microstructure-property relations for such steels, and their susceptibility to hydrogen embrittlement.

  19. Induction of Human iPSC-Derived Cardiomyocyte Proliferation Revealed by Combinatorial Screening in High Density Microbioreactor Arrays

    PubMed Central

    Titmarsh, Drew M.; Glass, Nick R.; Mills, Richard J.; Hidalgo, Alejandro; Wolvetang, Ernst J.; Porrello, Enzo R.; Hudson, James E.; Cooper-White, Justin J.

    2016-01-01

    Inducing cardiomyocyte proliferation in post-mitotic adult heart tissue is attracting significant attention as a therapeutic strategy to regenerate the heart after injury. Model animal screens have identified several candidate signalling pathways, however, it remains unclear as to what extent these pathways can be exploited, either individually or in combination, in the human system. The advent of human cardiac cells from directed differentiation of human pluripotent stem cells (hPSCs) now provides the ability to interrogate human cardiac biology in vitro, but it remains difficult with existing culture formats to simply and rapidly elucidate signalling pathway penetrance and interplay. To facilitate high-throughput combinatorial screening of candidate biologicals or factors driving relevant molecular pathways, we developed a high-density microbioreactor array (HDMA) – a microfluidic cell culture array containing 8100 culture chambers. We used HDMAs to combinatorially screen Wnt, Hedgehog, IGF and FGF pathway agonists. The Wnt activator CHIR99021 was identified as the most potent molecular inducer of human cardiomyocyte proliferation, inducing cell cycle activity marked by Ki67, and an increase in cardiomyocyte numbers compared to controls. The combination of human cardiomyocytes with the HDMA provides a versatile and rapid tool for stratifying combinations of factors for heart regeneration. PMID:27097795

  20. Induction of Human iPSC-Derived Cardiomyocyte Proliferation Revealed by Combinatorial Screening in High Density Microbioreactor Arrays.

    PubMed

    Titmarsh, Drew M; Glass, Nick R; Mills, Richard J; Hidalgo, Alejandro; Wolvetang, Ernst J; Porrello, Enzo R; Hudson, James E; Cooper-White, Justin J

    2016-01-01

    Inducing cardiomyocyte proliferation in post-mitotic adult heart tissue is attracting significant attention as a therapeutic strategy to regenerate the heart after injury. Model animal screens have identified several candidate signalling pathways, however, it remains unclear as to what extent these pathways can be exploited, either individually or in combination, in the human system. The advent of human cardiac cells from directed differentiation of human pluripotent stem cells (hPSCs) now provides the ability to interrogate human cardiac biology in vitro, but it remains difficult with existing culture formats to simply and rapidly elucidate signalling pathway penetrance and interplay. To facilitate high-throughput combinatorial screening of candidate biologicals or factors driving relevant molecular pathways, we developed a high-density microbioreactor array (HDMA) - a microfluidic cell culture array containing 8100 culture chambers. We used HDMAs to combinatorially screen Wnt, Hedgehog, IGF and FGF pathway agonists. The Wnt activator CHIR99021 was identified as the most potent molecular inducer of human cardiomyocyte proliferation, inducing cell cycle activity marked by Ki67, and an increase in cardiomyocyte numbers compared to controls. The combination of human cardiomyocytes with the HDMA provides a versatile and rapid tool for stratifying combinations of factors for heart regeneration. PMID:27097795

  1. Phenomenological Relativistic Energy Density Functionals

    SciTech Connect

    Lalazissis, G. A.; Kartzikos, S.; Niksic, T.; Paar, N.; Vretenar, D.; Ring, P.

    2009-08-26

    The framework of relativistic nuclear energy density functionals is applied to the description of a variety of nuclear structure phenomena, not only in spherical and deformed nuclei along the valley of beta-stability, but also in exotic systems with extreme isospin values and close to the particle drip-lines. Dynamical aspects of exotic nuclear structure is explored using the fully consistent quasiparticle random-phase approximation based on the relativistic Hartree-Bogoliubov model. Recent applications of energy density functionals with explicit density dependence of the meson-nucleon couplings are presented.

  2. Density Functionals of Chemical Bonding

    PubMed Central

    Putz, Mihai V.

    2008-01-01

    The behavior of electrons in general many-electronic systems throughout the density functionals of energy is reviewed. The basic physico-chemical concepts of density functional theory are employed to highlight the energy role in chemical structure while its extended influence in electronic localization function helps in chemical bonding understanding. In this context the energy functionals accompanied by electronic localization functions may provide a comprehensive description of the global-local levels electronic structures in general and of chemical bonds in special. Becke-Edgecombe and author’s Markovian electronic localization functions are discussed at atomic, molecular and solid state levels. Then, the analytical survey of the main workable kinetic, exchange, and correlation density functionals within local and gradient density approximations is undertaken. The hierarchy of various energy functionals is formulated by employing both the parabolic and statistical correlation degree of them with the electronegativity and chemical hardness indices by means of quantitative structure-property relationship (QSPR) analysis for basic atomic and molecular systems. PMID:19325846

  3. Insights into Hox Protein Function from a Large Scale Combinatorial Analysis of Protein Domains

    PubMed Central

    Karlsson, Daniel; Dixit, Richa; Saadaoui, Mehdi; Monier, Bruno; Brun, Christine; Thor, Stefan; Vijayraghavan, K.; Perrin, Laurent; Pradel, Jacques; Graba, Yacine

    2011-01-01

    Protein function is encoded within protein sequence and protein domains. However, how protein domains cooperate within a protein to modulate overall activity and how this impacts functional diversification at the molecular and organism levels remains largely unaddressed. Focusing on three domains of the central class Drosophila Hox transcription factor AbdominalA (AbdA), we used combinatorial domain mutations and most known AbdA developmental functions as biological readouts to investigate how protein domains collectively shape protein activity. The results uncover redundancy, interactivity, and multifunctionality of protein domains as salient features underlying overall AbdA protein activity, providing means to apprehend functional diversity and accounting for the robustness of Hox-controlled developmental programs. Importantly, the results highlight context-dependency in protein domain usage and interaction, allowing major modifications in domains to be tolerated without general functional loss. The non-pleoitropic effect of domain mutation suggests that protein modification may contribute more broadly to molecular changes underlying morphological diversification during evolution, so far thought to rely largely on modification in gene cis-regulatory sequences. PMID:22046139

  4. Insights into Hox protein function from a large scale combinatorial analysis of protein domains.

    PubMed

    Merabet, Samir; Litim-Mecheri, Isma; Karlsson, Daniel; Dixit, Richa; Saadaoui, Mehdi; Monier, Bruno; Brun, Christine; Thor, Stefan; Vijayraghavan, K; Perrin, Laurent; Pradel, Jacques; Graba, Yacine

    2011-10-01

    Protein function is encoded within protein sequence and protein domains. However, how protein domains cooperate within a protein to modulate overall activity and how this impacts functional diversification at the molecular and organism levels remains largely unaddressed. Focusing on three domains of the central class Drosophila Hox transcription factor AbdominalA (AbdA), we used combinatorial domain mutations and most known AbdA developmental functions as biological readouts to investigate how protein domains collectively shape protein activity. The results uncover redundancy, interactivity, and multifunctionality of protein domains as salient features underlying overall AbdA protein activity, providing means to apprehend functional diversity and accounting for the robustness of Hox-controlled developmental programs. Importantly, the results highlight context-dependency in protein domain usage and interaction, allowing major modifications in domains to be tolerated without general functional loss. The non-pleoitropic effect of domain mutation suggests that protein modification may contribute more broadly to molecular changes underlying morphological diversification during evolution, so far thought to rely largely on modification in gene cis-regulatory sequences. PMID:22046139

  5. Uniformly constructing combinatorial solutions, combining a rational function with hyperbolic or trigonometric functions, for the (2+1) dimensional Broer-Kaup-Kupershmidt equation

    NASA Astrophysics Data System (ADS)

    Liu, Qing; Wang, Zi-Hua

    2010-12-01

    According to two dependent rational solutions to a generalized Riccati equation together with the equation itself, a rational-exponent solution to a nonlinear partial differential equation can be constructed. By selecting different parameter values in the rational-exponent solution, many families of combinatorial solutions combined with a rational function such as hyperbolic functions or trigonometric functions, are rapidly derived. This method is applied to the Whitham-Broer-Kaup equation and a series of combinatorial solutions are obtained, showing that this method is a more concise and efficient approach and can uniformly construct many types of combined solutions to nonlinear partial differential equations.

  6. A combinatorial approach to determine functional group effects on lipidoid-mediated siRNA delivery

    PubMed Central

    Mahon, Kerry P.; Love, Kevin T.; Whitehead, Kathryn A.; Qin, June; Akinc, Akin; Leshchiner, Elizaveta; Leshchiner, Ignaty; Langer, Robert

    2010-01-01

    The application of RNA interference (RNAi), either in the clinic or laboratory, requires safe and effective delivery methods. Here we develop a combinatorial approach to synthesize a library of delivery vectors based on two lipid-like substrates with known siRNA delivery capabilities. Members of this library have a mixture of lipid-like tails and feature appendages containing hydroxyl, carbamate, ether or amine functional groups as well as variations in alkyl chain length and branching. Using a luciferase reporter system in HeLa cells, we study the relationship between lipid chemical modification and delivery performance in vitro. The impact of the functional group was shown to vary depending on the overall amine content and tail number of the delivery vector. Additionally, in vivo performance was evaluated using a Factor VII knockdown assay. Two library members, each containing ether groups, were found to knock down the target protein at levels comparable to the parent delivery vector. These results demonstrate that small chemical changes to the delivery vector impact knockdown efficiency and cell viability both in vitro and in vivo. The work described here identifies new materials for siRNA delivery, as well as provides new insight into the parameters for optimized chemical makeup of lipid-like siRNA delivery materials. PMID:20715849

  7. Semiclassical origins of density functionals

    NASA Astrophysics Data System (ADS)

    Burke, Kieron

    By careful numerical analysis of non-relativistic atomic correlation energies, we show that (a) the local density approximation becomes relatively exact for the correlation energy as the atomic number approaches infinity, (b) we find the leading correction, which is about 38.5 milliHartrees per atom, (c) show how this correction dominates for larger atoms and (d) how to construct a generalized gradient approximation that respects this limit (See KB, A. Cancio, T. Gould, S. Pittalis, arXiv:1409.4834). The relevance to density functional calculations will also be explained. Support provided by NSF CHE-1464795.

  8. Density functionals from deep learning

    NASA Astrophysics Data System (ADS)

    McMahon, Jeffrey

    Density-functional theory is a formally exact description of a many-body quantum system in terms of its density; in practice, however, approximations to the universal density functional (DF) are necessary. Machine learning has recently been proposed as a novel approach to discover such a DF (or components of it). Conventional machine learning algorithms, however, are limited in their ability to process data in their raw form, leading to invariance and/or sensitivity issues. In this presentation, an alternative approach based on deep learning will be demonstrated. Deep learning allows computational models that are capable of discovering intricate structure in large and/or high-dimensional data sets with multiple levels of abstraction, and do not suffer from the aforementioned issues. Results from the application of this approach to the prediction of the kinetic-energy DF of noninteracting electrons will be presented. Using theoretical results from computer science, a connection between the underlying model and the theorems of Hohenberg and Kohn will also be suggested.

  9. Density functional theory: Foundations reviewed

    NASA Astrophysics Data System (ADS)

    Kryachko, Eugene S.; Ludeña, Eduardo V.

    2014-11-01

    Guided by the above motto (quotation), we review a broad range of issues lying at the foundations of Density Functional Theory, DFT, a theory which is currently omnipresent in our everyday computational study of atoms and molecules, solids and nano-materials, and which lies at the heart of modern many-body computational technologies. The key goal is to demonstrate that there are definitely the ways to improve DFT. We start by considering DFT in the larger context provided by reduced density matrix theory (RDMT) and natural orbital functional theory (NOFT), and examine the implications that N-representability conditions on the second-order reduced density matrix (2-RDM) have not only on RDMT and NOFT but, also, by extension, on the functionals of DFT. This examination is timely in view of the fact that necessary and sufficient N-representability conditions on the 2-RDM have recently been attained. In the second place, we review some problems appearing in the original formulation of the first Hohenberg-Kohn theorem which is still a subject of some controversy. In this vein we recall Lieb's comment on this proof and the extension to this proof given by Pino et al. (2009), and in this context examine the conditions that must be met in order that the one-to-one correspondence between ground-state densities and external potentials remains valid for finite subspaces (namely, the subspaces where all Kohn-Sham solutions are obtained in practical applications). We also consider the issue of whether the Kohn-Sham equations can be derived from basic principles or whether they are postulated. We examine this problem in relation to ab initio DFT. The possibility of postulating arbitrary Kohn-Sham-type equations, where the effective potential is by definition some arbitrary mixture of local and non-local terms, is discussed. We also deal with the issue of whether there exists a universal functional, or whether one should advocate instead the construction of problem

  10. ]: a density functional theory investigation

    NASA Astrophysics Data System (ADS)

    Ulian, Gianfranco; Tosoni, Sergio; Valdrè, Giovanni

    2014-09-01

    In this work, we modelled the structure, the compressional behaviour and the physical properties of talc over a wide range of pressure using a quantum mechanical approach based on periodic boundary conditions. We adopted the density functional theory using the B3LYP-D* functional, which includes a correction for the dispersive forces and all-electron Gaussian-type orbitals basis sets. An atomic level description of the athermal pressure-induced structural modification of talc is provided. From the compression results, we obtained the athermal ( T = 0 K) bulk modulus ( K T0), its first derivative ( K') and the athermal volume at zero pressure ( V 0) by a third-order Birch-Murnaghan equation with parameters K T0 = 56.25 GPa, K' = 5.66 and V 0 = 450.34 Å3. The mechanical behaviour is highly anisotropic, as observed by the axial compressibility. The presented data are in very good agreement with recent experimental results obtained by single-crystal neutron and X-ray diffraction experiments.

  11. Gedanken densities and exact constraints in density functional theory

    SciTech Connect

    Perdew, John P.; Ruzsinszky, Adrienn; Sun, Jianwei; Burke, Kieron

    2014-05-14

    Approximations to the exact density functional for the exchange-correlation energy of a many-electron ground state can be constructed by satisfying constraints that are universal, i.e., valid for all electron densities. Gedanken densities are designed for the purpose of this construction, but need not be realistic. The uniform electron gas is an old gedanken density. Here, we propose a spherical two-electron gedanken density in which the dimensionless density gradient can be an arbitrary positive constant wherever the density is non-zero. The Lieb-Oxford lower bound on the exchange energy can be satisfied within a generalized gradient approximation (GGA) by bounding its enhancement factor or simplest GGA exchange-energy density. This enhancement-factor bound is well known to be sufficient, but our gedanken density shows that it is also necessary. The conventional exact exchange-energy density satisfies no such local bound, but energy densities are not unique, and the simplest GGA exchange-energy density is not an approximation to it. We further derive a strongly and optimally tightened bound on the exchange enhancement factor of a two-electron density, which is satisfied by the local density approximation but is violated by all published GGA's or meta-GGA’s. Finally, some consequences of the non-uniform density-scaling behavior for the asymptotics of the exchange enhancement factor of a GGA or meta-GGA are given.

  12. Nuclear level density and γ-ray strength function of 43Sc

    NASA Astrophysics Data System (ADS)

    Bürger, A.; Larsen, A. C.; Hilaire, S.; Guttormsen, M.; Harissopulos, S.; Kmiecik, M.; Konstantinopoulos, T.; Krtička, M.; Lagoyannis, A.; Lönnroth, T.; Mazurek, K.; Norrby, M.; Nyhus, H. T.; Perdikakis, G.; Siem, S.; Spyrou, A.; Syed, N. U. H.

    2012-06-01

    The nuclear level density and the γ-ray strength function have been determined for 43Sc in the energy range up to 2 MeV below the neutron separation energy using the Oslo method with the 46Ti(p,α)43Sc reaction. A comparison to 45Sc shows that the level density of 43Sc is smaller by an approximately constant factor of two. This behavior is well reproduced in a microscopic, combinatorial model calculation. The γ-ray strength function increases at low γ-ray energies, a feature which has been observed in several nuclei but which still awaits theoretical explanation.

  13. Combinatorial theory of the semiclassical evaluation of transport moments II: Algorithmic approach for moment generating functions

    SciTech Connect

    Berkolaiko, G.; Kuipers, J.

    2013-12-15

    Electronic transport through chaotic quantum dots exhibits universal behaviour which can be understood through the semiclassical approximation. Within the approximation, calculation of transport moments reduces to codifying classical correlations between scattering trajectories. These can be represented as ribbon graphs and we develop an algorithmic combinatorial method to generate all such graphs with a given genus. This provides an expansion of the linear transport moments for systems both with and without time reversal symmetry. The computational implementation is then able to progress several orders further than previous semiclassical formulae as well as those derived from an asymptotic expansion of random matrix results. The patterns observed also suggest a general form for the higher orders.

  14. Unexpected properties of a density functional

    SciTech Connect

    Karwowski, J.; Stanke, M.

    2005-02-01

    An observation on a pathological behavior of an exact density functional derived from either relativistic (Dirac) or nonrelativistic (Levy-Leblond) quantum-mechanical equation is reported. As expected, in the case of a one-electron atom the variational minimum of this functional is equal to the exact ground-state energy. However, apart from the correct density, this minimum is reached also by an infinite set of densities which do not correspond to the exact wave function. This paradoxical property of the functional is related to the multicomponent structure of both Dirac and Levy-Leblond wave functions. In particular, imposing the correct boundary conditions upon the trial densities removes only a part of the fake solutions. The results of this study demonstrate that in density-functional theories derived from models based on multicomponent wave functions, one should not expect any simple relation between the accuracy of the energy and the correctness of the corresponding density.

  15. A Safari Through Density Functional Theory

    NASA Astrophysics Data System (ADS)

    Dreizler, Reiner M.; Lüdde, Cora S.

    Density functional theory is widely used to treat quantum many body problems in many areas of physics and related fields. A brief survey of this method covering foundations, functionals and applications is presented here.

  16. Magnetic fields and density functional theory

    SciTech Connect

    Salsbury Jr., Freddie

    1999-02-01

    A major focus of this dissertation is the development of functionals for the magnetic susceptibility and the chemical shielding within the context of magnetic field density functional theory (BDFT). These functionals depend on the electron density in the absence of the field, which is unlike any other treatment of these responses. There have been several advances made within this theory. The first of which is the development of local density functionals for chemical shieldings and magnetic susceptibilities. There are the first such functionals ever proposed. These parameters have been studied by constructing functionals for the current density and then using the Biot-Savart equations to obtain the responses. In order to examine the advantages and disadvantages of the local functionals, they were tested numerically on some small molecules.

  17. Density functional theory for carbon dioxide crystal

    SciTech Connect

    Chang, Yiwen; Mi, Jianguo Zhong, Chongli

    2014-05-28

    We present a density functional approach to describe the solid−liquid phase transition, interfacial and crystal structure, and properties of polyatomic CO{sub 2}. Unlike previous phase field crystal model or density functional theory, which are derived from the second order direct correlation function, the present density functional approach is based on the fundamental measure theory for hard-sphere repulsion in solid. More importantly, the contributions of enthalpic interactions due to the dispersive attractions and of entropic interactions arising from the molecular architecture are integrated in the density functional model. Using the theoretical model, the predicted liquid and solid densities of CO{sub 2} at equilibrium triple point are in good agreement with the experimental values. Based on the structure of crystal-liquid interfaces in different planes, the corresponding interfacial tensions are predicted. Their respective accuracies need to be tested.

  18. Accurate density functional thermochemistry for larger molecules.

    SciTech Connect

    Raghavachari, K.; Stefanov, B. B.; Curtiss, L. A.; Lucent Tech.

    1997-06-20

    Density functional methods are combined with isodesmic bond separation reaction energies to yield accurate thermochemistry for larger molecules. Seven different density functionals are assessed for the evaluation of heats of formation, Delta H 0 (298 K), for a test set of 40 molecules composed of H, C, O and N. The use of bond separation energies results in a dramatic improvement in the accuracy of all the density functionals. The B3-LYP functional has the smallest mean absolute deviation from experiment (1.5 kcal mol/f).

  19. Locality of correlation in density functional theory

    NASA Astrophysics Data System (ADS)

    Burke, Kieron; Cancio, Antonio; Gould, Tim; Pittalis, Stefano

    2016-08-01

    The Hohenberg-Kohn density functional was long ago shown to reduce to the Thomas-Fermi (TF) approximation in the non-relativistic semiclassical (or large-Z) limit for all matter, i.e., the kinetic energy becomes local. Exchange also becomes local in this limit. Numerical data on the correlation energy of atoms support the conjecture that this is also true for correlation, but much less relevant to atoms. We illustrate how expansions around a large particle number are equivalent to local density approximations and their strong relevance to density functional approximations. Analyzing highly accurate atomic correlation energies, we show that EC → -AC ZlnZ + BCZ as Z → ∞, where Z is the atomic number, AC is known, and we estimate BC to be about 37 mhartree. The local density approximation yields AC exactly, but a very incorrect value for BC, showing that the local approximation is less relevant for the correlation alone. This limit is a benchmark for the non-empirical construction of density functional approximations. We conjecture that, beyond atoms, the leading correction to the local density approximation in the large-Z limit generally takes this form, but with BC a functional of the TF density for the system. The implications for the construction of approximate density functionals are discussed.

  20. Locality of correlation in density functional theory.

    PubMed

    Burke, Kieron; Cancio, Antonio; Gould, Tim; Pittalis, Stefano

    2016-08-01

    The Hohenberg-Kohn density functional was long ago shown to reduce to the Thomas-Fermi (TF) approximation in the non-relativistic semiclassical (or large-Z) limit for all matter, i.e., the kinetic energy becomes local. Exchange also becomes local in this limit. Numerical data on the correlation energy of atoms support the conjecture that this is also true for correlation, but much less relevant to atoms. We illustrate how expansions around a large particle number are equivalent to local density approximations and their strong relevance to density functional approximations. Analyzing highly accurate atomic correlation energies, we show that EC → -AC ZlnZ + BCZ as Z → ∞, where Z is the atomic number, AC is known, and we estimate BC to be about 37 mhartree. The local density approximation yields AC exactly, but a very incorrect value for BC, showing that the local approximation is less relevant for the correlation alone. This limit is a benchmark for the non-empirical construction of density functional approximations. We conjecture that, beyond atoms, the leading correction to the local density approximation in the large-Z limit generally takes this form, but with BC a functional of the TF density for the system. The implications for the construction of approximate density functionals are discussed. PMID:27497544

  1. Particle conservation in dynamical density functional theory.

    PubMed

    de Las Heras, Daniel; Brader, Joseph M; Fortini, Andrea; Schmidt, Matthias

    2016-06-22

    We present the exact adiabatic theory for the dynamics of the inhomogeneous density distribution of a classical fluid. Erroneous particle number fluctuations of dynamical density functional theory are absent, both for canonical and grand canonical initial conditions. We obtain the canonical free energy functional, which yields the adiabatic interparticle forces of overdamped Brownian motion. Using an exact and one of the most advanced approximate hard core free energy functionals, we obtain excellent agreement with simulations. The theory applies to finite systems in and out of equilibrium. PMID:27115673

  2. Particle conservation in dynamical density functional theory

    NASA Astrophysics Data System (ADS)

    de las Heras, Daniel; Brader, Joseph M.; Fortini, Andrea; Schmidt, Matthias

    2016-06-01

    We present the exact adiabatic theory for the dynamics of the inhomogeneous density distribution of a classical fluid. Erroneous particle number fluctuations of dynamical density functional theory are absent, both for canonical and grand canonical initial conditions. We obtain the canonical free energy functional, which yields the adiabatic interparticle forces of overdamped Brownian motion. Using an exact and one of the most advanced approximate hard core free energy functionals, we obtain excellent agreement with simulations. The theory applies to finite systems in and out of equilibrium.

  3. Combinatorial Motor Training Results in Functional Reorganization of Remaining Motor Cortex after Controlled Cortical Impact in Rats.

    PubMed

    Combs, Hannah L; Jones, Theresa A; Kozlowski, Dorothy A; Adkins, DeAnna L

    2016-04-15

    Cortical reorganization subsequent to post-stroke motor rehabilitative training (RT) has been extensively examined in animal models and humans. However, similar studies focused on the effects of motor training after traumatic brain injury (TBI) are lacking. We previously reported that after a moderate/severe TBI in adult male rats, functional improvements in forelimb use were accomplished only with a combination of skilled forelimb reach training and aerobic exercise, with or without nonimpaired forelimb constraint. Thus, the current study was designed to examine the relationship between functional motor cortical map reorganization after experimental TBI and the behavioral improvements resulting from this combinatorial rehabilitative regime. Adult male rats were trained to proficiency on a skilled reaching task, received a unilateral controlled cortical impact (CCI) over the forelimb area of the caudal motor cortex (CMC). Three days post-CCI, animals began RT (n = 13) or no rehabilitative training (NoRT) control procedures (n = 13). The RT group participated in daily skilled reach training, voluntary aerobic exercise, and nonimpaired forelimb constraint. This RT regimen significantly improved impaired forelimb reaching success and normalized reaching strategies, consistent with previous findings. RT also enlarged the area of motor cortical wrist representation, derived by intracortical microstimulation, compared to NoRT. These findings indicate that sufficient RT can greatly improve motor function and improve the functional integrity of remaining motor cortex after a moderate/severe CCI. When compared with findings from stroke models, these findings also suggest that more intense RT may be needed to improve motor function and remodel the injured cortex after TBI. PMID:26421759

  4. Nonlocal kinetic-energy-density functionals

    SciTech Connect

    Garcia-Gonzalez, P.; Alvarellos, J.E.; Chacon, E. |

    1996-04-01

    In this paper we present nonlocal kinetic-energy functionals {ital T}[{ital n}] within the average density approximation (ADA) framework, which do not require any extra input when applied to any electron system and recover the exact kinetic energy and the linear response function of a homogeneous system. In contrast with previous ADA functionals, these present good behavior of the long-range tail of the exact weight function. The averaging procedure for the kinetic functional (averaging the Fermi momentum of the electron gas, instead of averaging the electron density) leads to a functional without numerical difficulties in the calculation of extended systems, and it gives excellent results when applied to atoms and jellium surfaces. {copyright} {ital 1996 The American Physical Society.}

  5. Density Functionals with Broad Applicability in Chemistry

    SciTech Connect

    Zhao, Yan; Truhlar, Donald G.

    2008-02-01

    The research described in this product was performed in part in the Environmental Molecular Sciences Laboratory, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory. Although density functional theory is widely used in the computational chemistry community, the most popular density functional, B3LYP, has some serious shortcomings: (i) it is better for main-group chemistry than for transition metals; (ii) it systematically underestimates reaction barrier heights; (iii) it is inaccurate for interactions dominated by mediumrange correlation energy, such as van der Waals attraction, aromatic-aromatic stacking, and alkane isomerization energies. We have developed a variety of databases for testing and designing new density functionals. We used these data to design new density functionals, called M06-class (and, earlier, M05-class) functionals, for which we enforced some fundamental exact constraints such as the uniform-electron-gas limit and the absence of self-correlation energy. Our M06-class functionals depend on spin-up and spin-down electron densities (i.e., spin densities), spin density gradients, spin kinetic energy densities, and, for nonlocal (also called hybrid) functionals, Hartree-Fock exchange. We have developed four new functionals that overcome the above-mentioned difficulties: (a) M06, a hybrid meta functional, is a functional with good accuracy “across-theboard” for transition metals, main group thermochemistry, medium-range correlation energy, and barrier heights; (b) M06- 2X, another hybrid meta functional, is not good for transition metals but has excellent performance for main group chemistry, predicts accurate valence and Rydberg electronic excitation energies, and is an excellent functional for aromatic-aromatic stacking interactions; (c) M06-L is not as accurate as M06 for barrier heights but is the most accurate

  6. Periodic subsystem density-functional theory

    SciTech Connect

    Genova, Alessandro; Pavanello, Michele; Ceresoli, Davide

    2014-11-07

    By partitioning the electron density into subsystem contributions, the Frozen Density Embedding (FDE) formulation of subsystem Density Functional Theory (DFT) has recently emerged as a powerful tool for reducing the computational scaling of Kohn–Sham DFT. To date, however, FDE has been employed to molecular systems only. Periodic systems, such as metals, semiconductors, and other crystalline solids have been outside the applicability of FDE, mostly because of the lack of a periodic FDE implementation. To fill this gap, in this work we aim at extending FDE to treat subsystems of molecular and periodic character. This goal is achieved by a dual approach. On one side, the development of a theoretical framework for periodic subsystem DFT. On the other, the realization of the method into a parallel computer code. We find that periodic FDE is capable of reproducing total electron densities and (to a lesser extent) also interaction energies of molecular systems weakly interacting with metallic surfaces. In the pilot calculations considered, we find that FDE fails in those cases where there is appreciable density overlap between the subsystems. Conversely, we find FDE to be in semiquantitative agreement with Kohn–Sham DFT when the inter-subsystem density overlap is low. We also conclude that to make FDE a suitable method for describing molecular adsorption at surfaces, kinetic energy density functionals that go beyond the GGA level must be employed.

  7. Connection formulas for thermal density functional theory

    NASA Astrophysics Data System (ADS)

    Pribram-Jones, A.; Burke, K.

    2016-05-01

    The adiabatic connection formula of ground-state density functional theory relates the correlation energy to a coupling-constant integral over a purely potential contribution, and is widely used to understand and improve approximations. The corresponding formula for thermal density functional theory is cast as an integral over temperatures instead, ranging upward from the system's physical temperature. We also show how to relate different correlation components to each other, either in terms of temperature or coupling-constant integrations. We illustrate our results on the uniform electron gas.

  8. Combinatorial protein therapy of angiogenic and arteriogenic factors remarkably improves collaterogenesis and cardiac function in pigs

    PubMed Central

    Lu, Huixia; Xu, Xinsheng; Zhang, Mei; Cao, Renhai; Bråkenhielm, Ebba; Li, Changjiang; Lin, Huili; Yao, Guihua; Sun, Huiwen; Qi, Lihang; Tang, Mengxiong; Dai, Hongyan; Zhang, Yanen; Su, Runyi; Bi, Yanwen; Zhang, Yun; Cao, Yihai

    2007-01-01

    Establishment of functional and stable collaterals in the ischemic myocardium is crucial to restoring cardiac function after myocardial infarction. Here, we show that only dual delivery of a combination of angiogenic and arteriogenic factors to the ischemic myocardium could significantly reestablish stable collateral networks and improve myocardial perfusion and function. A combination of FGF-2 with PDGF-BB, two factors primarily targeting endothelial cells and vascular smooth muscle cells, remarkably promotes myocardial collateral growth and stabilizes the newly formed collateral networks, which significantly restore myocardial perfusion and function. Using various members of the PDGF family together with FGF-2 in an angiogenesis assay, we demonstrate that PDGFR-α is mainly involved in angiogenic synergism, whereas PDGFR-β mediates vessel stability signals. Our findings provide conceptual guidelines for the clinical development of proangiogenic/arteriogenic factors for the treatment of ischemic heart disease. PMID:17636133

  9. Koopmans' condition for density-functional theory

    SciTech Connect

    Dabo, Ismaila; Ferretti, Andrea; Poilvert, Nicolas; Marzari, Nicola; Li, Yanli; Cococcioni, Matteo

    2010-09-15

    In approximate Kohn-Sham density-functional theory, self-interaction manifests itself as the dependence of the energy of an orbital on its fractional occupation. This unphysical behavior translates into qualitative and quantitative errors that pervade many fundamental aspects of density-functional predictions. Here, we first examine self-interaction in terms of the discrepancy between total and partial electron removal energies, and then highlight the importance of imposing the generalized Koopmans' condition - that identifies orbital energies as opposite total electron removal energies - to resolve this discrepancy. In the process, we derive a correction to approximate functionals that, in the frozen-orbital approximation, eliminates the unphysical occupation dependence of orbital energies up to the third order in the single-particle densities. This non-Koopmans correction brings physical meaning to single-particle energies; when applied to common local or semilocal density functionals it provides results that are in excellent agreement with experimental data - with an accuracy comparable to that of GW many-body perturbation theory - while providing an explicit total energy functional that preserves or improves on the description of established structural properties.

  10. Density-functional expansion methods: Grand challenges

    PubMed Central

    Giese, Timothy J.; York, Darrin M.

    2016-01-01

    We discuss the source of errors in semiempirical density functional expansion (VE) methods. In particular, we show that VE methods are capable of well-reproducing their standard Kohn-Sham density functional method counterparts, but suffer from large errors upon using one or more of these approximations: the limited size of the atomic orbital basis, the Slater monopole auxiliary basis description of the response density, and the one- and two-body treatment of the core-Hamiltonian matrix elements. In the process of discussing these approximations and highlighting their symptoms, we introduce a new model that supplements the second-order density-functional tight-binding model with a self-consistent charge-dependent chemical potential equalization correction; we review our recently reported method for generalizing the auxiliary basis description of the atomic orbital response density; and we decompose the first-order potential into a summation of additive atomic components and many-body corrections, and from this examination, we provide new insights and preliminary results that motivate and inspire new approximate treatments of the core-Hamiltonian. PMID:27293378

  11. Generalized local-spin-density-functional theory

    NASA Astrophysics Data System (ADS)

    Guo, Yufei; Whitehead, M. A.

    1991-01-01

    An alternative density-functional theory, the generalized local-spin-density-functional (G-LSD) theory, is proposed based on the boundary conditions and sum rule of the Fermi-correlation factor in the Hartree-Fock (HF) limit. It avoids the physical restrictions to the boundary conditions and the sum rule used in the generalized exchange local-spin-density-functional (GX-LSD) theory completely, the homogeneous electron-density approximation in the Hartree-Fock-Slater (HFS) theory and in the Gáspár-Kohn-Sham (GKS) theory partly, and the time-consuming step to search for the optimal exchange parameter for each atom or ion in the Xα and Ξa theories. The alternative G-LSD theory generates the GX-LSD, HFS, GKS, Xα, and Ξa theories, when additional approximations or certain Fermi-hole shapes or high electron-density limit in a system are used. Theoretically, the G-LSD theory is more rigorous than the GX-LSD, HFS, GKS, and Ξa theories. Numerically, the statistical total energies for atoms in the G-LSD theory are in excellent agreement with the HF results, when the Gopinathan, Whitehead, and Bogdanovic [Phys. Rev. A 14, 1 (1976)] Fermi-hole parameters are used.

  12. Communication: Embedded fragment stochastic density functional theory

    SciTech Connect

    Neuhauser, Daniel; Baer, Roi; Rabani, Eran

    2014-07-28

    We develop a method in which the electronic densities of small fragments determined by Kohn-Sham density functional theory (DFT) are embedded using stochastic DFT to form the exact density of the full system. The new method preserves the scaling and the simplicity of the stochastic DFT but cures the slow convergence that occurs when weakly coupled subsystems are treated. It overcomes the spurious charge fluctuations that impair the applications of the original stochastic DFT approach. We demonstrate the new approach on a fullerene dimer and on clusters of water molecules and show that the density of states and the total energy can be accurately described with a relatively small number of stochastic orbitals.

  13. Density functional theory studies of etoricoxib

    NASA Astrophysics Data System (ADS)

    Sachdeva, Ritika; Kaur, Prabhjot; Singh, V. P.; Saini, G. S. S.

    2016-05-01

    Etoricoxib is a COX-2 selective inhibitor drug with molecular formula C18H15ClN2O2S. It is primarily used for the treatment of arthritis(rheumatoid, psoriatic, osteoarthritis), ankylosing spondylitis, gout and chronic low back pain. Theoretical studies of the molecule including geometry optimization and vibrational frequency calculations were carried out with the help of density functional theory calculations using 6-311++ g (d, p) basis set and B3LYP functional.

  14. Nuclear moments in covariant density functional theory

    NASA Astrophysics Data System (ADS)

    Meng, J.; Zhao, P. W.; Zhang, S. Q.; Hu, J. N.; Li, J.

    2014-05-01

    Recent progresses on microscopic and self-consistent description of the nuclear moments in covariant density functional theory based on a point-coupling interaction are briefly reviewed. In particular, the electric quadrupole moments of Cd isotopes and the magnetic moments of Pb isotopes are discussed.

  15. On an Integer Sequence Related to a Product of Trigonometric Functions, and Its Combinatorial Relevance

    NASA Astrophysics Data System (ADS)

    Andrica, Dorin; Tomescu, Ioan

    2002-11-01

    In this paper it is shown that for n == 0 or 3 (mod 4), the middle term S(n) in the expansion of the polynomial (1+x)(1+x^2)... (1+x^n) occurs naturally when one analyzes when a discontinuous product of trigonometric functions is a derivative of a function. This number also represents the number of partitions of T_n/2 = n(n+1)/4, (where T_n is the nth triangular number) into distinct parts less than or equal to n. It is proved in a constructive way that S(n)>= 6S(n-4) for every n >= 8, and an asymptotic evaluation of S(n)^{1/n} is obtained as a consequence of the unimodality of the coefficients of this polynomial. Also an integral expression of S(n) is deduced.

  16. Probability density function learning by unsupervised neurons.

    PubMed

    Fiori, S

    2001-10-01

    In a recent work, we introduced the concept of pseudo-polynomial adaptive activation function neuron (FAN) and presented an unsupervised information-theoretic learning theory for such structure. The learning model is based on entropy optimization and provides a way of learning probability distributions from incomplete data. The aim of the present paper is to illustrate some theoretical features of the FAN neuron, to extend its learning theory to asymmetrical density function approximation, and to provide an analytical and numerical comparison with other known density function estimation methods, with special emphasis to the universal approximation ability. The paper also provides a survey of PDF learning from incomplete data, as well as results of several experiments performed on real-world problems and signals. PMID:11709808

  17. Differentiability of Lieb functional in electronic density functional theory

    NASA Astrophysics Data System (ADS)

    Lammert, Paul E.

    A solid understanding of the Lieb functional FL is important because of its centrality in the foundations of electronic density functional theory. A basic question is whether directional derivatives of FL at an ensemble-V-representable density are given by (minus) the potential. A widely accepted purported proof that FL is Gâteaux differentiable at EV-representable densities would say, ?yes.? But that proof is fallacious, as shown here. FL is not Gâteaux differentiable in the normal sense, nor is it continuous. By means of a constructive approach, however, we are able to show that the derivative of FL at an EV-representable density ?0 in the direction of ?1 is given by the potential if ?0 and ?1 are everywhere strictly greater than zero, and they and the ground state wave function have square integrable derivatives through second order.

  18. Density gradient expansion of correlation functions

    NASA Astrophysics Data System (ADS)

    van Leeuwen, Robert

    2013-04-01

    We present a general scheme based on nonlinear response theory to calculate the expansion of correlation functions such as the pair-correlation function or the exchange-correlation hole of an inhomogeneous many-particle system in terms of density derivatives of arbitrary order. We further derive a consistency condition that is necessary for the existence of the gradient expansion. This condition is used to carry out an infinite summation of terms involving response functions up to infinite order from which it follows that the coefficient functions of the gradient expansion can be expressed in terms of the local density profile rather than the background density around which the expansion is carried out. We apply the method to the calculation of the gradient expansion of the one-particle density matrix to second order in the density gradients and recover in an alternative manner the result of Gross and Dreizler [Gross and Dreizler, Z. Phys. AZPAADB0340-219310.1007/BF01413038 302, 103 (1981)], which was derived using the Kirzhnits method. The nonlinear response method is more general and avoids the turning point problem of the Kirzhnits expansion. We further give a description of the exchange hole in momentum space and confirm the wave vector analysis of Langreth and Perdew [Langreth and Perdew, Phys. Rev. BPRBMDO1098-012110.1103/PhysRevB.21.5469 21, 5469 (1980)] for this case. This is used to derive that the second-order gradient expansion of the system averaged exchange hole satisfies the hole sum rule and to calculate the gradient coefficient of the exchange energy without the need to regularize divergent integrals.

  19. Density-Functional Theory of Thermal Transport

    NASA Astrophysics Data System (ADS)

    Eich, F. G.; Principi, A.; di Ventra, M.; Vignale, G.

    2014-03-01

    We have recently introduced a non-equilibrium density-functional theory of local temperature and associated energy density that is suitable for the study of thermoelectric phenomena from first principles. This theory rests on a local temperature field coupled to the energy-density operator. Here we apply the theory to a simple two-terminal setup, in which the terminals are held at different temperatures. We show that our treatment becomes equivalent to the standard Landauer-Büttiker formulation of thermal transport in the non-interacting limit. We gratefully acknowledge support from DOE under Grant No. DE-FG02-05ER46203 (FGE, AP, GV) and DE-FG02-05ER46204 (MD).

  20. Thermal density functional theory, ensemble density functional theory, and potential functional theory for warm dense matter

    NASA Astrophysics Data System (ADS)

    Pribram-Jones, Aurora

    Warm dense matter (WDM) is a high energy phase between solids and plasmas, with characteristics of both. It is present in the centers of giant planets, within the earth's core, and on the path to ignition of inertial confinement fusion. The high temperatures and pressures of warm dense matter lead to complications in its simulation, as both classical and quantum effects must be included. One of the most successful simulation methods is density functional theory-molecular dynamics (DFT-MD). Despite great success in a diverse array of applications, DFT-MD remains computationally expensive and it neglects the explicit temperature dependence of electron-electron interactions known to exist within exact DFT. Finite-temperature density functional theory (FT DFT) is an extension of the wildly successful ground-state DFT formalism via thermal ensembles, broadening its quantum mechanical treatment of electrons to include systems at non-zero temperatures. Exact mathematical conditions have been used to predict the behavior of approximations in limiting conditions and to connect FT DFT to the ground-state theory. An introduction to FT DFT is given within the context of ensemble DFT and the larger field of DFT is discussed for context. Ensemble DFT is used to describe ensembles of ground-state and excited systems. Exact conditions in ensemble DFT and the performance of approximations depend on ensemble weights. Using an inversion method, exact Kohn-Sham ensemble potentials are found and compared to approximations. The symmetry eigenstate Hartree-exchange approximation is in good agreement with exact calculations because of its inclusion of an ensemble derivative discontinuity. Since ensemble weights in FT DFT are temperature-dependent Fermi weights, this insight may help develop approximations well-suited to both ground-state and FT DFT. A novel, highly efficient approach to free energy calculations, finite-temperature potential functional theory, is derived, which has the

  1. Coverage of dynamic correlation effects by density functional theory functionals: density-based analysis for neon.

    PubMed

    Jankowski, K; Nowakowski, K; Grabowski, I; Wasilewski, J

    2009-04-28

    The problem of linking the dynamic electron correlation effects defined in traditional ab initio methods [or wave function theories (WFTs)] with the structure of the individual density functional theory (DFT) exchange and correlation functionals has been analyzed for the Ne atom, for which nondynamic correlation effects play a negligible role. A density-based approach directly hinged on difference radial-density (DRD) distributions defined with respect the Hartree-Fock radial density has been employed for analyzing the impact of dynamic correlation effects on the density. Attention has been paid to the elimination of basis-set incompleteness errors. The DRD distributions calculated by several ab initio methods have been compared to their DFT counterparts generated for representatives of several generations of broadly used exchange-correlation functionals and for the recently developed orbital-dependent OEP2 exchange-correlation functional [Bartlett et al., J. Chem. Phys. 122, 034104 (2005)]. For the local, generalized-gradient, and hybrid functionals it has been found that the dynamic correlation effects are to a large extend accounted for by densities resulting from exchange-only calculations. Additional calculations with self-interaction corrected exchange potentials indicate that this finding cannot be explained as an artifact caused by the self-interaction error. It has been demonstrated that the VWN5 and LYP correlation functionals do not represent any substantial dynamical correlation effects on the electron density, whereas these effects are well represented by the orbital-dependent OEP2 correlation functional. Critical comparison of the present results with their counterparts reported in literature has been made. Some attention has been paid to demonstrating the differences between the energy- and density-based perspectives. They indicate the usefulness of density-based criteria for developing new exchange-correlation functionals. PMID:19405556

  2. Coverage of dynamic correlation effects by density functional theory functionals: Density-based analysis for neon

    NASA Astrophysics Data System (ADS)

    Jankowski, K.; Nowakowski, K.; Grabowski, I.; Wasilewski, J.

    2009-04-01

    The problem of linking the dynamic electron correlation effects defined in traditional ab initio methods [or wave function theories (WFTs)] with the structure of the individual density functional theory (DFT) exchange and correlation functionals has been analyzed for the Ne atom, for which nondynamic correlation effects play a negligible role. A density-based approach directly hinged on difference radial-density (DRD) distributions defined with respect the Hartree-Fock radial density has been employed for analyzing the impact of dynamic correlation effects on the density. Attention has been paid to the elimination of basis-set incompleteness errors. The DRD distributions calculated by several ab initio methods have been compared to their DFT counterparts generated for representatives of several generations of broadly used exchange-correlation functionals and for the recently developed orbital-dependent OEP2 exchange-correlation functional [Bartlett et al., J. Chem. Phys. 122, 034104 (2005)]. For the local, generalized-gradient, and hybrid functionals it has been found that the dynamic correlation effects are to a large extend accounted for by densities resulting from exchange-only calculations. Additional calculations with self-interaction corrected exchange potentials indicate that this finding cannot be explained as an artifact caused by the self-interaction error. It has been demonstrated that the VWN5 and LYP correlation functionals do not represent any substantial dynamical correlation effects on the electron density, whereas these effects are well represented by the orbital-dependent OEP2 correlation functional. Critical comparison of the present results with their counterparts reported in literature has been made. Some attention has been paid to demonstrating the differences between the energy- and density-based perspectives. They indicate the usefulness of density-based criteria for developing new exchange-correlation functionals.

  3. Orbital-optimized density cumulant functional theory

    SciTech Connect

    Sokolov, Alexander Yu. Schaefer, Henry F.

    2013-11-28

    In density cumulant functional theory (DCFT) the electronic energy is evaluated from the one-particle density matrix and two-particle density cumulant, circumventing the computation of the wavefunction. To achieve this, the one-particle density matrix is decomposed exactly into the mean-field (idempotent) and correlation components. While the latter can be entirely derived from the density cumulant, the former must be obtained by choosing a specific set of orbitals. In the original DCFT formulation [W. Kutzelnigg, J. Chem. Phys. 125, 171101 (2006)] the orbitals were determined by diagonalizing the effective Fock operator, which introduces partial orbital relaxation. Here we present a new orbital-optimized formulation of DCFT where the energy is variationally minimized with respect to orbital rotations. This introduces important energy contributions and significantly improves the description of the dynamic correlation. In addition, it greatly simplifies the computation of analytic gradients, for which expressions are also presented. We offer a perturbative analysis of the new orbital stationarity conditions and benchmark their performance for a variety of chemical systems.

  4. Density functional calculations on model tyrosyl radicals.

    PubMed Central

    Himo, F; Gräslund, A; Eriksson, L A

    1997-01-01

    A gradient-corrected density functional theory approach (PWP86) has been applied, together with large basis sets (IGLO-III), to investigate the structure and hyperfine properties of model tyrosyl free radicals. In nature, these radicals are observed in, e.g., the charge transfer pathways in photosystem II (PSII) and in ribonucleotide reductases (RNRs). By comparing spin density distributions and proton hyperfine couplings with experimental data, it is confirmed that the tyrosyl radicals present in the proteins are neutral. It is shown that hydrogen bonding to the phenoxyl oxygen atom, when present, causes a reduction in spin density on O and a corresponding increase on C4. Calculated proton hyperfine coupling constants for the beta-protons show that the alpha-carbon is rotated 75-80 degrees out of the plane of the ring in PSII and Salmonella typhimurium RNR, but only 20-30 degrees in, e.g., Escherichia coli, mouse, herpes simplex, and bacteriophage T4-induced RNRs. Furthermore, based on the present calculations, we have revised the empirical parameters used in the experimental determination of the oxygen spin density in the tyrosyl radical in E. coli RNR and of the ring carbon spin densities, from measured hyperfine coupling constants. Images FIGURE 1 FIGURE 5 PMID:9083661

  5. Teaching Density Functional Theory Through Experiential Learning

    NASA Astrophysics Data System (ADS)

    Narasimhan, Shobhana

    2015-09-01

    Today, quantum mechanical density functional theory is often the method of choice for performing accurate calculations on atomic, molecular and condensed matter systems. Here, I share some of my experiences in teaching the necessary basics of solid state physics, as well as the theory and practice of density functional theory, in a number of workshops held in developing countries over the past two decades. I discuss the advantages of supplementing the usual mathematically formal teaching methods, characteristic of graduate courses, with the use of visual imagery and analogies. I also describe a successful experiment we carried out, which resulted in a joint publication co-authored by 67 lecturers and students participating in a summer school.

  6. Density Functional Theory Models for Radiation Damage

    NASA Astrophysics Data System (ADS)

    Dudarev, S. L.

    2013-07-01

    Density functional theory models developed over the past decade provide unique information about the structure of nanoscale defects produced by irradiation and about the nature of short-range interaction between radiation defects, clustering of defects, and their migration pathways. These ab initio models, involving no experimental input parameters, appear to be as quantitatively accurate and informative as the most advanced experimental techniques developed for the observation of radiation damage phenomena. Density functional theory models have effectively created a new paradigm for the scientific investigation and assessment of radiation damage effects, offering new insight into the origin of temperature- and dose-dependent response of materials to irradiation, a problem of pivotal significance for applications.

  7. Combinatorial Mechanical Metamaterials

    NASA Astrophysics Data System (ADS)

    van Hecke, Martin

    The structure of most mechanical metamaterials is periodic so that their design space is that of the unit cell. Here we introduce a combinatorial strategy to create a vast number of distinct mechanical metamaterials, each with a unique spatial texture and response. These are aperiodic stackings of anisotropic building blocks, and their functionality rests on both the block design and their stacking configuration which is governed by a tiling problem. We realize such metamaterials by 3D printing, and show that they act as soft machines, capable of pattern recognition and pattern analysis.

  8. The problem of the universal density functional and the density matrix functional theory

    SciTech Connect

    Bobrov, V. B. Trigger, S. A.

    2013-04-15

    The analysis in this paper shows that the Hohenberg-Kohn theorem is the constellation of two statements: (i) the mathematically rigorous Hohenberg-Kohn lemma, which demonstrates that the same ground-state density cannot correspond to two different potentials of an external field, and (ii) the hypothesis of the existence of the universal density functional. Based on the obtained explicit expression for the nonrel-ativistic particle energy in a local external field, we prove that the energy of the system of more than two non-interacting electrons cannot be a functional of the inhomogeneous density. This result is generalized to the system of interacting electrons. It means that the Hohenberg-Kohn lemma cannot provide justification of the universal density functional for fermions. At the same time, statements of the density functional theory remain valid when considering any number of noninteracting ground-state bosons due to the Bose condensation effect. In the framework of the density matrix functional theory, the hypothesis of the existence of the universal density matrix functional corresponds to the cases of noninteracting particles and to interaction in the Hartree-Fock approximation.

  9. Properties of the cosmological density distribution function

    NASA Astrophysics Data System (ADS)

    Bernardeau, Francis; Kofman, Lev

    1995-04-01

    The properties of the probability distribution function (PDF) of the cosmological continuous density field are studied. We focus our analysis on the quasi-linear regime where various calculations, based on dynamically motivated methods, have been presented: either by using the Zel'dovich approximation (ZA) or by using the perturbation theory to evaluate the behavior of the moments of the distribution function. We show how these two approaches are related to each other and that they can be used in a complementary way. For that respect, the one-dimensional dynamics, where the ZA is exact solution, has first been used as a testing ground. In particular, we show that, when the density PDF obtained with the ZA is regularized, its various moments exhibit the behavior expected by the perturbation theory applied to the ZA. We show that ZA approach can be used for arbitrary initial conditions (not only Gaussian) and that the nonlinear evolution of the moments can be obtained. The perturbation theory can be used for the exact dynamics. We take into account the final filtering of the density field both for ZA and perturbation theory. Applying these techniques, we got the generating function of the moments for the one-dimensional dynamics, the three-dimensional ZA, with and without smoothing effects. We also suggest methods to build PDFs. One is based on the Laplace inverse transform of the moment generating function. The other, the Edgeworth expansion, is obtained when the previous generating function is truncated at a given order and allows evaluation of the PDF out of limited number of moments. It provides insight on the relationship between the moments and the shape of the density PDF. In particular, it provides an alternative method to evaluate the skewness and kurtosis by measuring the PDF around its maximum. Eventually, results obtained from a numerical simulation with cold dark matter initial conditions have been used to validate the accuracy of the considered

  10. Pseudospectral time-dependent density functional theory

    NASA Astrophysics Data System (ADS)

    Ko, Chaehyuk; Malick, David K.; Braden, Dale A.; Friesner, Richard A.; Martínez, Todd J.

    2008-03-01

    Time-dependent density functional theory (TDDFT) is implemented within the Tamm-Dancoff approximation (TDA) using a pseudospectral approach to evaluate two-electron repulsion integrals. The pseudospectral approximation uses a split representation with both spectral basis functions and a physical space grid to achieve a reduction in the scaling behavior of electronic structure methods. We demonstrate here that exceptionally sparse grids may be used in the excitation energy calculation, following earlier work employing the pseudospectral approximation for determining correlation energies in wavefunction-based methods with similar conclusions. The pseudospectral TDA-TDDFT method is shown to be up to ten times faster than a conventional algorithm for hybrid functionals without sacrificing chemical accuracy.

  11. Optimization of constrained density functional theory

    NASA Astrophysics Data System (ADS)

    O'Regan, David D.; Teobaldi, Gilberto

    2016-07-01

    Constrained density functional theory (cDFT) is a versatile electronic structure method that enables ground-state calculations to be performed subject to physical constraints. It thereby broadens their applicability and utility. Automated Lagrange multiplier optimization is necessary for multiple constraints to be applied efficiently in cDFT, for it to be used in tandem with geometry optimization, or with molecular dynamics. In order to facilitate this, we comprehensively develop the connection between cDFT energy derivatives and response functions, providing a rigorous assessment of the uniqueness and character of cDFT stationary points while accounting for electronic interactions and screening. In particular, we provide a nonperturbative proof that stable stationary points of linear density constraints occur only at energy maxima with respect to their Lagrange multipliers. We show that multiple solutions, hysteresis, and energy discontinuities may occur in cDFT. Expressions are derived, in terms of convenient by-products of cDFT optimization, for quantities such as the dielectric function and a condition number quantifying ill definition in multiple constraint cDFT.

  12. Recent progress in density functional theory

    NASA Astrophysics Data System (ADS)

    Truhlar, Donald

    2014-03-01

    Ongoing work involves several areas of density functional theory: new methods for computing electronic excitation energies, including a new way to remove spin contamination in the spin-flip Tamm-Dancoff approximation and a configuration-interaction-corrected Tamm-Dancoff Approximation for treating conical intersections; new ways to treat open-shell states, including a reinterpreted broken-symmetry method and multi-configuration Kohn-Sham theory; a new exchange-correlation functional; new tests of density functional theory against databases for electronic transition energies and molecules and solids containing metal atoms; and applications. A selection of results will be presented. I am grateful to the following collaborators for contributions to the ongoing work: Boris Averkiev, Rebecca Carlson, Laura Fernandez, Laura Gagliardi, Chad Hoyer, Francesc Illas, Miho Isegawa, Shaohong Li, Giovanni Li Manni, Sijie Luo, Dongxia Ma, Remi Maurice, Rubén Means-Pañeda, Roberto Peverati, Nora Planas, Prasenjit Seal, Pragya Verma, Bo Wang, Xuefei Xu, Ke R. Yang, Haoyu Yu, Wenjing Zhang, and Jingjing Zheng. Supported in part by the AFOSR and U.S. DOE.

  13. Modulation Based on Probability Density Functions

    NASA Technical Reports Server (NTRS)

    Williams, Glenn L.

    2009-01-01

    A proposed method of modulating a sinusoidal carrier signal to convey digital information involves the use of histograms representing probability density functions (PDFs) that characterize samples of the signal waveform. The method is based partly on the observation that when a waveform is sampled (whether by analog or digital means) over a time interval at least as long as one half cycle of the waveform, the samples can be sorted by frequency of occurrence, thereby constructing a histogram representing a PDF of the waveform during that time interval.

  14. Adiabatic corrections to density functional theory energies and wave functions.

    PubMed

    Mohallem, José R; Coura, Thiago de O; Diniz, Leonardo G; de Castro, Gustavo; Assafrão, Denise; Heine, Thomas

    2008-09-25

    The adiabatic finite-nuclear-mass-correction (FNMC) to the electronic energies and wave functions of atoms and molecules is formulated for density-functional theory and implemented in the deMon code. The approach is tested for a series of local and gradient corrected density functionals, using MP2 results and diagonal-Born-Oppenheimer corrections from the literature for comparison. In the evaluation of absolute energy corrections of nonorganic molecules the LDA PZ81 functional works surprisingly better than the others. For organic molecules the GGA BLYP functional has the best performance. FNMC with GGA functionals, mainly BLYP, show a good performance in the evaluation of relative corrections, except for nonorganic molecules containing H atoms. The PW86 functional stands out with the best evaluation of the barrier of linearity of H2O and the isotopic dipole moment of HDO. In general, DFT functionals display an accuracy superior than the common belief and because the corrections are based on a change of the electronic kinetic energy they are here ranked in a new appropriate way. The approach is applied to obtain the adiabatic correction for full atomization of alcanes C(n)H(2n+2), n = 4-10. The barrier of 1 mHartree is approached for adiabatic corrections, justifying its insertion into DFT. PMID:18537228

  15. Level densities and γ-ray strength functions in Sn isotopes

    NASA Astrophysics Data System (ADS)

    Toft, H. K.; Larsen, A. C.; Agvaanluvsan, U.; Bürger, A.; Guttormsen, M.; Mitchell, G. E.; Nyhus, H. T.; Schiller, A.; Siem, S.; Syed, N. U. H.; Voinov, A.

    2010-06-01

    The nuclear level densities of Sn118,119 and the γ-ray strength functions of Sn116,118,119 below the neutron separation energy are extracted with the Oslo method using the (He3,αγ) and (He3,He3'γ) reactions. The level-density function of Sn119 displays steplike structures. The microcanonical entropies are deduced from the level densities, and the single neutron entropy of Sn119 is determined to be 1.7 ± 0.2 kB. Results from a combinatorial model support the interpretation that some of the low-energy steps in the level density function are caused by neutron pair breaking. An enhancement in all the γ-ray strength functions of Sn116-119, compared to standard models for radiative strength, is observed for the γ-ray energy region of ≃4-11 MeV. These small resonances all have a centroid energy of 8.0(1) MeV and an integrated strength corresponding to 1.7(9)% of the classical Thomas-Reiche-Kuhn sum rule. The Sn resonances may be due to electric dipole neutron skin oscillations or to an enhancement of the giant magnetic dipole resonance.

  16. Density-functional theory of superconductivity

    NASA Astrophysics Data System (ADS)

    Gross, E. K. U.

    2008-03-01

    A prominent challenge of modern condensed-matter theory is to predict reliably material-specific properties of superconductors, such as the critical temperature. The traditional model of Bardeen, Cooper and Schrieffer (BCS) properly describes the universal features that all conventional superconductors have in common, but it is not able to make accurate predictions of material-specific properties. To tackle this problem, a density-functional formalism has been developed [1] which describes superconductors in thermal equilibrium in terms of three quantities: the ordinary density, the superconducting order parameter, and the nuclear N-body density. These three ``densities'' are determined self-consistently through a set of Kohn-Sham equations. Approximations of the universal exchange-correlation functional are derived on the basis of many-body perturbation theory. In this way, a true ab-initio description is achieved which does not contain any adjustable parameters such as the μ* of Eliashberg theory. Numerical results for the critical temperature, the isotope effect, the gap function and the jump of the specific heat will be presented for simple metals, for MgB2 [2] and CaBeSi, and for calcium intercalated graphite (CaC6) [3]. Furthermore, results for Li, Al, K, and H under pressure will be discussed. The calculations explain why Li and Al behave very differently, leading to a strong enhancement of superconductivity for Li and to a clear suppression for Al with increasing pressure [4]. For K we predict a behavior similar to Li, i.e. a strong increase of Tc with increasing pressure. Finally, hydrogen is found to be a three-gap superconductor whose critical temperature increases with increasing pressure until about 100K (at 500 GPa). [1] M. Lüders, M.A.L. Marques, N.N. Lathiotakis, A. Floris,G. Profeta, L. Fast, A.Continenza, S. Massidda, E.K.U. Gross, PRB 72, 024545 (2005). [2] A. Floris, G. Profeta, N.N. Lathiotakis, M. Lüders, M.A.L. Marques, C. Franchini, E

  17. Local spin analyses using density functional theory

    NASA Astrophysics Data System (ADS)

    Abate, Bayileyegn; Peralta, Juan

    Local spin analysis is a valuable technique in computational investigations magnetic interactions on mono- and polynuclear transition metal complexes, which play vital roles in catalysis, molecular magnetism, artificial photosynthesis, and several other commercially important materials. The relative size and complex electronic structure of transition metal complexes often prohibits the use of multi-determinant approaches, and hence, practical calculations are often limited to single-determinant methods. Density functional theory (DFT) has become one of the most successful and widely used computational tools for the electronic structure study of complex chemical systems; transition metal complexes in particular. Within the DFT formalism, a more flexible and complete theoretical modeling of transition metal complexes can be achieved by considering noncollinear spins, in which the spin density is 'allowed to' adopt noncollinear structures in stead of being constrained to align parallel/antiparallel to a universal axis of magnetization. In this meeting, I will present local spin analyses results obtained using different DFT functionals. Local projection operators are used to decompose the expectation value of the total spin operator; first introduced by Clark and Davidson.

  18. COMBINATORIAL MATERIALS SYNTHESIS AND SCREENING: An Integrated Materials Chip Approach to Discovery and Optimization of Functional Materials

    NASA Astrophysics Data System (ADS)

    Xiang, X.-D.

    1999-08-01

    Combinatorial materials synthesis methods and high throughput evaluation techniques have been developed to accelerate the process of materials discovery and optimization. Analogous to integrated circuit chips, integrated materials chips containing thousands, possibly millions, of different compounds/materials, often in the form of high-quality epitaxial thin film can be fabricated and screened for interesting physical or chemical properties. Microspot X-ray methods, various optical measurement techniques, and a novel evanescent microwave microscope have been used to characterize the structural, optical, magnetic, and electrical properties of samples on materials chips. These techniques are routinely used to discover and optimize luminescent, ferroelectric, dielectric, and magnetic materials.

  19. Chemistry by Way of Density Functional Theory

    NASA Technical Reports Server (NTRS)

    Bauschlicher, Charles W., Jr.; Ricca, Alessandra; Partridge, Harry; Langohff, Stephen R.; Arnold, James O. (Technical Monitor)

    1996-01-01

    In this work we demonstrate that density functional theory (DFT) methods make an important contribution to understanding chemical systems and are an important additional method for the computational chemist. We report calibration calculations obtained with different functionals for the 55 G2 molecules to justify our selection of the B3LYP functional. We show that accurate geometries and vibrational frequencies obtained at the B3LYP level can be combined with traditional methods to simplify the calculation of accurate heats of formation. We illustrate the application of the B3LYP approach to a variety of chemical problems from the vibrational frequencies of polycyclic aromatic hydrocarbons to transition metal systems. We show that the B3LYP method typically performs better than the MP2 method at a significantly lower computational cost. Thus the B3LYP method allows us to extend our studies to much larger systems while maintaining a high degree of accuracy. We show that for transition metal systems, the B3LYP bond energies are typically of sufficient accuracy that they can be used to explain experimental trends and even differentiate between different experimental values. We show that for boron clusters the B3LYP energetics are not as good as for many of the other systems presented, but even in this case the B3LYP approach is able to help understand the experimental trends.

  20. Semiconductor Thermochemistry in Density Functional Calculations

    SciTech Connect

    Lany, S.

    2008-01-01

    The local-density and generalized gradient approximations (LDA and GGA) to density functional theory (DFT) exhibit incomplete error cancellation when energy differences are taken between chemically dissimilar systems. This energy inconsistency is manifested, e.g., in the tendency to underestimate the heat (enthalpy) of formation of semiconducting and insulating compounds in LDA and, even more so, in GGA. Considering a set of 61 compounds that can be formed from 14 elements (cations: Cu, Mg, Ca, Zn, Cd, Al, Ga, and In; anions: N, P, As, O, S, and Se), optimized elemental reference energies are determined by least-squares error minimization of an overdetermined set of linear equations. These elemental energies are 'optimally consistent' with the DFT energies of the semiconductor compounds and imply corrections of up to 1 eV compared to the respective LDA or GGA energies. While these 'corrections' are not to be understood to yield the correct absolute total energies of the elements, they are proposed to give appropriate bounds for the chemical potentials for thermodynamic processes in semiconductors and insulators, such as, e.g., defect formation, surface reconstruction, or catalytic processes. The present model allows to evaluate thermodynamic processes using DFT energy differences taken only between systems that are expected to show good error cancellation.

  1. Density functional studies of representative pericyclic reactions

    NASA Astrophysics Data System (ADS)

    Carpenter, John E.; Sosa, Carlos P.

    1994-07-01

    Density functional theory (DFT) has traditionally been shunned by computational chemists, but has long seen widespread use in the physics community. Recently, however, DFT has been adopted by the ab initio quantum chemistry community and much activity has been devoted to refining the methodology and exploring the range of its applicability. We investigate the use of DFT (both local (LDF) and non-local (NLDF) spin density approximations) to calculate transition and equilibrium structures for three representative pericyclic reactions: the electrocyclic ring opening reaction of cyclobutene, the [1,5] sigmatropic hydrogen shift reaction in (Z)-1,3-pentadiene, and the Diels-Alder cycloaddition reaction between ethylene and butadiene. LDF theory tends to overemphasize the stability of the ringed structures in each of these reactions. For example, LDF predicts a very low (6 kcal mol -1) barrier to reaction for the Deils-Alder reaction. NLDF theory substantially improves the calculated reaction barrier (20 kcal mol -1), but it is still low with respect to experiment.

  2. Phases of Polonium via Density Functional Theory

    NASA Astrophysics Data System (ADS)

    Verstraete, Matthieu J.

    2010-01-01

    The thermodynamical properties of the main phases of metallic polonium are examined using density functional theory. The exceptional nature of the solid-solid phase transition of α to β Po is underlined: it induces a lowering in symmetry, from cubic to rhombohedral, with increasing temperature. This is explained as the result of a delicate balance between bonding and entropic effects. Overall agreement with existing experimental data is good by state-of-the-art standards. The phonons of Po present Kohn anomalies, and it is shown that the effect of spin-orbit interactions is the inverse of that in normal metals: due to the nonspherical nature of the Fermi Surface, spin-orbit effects reduce nesting and harden most phonon frequencies.

  3. Carrier Modulation Via Waveform Probability Density Function

    NASA Technical Reports Server (NTRS)

    Williams, Glenn L.

    2004-01-01

    Beyond the classic modes of carrier modulation by varying amplitude (AM), phase (PM), or frequency (FM), we extend the modulation domain of an analog carrier signal to include a class of general modulations which are distinguished by their probability density function histogram. Separate waveform states are easily created by varying the pdf of the transmitted waveform. Individual waveform states are assignable as proxies for digital ONEs or ZEROs. At the receiver, these states are easily detected by accumulating sampled waveform statistics and performing periodic pattern matching, correlation, or statistical filtering. No fundamental natural laws are broken in the detection process. We show how a typical modulation scheme would work in the digital domain and suggest how to build an analog version. We propose that clever variations of the modulating waveform (and thus the histogram) can provide simple steganographic encoding.

  4. Carrier Modulation Via Waveform Probability Density Function

    NASA Technical Reports Server (NTRS)

    Williams, Glenn L.

    2006-01-01

    Beyond the classic modes of carrier modulation by varying amplitude (AM), phase (PM), or frequency (FM), we extend the modulation domain of an analog carrier signal to include a class of general modulations which are distinguished by their probability density function histogram. Separate waveform states are easily created by varying the pdf of the transmitted waveform. Individual waveform states are assignable as proxies for digital one's or zero's. At the receiver, these states are easily detected by accumulating sampled waveform statistics and performing periodic pattern matching, correlation, or statistical filtering. No fundamental physical laws are broken in the detection process. We show how a typical modulation scheme would work in the digital domain and suggest how to build an analog version. We propose that clever variations of the modulating waveform (and thus the histogram) can provide simple steganographic encoding.

  5. Effective potential in density matrix functional theory.

    PubMed

    Nagy, A; Amovilli, C

    2004-10-01

    In the previous paper it was shown that in the ground state the diagonal of the spin independent second-order density matrix n can be determined by solving a single auxiliary equation of a two-particle problem. Thus the problem of an arbitrary system with even electrons can be reduced to a two-particle problem. The effective potential of the two-particle equation contains a term v(p) of completely kinetic origin. Virial theorem and hierarchy of equations are derived for v(p) and simple approximations are proposed. A relationship between the effective potential u(p) of the shape function equation and the potential v(p) is established. PMID:15473719

  6. Building a Universal Nuclear Energy Density Functional

    SciTech Connect

    Carlson, Joe A.; Furnstahl, Dick; Horoi, Mihai; Lust, Rusty; Nazaewicc, Witek; Ng, Esmond; Thompson, Ian; Vary, James

    2012-12-30

    During the period of Dec. 1 2006 – Jun. 30, 2012, the UNEDF collaboration carried out a comprehensive study of all nuclei, based on the most accurate knowledge of the strong nuclear interaction, the most reliable theoretical approaches, the most advanced algorithms, and extensive computational resources, with a view towards scaling to the petaflop platforms and beyond. The long-term vision initiated with UNEDF is to arrive at a comprehensive, quantitative, and unified description of nuclei and their reactions, grounded in the fundamental interactions between the constituent nucleons. We seek to replace current phenomenological models of nuclear structure and reactions with a well-founded microscopic theory that delivers maximum predictive power with well-quantified uncertainties. Specifically, the mission of this project has been three-fold:  First, to find an optimal energy density functional (EDF) using all our knowledge of the nucleonic Hamiltonian and basic nuclear properties;  Second, to apply the EDF theory and its extensions to validate the functional using all the available relevant nuclear structure and reaction data;  Third, to apply the validated theory to properties of interest that cannot be measured, in particular the properties needed for reaction theory.

  7. Improving Density Functionals with Quantum Harmonic Oscillators

    NASA Astrophysics Data System (ADS)

    Tkatchenko, Alexandre

    2013-03-01

    Density functional theory (DFT) is the most widely used and successful approach for electronic structure calculations. However, one of the pressing challenges for DFT is developing efficient functionals that can accurately capture the omnipresent long-range electron correlations, which determine the structure and stability of many molecules and materials. Here we show that, under certain conditions, the problem of computing the long-range correlation energy of interacting electrons can be mapped to a system of coupled quantum harmonic oscillators (QHOs). The proposed model allows us to synergistically combine concepts from DFT, quantum chemistry, and the widely discussed random-phase approximation for the correlation energy. In the dipole limit, the interaction energy for a system of coupled QHOs can be calculated exactly, thereby leading to an efficient and accurate model for the many-body dispersion energy of complex molecules and materials. The studied examples include intermolecular binding energies, the conformational hierarchy of DNA structures, the geometry and stability of molecular crystals, and supramolecular host-guest complexes (A. Tkatchenko, R. A. DiStasio Jr., R. Car, M. Scheffler, Phys. Rev. Lett. 108, 236402 (2012); R. A. DiStasio Jr., A. von Lilienfeld, A. Tkatchenko, PNAS 109, 14791 (2012); A. Tkatchenko, D. Alfe, K. S. Kim, J. Chem. Theory and Comp. (2012), doi: 10.1021/ct300711r; A. Tkatchenko, A. Ambrosetti, R. A. DiStasio Jr., arXiv:1210.8343v1).

  8. A density functional for sparse matter

    NASA Astrophysics Data System (ADS)

    Langreth, D. C.; Lundqvist, B. I.; Chakarova-Käck, S. D.; Cooper, V. R.; Dion, M.; Hyldgaard, P.; Kelkkanen, A.; Kleis, J.; Kong, Lingzhu; Li, Shen; Moses, P. G.; Murray, E.; Puzder, A.; Rydberg, H.; Schröder, E.; Thonhauser, T.

    2009-02-01

    Sparse matter is abundant and has both strong local bonds and weak nonbonding forces, in particular nonlocal van der Waals (vdW) forces between atoms separated by empty space. It encompasses a broad spectrum of systems, like soft matter, adsorption systems and biostructures. Density-functional theory (DFT), long since proven successful for dense matter, seems now to have come to a point, where useful extensions to sparse matter are available. In particular, a functional form, vdW-DF (Dion et al 2004 Phys. Rev. Lett. 92 246401; Thonhauser et al 2007 Phys. Rev. B 76 125112), has been proposed for the nonlocal correlations between electrons and applied to various relevant molecules and materials, including to those layered systems like graphite, boron nitride and molybdenum sulfide, to dimers of benzene, polycyclic aromatic hydrocarbons (PAHs), doped benzene, cytosine and DNA base pairs, to nonbonding forces in molecules, to adsorbed molecules, like benzene, naphthalene, phenol and adenine on graphite, alumina and metals, to polymer and carbon nanotube (CNT) crystals, and hydrogen storage in graphite and metal-organic frameworks (MOFs), and to the structure of DNA and of DNA with intercalators. Comparison with results from wavefunction calculations for the smaller systems and with experimental data for the extended ones show the vdW-DF path to be promising. This could have great ramifications.

  9. Semiclassical origins of density functional theory

    NASA Astrophysics Data System (ADS)

    Burke, Kieron

    2009-03-01

    Until the seminal work of Hohenberg, Kohn, and Sham of the mid 60's, most density functional theory (DFT) was derived from semiclassical approximations. This non-empirical approach shows an intrinsic difference between solids (for which DFT was originally developed) and molecules, and explains many of its more mysterious manifestations. For example, the success of DFT for molecules has nothing to do with the uniform gas. Results include [1] a derivation of the empirical parameter in the B88 exchange functional, [2] PBEsol, a new GGA that restores the exchange gradient expansion and improves lattice constants in solids, [3] a novel approach to ``orbital-free'' DFT that, in preliminary tests, is 40 times more accurate than its DFT counterpart. The talk is aimed at a general theoretical audience. Detailed technical knowledge of DFT is neither needed, nor desirable. [4pt] [1] J.P. Perdew, L.A. Constantin, E. Sagvolden, and KB, Phys. Rev. Lett. 97, 223002 (2006). [0pt] [2] J.P. Perdew, A. Ruzsinszky, G.I. Csonka, O.A. Vydrov, G.E. Scuseria, L.A. Constantin, X. Zhou, and KB, Phys. Rev. Lett. 100, 136406 (2008). [0pt] [3] Peter Elliott, Donghyung Lee, Attila Cangi, KB, Phys. Rev. Lett. 100, 256406 (2008).

  10. Rationale for switching to nonlocal functionals in density functional theory

    NASA Astrophysics Data System (ADS)

    Lazić, P.; Atodiresei, N.; Caciuc, V.; Brako, R.; Gumhalter, B.; Blügel, S.

    2012-10-01

    Density functional theory (DFT) has been steadily improving over the past few decades, becoming the standard tool for electronic structure calculations. The early local functionals (LDA) were eventually replaced by more accurate semilocal functionals (GGA) which are in use today. A major persisting drawback is the lack of the nonlocal correlation which is at the core of dispersive (van der Waals) forces, so that a large and important class of systems remains outside the scope of DFT. The vdW-DF correlation functional of Langreth and Lundqvist, published in 2004, was the first nonlocal functional which could be easily implemented. Beyond expectations, the nonlocal functional has brought significant improvement to systems that were believed not to be sensitive to nonlocal correlations. In this paper, we use the example of graphene nanodomes growing on the Ir(111) surface, where with an increase of the size of the graphene islands the character of the bonding changes from strong chemisorption towards almost pure physisorption. We demonstrate how the seamless character of the vdW-DF functionals makes it possible to treat all regimes self-consistently, proving to be a systematic and consistent improvement of DFT regardless of the nature of bonding. We also discuss the typical surface science example of CO adsorption on (111) surfaces of metals, which shows that the nonlocal correlation may also be crucial for strongly chemisorbed systems. We briefly discuss open questions, in particular the choice of the most appropriate exchange part of the functional. As the vdW-DF begins to appear implemented self-consistently in a number of popular DFT codes, with numerical costs close to the GGA calculations, we draw the attention of the DFT community to the advantages and benefits of the adoption of this new class of functionals.

  11. Density Functional Theory study of the equilibrium density of water at normal conditions

    NASA Astrophysics Data System (ADS)

    Wang, Jue; Roman-Perez, Guillermo; Soler, Jose M.; Artacho, Emilio; Fernandez-Serra, Marivi

    2010-03-01

    Ab initio molecular dynamics of liquid water with the use of density functional theory (DFT) currently underperform experimental equilibrium density 1g/cm^3 under room temperature. At constant density, not much is known about the equilibrium density of commonly used GGA functionals in liquid water simulations. We present a DFT-based AIMD study of liquid water at different densities and analyze the structure and diffusivity of water with different exchange and correlation functionals. We show that all current GGA functionals fail to reproduce experimental density, however, the explicit description of long range correlations through a Van der Walls density functional (DRSLL)footnotetextM. Dion, H. Rydberg, E. Schr"oder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004) can potentially transform our current understanding of the structure of liquid water. Our results shows that this new functional improves density, with only 2% error to experiment. But it underperforms GGA functionals in terms of structure.

  12. Flexoelectricity from density-functional perturbation theory

    NASA Astrophysics Data System (ADS)

    Stengel, Massimiliano

    2013-11-01

    We derive the complete flexoelectric tensor, including electronic and lattice-mediated effects, of an arbitrary insulator in terms of the microscopic linear response of the crystal to atomic displacements. The basic ingredient, which can be readily calculated from first principles in the framework of density-functional perturbation theory, is the quantum-mechanical probability current response to a long-wavelength acoustic phonon. Its second-order Taylor expansion in the wave vector q around the Γ (q=0) point in the Brillouin zone naturally yields the flexoelectric tensor. At order one in q we recover Martin's theory of piezoelectricity [Martin, Phys. Rev. B 5, 1607 (1972)], thus providing an alternative derivation thereof. To put our derivations on firm theoretical grounds, we perform a thorough analysis of the nonanalytic behavior of the dynamical matrix and other response functions in a vicinity of Γ. Based on this analysis, we find that there is an ambiguity in the specification of the “zero macroscopic field” condition in the flexoelectric case; such arbitrariness can be related to an analytic band-structure term, in close analogy to the theory of deformation potentials. As a by-product, we derive a rigorous generalization of the Cochran-Cowley formula [Cochran and Cowley, J. Phys. Chem. Solids 23, 447 (1962)] to higher orders in q. This can be of great utility in building reliable atomistic models of electromechanical phenomena, as well as for improving the accuracy of the calculation of phonon dispersion curves. Finally, we discuss the physical interpretation of the various contributions to the flexoelectric response, either in the static or dynamic regime, and we relate our findings to earlier theoretical works on the subject.

  13. Dynamical density functional theory for microswimmers

    NASA Astrophysics Data System (ADS)

    Menzel, Andreas M.; Saha, Arnab; Hoell, Christian; Löwen, Hartmut

    2016-01-01

    Dynamical density functional theory (DDFT) has been successfully derived and applied to describe on one hand passive colloidal suspensions, including hydrodynamic interactions between individual particles. On the other hand, active "dry" crowds of self-propelled particles have been characterized using DDFT. Here, we go one essential step further and combine these two approaches. We establish a DDFT for active microswimmer suspensions. For this purpose, simple minimal model microswimmers are introduced. These microswimmers self-propel by setting the surrounding fluid into motion. They hydrodynamically interact with each other through their actively self-induced fluid flows and via the common "passive" hydrodynamic interactions. An effective soft steric repulsion is also taken into account. We derive the DDFT starting from common statistical approaches. Our DDFT is then tested and applied by characterizing a suspension of microswimmers, the motion of which is restricted to a plane within a three-dimensional bulk fluid. Moreover, the swimmers are confined by a radially symmetric trapping potential. In certain parameter ranges, we find rotational symmetry breaking in combination with the formation of a "hydrodynamic pumping state," which has previously been observed in the literature as a result of particle-based simulations. An additional instability of this pumping state is revealed.

  14. Density functional theory in the solid state

    PubMed Central

    Hasnip, Philip J.; Refson, Keith; Probert, Matt I. J.; Yates, Jonathan R.; Clark, Stewart J.; Pickard, Chris J.

    2014-01-01

    Density functional theory (DFT) has been used in many fields of the physical sciences, but none so successfully as in the solid state. From its origins in condensed matter physics, it has expanded into materials science, high-pressure physics and mineralogy, solid-state chemistry and more, powering entire computational subdisciplines. Modern DFT simulation codes can calculate a vast range of structural, chemical, optical, spectroscopic, elastic, vibrational and thermodynamic phenomena. The ability to predict structure–property relationships has revolutionized experimental fields, such as vibrational and solid-state NMR spectroscopy, where it is the primary method to analyse and interpret experimental spectra. In semiconductor physics, great progress has been made in the electronic structure of bulk and defect states despite the severe challenges presented by the description of excited states. Studies are no longer restricted to known crystallographic structures. DFT is increasingly used as an exploratory tool for materials discovery and computational experiments, culminating in ex nihilo crystal structure prediction, which addresses the long-standing difficult problem of how to predict crystal structure polymorphs from nothing but a specified chemical composition. We present an overview of the capabilities of solid-state DFT simulations in all of these topics, illustrated with recent examples using the CASTEP computer program. PMID:24516184

  15. Bone mineral density, adiposity, and cognitive functions

    PubMed Central

    Sohrabi, Hamid R.; Bates, Kristyn A.; Weinborn, Michael; Bucks, Romola S.; Rainey-Smith, Stephanie R.; Rodrigues, Mark A.; Bird, Sabine M.; Brown, Belinda M.; Beilby, John; Howard, Matthew; Criddle, Arthur; Wraith, Megan; Taddei, Kevin; Martins, Georgia; Paton, Athena; Shah, Tejal; Dhaliwal, Satvinder S.; Mehta, Pankaj D.; Foster, Jonathan K.; Martins, Ian J.; Lautenschlager, Nicola T.; Mastaglia, Francis; Laws, Simon M.; Martins, Ralph N.

    2015-01-01

    Cognitive decline and dementia due to Alzheimer's disease (AD) have been associated with genetic, lifestyle, and environmental factors. A number of potentially modifiable risk factors should be taken into account when preventive or ameliorative interventions targeting dementia and its preclinical stages are investigated. Bone mineral density (BMD) and body composition are two such potentially modifiable risk factors, and their association with cognitive decline was investigated in this study. 164 participants, aged 34–87 years old (62.78 ± 9.27), were recruited for this longitudinal study and underwent cognitive and clinical examinations at baseline and after 3 years. Blood samples were collected for apolipoprotein E (APOE) genotyping and dual energy x-ray absorptiometry (DXA) was conducted at the same day as cognitive assessment. Using hierarchical regression analysis, we found that BMD and lean body mass, as measured using DXA were significant predictors of episodic memory. Age, gender, APOE status, and premorbid IQ were controlled for. Specifically, the List A learning from California Verbal Learning Test was significantly associated with BMD and lean mass both at baseline and at follow up assessment. Our findings indicate that there is a significant association between BMD and lean body mass and episodic verbal learning. While the involvement of modifiable lifestyle factors in human cognitive function has been examined in different studies, there is a need for further research to understand the potential underlying mechanisms. PMID:25741279

  16. Bone mineral density, adiposity, and cognitive functions.

    PubMed

    Sohrabi, Hamid R; Bates, Kristyn A; Weinborn, Michael; Bucks, Romola S; Rainey-Smith, Stephanie R; Rodrigues, Mark A; Bird, Sabine M; Brown, Belinda M; Beilby, John; Howard, Matthew; Criddle, Arthur; Wraith, Megan; Taddei, Kevin; Martins, Georgia; Paton, Athena; Shah, Tejal; Dhaliwal, Satvinder S; Mehta, Pankaj D; Foster, Jonathan K; Martins, Ian J; Lautenschlager, Nicola T; Mastaglia, Francis; Laws, Simon M; Martins, Ralph N

    2015-01-01

    Cognitive decline and dementia due to Alzheimer's disease (AD) have been associated with genetic, lifestyle, and environmental factors. A number of potentially modifiable risk factors should be taken into account when preventive or ameliorative interventions targeting dementia and its preclinical stages are investigated. Bone mineral density (BMD) and body composition are two such potentially modifiable risk factors, and their association with cognitive decline was investigated in this study. 164 participants, aged 34-87 years old (62.78 ± 9.27), were recruited for this longitudinal study and underwent cognitive and clinical examinations at baseline and after 3 years. Blood samples were collected for apolipoprotein E (APOE) genotyping and dual energy x-ray absorptiometry (DXA) was conducted at the same day as cognitive assessment. Using hierarchical regression analysis, we found that BMD and lean body mass, as measured using DXA were significant predictors of episodic memory. Age, gender, APOE status, and premorbid IQ were controlled for. Specifically, the List A learning from California Verbal Learning Test was significantly associated with BMD and lean mass both at baseline and at follow up assessment. Our findings indicate that there is a significant association between BMD and lean body mass and episodic verbal learning. While the involvement of modifiable lifestyle factors in human cognitive function has been examined in different studies, there is a need for further research to understand the potential underlying mechanisms. PMID:25741279

  17. Photostriction in Ferroelectrics from Density Functional Theory

    NASA Astrophysics Data System (ADS)

    Paillard, Charles; Xu, Bin; Dkhil, Brahim; Geneste, Grégory; Bellaiche, L.

    2016-06-01

    An ab initio procedure allowing the computation of the deformation of ferroelectric-based materials under light is presented. This numerical scheme consists in structurally relaxing the system under the constraint of a fixed ne concentration of electrons photoexcited into a specific conduction band edge state from a chosen valence band state, via the use of a constrained density functional theory method. The resulting change in lattice constant along a selected crystallographic direction is then calculated for a reasonable estimate of ne. This method is applied to bulk multiferroic BiFeO3 and predicts a photostriction effect of the same order of magnitude than the ones recently observed. A strong dependence of photostrictive response on both the reached conduction state and the crystallographic direction (along which this effect is determined) is also revealed. Furthermore, analysis of the results demonstrates that the photostriction mechanism mostly originates from the screening of the spontaneous polarization by the photoexcited electrons in combination with the inverse piezoelectric effect.

  18. Photostriction in Ferroelectrics from Density Functional Theory.

    PubMed

    Paillard, Charles; Xu, Bin; Dkhil, Brahim; Geneste, Grégory; Bellaiche, L

    2016-06-17

    An ab initio procedure allowing the computation of the deformation of ferroelectric-based materials under light is presented. This numerical scheme consists in structurally relaxing the system under the constraint of a fixed n_{e} concentration of electrons photoexcited into a specific conduction band edge state from a chosen valence band state, via the use of a constrained density functional theory method. The resulting change in lattice constant along a selected crystallographic direction is then calculated for a reasonable estimate of n_{e}. This method is applied to bulk multiferroic BiFeO_{3} and predicts a photostriction effect of the same order of magnitude than the ones recently observed. A strong dependence of photostrictive response on both the reached conduction state and the crystallographic direction (along which this effect is determined) is also revealed. Furthermore, analysis of the results demonstrates that the photostriction mechanism mostly originates from the screening of the spontaneous polarization by the photoexcited electrons in combination with the inverse piezoelectric effect. PMID:27367406

  19. The Nuclear Energy Density Functional Formalism

    NASA Astrophysics Data System (ADS)

    Duguet, T.

    The present document focuses on the theoretical foundations of the nuclear energy density functional (EDF) method. As such, it does not aim at reviewing the status of the field, at covering all possible ramifications of the approach or at presenting recent achievements and applications. The objective is to provide a modern account of the nuclear EDF formalism that is at variance with traditional presentations that rely, at one point or another, on a Hamiltonian-based picture. The latter is not general enough to encompass what the nuclear EDF method represents as of today. Specifically, the traditional Hamiltonian-based picture does not allow one to grasp the difficulties associated with the fact that currently available parametrizations of the energy kernel E[g',g] at play in the method do not derive from a genuine Hamilton operator, would the latter be effective. The method is formulated from the outset through the most general multi-reference, i.e. beyond mean-field, implementation such that the single-reference, i.e. "mean-field", derives as a particular case. As such, a key point of the presentation provided here is to demonstrate that the multi-reference EDF method can indeed be formulated in a mathematically meaningful fashion even if E[g',g] does not derive from a genuine Hamilton operator. In particular, the restoration of symmetries can be entirely formulated without making any reference to a projected state, i.e. within a genuine EDF framework. However, and as is illustrated in the present document, a mathematically meaningful formulation does not guarantee that the formalism is sound from a physical standpoint. The price at which the latter can be enforced as well in the future is eventually alluded to.

  20. Combinatorial 3D Mechanical Metamaterials

    NASA Astrophysics Data System (ADS)

    Coulais, Corentin; Teomy, Eial; de Reus, Koen; Shokef, Yair; van Hecke, Martin

    2015-03-01

    We present a class of elastic structures which exhibit 3D-folding motion. Our structures consist of cubic lattices of anisotropic unit cells that can be tiled in a complex combinatorial fashion. We design and 3d-print this complex ordered mechanism, in which we combine elastic hinges and defects to tailor the mechanics of the material. Finally, we use this large design space to encode smart functionalities such as surface patterning and multistability.

  1. A half century of density functional theory

    SciTech Connect

    Zangwill, Andrew

    2015-07-15

    Today’s most popular method for calculating the electronic structure of atoms, molecules, liquids, solids, and plasmas began as a bold hypothesis: The electron density distribution completely characterizes the ground state of a many-electron system.

  2. Ions in solution: Density corrected density functional theory (DC-DFT)

    SciTech Connect

    Kim, Min-Cheol; Sim, Eunji; Burke, Kieron

    2014-05-14

    Standard density functional approximations often give questionable results for odd-electron radical complexes, with the error typically attributed to self-interaction. In density corrected density functional theory (DC-DFT), certain classes of density functional theory calculations are significantly improved by using densities more accurate than the self-consistent densities. We discuss how to identify such cases, and how DC-DFT applies more generally. To illustrate, we calculate potential energy surfaces of HO·Cl{sup −} and HO·H{sub 2}O complexes using various common approximate functionals, with and without this density correction. Commonly used approximations yield wrongly shaped surfaces and/or incorrect minima when calculated self consistently, while yielding almost identical shapes and minima when density corrected. This improvement is retained even in the presence of implicit solvent.

  3. Combinatorial Geometry Printer Plotting.

    Energy Science and Technology Software Center (ESTSC)

    1987-01-05

    Picture generates plots of two-dimensional slices through the three-dimensional geometry described by the combinatorial geometry (CG) package used in such codes as MORSE and QAD-CG. These plots are printed on a standard line printer.

  4. Combinatorial Interdependence in Lottery

    ERIC Educational Resources Information Center

    Helman, Danny

    2005-01-01

    This paper examines a real life question of gamble facing lottery players. Combinatorial dependence plays a central role in shaping the game probabilistic structure, but might not carry the merited weight in punters' considerations.

  5. a Renormalization Group Calculation of the Velocity - and Density-Density Correlation Functions.

    NASA Astrophysics Data System (ADS)

    Cowan, Mark Timothy

    The velocity-velocity correlation function of a free field theory is obtained. The renormalization group, along with a 4-varepsilon expansion, is then used to find the leading order behavior of the velocity-velocity correlation function for an interacting field theory in the high temperature phase near the critical point. The details of the calculation of the density-density correlation function for Hedgehogs, in the context of a free field theory, is presented next. Finally the renormalization group, along with a 4-varepsilon expansion, is used to find the leading order behavior of the density-density correlation function for Hedgehogs in an interacting field theory near the critical point.

  6. Investigating the Mg-Si Binary System via Combinatorial Sputter Deposition As High Energy Density Anodes for Lithium-Ion Batteries.

    PubMed

    Schmuelling, Guido; Winter, Martin; Placke, Tobias

    2015-09-16

    Mg-Si thin films with various elemental compositions ranging from 0≤x≤1 in MgxSi(1-x) were obtained via combinatorial magnetron sputter deposition of Si and Mg in order to improve the electrochemical lithiation/delithiation process of pure Si by embedding Si in an active Mg-Si matrix. Scanning electron microscopy, energy-dispersive X-ray spectroscopy, X-ray diffraction, and Raman spectroscopy methods were used to investigate the morphology, stoichiometry, and structure of the different thin film samples. Constant current charge/discharge cycling revealed significant electrochemical changes depending on the Mg content in comparison to the pure Si active material improving the capacity retention to 96% over 400 cycles. PMID:26313948

  7. Density functional theory for Yukawa fluids

    NASA Astrophysics Data System (ADS)

    Hatlo, Marius M.; Banerjee, Priyanka; Forsman, Jan; Lue, Leo

    2012-08-01

    We develop an approximate field theory for particles interacting with a generalized Yukawa potential. This theory improves and extends a previous splitting field theory, originally developed for counterions around a fixed charge distribution. The resulting theory bridges between the second virial approximation, which is accurate at low particle densities, and the mean-field approximation, accurate at high densities. We apply this theory to charged, screened ions in bulk solution, modeled to interact with a Yukawa potential; the theory is able to accurately reproduce the thermodynamic properties of the system over a broad range of conditions. The theory is also applied to "dressed counterions," interacting with a screened electrostatic potential, contained between charged plates. It is found to work well from the weak coupling to the strong coupling limits. The theory is able to reproduce the counterion profiles and force curves for closed and open systems obtained from Monte Carlo simulations.

  8. Fermionic density functional at a Feshbach resonance

    SciTech Connect

    Seidl, Michael; Bhaduri, Rajat K.

    2007-05-15

    We consider a dilute gas of neutral unpolarized fermionic atoms at zero temperature. The atoms interact via a short-range (tunable) attractive interaction. We demonstrate analytically a curious property of the gas at unitarity. Namely, the correlation energy of the gas, evaluated by second-order perturbation theory, has the same density dependence as the first-order exchange energy, and the two almost exactly cancel each other at a Feshbach resonance irrespective of the shape of the potential, provided ({mu}r{sub s})>>1. Here ({mu}){sup -1} is the range of the two-body potential, and r{sub s} is defined through the number density, n=3/(4{pi}r{sub s}{sup 3}). The implications of this result for universality are discussed.

  9. Molecular dynamics simulation of liquid water: Hybrid density functionals

    SciTech Connect

    Todorova, T; Seitsonen, A; Hutter, J; Kuo, W; Mundy, C

    2005-09-12

    The structure, dynamical and electronic properties of liquid water utilizing different hybrid density functionals were tested within the plane wave framework of first principles molecular dynamics simulations. The computational approach, which employs modified functionals with short-ranged Hartree-Fock exchange, was first tested in calculations of the structural and bonding properties of the water dimer and cyclic water trimer. Liquid water simulations were performed at the state point of 350 K at the experimental density. Simulations included three different hybrid functionals, a meta functional, four gradient corrected functionals, the local density and Hartree-Fock approximation. It is found that hybrid functionals are superior in reproducing the experimental structure and dynamical properties as measured by the radial distribution function and self diffusion constant when compared to the pure density functionals. The local density and Hartree-Fock approximations show strongly over- and under-structured liquids, respectively. Hydrogen bond analysis shows that the hybrid functionals give slightly smaller averaged numbers of hydrogen bonds and similar hydrogen bond populations as pure density functionals. The average molecular dipole moments in the liquid from the three hybrid functionals are lower than from the corresponding pure density functionals.

  10. Exact conditions on the temperature dependence of density functionals

    NASA Astrophysics Data System (ADS)

    Burke, K.; Smith, J. C.; Grabowski, P. E.; Pribram-Jones, A.

    2016-05-01

    Universal exact conditions guided the construction of most ground-state density functional approximations in use today. We derive the relation between the entropy and Mermin free energy density functionals for thermal density functional theory. Both the entropy and sum of kinetic and electron-electron repulsion functionals are shown to be monotonically increasing with temperature, while the Mermin functional is concave downwards. Analogous relations are found for both exchange and correlation. The importance of these conditions is illustrated in two extremes: the Hubbard dimer and the uniform gas.

  11. Phonons in nonlocal van der Waals density functional theory

    NASA Astrophysics Data System (ADS)

    Sabatini, Riccardo; Küçükbenli, Emine; Pham, Cong Huy; de Gironcoli, Stefano

    2016-06-01

    We extend the formulation of density functional perturbation theory to treat nonlocal density functionals, accounting for van der Waals interactions, in a rigorous and efficient way. We provide a general formalism, suitable for any functional in this family, and give specific equations for the most widely used ones. We then study the lattice dynamics of graphite, comparing several nonlocal functionals and the local density approximation, showing that our recent revision of the VV10 functional [R. Sabatini et al., Phys. Rev. B 87, 041108(R) (2013), 10.1103/PhysRevB.87.041108] gives the best comparison with experiments.

  12. Density-potential mapping in time-dependent density-functional theory

    SciTech Connect

    Maitra, N. T.; Todorov, T. N.; Woodward, C.; Burke, K.

    2010-04-15

    The key questions of uniqueness and existence in time-dependent density-functional theory are usually formulated only for potentials and densities that are analytic in time. Simple examples, standard in quantum mechanics, lead, however, to nonanalyticities. We reformulate these questions in terms of a nonlinear Schroedinger equation with a potential that depends nonlocally on the wave function.

  13. Spin constraints on nuclear energy density functionals

    NASA Astrophysics Data System (ADS)

    Robledo, L. M.; Bernard, R. N.; Bertsch, G. F.

    2014-02-01

    The Gallagher-Moszkowski rule in the spectroscopy of odd-odd nuclei imposes a new spin constraint on the energy functionals for self-consistent mean field theory. The commonly used parametrization of the effective three-body interaction in the Gogny and Skyrme families of energy functionals is ill suited to satisfy the spin constraint. In particular, the Gogny parametrization of the three-body interaction has the spin dependence opposite to that required by the observed spectra. The two-body part has a correct sign, but in combination the rule is violated as often as not. We conclude that a new functional form is needed for the effective three-body interaction that can take into better account the different spin-isospin channels of the interaction.

  14. Introduction to Classical Density Functional Theory by a Computational Experiment

    ERIC Educational Resources Information Center

    Jeanmairet, Guillaume; Levy, Nicolas; Levesque, Maximilien; Borgis, Daniel

    2014-01-01

    We propose an in silico experiment to introduce the classical density functional theory (cDFT). Density functional theories, whether quantum or classical, rely on abstract concepts that are nonintuitive; however, they are at the heart of powerful tools and active fields of research in both physics and chemistry. They led to the 1998 Nobel Prize in…

  15. Nuclear collective excitations: A relativistic density functional approach

    NASA Astrophysics Data System (ADS)

    Piekarewicz, J.

    2015-08-01

    Density functional theory provides the most promising, and likely unique, microscopic framework to describe nuclear systems ranging from finite nuclei to neutron stars. Properly optimized energy density functionals define a new paradigm in nuclear theory where predictive capability is possible and uncertainty quantification is demanded. Moreover, density functional theory offers a consistent approach to the linear response of the nuclear ground state. In this paper, we review the fundamental role played by nuclear collective modes in uncovering novel excitations and in guiding the optimization of the density functional. Indeed, without collective excitations the determination of the density functional remains incomplete. Without collective excitations, the equation of state of neutron-rich matter continues to be poorly constrained. We conclude with a discussion of some of the remaining challenges in this field and propose a path forward to address these challenges.

  16. Single and combinatorial chromatin coupling events underlies the function of transcript factor krüppel-like factor 11 in the regulation of gene networks

    PubMed Central

    2014-01-01

    Background Krüppel-like factors (KLFs) are a group of master regulators of gene expression conserved from flies to human. However, scant information is available on either the mechanisms or functional impact of the coupling of KLF proteins to chromatin remodeling machines, a deterministic step in transcriptional regulation. Results and discussion In the current study, we use genome-wide analyses of chromatin immunoprecipitation (ChIP-on-Chip) and Affymetrix-based expression profiling to gain insight into how KLF11, a human transcription factor involved in tumor suppression and metabolic diseases, works by coupling to three co-factor groups: the Sin3-histone deacetylase system, WD40-domain containing proteins, and the HP1-histone methyltransferase system. Our results reveal that KLF11 regulates distinct gene networks involved in metabolism and growth by using single or combinatorial coupling events. Conclusion This study, the first of its type for any KLF protein, reveals that interactions with multiple chromatin systems are required for the full gene regulatory function of these proteins. PMID:24885560

  17. Efficient time-dependent density functional theory approximations for hybrid density functionals: analytical gradients and parallelization.

    PubMed

    Petrenko, Taras; Kossmann, Simone; Neese, Frank

    2011-02-01

    In this paper, we present the implementation of efficient approximations to time-dependent density functional theory (TDDFT) within the Tamm-Dancoff approximation (TDA) for hybrid density functionals. For the calculation of the TDDFT/TDA excitation energies and analytical gradients, we combine the resolution of identity (RI-J) algorithm for the computation of the Coulomb terms and the recently introduced "chain of spheres exchange" (COSX) algorithm for the calculation of the exchange terms. It is shown that for extended basis sets, the RIJCOSX approximation leads to speedups of up to 2 orders of magnitude compared to traditional methods, as demonstrated for hydrocarbon chains. The accuracy of the adiabatic transition energies, excited state structures, and vibrational frequencies is assessed on a set of 27 excited states for 25 molecules with the configuration interaction singles and hybrid TDDFT/TDA methods using various basis sets. Compared to the canonical values, the typical error in transition energies is of the order of 0.01 eV. Similar to the ground-state results, excited state equilibrium geometries differ by less than 0.3 pm in the bond distances and 0.5° in the bond angles from the canonical values. The typical error in the calculated excited state normal coordinate displacements is of the order of 0.01, and relative error in the calculated excited state vibrational frequencies is less than 1%. The errors introduced by the RIJCOSX approximation are, thus, insignificant compared to the errors related to the approximate nature of the TDDFT methods and basis set truncation. For TDDFT/TDA energy and gradient calculations on Ag-TB2-helicate (156 atoms, 2732 basis functions), it is demonstrated that the COSX algorithm parallelizes almost perfectly (speedup ~26-29 for 30 processors). The exchange-correlation terms also parallelize well (speedup ~27-29 for 30 processors). The solution of the Z-vector equations shows a speedup of ~24 on 30 processors. The

  18. Efficient time-dependent density functional theory approximations for hybrid density functionals: Analytical gradients and parallelization

    NASA Astrophysics Data System (ADS)

    Petrenko, Taras; Kossmann, Simone; Neese, Frank

    2011-02-01

    In this paper, we present the implementation of efficient approximations to time-dependent density functional theory (TDDFT) within the Tamm-Dancoff approximation (TDA) for hybrid density functionals. For the calculation of the TDDFT/TDA excitation energies and analytical gradients, we combine the resolution of identity (RI-J) algorithm for the computation of the Coulomb terms and the recently introduced "chain of spheres exchange" (COSX) algorithm for the calculation of the exchange terms. It is shown that for extended basis sets, the RIJCOSX approximation leads to speedups of up to 2 orders of magnitude compared to traditional methods, as demonstrated for hydrocarbon chains. The accuracy of the adiabatic transition energies, excited state structures, and vibrational frequencies is assessed on a set of 27 excited states for 25 molecules with the configuration interaction singles and hybrid TDDFT/TDA methods using various basis sets. Compared to the canonical values, the typical error in transition energies is of the order of 0.01 eV. Similar to the ground-state results, excited state equilibrium geometries differ by less than 0.3 pm in the bond distances and 0.5° in the bond angles from the canonical values. The typical error in the calculated excited state normal coordinate displacements is of the order of 0.01, and relative error in the calculated excited state vibrational frequencies is less than 1%. The errors introduced by the RIJCOSX approximation are, thus, insignificant compared to the errors related to the approximate nature of the TDDFT methods and basis set truncation. For TDDFT/TDA energy and gradient calculations on Ag-TB2-helicate (156 atoms, 2732 basis functions), it is demonstrated that the COSX algorithm parallelizes almost perfectly (speedup ˜26-29 for 30 processors). The exchange-correlation terms also parallelize well (speedup ˜27-29 for 30 processors). The solution of the Z-vector equations shows a speedup of ˜24 on 30 processors. The

  19. The densities produced by the density functional theory: Comparison to full configuration interaction

    NASA Astrophysics Data System (ADS)

    Bochevarov, Arteum D.; Friesner, Richard A.

    2008-01-01

    We investigate one of the fundamental observables, electronic charge density, as produced by a number of popular functionals of the density functional theory (DFT): SVWN5, B3LYP, B3LYP, OLYP, O3LYP, BP86, B3P86, O3P86, and PBE using restricted and unrestricted orbitals. Measuring and comparing the quality of the densities could tell us more about the physical soundness of the functional models. The study is performed on the small molecules He, H2, LiH, H4 in an extensive range of correlation-consistent basis sets. We compare DFT densities to those of full configuration interaction (FCI) under the assumption that the FCI density in the largest employed basis set is sufficiently close to the exact one. For LiH and H4, we also compare the DFT densities to those of CCSD. The SVWN5 functional consistently shows the worst performance. The OPTX exchange functional regularly beats the Becke exchange. Among the best performers are all the hybrid functionals, the novel O3P86 being the most accurate in most cases. The popular functional B3LYP was consistently outmatched by O3LYP, and produced, in fact, some of the poorest densities among the hybrids. CCSD was found to produce much more accurate densities than any DFT functional in the case of LiH in equilibrium geometry, but was sometimes outperformed by DFT in the case of slightly stretched H4, where CCSD theory itself starts to break down. Surprisingly, as one stretches the H2 molecule, BP86 and PBE improve the description of density although such behavior is not observed in other systems. We conclude by reasoning how functionals such as B3LYP, despite being quite average for density, could still be very successful in predicting thermodynamic properties.

  20. Molecular density functional theory of water including density-polarization coupling.

    PubMed

    Jeanmairet, Guillaume; Levy, Nicolas; Levesque, Maximilien; Borgis, Daniel

    2016-06-22

    We present a three-dimensional molecular density functional theory of water derived from first-principles that relies on the particle's density and multipolar polarization density and includes the density-polarization coupling. This brings two main benefits: (i) scalar density and vectorial multipolar polarization density fields are much more tractable and give more physical insight than the full position and orientation densities, and (ii) it includes the full density-polarization coupling of water, that is known to be non-vanishing but has never been taken into account. Furthermore, the theory requires only the partial charge distribution of a water molecule and three measurable bulk properties, namely the structure factor and the Fourier components of the longitudinal and transverse dielectric susceptibilities. PMID:27116250

  1. The benchmark of gutzwiller density functional theory in hydrogen systems

    SciTech Connect

    Yao, Y.; Wang, Cai-Zhuang; Ho, Kai-Ming

    2012-02-23

    We propose an approximate form of the exchange-correlation energy functional for the Gutzwiller density functional theory. It satisfies certain physical constraints in both weak and strong electron correlation limits. We benchmark the Gutzwiller density functional approximation in the hydrogen systems, where the static correlation error is shown to be negligible. The good transferability is demonstrated by applications to the hydrogen molecule and some crystal structures.

  2. The Benchmark of Gutzwiller Density Functional Theory in Hydrogen Systems

    SciTech Connect

    Yao, Yongxin; Wang, Cai-Zhuang; Ho, Kai-Ming

    2011-01-13

    We propose an approximate form of the exchange-correlation energy functional for the Gutzwiller density functional theory. It satisfies certain physical constraints in both weak and strong electron correlation limits. We benchmark the Gutzwiller density functional approximation in the hydrogen systems, where the static correlation error is shown to be negligible. The good transferability is demonstrated by applications to the hydrogen molecule and some crystal structures. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012

  3. Generalized van der Waals density functional theory for nonuniform polymers

    SciTech Connect

    Patra, Chandra N.; Yethiraj, Arun

    2000-01-15

    A density functional theory is presented for the effect of attractions on the structure of polymers at surfaces. The theory treats the ideal gas functional exactly, and uses a weighted density approximation for the hard chain contribution to the excess free energy functional. The attractive interactions are treated using a van der Waals approximation. The theory is in good agreement with computer simulations for the density profiles at surfaces for a wide range of densities and temperatures, except for low polymer densities at low temperatures where it overestimates the depletion of chains from the surface. This deficiency is attributed to the neglect of liquid state correlations in the van der Waals term of the free energy functional. (c) 2000 American Institute of Physics.

  4. Exact maps in density functional theory for lattice models

    NASA Astrophysics Data System (ADS)

    Dimitrov, Tanja; Appel, Heiko; Fuks, Johanna I.; Rubio, Angel

    2016-08-01

    In the present work, we employ exact diagonalization for model systems on a real-space lattice to explicitly construct the exact density-to-potential and graphically illustrate the complete exact density-to-wavefunction map that underly the Hohenberg–Kohn theorem in density functional theory. Having the explicit wavefunction-to-density map at hand, we are able to construct arbitrary observables as functionals of the ground-state density. We analyze the density-to-potential map as the distance between the fragments of a system increases and the correlation in the system grows. We observe a feature that gradually develops in the density-to-potential map as well as in the density-to-wavefunction map. This feature is inherited by arbitrary expectation values as functional of the ground-state density. We explicitly show the excited-state energies, the excited-state densities, and the correlation entropy as functionals of the ground-state density. All of them show this exact feature that sharpens as the coupling of the fragments decreases and the correlation grows. We denominate this feature as intra-system steepening and discuss how it relates to the well-known inter-system derivative discontinuity. The inter-system derivative discontinuity is an exact concept for coupled subsystems with degenerate ground state. However, the coupling between subsystems as in charge transfer processes can lift the degeneracy. An important conclusion is that for such systems with a near-degenerate ground state, the corresponding cut along the particle number N of the exact density functionals is differentiable with a well-defined gradient near integer particle number.

  5. Analysis of Combinatorial Epigenomic States.

    PubMed

    Soloway, Paul D

    2016-03-18

    Hundreds of distinct chemical modifications to DNA and histone amino acids have been described. Regulation exerted by these so-called epigenetic marks is vital to normal development, stability of cell identity through mitosis, and nongenetic transmission of traits between generations through meiosis. Loss of this regulation contributes to many diseases. Evidence indicates epigenetic marks function in combinations, whereby a given modification has distinct effects on local genome control, depending on which additional modifications are locally present. This review summarizes emerging methods for assessing combinatorial epigenomic states, as well as challenges and opportunities for their refinement. PMID:26555135

  6. Postfragmentation density function for bacterial aggregates in laminar flow

    NASA Astrophysics Data System (ADS)

    Byrne, Erin; Bortz, David M.; Dzul, Steve; Solomon, Michael; Younger, John

    2011-04-01

    The postfragmentation probability density of daughter flocs is one of the least well-understood aspects of modeling flocculation. We use three-dimensional positional data of Klebsiella pneumoniae bacterial flocs in suspension and the knowledge of hydrodynamic properties of a laminar flow field to construct a probability density function of floc volumes after a fragmentation event. We provide computational results which predict that the primary fragmentation mechanism for large flocs is erosion. The postfragmentation probability density function has a strong dependence on the size of the original floc and indicates that most fragmentation events result in clumps of one to three bacteria eroding from the original floc. We also provide numerical evidence that exhaustive fragmentation yields a limiting density inconsistent with the log-normal density predicted in the literature, most likely due to the heterogeneous nature of K. pneumoniae flocs. To support our conclusions, artificial flocs were generated and display similar postfragmentation density and exhaustive fragmentation.

  7. Integration, Continuity and a Connection with Probability Density Functions

    ERIC Educational Resources Information Center

    Samuels, M.

    2006-01-01

    This note considers functions of two variables which are continuous on a possibly unbounded closed region in [vertical bar]R[squared], and the functions of one variable obtained by integrating out the other variable over this region. The question of continuity of these functions is investigated, as are connections with joint density and marginal…

  8. Testing for parity violation in nuclei using spin density matrices for nuclear density functionals

    NASA Astrophysics Data System (ADS)

    Barrett, B. R.; Giraud, B. G.

    2015-06-01

    The spin density matrix (SDM) used in atomic and molecular physics is revisited for nuclear physics, in the context of the radial density functional theory. The vector part of the SDM defines a ‘hedgehog’ situation, which exists only if nuclear states contain some amount of parity violation. A toy model is given as an illustrative example.

  9. Fluctuation-dissipation theorem density-functional theory

    NASA Astrophysics Data System (ADS)

    Furche, Filipp; Van Voorhis, Troy

    2005-04-01

    Using the fluctuation-dissipation theorem (FDT) in the context of density-functional theory (DFT), one can derive an exact expression for the ground-state correlation energy in terms of the frequency-dependent density response function. When combined with time-dependent density-functional theory, a new class of density functionals results that use approximations to the exchange-correlation kernel fxc as input. This FDT-DFT scheme holds promise to solve two of the most distressing problems of conventional Kohn-Sham DFT: (i) It leads to correlation energy functionals compatible with exact exchange, and (ii) it naturally includes dispersion. The price is a moderately expensive O(N6) scaling of computational cost and a slower basis set convergence. These general features of FDT-DFT have all been recognized previously. In this paper, we present the first benchmark results for a set of molecules using FDT-DFT beyond the random-phase approximation (RPA)—that is, the first such results with fxc≠0. We show that kernels derived from the adiabatic local-density approximation and other semilocal functionals suffer from an "ultraviolet catastrophe," producing a pair density that diverges at small interparticle distance. Nevertheless, dispersion interactions can be treated accurately if hybrid functionals are employed, as is demonstrated for He2 and HeNe. We outline constraints that future approximations to fxc should satisfy and discuss the prospects of FDT-DFT.

  10. A 3-D microfluidic combinatorial cell array.

    PubMed

    Liu, Mike C; Tai, Yu-Chong

    2011-02-01

    We present the development of a three-dimensional (3-D) combinatorial cell culture array device featured with integrated three-input, eight-output combinatorial mixer and cell culture chambers. The device is designed for cell-based screening of multiple compounds simultaneously on a microfluidic platform. The final assembled device is composed of a porous membrane integrated in between a Parylene 3-D microfluidic chip and a PDMS microfluidic chip. The membrane turned the cell culture chambers into two-level configuration to facilitate cell loading and to maintain cells in a diffusion dominated space during device operation. Experimentally, we first characterized the combined compound concentration profile at each chamber using a fluorescence method. We then successfully demonstrated the functionality of the quantitative cell-based assay by culturing B35 rat neuronal cells on this device and screening the ability of three compounds (1,5-dihydroxyisoquinoline, deferoxamine, and 3-aminobenzoic acid) to attenuate cell death caused by cytotoxic hydrogen peroxide. In another experiment, we assayed for the combinatorial effects of three chemotherapeutic compound exposures (vinorelbine, paclitaxel, and γ-linolenic acid) on human breast cancer cells, MDA-MB-231. The same technology will enable the construction of inexpensive lab-on-a-chip devices with high-input combinatorial mixer for performing high-throughput cell-based assay and highly parallel and combinatorial chemical or biochemical reactions. PMID:21063783

  11. Charge and spin fluctuations in the density functional theory

    SciTech Connect

    Gyoerffy, B.L.; Barbieri, A. . H.H. Wills Physics Lab.); Staunton, J.B. . Dept. of Physics); Shelton, W.A.; Stocks, G.M. )

    1990-01-01

    We introduce a conceptual framework which allow us to treat charge and spin fluctuations about the Local density Approximation (LDA) to the Density Functional Theory (DFT). We illustrate the approach by explicit study of the Disordered Local Moment (DLM) state in Fe above the Curie Temperature {Tc} and the Mott insulating state in MnO. 27 refs., 6 figs.

  12. Density Functional Model for Nondynamic and Strong Correlation.

    PubMed

    Kong, Jing; Proynov, Emil

    2016-01-12

    A single-term density functional model for the left-right nondynamic/strong electron correlation is presented based on single-determinant Kohn-Sham density functional theory. It is derived from modeling the adiabatic connection for kinetic correlation energy based on physical arguments, with the correlation potential energy based on the Becke'13 model ( Becke, A.D. J. Chem. Phys . 2013 , 138 , 074109 ). This functional satisfies some known scaling relationships for correlation functionals. The fractional spin error is further reduced substantially with a new density-functional correction. Preliminary tests with self-consistent-field implementation show that the model, with only three empirical parameters, recovers the majority of left-right nondynamic/strong correlation upon bond dissociation and performs reasonably well for atomization energies and singlet-triplet energy splittings. This study also demonstrates the feasibility of developing DFT functionals for nondynamic and strong correlation within the single-determinant KS scheme. PMID:26636190

  13. Density and Spin Response Functions in Ultracold Fermionic Atom Gases

    SciTech Connect

    Mihaila, Bogdan; Blagoev, Krastan B.; Balatsky, Alexander V.; Smith, Darryl L.; Gaudio, Sergio; Littlewood, Peter B.

    2005-08-26

    We propose a new method of detecting the onset of superfluidity in a two-component ultracold fermionic gas of atoms governed by an attractive short-range interaction. By studying the two-body correlation functions we find that a measurement of the momentum distribution of the density and spin-response functions allows one to access separately the normal and anomalous densities. The change in sign at low momentum transfer of the normal-ordered part of the density response function signals the transition between a BEC and a BCS regime, characterized by small and large pairs, respectively. This change in sign of the density response function represents an unambiguous signature of the BEC-to-BCS crossover. Spin rotational symmetry breaking due to the magnetic field, if observed, can be used to validate the one-channel model.

  14. Manipulating Combinatorial Structures.

    ERIC Educational Resources Information Center

    Labelle, Gilbert

    This set of transparencies shows how the manipulation of combinatorial structures in the context of modern combinatorics can easily lead to interesting teaching and learning activities at every level of education from elementary school to university. The transparencies describe: (1) the importance and relations of combinatorics to science and…

  15. Reflection-Asymmetric Nuclear Deformations within the Density Functional Theory

    SciTech Connect

    Olsen, E; Erler, J; Nazarewicz, W.; Stoitsov, M

    2012-01-01

    Within the nuclear density functional theory (DFT) we study the effect of reflection- asymmetric shapes on ground-state binding energies and binding energy differences. To this end, we developed the new DFT solver axialhfb that uses an approximate second-order gradient to solve the Hartree-Fock-Bogoliubov equations of superconducting DFT with the quasi-local Skyrme energy density functionals. Illustrative calculations are carried out for even- even isotopes of radium and thorium.

  16. Gutzwiller density functional theory for correlated electron systems

    SciTech Connect

    Ho, K. M.; Schmalian, J.; Wang, C. Z.

    2008-02-04

    We develop a density functional theory (DFT) and formalism for correlated electron systems by taking as reference an interacting electron system that has a ground state wave function which exactly obeys the Gutzwiller approximation for all one-particle operators. The solution of the many-electron problem is mapped onto the self-consistent solution of a set of single-particle Schroedinger equations, analogously to standard DFT-local density approximation calculations.

  17. Preface: Special Topic on Advances in Density Functional Theory

    SciTech Connect

    Yang, Weitao

    2014-05-14

    This Special Topic Issue on the Advances in Density Functional Theory, published as a celebration of the fifty years of density functional theory, contains a retrospective article, a perspective article, and a collection of original research articles that showcase recent theoretical advances in the field. It provides a timely discussion reflecting a cross section of our understanding, and the theoretical and computational developments, which have significant implications in broad areas of sciences and engineering.

  18. Preface: Special topic on advances in density functional theory.

    PubMed

    Yang, Weitao

    2014-05-14

    This Special Topic Issue on the Advances in Density Functional Theory, published as a celebration of the fifty years of density functional theory, contains a retrospective article, a perspective article, and a collection of original research articles that showcase recent theoretical advances in the field. It provides a timely discussion reflecting a cross section of our understanding, and the theoretical and computational developments, which have significant implications in broad areas of sciences and engineering. PMID:24832306

  19. Density Functional Theory with Dissipation: Transport through Single Molecules

    SciTech Connect

    Kieron Burke

    2012-04-30

    A huge amount of fundamental research was performed on this grant. Most of it focussed on fundamental issues of electronic structure calculations of transport through single molecules, using density functional theory. Achievements were: (1) First density functional theory with dissipation; (2) Pseudopotential plane wave calculations with master equation; (3) Weak bias limit; (4) Long-chain conductance; and (5) Self-interaction effects in tunneling.

  20. Energy Densities in the Strong-Interaction Limit of Density Functional Theory.

    PubMed

    Mirtschink, André; Seidl, Michael; Gori-Giorgi, Paola

    2012-09-11

    We discuss energy densities in the strong-interaction limit of density functional theory, deriving an exact expression within the definition (gauge) of the electrostatic potential of the exchange-correlation hole. Exact results for small atoms and small model quantum dots (Hooke's atoms) are compared with available approximations defined in the same gauge. The idea of a local interpolation along the adiabatic connection is discussed, comparing the energy densities of the Kohn-Sham, the physical, and the strong-interacting systems. We also use our results to analyze the local version of the Lieb-Oxford bound, widely used in the construction of approximate exchange-correlation functionals. PMID:26605721

  1. Comments on the locality in density-functional theory

    SciTech Connect

    Lindgren, Ingvar; Salomonson, Sten

    2003-05-01

    The 'locality hypothesis' in density-functional theory (DFT), implying that the functional derivative is equivalent to a multiplicative local function, forms the basis of models of Kohn-Sham type. This has been generally accepted by the community since the advent of the model, and has later been formally proved for a large class of functionals. The hypothesis has recently been questioned by Nesbet [Phys. Rev. A 58, R12 (1998) and Phys. Rev. A 65, 010502 (2001)], who claims that it fails for the kinetic-energy functional for a system with more than two noninteracting electrons with a nondegenerate ground state. This conclusion has been questioned by Gal [Phys. Rev. A 62, 044501 (2000)] and by Holas and March [Phys. Rev. A 64, 016501 (2001)]. We claim that the arguments of Nesbet are incorrect, since the orbital functional used for the kinetic energy is not a unique functional of the total density in the domain of unnormalized orbitals. We have demonstrated that with a proper definition of the kinetic energy, which is a unique density functional also in the unnormalized region, the derivative can be represented by a single local multiplicative function for all v-representable densities. Therefore, we consider the controversy connected with the issue raised by Nesbet as resolved. We believe that the proof of the differentiability given here can be extended to larger groups of DFT functionals, and works along these lines are in progress.

  2. A defined synthetic substrate for serum-free culture of human stem cell derived cardiomyocytes with improved functional maturity identified using combinatorial materials microarrays

    PubMed Central

    Patel, Asha K.; Celiz, Adam D.; Rajamohan, Divya; Anderson, Daniel G.; Langer, Robert; Davies, Martyn C.

    2016-01-01

    Cardiomyocytes from human stem cells have applications in regenerative medicine and can provide models for heart disease and toxicity screening. Soluble components of the culture system such as growth factors within serum and insoluble components such as the substrate on which cells adhere to are important variables controlling the biological activity of cells. Using a combinatorial materials approach we develop a synthetic, chemically defined cellular niche for the support of functional cardiomyocytes derived from human embryonic stem cells (hESC-CMs) in a serum-free fully defined culture system. Almost 700 polymers were synthesized and evaluated for their utility as growth substrates. From this group, 20 polymers were identified that supported cardiomyocyte adhesion and spreading. The most promising 3 polymers were scaled up for extended culture of hESC-CMs for 15 days and were characterized using patch clamp electrophysiology and myofibril analysis to find that functional and structural phenotype was maintained on these synthetic substrates without the need for coating with extracellular matrix protein. In addition, we found that hESC-CMs cultured on a co-polymer of isobornyl methacrylate and tert-butylamino-ethyl methacrylate exhibited significantly longer sarcomeres relative to gelatin control. The potential utility of increased structural integrity was demonstrated in an in vitro toxicity assay that found an increase in detection sensitivity of myofibril disruption by the anti-cancer drug doxorubicin at a concentration of 0.05 µM in cardiomyocytes cultured on the co-polymer compared to 0.5 µM on gelatin. The chemical moieties identified in this large-scale screen provide chemically defined conditions for the culture and manipulation of hESC-CMs, as well as a framework for the rational design of superior biomaterials. PMID:26005764

  3. High-Density IgE Recognition of the Major Grass Pollen Allergen Phl p 1 Revealed with Single-Chain IgE Antibody Fragments Obtained by Combinatorial Cloning

    PubMed Central

    Madritsch, Christoph; Gadermaier, Elisabeth; Roder, Uwe W.; Lupinek, Christian; Valenta, Rudolf; Flicker, Sabine

    2015-01-01

    The timothy grass pollen allergen Phl p 1 belongs to the group 1 of highly cross-reactive grass pollen allergens with a molecular mass of ~25–30 kDa. Group 1 allergens are recognized by >95% of grass pollen allergic patients. We investigated the IgE recognition of Phl p 1 using allergen-specific IgE-derived single-chain variable Ab fragments (IgE-ScFvs) isolated from a combinatorial library constructed from PBMCs of a grass pollen–allergic patient. IgE-ScFvs reacted with recombinant Phl p 1 and natural group 1 grass pollen allergens. Using synthetic Phl p 1–derived peptides, the binding sites of two ScFvs were mapped to the N terminus of the allergen. In surface plasmon resonance experiments they showed comparable high-affinity binding to Phl p 1 as a complete human IgE-derived Ab recognizing the allergens’ C terminus. In a set of surface plasmon resonance experiments simultaneous allergen recognition of all three binders was demonstrated. Even in the presence of the three binders, allergic patients’ polyclonal IgE reacted with Phl p 1, indicating high-density IgE recognition of the Phl p 1 allergen. Our results show that multiple IgE Abs can bind with high density to Phl p 1, which may explain the high allergenic activity and sensitizing capacity of this allergen. PMID:25637023

  4. High-density IgE recognition of the major grass pollen allergen Phl p 1 revealed with single-chain IgE antibody fragments obtained by combinatorial cloning.

    PubMed

    Madritsch, Christoph; Gadermaier, Elisabeth; Roder, Uwe W; Lupinek, Christian; Valenta, Rudolf; Flicker, Sabine

    2015-03-01

    The timothy grass pollen allergen Phl p 1 belongs to the group 1 of highly cross-reactive grass pollen allergens with a molecular mass of ∼25-30 kDa. Group 1 allergens are recognized by >95% of grass pollen allergic patients. We investigated the IgE recognition of Phl p 1 using allergen-specific IgE-derived single-chain variable Ab fragments (IgE-ScFvs) isolated from a combinatorial library constructed from PBMCs of a grass pollen-allergic patient. IgE-ScFvs reacted with recombinant Phl p 1 and natural group 1 grass pollen allergens. Using synthetic Phl p 1-derived peptides, the binding sites of two ScFvs were mapped to the N terminus of the allergen. In surface plasmon resonance experiments they showed comparable high-affinity binding to Phl p 1 as a complete human IgE-derived Ab recognizing the allergens' C terminus. In a set of surface plasmon resonance experiments simultaneous allergen recognition of all three binders was demonstrated. Even in the presence of the three binders, allergic patients' polyclonal IgE reacted with Phl p 1, indicating high-density IgE recognition of the Phl p 1 allergen. Our results show that multiple IgE Abs can bind with high density to Phl p 1, which may explain the high allergenic activity and sensitizing capacity of this allergen. PMID:25637023

  5. Density functional theory for polymeric systems in 2D

    NASA Astrophysics Data System (ADS)

    Słyk, Edyta; Roth, Roland; Bryk, Paweł

    2016-06-01

    We propose density functional theory for polymeric fluids in two dimensions. The approach is based on Wertheim’s first order thermodynamic perturbation theory (TPT) and closely follows density functional theory for polymers proposed by Yu and Wu (2002 J. Chem. Phys. 117 2368). As a simple application we evaluate the density profiles of tangent hard-disk polymers at hard walls. The theoretical predictions are compared against the results of the Monte Carlo simulations. We find that for short chain lengths the theoretical density profiles are in an excellent agreement with the Monte Carlo data. The agreement is less satisfactory for longer chains. The performance of the theory can be improved by recasting the approach using the self-consistent field theory formalism. When the self-avoiding chain statistics is used, the theory yields a marked improvement in the low density limit. Further improvements for long chains could be reached by going beyond the first order of TPT.

  6. Density functional theory for polymeric systems in 2D.

    PubMed

    Słyk, Edyta; Roth, Roland; Bryk, Paweł

    2016-06-22

    We propose density functional theory for polymeric fluids in two dimensions. The approach is based on Wertheim's first order thermodynamic perturbation theory (TPT) and closely follows density functional theory for polymers proposed by Yu and Wu (2002 J. Chem. Phys. 117 2368). As a simple application we evaluate the density profiles of tangent hard-disk polymers at hard walls. The theoretical predictions are compared against the results of the Monte Carlo simulations. We find that for short chain lengths the theoretical density profiles are in an excellent agreement with the Monte Carlo data. The agreement is less satisfactory for longer chains. The performance of the theory can be improved by recasting the approach using the self-consistent field theory formalism. When the self-avoiding chain statistics is used, the theory yields a marked improvement in the low density limit. Further improvements for long chains could be reached by going beyond the first order of TPT. PMID:27115343

  7. Research on universal combinatorial coding.

    PubMed

    Lu, Jun; Zhang, Zhuo; Mo, Juan

    2014-01-01

    The conception of universal combinatorial coding is proposed. Relations exist more or less in many coding methods. It means that a kind of universal coding method is objectively existent. It can be a bridge connecting many coding methods. Universal combinatorial coding is lossless and it is based on the combinatorics theory. The combinational and exhaustive property make it closely related with the existing code methods. Universal combinatorial coding does not depend on the probability statistic characteristic of information source, and it has the characteristics across three coding branches. It has analyzed the relationship between the universal combinatorial coding and the variety of coding method and has researched many applications technologies of this coding method. In addition, the efficiency of universal combinatorial coding is analyzed theoretically. The multicharacteristic and multiapplication of universal combinatorial coding are unique in the existing coding methods. Universal combinatorial coding has theoretical research and practical application value. PMID:24772019

  8. Research on Universal Combinatorial Coding

    PubMed Central

    Lu, Jun; Zhang, Zhuo; Mo, Juan

    2014-01-01

    The conception of universal combinatorial coding is proposed. Relations exist more or less in many coding methods. It means that a kind of universal coding method is objectively existent. It can be a bridge connecting many coding methods. Universal combinatorial coding is lossless and it is based on the combinatorics theory. The combinational and exhaustive property make it closely related with the existing code methods. Universal combinatorial coding does not depend on the probability statistic characteristic of information source, and it has the characteristics across three coding branches. It has analyzed the relationship between the universal combinatorial coding and the variety of coding method and has researched many applications technologies of this coding method. In addition, the efficiency of universal combinatorial coding is analyzed theoretically. The multicharacteristic and multiapplication of universal combinatorial coding are unique in the existing coding methods. Universal combinatorial coding has theoretical research and practical application value. PMID:24772019

  9. Entanglement from density measurements: Analytical density functional for the entanglement of strongly correlated fermions

    SciTech Connect

    Franca, Vivian V.; D'Amico, Irene

    2011-04-15

    We derive an analytical density functional for the single-site entanglement of the one-dimensional homogeneous Hubbard model by means of an approximation to the linear entropy. We show that this very simple density functional reproduces quantitatively the exact results. We then use this functional as input for a local-density approximation to the single-site entanglement of inhomogeneous systems. We illustrate the power of this approach in a harmonically confined system, which could simulate recent experiments with ultracold atoms in optical lattices as well as in a superlattice and in an impurity system. The impressive quantitative agreement with numerical calculations--which includes reproducing subtle signatures of the particle density stages--shows that our density functional can provide entanglement calculations for actual experiments via density measurements. Next we use our functional to calculate the entanglement in disordered systems. We find that, in contrast with the expectation that disorder destroys the entanglement, there exist regimes for which the entanglement remains almost unaffected by the presence of disordered impurities.

  10. Reactivity of Graphene Investigated by Density-Functional Theory

    NASA Astrophysics Data System (ADS)

    Soni, Himadri; Gebhardt, Julian; Görling, Andreas; Chair of Theoretical Chemistry Team

    Using spin-polarized density-functional theory, we study the adsorption and reaction of hydrogen and fluorine with graphene. Graphene has a bipartite lattice with two different sublattices and hence, due to Lieb's theorem, the inequality between two sublattices should lead to a net magnetic moment upon adsorption of hydrogen or fluorine. Our calculations using density-functional theory with the generalized gradient approximation predict a magnetic moment of 1 µB for a single hydrogen adsorbed on graphene but not for a single fluorine atom adsorbed on graphene. Switching to hybrid density-functional theory with the HSE functional, we obtain a magnetic moment of 1 µB for of a single fluorine atom adsorption on graphene. This is in line with work of Kim et al., who also found in density-functional theory calculations with the HSE exchange-correlation functional spin-polarization for a fluorine adatom on graphene. Here, we present a systematic study of the reactivity and relevant adsorption mechanism for single-sided graphene, i.e., a graphene sheet which is accessible by an adsorbate from only one side with hydrogen and fluorine using hybrid density-functional theory. German Research Council (DFG) by the Collaborative Research Center 953.

  11. Range Separation and Local Hybridization in Density Functional Theory†

    PubMed Central

    Henderson, Thomas M.; Janesko, Benjamin G.; Scuseria, Gustavo E.

    2016-01-01

    Kohn–Sham density functional theory has become a standard method for modeling energetic, spectroscopic, and chemical reactivity properties of large molecules and solids. Density functional theory provides a rigorous theoretical framework for modeling the many-body exchange-correlation effects that dominate the computational cost of traditional wave function approaches. The advent of hybrid exchange-correlation functionals which incorporate a fraction of nonlocal exact exchange has solidified the prominence of density functional theory within computational chemistry. Hybrids provide accurate treatments of properties such as thermochemistry and molecular geometry. But they also exhibit some rather spectacular failures, and often contain multiple empirical parameters. This article reviews our work on developing novel exchange-correlation functionals that build upon the successes of global hybrids. We focus on more flexible functional forms, including local and range-separated hybrid functionals, constructed to obey known exact constraints and (ideally) to incorporate a minimum of empirical parametrization. The article places our work within the context of some other new approximate density functionals and discusses prospects for future work. PMID:19006280

  12. Optimal-transport formulation of electronic density-functional theory

    NASA Astrophysics Data System (ADS)

    Buttazzo, Giuseppe; De Pascale, Luigi; Gori-Giorgi, Paola

    2012-06-01

    The most challenging scenario for Kohn-Sham density-functional theory, that is, when the electrons move relatively slowly trying to avoid each other as much as possible because of their repulsion (strong-interaction limit), is reformulated here as an optimal transport (or mass transportation theory) problem, a well-established field of mathematics and economics. In practice, we show that to solve the problem of finding the minimum possible internal repulsion energy for N electrons in a given density ρ(r) is equivalent to find the optimal way of transporting N-1 times the density ρ into itself, with the cost function given by the Coulomb repulsion. We use this link to set the strong-interaction limit of density-functional theory on firm ground and to discuss the potential practical aspects of this reformulation.

  13. Optimization of an exchange-correlation density functional for water.

    PubMed

    Fritz, Michelle; Fernández-Serra, Marivi; Soler, José M

    2016-06-14

    We describe a method, that we call data projection onto parameter space (DPPS), to optimize an energy functional of the electron density, so that it reproduces a dataset of experimental magnitudes. Our scheme, based on Bayes theorem, constrains the optimized functional not to depart unphysically from existing ab initio functionals. The resulting functional maximizes the probability of being the "correct" parameterization of a given functional form, in the sense of Bayes theory. The application of DPPS to water sheds new light on why density functional theory has performed rather poorly for liquid water, on what improvements are needed, and on the intrinsic limitations of the generalized gradient approximation to electron exchange and correlation. Finally, we present tests of our water-optimized functional, that we call vdW-DF-w, showing that it performs very well for a variety of condensed water systems. PMID:27305990

  14. Optimization of an exchange-correlation density functional for water

    NASA Astrophysics Data System (ADS)

    Fritz, Michelle; Fernández-Serra, Marivi; Soler, José M.

    2016-06-01

    We describe a method, that we call data projection onto parameter space (DPPS), to optimize an energy functional of the electron density, so that it reproduces a dataset of experimental magnitudes. Our scheme, based on Bayes theorem, constrains the optimized functional not to depart unphysically from existing ab initio functionals. The resulting functional maximizes the probability of being the "correct" parameterization of a given functional form, in the sense of Bayes theory. The application of DPPS to water sheds new light on why density functional theory has performed rather poorly for liquid water, on what improvements are needed, and on the intrinsic limitations of the generalized gradient approximation to electron exchange and correlation. Finally, we present tests of our water-optimized functional, that we call vdW-DF-w, showing that it performs very well for a variety of condensed water systems.

  15. Density functional for ternary non-additive hard sphere mixtures.

    PubMed

    Schmidt, Matthias

    2011-10-19

    Based on fundamental measure theory, a Helmholtz free energy density functional for three-component mixtures of hard spheres with general, non-additive interaction distances is constructed. The functional constitutes a generalization of the previously given theory for binary non-additive mixtures. The diagrammatic structure of the spatial integrals in both functionals is of star-like (or tree-like) topology. The ternary diagrams possess a higher degree of complexity than the binary diagrams. Results for partial pair correlation functions, obtained via the Ornstein-Zernike route from the second functional derivatives of the excess free energy functional, agree well with Monte Carlo simulation data. PMID:21946780

  16. Osteoblast response to dimethacrylate composites varying in composition, conversion and roughness using a combinatorial approach

    PubMed Central

    Lin-Gibson, Sheng

    2009-01-01

    Dimethacrylate polymers and composites are seeing increased usage in orthopedics. As these applications require the material to integrate with the surrounding tissues, direct contact cytotoxicity assays should be used to assess the biocompatibility. This study utilized a combinatorial testing platform to evaluate the cell response to dimethacrylate composites with a variety of properties on a single sample. MC3T3-E1 pre-osteoblasts were cultured directly on composites with varying filler content, filler type, degree of conversion (DC), and surface topography. Cell viability, density, and area depended on an interplay of the material properties, with low DC causing a reduction in cell area but having minimal effect on cell viability, high filler content causing an increase in cell density, and filler content/type altering the surface roughness as a function of DC. The combinatorial testing platform successfully quantified the effects of numerous material properties on several aspects of the osteoblast response. PMID:19520423

  17. Growth, differentiation capacity, and function of mesenchymal stem cells expanded in serum-free medium developed via combinatorial screening.

    PubMed

    Crapnell, Kirsten; Blaesius, Rainer; Hastings, Abel; Lennon, Donald P; Caplan, Arnold I; Bruder, Scott P

    2013-06-10

    The presence of serum in cell culture medium presents an obstacle to safe and efficient production of hMSCs for therapeutic purposes. Availability of defined medium will be crucial to elucidating the mechanism of action of hMSCs in many indications as well as a prerequisite to consistently produce cells with predictable performance characteristics. Using a bioinformatics driven approach, which we call the BD Discovery Platform, we have developed a novel serum-free medium that supports highly efficient growth while maintaining the surface markers and functional characteristics defining hMSCs. In a comparison with serum-containing and other commercially available serum-free formulations, all conditions led to expansion of cells that meet the minimal criteria for hMSCs as set by the International Society for Cellular Therapy (ISCT). However, differences in growth characteristics and gene expression patterns suggest that expansion in serum-free growth conditions can provide greater yields in a shorter time. The mRNA expression profile observed in cells grown without serum suggests upregulation of several genes implicated in hMSC function as well as downregulation of the proinflammatory cytokine IL6. PMID:23597555

  18. Systematic control of experimental inconsistency in combinatorial materials science.

    PubMed

    Sharma, Asish Kumar; Kulshreshtha, Chandramouli; Sohn, Keemin; Sohn, Kee-Sun

    2009-01-01

    We developed a method to systematically control experimental inconsistency, which is one of the most troublesome and difficult problems in high-throughput combinatorial experiments. The topic of experimental inconsistency is never addressed, even though all scientists in the field of combinatorial materials science face this very serious problem. Experimental inconsistency and material property were selected as dual objective functions that were simultaneously optimized. Specifically, in an attempt to search for promising phosphors with high reproducibility, photoluminescence (PL) intensity was maximized, and experimental inconsistency was minimized by employing a multiobjective evolutionary optimization-assisted combinatorial materials search (MOEO combinatorial material search) strategy. A tetravalent manganese-doped alkali earth germanium/titanium oxide system was used as a model system to be screened using MOEO combinatorial materials search. As a result of MOEO reiteration, we identified a halide-detached deep red phosphor with improved PL intensity and reliable reproducibility. PMID:19061418

  19. Modal density function and number of propagating modes in ducts

    NASA Technical Reports Server (NTRS)

    Rice, E. J.

    1976-01-01

    Often raised questions in duct sound propagation studies involve the total number of propagating modes, the number of propagating radial modes for a particular spinning lobe number, and the number of modes possible between two given values of cutoff ratio or eigenvalue. These questions can be answered approximately by using the modal distribution function which is the integral of the modal density function for ducts in a manner similar to that previously published for architectural acoustics. The modal density functions are derived for rectangular and circular ducts with a uniform steady flow. Results from this continuous theory are compared to the actual (discrete) modal distributions.

  20. Van der Waals density functional applied to adsorption systems

    NASA Astrophysics Data System (ADS)

    Hamada, Ikutaro

    2013-03-01

    The van der Waals density functional (vdW-DF) is a promising density functional to describe the van der Waals forces within density functional theory. However, despite the recent efforts, there is still room for further improvement, especially for describing molecular adsorption on metal surfaces. I will show that by choosing appropriate exchange and nonlocal correlation functionals, it is possible to calculate geometries and electronic structures for adsorption systems accurately within the framework of vdW-DF. Applicability of the present approach will be illustrated with its applications to graphene/metal, fullerene/metal, and water/graphene interfaces. This work is partly supported by a Grant-in-Aid for Scientific Research on Innovative Area (No. 23104501). AIMR was established by the World Premier International Research Center Initiative (WPI), MEXT, Japan.

  1. Force Density Function Relationships in 2-D Granular Media

    NASA Technical Reports Server (NTRS)

    Youngquist, Robert C.; Metzger, Philip T.; Kilts, Kelly N.

    2004-01-01

    An integral transform relationship is developed to convert between two important probability density functions (distributions) used in the study of contact forces in granular physics. Developing this transform has now made it possible to compare and relate various theoretical approaches with one another and with the experimental data despite the fact that one may predict the Cartesian probability density and another the force magnitude probability density. Also, the transforms identify which functional forms are relevant to describe the probability density observed in nature, and so the modified Bessel function of the second kind has been identified as the relevant form for the Cartesian probability density corresponding to exponential forms in the force magnitude distribution. Furthermore, it is shown that this transform pair supplies a sufficient mathematical framework to describe the evolution of the force magnitude distribution under shearing. Apart from the choice of several coefficients, whose evolution of values must be explained in the physics, this framework successfully reproduces the features of the distribution that are taken to be an indicator of jamming and unjamming in a granular packing. Key words. Granular Physics, Probability Density Functions, Fourier Transforms

  2. Improved association in a classical density functional theory for water.

    PubMed

    Krebs, Eric J; Schulte, Jeff B; Roundy, David

    2014-03-28

    We present a modification to our recently published statistical associating fluid theory-based classical density functional theory for water. We have recently developed and tested a functional for the averaged radial distribution function at contact of the hard-sphere fluid that is dramatically more accurate at interfaces than earlier approximations. We now incorporate this improved functional into the association term of our free energy functional for water, improving its description of hydrogen bonding. We examine the effect of this improvement by studying two hard solutes (a hard hydrophobic rod and a hard sphere) and a Lennard-Jones approximation of a krypton atom solute. The improved functional leads to a moderate change in the density profile and a large decrease in the number of hydrogen bonds broken in the vicinity of the hard solutes. We find an improvement of the partial radial distribution for a krypton atom in water when compared with experiment. PMID:24697459

  3. Improved association in a classical density functional theory for water

    SciTech Connect

    Krebs, Eric J.; Schulte, Jeff B.; Roundy, David

    2014-03-28

    We present a modification to our recently published statistical associating fluid theory-based classical density functional theory for water. We have recently developed and tested a functional for the averaged radial distribution function at contact of the hard-sphere fluid that is dramatically more accurate at interfaces than earlier approximations. We now incorporate this improved functional into the association term of our free energy functional for water, improving its description of hydrogen bonding. We examine the effect of this improvement by studying two hard solutes (a hard hydrophobic rod and a hard sphere) and a Lennard-Jones approximation of a krypton atom solute. The improved functional leads to a moderate change in the density profile and a large decrease in the number of hydrogen bonds broken in the vicinity of the hard solutes. We find an improvement of the partial radial distribution for a krypton atom in water when compared with experiment.

  4. Optimization of Combinatorial Mutagenesis

    NASA Astrophysics Data System (ADS)

    Parker, Andrew S.; Griswold, Karl E.; Bailey-Kellogg, Chris

    Protein engineering by combinatorial site-directed mutagenesis evaluates a portion of the sequence space near a target protein, seeking variants with improved properties (stability, activity, immunogenicity, etc.). In order to improve the hit-rate of beneficial variants in such mutagenesis libraries, we develop methods to select optimal positions and corresponding sets of the mutations that will be used, in all combinations, in constructing a library for experimental evaluation. Our approach, OCoM (Optimization of Combinatorial Mutagenesis), encompasses both degenerate oligonucleotides and specified point mutations, and can be directed accordingly by requirements of experimental cost and library size. It evaluates the quality of the resulting library by one- and two-body sequence potentials, averaged over the variants. To ensure that it is not simply recapitulating extant sequences, it balances the quality of a library with an explicit evaluation of the novelty of its members. We show that, despite dealing with a combinatorial set of variants, in our approach the resulting library optimization problem is actually isomorphic to single-variant optimization. By the same token, this means that the two-body sequence potential results in an NP-hard optimization problem. We present an efficient dynamic programming algorithm for the one-body case and a practically-efficient integer programming approach for the general two-body case. We demonstrate the effectiveness of our approach in designing libraries for three different case study proteins targeted by previous combinatorial libraries - a green fluorescent protein, a cytochrome P450, and a beta lactamase. We found that OCoM worked quite efficiently in practice, requiring only 1 hour even for the massive design problem of selecting 18 mutations to generate 107 variants of a 443-residue P450. We demonstrate the general ability of OCoM in enabling the protein engineer to explore and evaluate trade-offs between quality and

  5. Addressing spectroscopic quality of covariant density functional theory

    NASA Astrophysics Data System (ADS)

    Afanasjev, A. V.

    2015-03-01

    The spectroscopic quality of covariant density functional theory has been accessed by analyzing the accuracy and theoretical uncertainties in the description of spectroscopic observables. Such analysis is first presented for the energies of the single-particle states in spherical and deformed nuclei. It is also shown that the inclusion of particle-vibration coupling improves the description of the energies of predominantly single-particle states in medium and heavy-mass spherical nuclei. However, the remaining differences between theory and experiment clearly indicate missing physics and missing terms in covariant energy density functionals. The uncertainties in the predictions of the position of two-neutron drip line sensitively depend on the uncertainties in the prediction of the energies of the single-particle states. On the other hand, many spectroscopic observables in well deformed nuclei at ground state and finite spin only weakly depend on the choice of covariant energy density functional.

  6. Exploration of a modified density dependence in the Skyrme functional

    SciTech Connect

    Erler, J.; Reinhard, P.-G.; Kluepfel, P.

    2010-10-15

    A variant of the basic Skyrme-Hartree-Fock functional is considered dealing with a new form of density dependence. It employs only integer powers and thus will allow a more sound basis for projection schemes (particle number, angular momentum). We optimize the new functional with exactly the same adjustment strategy as used in an earlier study with a standard Skyrme functional. This allows direct comparisons of the performance of the new functional relative to the standard one. We discuss various observables: bulk properties of finite nuclei, nuclear matter, giant resonances, superheavy elements, and energy systematics. The new functional performs at least as well as the standard one, but offers a wider range of applicability (e.g., for projection) and more flexibility in the regime of high densities.

  7. Kinetic-energy density functional: Atoms and shell structure

    SciTech Connect

    Garcia-Gonzalez, P.; Alvarellos, J.E.; Chacon, E. |

    1996-09-01

    We present a nonlocal kinetic-energy functional which includes an anisotropic average of the density through a symmetrization procedure. This functional allows a better description of the nonlocal effects of the electron system. The main consequence of the symmetrization is the appearance of a clear shell structure in the atomic density profiles, obtained after the minimization of the total energy. Although previous results with some of the nonlocal kinetic functionals have given incipient structures for heavy atoms, only our functional shows a clear shell structure for most of the atoms. The atomic total energies have a good agreement with the exact calculations. Discussion of the chemical potential and the first ionization potential in atoms is included. The functional is also extended to spin-polarized systems. {copyright} {ital 1996 The American Physical Society.}

  8. An Accurate Density Functional from Exchange-Correlation Hole

    NASA Astrophysics Data System (ADS)

    Tao, Jianmin; Mo, Yuxiang

    The exchange-correlation hole is most fundamentally important in the development and understanding of density functional theory (DFT). However, due to the nonlocal nature of the exchange-correlation hole, development of DFT from the underlying hole presents a great challenge, and the works along this direction are limited. Here I will discuss a novel nonempirical DFT based on a semilocal hole, which is obtained from the density matrix expansion. Extensive tests on molecules and solids show that this functional can achieve remarkable accuracy for wide-ranging properties in condensed matter physics and quantum chemistry. This work was supported by NSF under Grant No. CHE-1261918.

  9. Nonlinear eigenvalue problems in Density Functional Theory calculations

    SciTech Connect

    Fattebert, J

    2009-08-28

    Developed in the 1960's by W. Kohn and coauthors, Density Functional Theory (DFT) is a very popular quantum model for First-Principles simulations in chemistry and material sciences. It allows calculations of systems made of hundreds of atoms. Indeed DFT reduces the 3N-dimensional Schroedinger electronic structure problem to the search for a ground state electronic density in 3D. In practice it leads to the search for N electronic wave functions solutions of an energy minimization problem in 3D, or equivalently the solution of an eigenvalue problem with a non-linear operator.

  10. Constrained Density Functional Calculations of alpha and delta Pu

    NASA Astrophysics Data System (ADS)

    Eriksson, Olle

    2002-03-01

    The electronic structure of α and δ Pu are described using a modified density functional theory that incorporates localization effects of the 5f shell. It is argued that a Russel-Saunders coupled state involving a 5f^4 multiplet, together with one itinerant 5f electron explains most of the observed ground state properties of δ Pu (equilibrium volume, elastic constants, near degeneracy with the alpha phase). This 5f electrons in the α phase are argued to form itinerant states, that are well described in density functional theory. The two distinctly different electronic ground states give rise to different excitation spectra and a comparison with experimental data is made.

  11. Direct propagation of probability density functions in hydrological equations

    NASA Astrophysics Data System (ADS)

    Kunstmann, Harald; Kastens, Marko

    2006-06-01

    Sustainable decisions in hydrological risk management require detailed information on the probability density function ( pdf) of the model output. Only then probabilities for the failure of a specific management option or the exceedance of critical thresholds (e.g. of pollutants) can be derived. A new approach of uncertainty propagation in hydrological equations is developed that directly propagates the probability density functions of uncertain model input parameters into the corresponding probability density functions of model output. The basics of the methodology are presented and central applications to different disciplines in hydrology are shown. This work focuses on the following basic hydrological equations: (1) pumping test analysis (Theis-equation, propagation of uncertainties in recharge and transmissivity), (2) 1-dim groundwater contaminant transport equation (Gauss-equation, propagation of uncertainties in decay constant and dispersivity), (3) evapotranspiration estimation (Penman-Monteith-equation, propagation of uncertainty in roughness length). The direct propagation of probability densities is restricted to functions that are monotonically increasing or decreasing or that can be separated in corresponding monotonic branches so that inverse functions can be derived. In case no analytic solutions for inverse functions could be derived, semi-analytical approximations were used. It is shown that the results of direct probability density function propagation are in perfect agreement with results obtained from corresponding Monte Carlo derived frequency distributions. Direct pdf propagation, however, has the advantage that is yields exact solutions for the resulting hydrological pdfs rather than approximating discontinuous frequency distributions. It is additionally shown that the type of the resulting pdf depends on the specific values (order of magnitude, respectively) of the standard deviation of the input pdf. The dependency of skewness and kurtosis

  12. The exact density functional for two electrons in one dimension

    NASA Astrophysics Data System (ADS)

    Cohen, Aron; Mori-Sanchez, Paula

    The exact universal density functional F [ ρ ] is calculated for real space two-electron densities in one dimension ρ (x) with a soft-Coulomb interaction. It is calculated by the Levy constrained search F [ ρ ] =minΨ-->ρ < Ψ | \\Tcirc +\\Vcircee | Ψ > over wavefunctions of a two-dimensional Hilbert space Ψ (x1 ,x2) --> ρ (x1) and can be directly visualized. We do an approximate constrained search via density matrices and a direct approximation to natural orbitals. This allows us to make an accurate approximation to the exact functional that is calculated using a search over potentials. We investigate the exact functional and the performance of many approximations on some of the most challenging electronic structure in two-electron systems, from strongly-correlated electron transfer to the description of a localized-delocalized transition. The exact Kohn-Sham potential, vs (x) , and exact Kohn-Sham eigenvalues, ɛi, are calculated and this allows us to discuss the band-gap problem versus the perspective of the exact density functional F [ ρ ] for all numbers of electrons. We calculate the derivative discontinuity of the exact functional in an example of a Mott-Insulator, one-dimensional stretched H2.

  13. Probability density function method for variable-density pressure-gradient-driven turbulence and mixing

    SciTech Connect

    Bakosi, Jozsef; Ristorcelli, Raymond J

    2010-01-01

    Probability density function (PDF) methods are extended to variable-density pressure-gradient-driven turbulence. We apply the new method to compute the joint PDF of density and velocity in a non-premixed binary mixture of different-density molecularly mixing fluids under gravity. The full time-evolution of the joint PDF is captured in the highly non-equilibrium flow: starting from a quiescent state, transitioning to fully developed turbulence and finally dissipated by molecular diffusion. High-Atwood-number effects (as distinguished from the Boussinesq case) are accounted for: both hydrodynamic turbulence and material mixing are treated at arbitrary density ratios, with the specific volume, mass flux and all their correlations in closed form. An extension of the generalized Langevin model, originally developed for the Lagrangian fluid particle velocity in constant-density shear-driven turbulence, is constructed for variable-density pressure-gradient-driven flows. The persistent small-scale anisotropy, a fundamentally 'non-Kolmogorovian' feature of flows under external acceleration forces, is captured by a tensorial diffusion term based on the external body force. The material mixing model for the fluid density, an active scalar, is developed based on the beta distribution. The beta-PDF is shown to be capable of capturing the mixing asymmetry and that it can accurately represent the density through transition, in fully developed turbulence and in the decay process. The joint model for hydrodynamics and active material mixing yields a time-accurate evolution of the turbulent kinetic energy and Reynolds stress anisotropy without resorting to gradient diffusion hypotheses, and represents the mixing state by the density PDF itself, eliminating the need for dubious mixing measures. Direct numerical simulations of the homogeneous Rayleigh-Taylor instability are used for model validation.

  14. On the Evolution of the Density Probability Density Function in Strongly Self-gravitating Systems

    NASA Astrophysics Data System (ADS)

    Girichidis, Philipp; Konstandin, Lukas; Whitworth, Anthony P.; Klessen, Ralf S.

    2014-02-01

    The time evolution of the probability density function (PDF) of the mass density is formulated and solved for systems in free-fall using a simple approximate function for the collapse of a sphere. We demonstrate that a pressure-free collapse results in a power-law tail on the high-density side of the PDF. The slope quickly asymptotes to the functional form PV (ρ)vpropρ-1.54 for the (volume-weighted) PDF and PM (ρ)vpropρ-0.54 for the corresponding mass-weighted distribution. From the simple approximation of the PDF we derive analytic descriptions for mass accretion, finding that dynamically quiet systems with narrow density PDFs lead to retarded star formation and low star formation rates (SFRs). Conversely, strong turbulent motions that broaden the PDF accelerate the collapse causing a bursting mode of star formation. Finally, we compare our theoretical work with observations. The measured SFRs are consistent with our model during the early phases of the collapse. Comparison of observed column density PDFs with those derived from our model suggests that observed star-forming cores are roughly in free-fall.

  15. On the evolution of the density probability density function in strongly self-gravitating systems

    SciTech Connect

    Girichidis, Philipp; Konstandin, Lukas; Klessen, Ralf S.; Whitworth, Anthony P.

    2014-02-01

    The time evolution of the probability density function (PDF) of the mass density is formulated and solved for systems in free-fall using a simple approximate function for the collapse of a sphere. We demonstrate that a pressure-free collapse results in a power-law tail on the high-density side of the PDF. The slope quickly asymptotes to the functional form P{sub V} (ρ)∝ρ{sup –1.54} for the (volume-weighted) PDF and P{sub M} (ρ)∝ρ{sup –0.54} for the corresponding mass-weighted distribution. From the simple approximation of the PDF we derive analytic descriptions for mass accretion, finding that dynamically quiet systems with narrow density PDFs lead to retarded star formation and low star formation rates (SFRs). Conversely, strong turbulent motions that broaden the PDF accelerate the collapse causing a bursting mode of star formation. Finally, we compare our theoretical work with observations. The measured SFRs are consistent with our model during the early phases of the collapse. Comparison of observed column density PDFs with those derived from our model suggests that observed star-forming cores are roughly in free-fall.

  16. Subsystem real-time time dependent density functional theory

    NASA Astrophysics Data System (ADS)

    Krishtal, Alisa; Ceresoli, Davide; Pavanello, Michele

    2015-04-01

    We present the extension of Frozen Density Embedding (FDE) formulation of subsystem Density Functional Theory (DFT) to real-time Time Dependent Density Functional Theory (rt-TDDFT). FDE is a DFT-in-DFT embedding method that allows to partition a larger Kohn-Sham system into a set of smaller, coupled Kohn-Sham systems. Additional to the computational advantage, FDE provides physical insight into the properties of embedded systems and the coupling interactions between them. The extension to rt-TDDFT is done straightforwardly by evolving the Kohn-Sham subsystems in time simultaneously, while updating the embedding potential between the systems at every time step. Two main applications are presented: the explicit excitation energy transfer in real time between subsystems is demonstrated for the case of the Na4 cluster and the effect of the embedding on optical spectra of coupled chromophores. In particular, the importance of including the full dynamic response in the embedding potential is demonstrated.

  17. External orthogonality in subsystem time-dependent density functional theory.

    PubMed

    Chulhai, Dhabih V; Jensen, Lasse

    2016-08-01

    Subsystem density functional theory (subsystem DFT) is a DFT partitioning method that is exact in principle, but depends on approximations to the kinetic energy density functional (KEDF). One may avoid the use of approximate KEDFs by ensuring that the inter-subsystem molecular orbitals are orthogonal, termed external orthogonality (EO). We present a method that extends a subsystem DFT method, that includes EO, into the time-dependent DFT (TDDFT) regime. This method therefore removes the need for approximations to the kinetic energy potential and kernel, and we show that it can accurately reproduce the supermolecular TDDFT results for weakly and strongly coupled subsystems, and for systems with strongly overlapping densities (where KEDF approximations traditionally fail). PMID:26932176

  18. Applying Density Functional Theory for Atomic Vacancies in Solids

    NASA Astrophysics Data System (ADS)

    Zhou, Xiaolan; Perdew, John P.

    2008-03-01

    We use a new generalized gradient approximation of density functional theory -- PBEsol, a revised Perdew-Burke-Ernzerhof GGA, to calculate the vacancy formation energies and other properties of metals and semiconductors: Al, Pt, Pd, GaN [1], etc. By restoring the gradient expansion over a wide range of density gradients, PBEsol [2] yields excellent jellium exchange and correlation surface energies. We expect that this new functional will improve the description of vacancies in real materials, since the vacancy formation energy is essentially the work needed to create an interior surface. [1] Thomas R.Mattsson and Ann E. Mattson. Phys. Rev. B 66, 214410 (2002). [2] John P. Perdew, Adrienn Ruzsinszky, Gabor I. Csonka, Oleg A. Vydrov, Gustavo E. Scuseria, Lucian A. Constantin, Xiaolan Zhou, and Kieron Burke, Restoring the density-gradient expansion for exchange in solids and surfaces, http://arxiv.org/abs/0711.0156

  19. Catalytic activities of platinum nanotubes: a density functional study

    NASA Astrophysics Data System (ADS)

    Mukherjee, Prajna; Gupta, Bikash C.; Jena, Puru

    2015-10-01

    In this work we investigate the catalytic properties of platinum nanotubes using density functional theory based calculations. In particular, we study the dissociation of hydrogen and oxygen molecules as well as oxidation of CO molecules. The results indicate that platinum nanotubes have good catalytic properties and can be effectively used in converting CO molecule to CO2.

  20. Comparative density functional theory and density functional tight binding study of 2-anthroic acid on TiO2

    NASA Astrophysics Data System (ADS)

    Manzhos, Sergei

    2016-01-01

    A comparative DFTB (density functional tight binding)-DFT (density functional theory) study of several adsorption modes of 2-anthroic acid on titania is presented. Two parameterizations of DFTB previously used for dye-TiO2 interfaces are tested. DFTB predicts adsorption energies which differ from those computed by DFT not only in magnitude (by up to 0.5 eV) but also in the order among different configurations. The band alignment computed with DFTB is not consistent with DFT results and with experimental data. The strategy of geometry optimization with DFTB followed by single-point DFT calculations also does not necessarily result in plausible adsorption energies.

  1. Combinatorial protein design strategies using computational methods.

    PubMed

    Kono, Hidetoshi; Wang, Wei; Saven, Jeffery G

    2007-01-01

    Computational methods continue to facilitate efforts in protein design. Most of this work has focused on searching sequence space to identify one or a few sequences compatible with a given structure and functionality. Probabilistic computational methods provide information regarding the range of amino acid variability permitted by desired functional and structural constraints. Such methods may be used to guide the construction of both individual sequences and combinatorial libraries of proteins. PMID:17041256

  2. Optimization of combinatorial mutagenesis.

    PubMed

    Parker, Andrew S; Griswold, Karl E; Bailey-Kellogg, Chris

    2011-11-01

    Protein engineering by combinatorial site-directed mutagenesis evaluates a portion of the sequence space near a target protein, seeking variants with improved properties (e.g., stability, activity, immunogenicity). In order to improve the hit-rate of beneficial variants in such mutagenesis libraries, we develop methods to select optimal positions and corresponding sets of the mutations that will be used, in all combinations, in constructing a library for experimental evaluation. Our approach, OCoM (Optimization of Combinatorial Mutagenesis), encompasses both degenerate oligonucleotides and specified point mutations, and can be directed accordingly by requirements of experimental cost and library size. It evaluates the quality of the resulting library by one- and two-body sequence potentials, averaged over the variants. To ensure that it is not simply recapitulating extant sequences, it balances the quality of a library with an explicit evaluation of the novelty of its members. We show that, despite dealing with a combinatorial set of variants, in our approach the resulting library optimization problem is actually isomorphic to single-variant optimization. By the same token, this means that the two-body sequence potential results in an NP-hard optimization problem. We present an efficient dynamic programming algorithm for the one-body case and a practically-efficient integer programming approach for the general two-body case. We demonstrate the effectiveness of our approach in designing libraries for three different case study proteins targeted by previous combinatorial libraries--a green fluorescent protein, a cytochrome P450, and a beta lactamase. We found that OCoM worked quite efficiently in practice, requiring only 1 hour even for the massive design problem of selecting 18 mutations to generate 10⁷ variants of a 443-residue P450. We demonstrate the general ability of OCoM in enabling the protein engineer to explore and evaluate trade-offs between quality and

  3. Reduced density-matrix functionals applied to the Hubbard dimer

    NASA Astrophysics Data System (ADS)

    Kamil, Ebad; Schade, Robert; Pruschke, Thomas; Blöchl, Peter E.

    2016-02-01

    Common density-matrix functionals, the Müller and the power functional, have been benchmarked for the half-filled Hubbard dimer, which allows us to model the bond dissociation problem and the transition from the weakly to the strongly correlated limit. Unbiased numerical calculations are combined with analytical results. Despite the well known successes of the Müller functional, the ground state is degenerate with a one-dimensional manifold of ferromagnetic solutions. The resulting infinite magnetic susceptibility indicates another qualitative flaw of the Müller functional. The derivative discontinuity with respect to particle number is not present indicating an incorrect metal-like behavior. The power functional actually favors the ferromagnetic state for weak interaction. Analogous to the Hartree-Fock approximation, the power functional undergoes a transition beyond a critical interaction strength, in this case, however, to a noncollinear antiferromagnetic state.

  4. Density Functional Study of Perovskite Superconductor MgCNi3

    NASA Astrophysics Data System (ADS)

    Kumar, Jagdish; Sharma, Devina; Kumar, Ranjan; Awana, V. P. S.; Ahluwalia, P. K.

    2011-12-01

    We here report the first principle density functional study of MgCNi3 which crystallize in cubic perovskite structure having critical transition temperature of 8 K. The interesting aspect of this compound is that in normal state it is non magnetic in nature despite conduction electrons in it are derived from partially filled Ni d states, which typically lead to ferromagnetism in metallic Ni and many Ni-based binary alloys. To investigate the detailed microscopic origin of the non magnetic nature we have done density functional based calculations on this compound. The lattice constant is calculated using minimum energy criteria from total energy versus lattice constant plot. By taking the calculated values of lattice constant we have done the precise calculations on the compound using Full Potential Linear Augmented Plane Wave (FP-LAPW) method implemented in ELK code. The electronic density of states is found spin degenerate that corresponds to a non-magnetic ground state. The density of states (DOS) at Fermi level, N(EF) is dominated by Ni-d states. The sharp peak observed just below Fermi level corresponds to van Hove singularity (vHs). The projected density of states (PDOS) suggests a strong hybridization of Ni-3d and C-2p states which is responsible for the observed non magnetic nature of MgCNi3.

  5. Density functional calculations of point defects in InAs

    NASA Astrophysics Data System (ADS)

    Moussa, Jonathan; Schultz, Peter

    2013-03-01

    Standard semilocal density functionals do not generate a gap in the Kohn-Sham eigenvalues for InAs, a semiconductor with an experimental gap of 0.4 eV. Without a theoretical band gap, it becomes difficult to identify, specify, and characterize pure localized states of point defects with energy levels within the experimental band gap. The bulk band gap problem can be alleviated with screened hybrid density functionals, such as the Heyd-Scuseria-Ernzerhof (HSE) functional, that open the generalized Kohn-Sham eigenvalue gap of InAs to near the experimental value. However, even without a Kohn-Sham gap, the local moment countercharge (LMCC) method [Phys. Rev. Lett. 96, 246401 (2006)] is able to predict charge transition energy levels of localized defect states, using standard semi-local functionals. We present an LMCC-based study of standard point defects in InAs using semilocal density functionals and compare the results to HSE-based calculations to assess the validity of LMCC calculations in this situation. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. DOE's National Nuclear Security Administration under contract DE-AC04-94AL85000.

  6. Density Functional Theory of Biologically Relevant Metal Centers

    NASA Astrophysics Data System (ADS)

    Siegbahn, Per E. M.; Blomberg, Margareta R. A.

    1999-10-01

    Recent applications of density functional theory to biologically relevant metal centers are reviewed. The emphasis is on reaction mechanisms, structures, and modeling. The accuracy of different functionals is discussed for standard benchmark tests of first- and second-row molecules and for transition metal systems. Modeling aspects of the protein metal complexes are discussed regarding both the size of the model being treated quantum mechanically and the treatment of the protein surrounding it. To illustrate the effects, structures computed without the effects of the protein are compared with experimental structures from enzymes, and results from simple dielectric models of the protein for electron transfer processes are described. The choice of spin state is discussed for multimetal complexes. Examples of mechanisms studied recently by density functional theory are described, such as O2 and methane activation in methane monooxygenase and O2 formation in photosystem II.

  7. Ab initio derivation of model energy density functionals

    NASA Astrophysics Data System (ADS)

    Dobaczewski, Jacek

    2016-08-01

    I propose a simple and manageable method that allows for deriving coupling constants of model energy density functionals (EDFs) directly from ab initio calculations performed for finite fermion systems. A proof-of-principle application allows for linking properties of finite nuclei, determined by using the nuclear nonlocal Gogny functional, to the coupling constants of the quasilocal Skyrme functional. The method does not rely on properties of infinite fermion systems but on the ab initio calculations in finite systems. It also allows for quantifying merits of different model EDFs in describing the ab initio results.

  8. Perspective: Fundamental aspects of time-dependent density functional theory

    NASA Astrophysics Data System (ADS)

    Maitra, Neepa T.

    2016-06-01

    In the thirty-two years since the birth of the foundational theorems, time-dependent density functional theory has had a tremendous impact on calculations of electronic spectra and dynamics in chemistry, biology, solid-state physics, and materials science. Alongside the wide-ranging applications, there has been much progress in understanding fundamental aspects of the functionals and the theory itself. This Perspective looks back to some of these developments, reports on some recent progress and current challenges for functionals, and speculates on future directions to improve the accuracy of approximations used in this relatively young theory.

  9. Reformulation of Density Functional Theory for N-Representable Densities and the Resolution of the v-Representability Problem

    SciTech Connect

    Gonis, A.; Zhang, X. G.; Stocks, G. M.; Nicholson, D. M.

    2015-10-23

    Density functional theory for the case of general, N-representable densities is reformulated in terms of density functional derivatives of expectation values of operators evaluated with wave functions leading to a density, making no reference to the concept of potential. The developments provide a complete solution of the v-representability problem by establishing a mathematical procedure that determines whether a density is v-representable and in the case of an affirmative answer determines the potential (within an additive constant) as a derivative with respect to the density of a constrained search functional. It also establishes the existence of an energy functional of the density that, for v-representable densities, assumes its minimum value at the density describing the ground state of an interacting many-particle system. The theorems of Hohenberg and Kohn emerge as special cases of the formalism.

  10. Reformulation of Density Functional Theory for N-Representable Densities and the Resolution of the v-Representability Problem

    DOE PAGESBeta

    Gonis, A.; Zhang, X. G.; Stocks, G. M.; Nicholson, D. M.

    2015-10-23

    Density functional theory for the case of general, N-representable densities is reformulated in terms of density functional derivatives of expectation values of operators evaluated with wave functions leading to a density, making no reference to the concept of potential. The developments provide a complete solution of the v-representability problem by establishing a mathematical procedure that determines whether a density is v-representable and in the case of an affirmative answer determines the potential (within an additive constant) as a derivative with respect to the density of a constrained search functional. It also establishes the existence of an energy functional of themore » density that, for v-representable densities, assumes its minimum value at the density describing the ground state of an interacting many-particle system. The theorems of Hohenberg and Kohn emerge as special cases of the formalism.« less

  11. Modal density function and number of propagating modes in ducts

    NASA Technical Reports Server (NTRS)

    Rice, E. J.

    1976-01-01

    The question of the number of propagating modes within a small range of mode cut off ratio was raised. The population density of modes were shown to be greatest near cut off and least for the well propagating modes. It was shown that modes of nearly the same cut off ratio behave nearly the same in a sound absorbing duct as well as in the way they propagate to the far. Handling all of the propagating modes individually, they can be grouped into several cut off ratio ranges. It is important to know the modal density function to estimate acoustic power distribution.

  12. Cryptographic Combinatorial Securities Exchanges

    NASA Astrophysics Data System (ADS)

    Thorpe, Christopher; Parkes, David C.

    We present a useful new mechanism that facilitates the atomic exchange of many large baskets of securities in a combinatorial exchange. Cryptography prevents information about the securities in the baskets from being exploited, enhancing trust. Our exchange offers institutions who wish to trade large positions a new alternative to existing methods of block trading: they can reduce transaction costs by taking advantage of other institutions’ available liquidity, while third party liquidity providers guarantee execution—preserving their desired portfolio composition at all times. In our exchange, institutions submit encrypted orders which are crossed, leaving a “remainder”. The exchange proves facts about the portfolio risk of this remainder to third party liquidity providers without revealing the securities in the remainder, the knowledge of which could also be exploited. The third parties learn either (depending on the setting) the portfolio risk parameters of the remainder itself, or how their own portfolio risk would change if they were to incorporate the remainder into a portfolio they submit. In one setting, these third parties submit bids on the commission, and the winner supplies necessary liquidity for the entire exchange to clear. This guaranteed clearing, coupled with external price discovery from the primary markets for the securities, sidesteps difficult combinatorial optimization problems. This latter method of proving how taking on the remainder would change risk parameters of one’s own portfolio, without revealing the remainder’s contents or its own risk parameters, is a useful protocol of independent interest.

  13. Numerical density-to-potential inversions in time-dependent density functional theory.

    PubMed

    Jensen, Daniel S; Wasserman, Adam

    2016-08-01

    We treat the density-to-potential inverse problem of time-dependent density functional theory as an optimization problem with a partial differential equation constraint. The unknown potential is recovered from a target density by applying a multilevel optimization method controlled by error estimates. We employ a classical optimization routine using gradients efficiently computed by the discrete adjoint method. The inverted potential has both a real and imaginary part to reduce reflections at the boundaries and other numerical artifacts. We demonstrate this method on model one-dimensional systems. The method can be straightforwardly extended to a variety of numerical solvers of the time-dependent Kohn-Sham equations and to systems in higher dimensions. PMID:26984427

  14. SURFACE SYMMETRY ENERGY OF NUCLEAR ENERGY DENSITY FUNCTIONALS

    SciTech Connect

    Nikolov, N; Schunck, N; Nazarewicz, W; Bender, M; Pei, J

    2010-12-20

    We study the bulk deformation properties of the Skyrme nuclear energy density functionals. Following simple arguments based on the leptodermous expansion and liquid drop model, we apply the nuclear density functional theory to assess the role of the surface symmetry energy in nuclei. To this end, we validate the commonly used functional parametrizations against the data on excitation energies of superdeformed band-heads in Hg and Pb isotopes, and fission isomers in actinide nuclei. After subtracting shell effects, the results of our self-consistent calculations are consistent with macroscopic arguments and indicate that experimental data on strongly deformed configurations in neutron-rich nuclei are essential for optimizing future nuclear energy density functionals. The resulting survey provides a useful benchmark for further theoretical improvements. Unlike in nuclei close to the stability valley, whose macroscopic deformability hangs on the balance of surface and Coulomb terms, the deformability of neutron-rich nuclei strongly depends on the surface-symmetry energy; hence, its proper determination is crucial for the stability of deformed phases of the neutron-rich matter and description of fission rates for r-process nucleosynthesis.

  15. Computational complexity of time-dependent density functional theory

    NASA Astrophysics Data System (ADS)

    Whitfield, J. D.; Yung, M.-H.; Tempel, D. G.; Boixo, S.; Aspuru-Guzik, A.

    2014-08-01

    Time-dependent density functional theory (TDDFT) is rapidly emerging as a premier method for solving dynamical many-body problems in physics and chemistry. The mathematical foundations of TDDFT are established through the formal existence of a fictitious non-interacting system (known as the Kohn-Sham system), which can reproduce the one-electron reduced probability density of the actual system. We build upon these works and show that on the interior of the domain of existence, the Kohn-Sham system can be efficiently obtained given the time-dependent density. We introduce a V-representability parameter which diverges at the boundary of the existence domain and serves to quantify the numerical difficulty of constructing the Kohn-Sham potential. For bounded values of V-representability, we present a polynomial time quantum algorithm to generate the time-dependent Kohn-Sham potential with controllable error bounds.

  16. Steady-State Density Functional Theory for Finite Bias Conductances.

    PubMed

    Stefanucci, G; Kurth, S

    2015-12-01

    In the framework of density functional theory, a formalism to describe electronic transport in the steady state is proposed which uses the density on the junction and the steady current as basic variables. We prove that, in a finite window around zero bias, there is a one-to-one map between the basic variables and both local potential on as well as bias across the junction. The resulting Kohn-Sham system features two exchange-correlation (xc) potentials, a local xc potential, and an xc contribution to the bias. For weakly coupled junctions the xc potentials exhibit steps in the density-current plane which are shown to be crucial to describe the Coulomb blockade diamonds. At small currents these steps emerge as the equilibrium xc discontinuity bifurcates. The formalism is applied to a model benzene junction, finding perfect agreement with the orthodox theory of Coulomb blockade. PMID:26571349

  17. Scrutinizing "Invisible" astatine: A challenge for modern density functionals.

    PubMed

    Sergentu, Dumitru-Claudiu; David, Grégoire; Montavon, Gilles; Maurice, Rémi; Galland, Nicolas

    2016-06-01

    The main-group 6p elements did not receive much attention in the development of recent density functionals. In many cases it is still difficult to choose among the modern ones a relevant functional for various applications. Here, we illustrate the case of astatine species (At, Z = 85) and we report the first, and quite complete, benchmark study on several properties concerning such species. Insights on geometries, transition energies and thermodynamic properties of a set of 19 astatine species, for which reference experimental or theoretical data has been reported, are obtained with relativistic (two-component) density functional theory calculations. An extensive set of widely used functionals is employed. The hybrid meta-generalized gradient approximation (meta-GGA) PW6B95 functional is overall the best choice. It is worth noting that the range-separated HSE06 functional as well as the old and very popular B3LYP and PBE0 hybrid-GGAs appear to perform quite well too. Moreover, we found that astatine chemistry in solution can accurately be predicted using implicit solvent models, provided that specific parameters are used to build At cavities. © 2016 Wiley Periodicals, Inc. PMID:27059181

  18. What Density Functional Theory could do for Quantum Information

    NASA Astrophysics Data System (ADS)

    Mattsson, Ann

    2015-03-01

    The Hohenberg-Kohn theorem of Density Functional Theory (DFT), and extensions thereof, tells us that all properties of a system of electrons can be determined through their density, which uniquely determines the many-body wave-function. Given access to the appropriate, universal, functionals of the density we would, in theory, be able to determine all observables of any electronic system, without explicit reference to the wave-function. On the other hand, the wave-function is at the core of Quantum Information (QI), with the wave-function of a set of qubits being the central computational resource in a quantum computer. While there is seemingly little overlap between DFT and QI, reliance upon observables form a key connection. Though the time-evolution of the wave-function and associated phase information is fundamental to quantum computation, the initial and final states of a quantum computer are characterized by observables of the system. While observables can be extracted directly from a system's wave-function, DFT tells us that we may be able to intuit a method for extracting them from its density. In this talk, I will review the fundamentals of DFT and how these principles connect to the world of QI. This will range from DFT's utility in the engineering of physical qubits, to the possibility of using it to efficiently (but approximately) simulate Hamiltonians at the logical level. The apparent paradox of describing algorithms based on the quantum mechanical many-body wave-function with a DFT-like theory based on observables will remain a focus throughout. The ultimate goal of this talk is to initiate a dialog about what DFT could do for QI, in theory and in practice. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

  19. Dynamical density functional theory for colloidal dispersions including hydrodynamic interactions

    NASA Astrophysics Data System (ADS)

    Rex, M.; Löwen, H.

    2009-02-01

    A dynamical density functional theory (DDFT) for translational Brownian dynamics is derived which includes hydrodynamic interactions. The theory reduces to the simple Brownian DDFT proposed by Marconi and Tarazona (U. Marini Bettolo Marconi and P. Tarazona, J. Chem. Phys. 110, 8032 (1999); J. Phys.: Condens. Matter 12, A413 (2000)) when hydrodynamic interactions are neglected. The derivation is based on Smoluchowski’s equation for the time evolution of the probability density with pairwise hydrodynamic interactions. The theory is applied to hard-sphere colloids in an oscillating spherical optical trap which switches periodically in time from a stable confining to an unstable potential. Rosenfeld’s fundamental measure theory for the equilibrium density functional is used and hydrodynamics are incorporated on the Rotne-Prager level. The results for the time-dependent density profiles are compared to extensive Brownian dynamics simulations which are performed on the same Rotne-Prager level and excellent agreement is obtained. It is further found that hydrodynamic interactions damp and slow the dynamics of the confined colloid cluster in comparison to the same situation with neglected hydrodynamic interactions.

  20. Precise effective masses from density functional perturbation theory

    NASA Astrophysics Data System (ADS)

    Laflamme Janssen, J.; Gillet, Y.; Poncé, S.; Martin, A.; Torrent, M.; Gonze, X.

    2016-05-01

    The knowledge of effective masses is a key ingredient to analyze numerous properties of semiconductors, like carrier mobilities, (magneto)transport properties, or band extrema characteristics yielding carrier densities and density of states. Currently, these masses are usually calculated using finite-difference estimation of density functional theory (DFT) electronic band curvatures. However, finite differences require an additional convergence study and are prone to numerical noise. Moreover, the concept of effective mass breaks down at degenerate band extrema. We assess the former limitation by developing a method that allows to obtain the Hessian of DFT bands directly, using density functional perturbation theory. Then, we solve the latter issue by adapting the concept of "transport equivalent effective mass" to the k .p ̂ framework. The numerical noise inherent to finite-difference methods is thus eliminated, along with the associated convergence study. The resulting method is therefore more general, more robust, and simpler to use, which makes it especially appropriate for high-throughput computing. After validating the developed techniques, we apply them to the study of silicon, graphane, and arsenic. The formalism is implemented into the abinit software and supports the norm-conserving pseudopotential approach, the projector augmented-wave method, and the inclusion of spin-orbit coupling. The derived expressions also apply to the ultrasoft pseudopotential method.

  1. Density Functional Theory for Steady-State Nonequilibrium Molecular Junctions

    NASA Astrophysics Data System (ADS)

    Liu, Shuanglong; Nurbawono, Argo; Zhang, Chun

    2015-10-01

    We present a density functional theory (DFT) for steady-state nonequilibrium quantum systems such as molecular junctions under a finite bias. Based on the steady-state nonequilibrium statistics that maps nonequilibrium to an effective equilibrium, we show that ground-state DFT (GS-DFT) is not applicable in this case and two densities, the total electron density and the density of current-carrying electrons, are needed to uniquely determine the properties of the corresponding nonequilibrium system. A self-consistent mean-field approach based on two densities is then derived. The theory is implemented into SIESTA computational package and applied to study nonequilibrium electronic/transport properties of a realistic carbon-nanotube (CNT)/Benzene junction. Results obtained from our steady-state DFT (SS-DFT) are compared with those of conventional GS-DFT based transport calculations. We show that SS-DFT yields energetically more stable nonequilibrium steady state, predicts significantly lower electric current, and is able to produce correct electronic structures in local equilibrium under a limiting case.

  2. Analyzing the financial crisis using the entropy density function

    NASA Astrophysics Data System (ADS)

    Oh, Gabjin; Kim, Ho-yong; Ahn, Seok-Won; Kwak, Wooseop

    2015-02-01

    The risk that is created by nonlinear interactions among subjects in economic systems is assumed to increase during an abnormal state of a financial market. Nevertheless, investigating the systemic risk in financial markets following the global financial crisis is not sufficient. In this paper, we analyze the entropy density function in the return time series for several financial markets, such as the S&P500, KOSPI, and DAX indices, from October 2002 to December 2011 and analyze the variability in the entropy value over time. We find that the entropy density function of the S&P500 index during the subprime crisis exhibits a significant decrease compared to that in other periods, whereas the other markets, such as those in Germany and Korea, exhibit no significant decrease during the market crisis. These findings demonstrate that the S&P500 index generated a regular pattern in the return time series during the financial crisis.

  3. Covariant density functional theory: The role of the pion

    SciTech Connect

    Lalazissis, G. A.; Karatzikos, S.; Serra, M.; Otsuka, T.; Ring, P.

    2009-10-15

    We investigate the role of the pion in covariant density functional theory. Starting from conventional relativistic mean field (RMF) theory with a nonlinear coupling of the {sigma} meson and without exchange terms we add pions with a pseudovector coupling to the nucleons in relativistic Hartree-Fock approximation. In order to take into account the change of the pion field in the nuclear medium the effective coupling constant of the pion is treated as a free parameter. It is found that the inclusion of the pion to this sort of density functionals does not destroy the overall description of the bulk properties by RMF. On the other hand, the noncentral contribution of the pion (tensor coupling) does have effects on single particle energies and on binding energies of certain nuclei.

  4. Semilocal density functional theory with correct surface asymptotics

    NASA Astrophysics Data System (ADS)

    Constantin, Lucian A.; Fabiano, Eduardo; Pitarke, J. M.; Della Sala, Fabio

    2016-03-01

    Semilocal density functional theory is the most used computational method for electronic structure calculations in theoretical solid-state physics and quantum chemistry of large systems, providing good accuracy with a very attractive computational cost. Nevertheless, because of the nonlocality of the exchange-correlation hole outside a metal surface, it was always considered inappropriate to describe the correct surface asymptotics. Here, we derive, within the semilocal density functional theory formalism, an exact condition for the imagelike surface asymptotics of both the exchange-correlation energy per particle and potential. We show that this condition can be easily incorporated into a practical computational tool, at the simple meta-generalized-gradient approximation level of theory. Using this tool, we also show that the Airy-gas model exhibits asymptotic properties that are closely related to those at metal surfaces. This result highlights the relevance of the linear effective potential model to the metal surface asymptotics.

  5. Excitations and benchmark ensemble density functional theory for two electrons

    SciTech Connect

    Pribram-Jones, Aurora; Burke, Kieron; Yang, Zeng-hui; Ullrich, Carsten A.; Trail, John R.; Needs, Richard J.

    2014-05-14

    A new method for extracting ensemble Kohn-Sham potentials from accurate excited state densities is applied to a variety of two-electron systems, exploring the behavior of exact ensemble density functional theory. The issue of separating the Hartree energy and the choice of degenerate eigenstates is explored. A new approximation, spin eigenstate Hartree-exchange, is derived. Exact conditions that are proven include the signs of the correlation energy components and the asymptotic behavior of the potential for small weights of the excited states. Many energy components are given as a function of the weights for two electrons in a one-dimensional flat box, in a box with a large barrier to create charge transfer excitations, in a three-dimensional harmonic well (Hooke's atom), and for the He atom singlet-triplet ensemble, singlet-triplet-singlet ensemble, and triplet bi-ensemble.

  6. Density-functional theory of inhomogeneous systems of hard spherocylinders

    PubMed

    Velasco; Mederos; Sullivan

    2000-09-01

    The smectic-A phase boundaries of a hard-spherocylinder fluid are calculated using a density-functional theory based on one proposed earlier by Somoza and Tarazona [Phys. Rev. A 41, 965 (1990)]. Our calculations do not employ the translation-rotation decoupling approximation used in previous density-functional theories. The calculated phase boundaries agree well with computer simulation results up to aspect ratios L/D approximately 5 and are in better agreement with the simulations than are previous theories. We generalize the model fluid by including long-range interactions with quadrupolar orientational symmetry, which are taken into account by mean-field approximation. For sufficiently large strength, these interactions produce a smectic-C phase, which undergoes either a continuous or weakly first-order transition to the smectic-A phase. The theory and numerical methods discussed here can be applied to the analysis of interfacial phenomena. PMID:11088887

  7. Neutrinoless double-beta decay in covariant density functional theory

    NASA Astrophysics Data System (ADS)

    Ring, P.; Yao, J. M.; Song, L. S.; Hagino, K.; Meng, J.

    2015-10-01

    We use covariant density functional theory beyond mean field in order to describe neutrinoless double-beta decay in a fully relativistic way. The dynamic effects of particle-number and angular-momentum conservations as well as shape fluctuations of quadrupole character are taken into account within the generator coordinate method for both initial and final nuclei. The calculations are based on the full relativistic transition operator. The nuclear matrix elements (NME's) for a large number of possible transitions are investigated. The results are compared with various non-relativistic calculations, in particular also with the density functional theory based on the Gogny force. We find that the non-relativistic approximation is justified and that the total NME's can be well approximated by the pure axial-vector coupling term. This corresponds to a considerable reduction of the computational effort.

  8. Neutrinoless double-beta decay in covariant density functional theory

    SciTech Connect

    Ring, P.; Yao, J. M.; Song, L. S.; Hagino, K.; Meng, J.

    2015-10-15

    We use covariant density functional theory beyond mean field in order to describe neutrinoless double-beta decay in a fully relativistic way. The dynamic effects of particle-number and angular-momentum conservations as well as shape fluctuations of quadrupole character are taken into account within the generator coordinate method for both initial and final nuclei. The calculations are based on the full relativistic transition operator. The nuclear matrix elements (NME’s) for a large number of possible transitions are investigated. The results are compared with various non-relativistic calculations, in particular also with the density functional theory based on the Gogny force. We find that the non-relativistic approximation is justified and that the total NME’s can be well approximated by the pure axial-vector coupling term. This corresponds to a considerable reduction of the computational effort.

  9. Density versus spin-density functional in DFT+U and DFT+DMFT

    NASA Astrophysics Data System (ADS)

    Park, Hyowon; Millis, Andrew; Marianetti, Chris

    2015-03-01

    The construction of multi-variable effective action theories such as DFT+U and DFT+DMFT requires the choice of a local subspace of correlated orbitals and an additional variable being either the charge density or spin density. This talk examines the differences between using charge-only and spin-dependent exchange-correlation functionals with the aim of providing guidance for constructing more sophisticated beyond-density functional theories. The widely used spin-dependent approximations to the exchange-correlation functional are found to lead to a large and in some cases unphysical effective exchange coupling within the correlated subspace. Additionally, the differences between Wannier and Projector based definitions of the correlated orbitals are examined, and only small differences are found provided that the orbitals are orthonormal and strongly localized. These results are documented in the context of the rare earth nickelates. This work is supported under the Grant DOE-ER-046169 and under the FAME grant, one of six centers of STARnet, a Semiconductor Research Corporation program sponsored by MARCO and DARPA.

  10. Subsystem density functional theory with meta-generalized gradient approximation exchange-correlation functionals.

    PubMed

    Śmiga, Szymon; Fabiano, Eduardo; Laricchia, Savio; Constantin, Lucian A; Della Sala, Fabio

    2015-04-21

    We analyze the methodology and the performance of subsystem density functional theory (DFT) with meta-generalized gradient approximation (meta-GGA) exchange-correlation functionals for non-bonded molecular systems. Meta-GGA functionals depend on the Kohn-Sham kinetic energy density (KED), which is not known as an explicit functional of the density. Therefore, they cannot be directly applied in subsystem DFT calculations. We propose a Laplacian-level approximation to the KED which overcomes this limitation and provides a simple and accurate way to apply meta-GGA exchange-correlation functionals in subsystem DFT calculations. The so obtained density and energy errors, with respect to the corresponding supermolecular calculations, are comparable with conventional approaches, depending almost exclusively on the approximations in the non-additive kinetic embedding term. An embedding energy error decomposition explains the accuracy of our method. PMID:25903880

  11. Probability density function modeling for sub-powered interconnects

    NASA Astrophysics Data System (ADS)

    Pater, Flavius; Amaricǎi, Alexandru

    2016-06-01

    This paper proposes three mathematical models for reliability probability density function modeling the interconnect supplied at sub-threshold voltages: spline curve approximations, Gaussian models,and sine interpolation. The proposed analysis aims at determining the most appropriate fitting for the switching delay - probability of correct switching for sub-powered interconnects. We compare the three mathematical models with the Monte-Carlo simulations of interconnects for 45 nm CMOS technology supplied at 0.25V.

  12. Autoionization in time-dependent density-functional theory

    NASA Astrophysics Data System (ADS)

    Kapoor, V.

    2016-06-01

    We compute the exact exchange-correlation potential of the time-dependent density-functional theory (TDDFT) for the correlated process of autoionization. The potential develops barriers which regulate the autoionization rate. TDDFT employing known and practicable exchange-correlation potentials does not capture any autoionization dynamics. Approximate exchange-correlation potentials capturing such dynamics would necessarily require memory effects and are unlikely to be developed, as will be illustrated.

  13. Density functional theory across chemistry, physics and biology

    PubMed Central

    van Mourik, Tanja; Bühl, Michael; Gaigeot, Marie-Pierre

    2014-01-01

    The past decades have seen density functional theory (DFT) evolve from a rising star in computational quantum chemistry to one of its major players. This Theme Issue, which comes half a century after the publication of the Hohenberg–Kohn theorems that laid the foundations of modern DFT, reviews progress and challenges in present-day DFT research. Rather than trying to be comprehensive, this Theme Issue attempts to give a flavour of selected aspects of DFT. PMID:24516181

  14. Association between lung function and airway wall density

    NASA Astrophysics Data System (ADS)

    Leader, J. Ken; Zheng, Bin; Fuhrman, Carl R.; Tedrow, John; Park, Sang C.; Tan, Jun; Pu, Jiantao; Drescher, John M.; Gur, David; Sciurba, Frank C.

    2009-02-01

    Computed tomography (CT) examination is often used to quantify the relation between lung function and airway remodeling in chronic obstructive pulmonary disease (COPD). In this preliminary study, we examined the association between lung function and airway wall computed attenuation ("density") in 200 COPD screening subjects. Percent predicted FVC (FVC%), percent predicted FEV1 (FEV1%), and the ratio of FEV1 to FVC as a percentage (FEV1/FVC%) were measured post-bronchodilator. The apical bronchus of the right upper lobe was manually selected from CT examinations for evaluation. Total airway area, lumen area, wall area, lumen perimeter and wall area as fraction of the total airway area were computed. Mean HU (meanHU) and maximum HU (maxHU) values were computed across pixels assigned membership in the wall and with a HU value greater than -550. The Pearson correlation coefficients (PCC) between FVC%, FEV1%, and FEV1/FVC% and meanHU were -0.221 (p = 0.002), -0.175 (p = 0.014), and -0.110 (p = 0.123), respectively. The PCCs for maxHU were only significant for FVC%. The correlations between lung function and the airway morphometry parameters were slightly stronger compared to airway wall density. MeanHU was significantly correlated with wall area (PCC = 0.720), airway area (0.498) and wall area percent (0.611). This preliminary work demonstrates that airway wall density is associated with lung function. Although the correlations in our study were weaker than a recent study, airway wall density initially appears to be an important parameter in quantitative CT analysis of COPD.

  15. Next generation high density self assembling functional protein arrays

    PubMed Central

    Ramachandran, Niroshan; Raphael, Jacob V.; Hainsworth, Eugenie; Demirkan, Gokhan; Fuentes, Manuel G.; Rolfs, Andreas; Hu, Yanhui; LaBaer, Joshua

    2009-01-01

    We report a high-density self assembling protein microarray that displays thousands of proteins, produced and captured in situ from immobilized cDNA templates. Over 1500 unique cDNAs were tested with > 90% success with nearly all proteins displaying yields within 2 fold of the mean, minimal sample variation and good day to day reproducibility. The displayed proteins revealed selective protein interactions. This method will enable various experimental approaches to study protein function in high throughput. PMID:18469824

  16. Probability Density Function Method for Langevin Equations with Colored Noise

    SciTech Connect

    Wang, Peng; Tartakovsky, Alexandre M.; Tartakovsky, Daniel M.

    2013-04-05

    We present a novel method to derive closed-form, computable PDF equations for Langevin systems with colored noise. The derived equations govern the dynamics of joint or marginal probability density functions (PDFs) of state variables, and rely on a so-called Large-Eddy-Diffusivity (LED) closure. We demonstrate the accuracy of the proposed PDF method for linear and nonlinear Langevin equations, describing the classical Brownian displacement and dispersion in porous media.

  17. Density functional theory across chemistry, physics and biology.

    PubMed

    van Mourik, Tanja; Bühl, Michael; Gaigeot, Marie-Pierre

    2014-03-13

    The past decades have seen density functional theory (DFT) evolve from a rising star in computational quantum chemistry to one of its major players. This Theme Issue, which comes half a century after the publication of the Hohenberg-Kohn theorems that laid the foundations of modern DFT, reviews progress and challenges in present-day DFT research. Rather than trying to be comprehensive, this Theme Issue attempts to give a flavour of selected aspects of DFT. PMID:24516181

  18. Plato: A localised orbital based density functional theory code

    NASA Astrophysics Data System (ADS)

    Kenny, S. D.; Horsfield, A. P.

    2009-12-01

    The Plato package allows both orthogonal and non-orthogonal tight-binding as well as density functional theory (DFT) calculations to be performed within a single framework. The package also provides extensive tools for analysing the results of simulations as well as a number of tools for creating input files. The code is based upon the ideas first discussed in Sankey and Niklewski (1989) [1] with extensions to allow high-quality DFT calculations to be performed. DFT calculations can utilise either the local density approximation or the generalised gradient approximation. Basis sets from minimal basis through to ones containing multiple radial functions per angular momenta and polarisation functions can be used. Illustrations of how the package has been employed are given along with instructions for its utilisation. Program summaryProgram title: Plato Catalogue identifier: AEFC_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFC_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 219 974 No. of bytes in distributed program, including test data, etc.: 1 821 493 Distribution format: tar.gz Programming language: C/MPI and PERL Computer: Apple Macintosh, PC, Unix machines Operating system: Unix, Linux and Mac OS X Has the code been vectorised or parallelised?: Yes, up to 256 processors tested RAM: Up to 2 Gbytes per processor Classification: 7.3 External routines: LAPACK, BLAS and optionally ScaLAPACK, BLACS, PBLAS, FFTW Nature of problem: Density functional theory study of electronic structure and total energies of molecules, crystals and surfaces. Solution method: Localised orbital based density functional theory. Restrictions: Tight-binding and density functional theory only, no exact exchange. Unusual features: Both atom centred and uniform meshes available

  19. Local and linear chemical reactivity response functions at finite temperature in density functional theory.

    PubMed

    Franco-Pérez, Marco; Ayers, Paul W; Gázquez, José L; Vela, Alberto

    2015-12-28

    We explore the local and nonlocal response functions of the grand canonical potential density functional at nonzero temperature. In analogy to the zero-temperature treatment, local (e.g., the average electron density and the local softness) and nonlocal (e.g., the softness kernel) intrinsic response functions are defined as partial derivatives of the grand canonical potential with respect to its thermodynamic variables (i.e., the chemical potential of the electron reservoir and the external potential generated by the atomic nuclei). To define the local and nonlocal response functions of the electron density (e.g., the Fukui function, the linear density response function, and the dual descriptor), we differentiate with respect to the average electron number and the external potential. The well-known mathematical relationships between the intrinsic response functions and the electron-density responses are generalized to nonzero temperature, and we prove that in the zero-temperature limit, our results recover well-known identities from the density functional theory of chemical reactivity. Specific working equations and numerical results are provided for the 3-state ensemble model. PMID:26723661

  20. Local and linear chemical reactivity response functions at finite temperature in density functional theory

    SciTech Connect

    Franco-Pérez, Marco E-mail: ayers@mcmaster.ca E-mail: avela@cinvestav.mx; Ayers, Paul W. E-mail: ayers@mcmaster.ca E-mail: avela@cinvestav.mx; Gázquez, José L. E-mail: ayers@mcmaster.ca E-mail: avela@cinvestav.mx; Vela, Alberto E-mail: ayers@mcmaster.ca E-mail: avela@cinvestav.mx

    2015-12-28

    We explore the local and nonlocal response functions of the grand canonical potential density functional at nonzero temperature. In analogy to the zero-temperature treatment, local (e.g., the average electron density and the local softness) and nonlocal (e.g., the softness kernel) intrinsic response functions are defined as partial derivatives of the grand canonical potential with respect to its thermodynamic variables (i.e., the chemical potential of the electron reservoir and the external potential generated by the atomic nuclei). To define the local and nonlocal response functions of the electron density (e.g., the Fukui function, the linear density response function, and the dual descriptor), we differentiate with respect to the average electron number and the external potential. The well-known mathematical relationships between the intrinsic response functions and the electron-density responses are generalized to nonzero temperature, and we prove that in the zero-temperature limit, our results recover well-known identities from the density functional theory of chemical reactivity. Specific working equations and numerical results are provided for the 3-state ensemble model.

  1. Advances in time-dependent current-density functional theory

    NASA Astrophysics Data System (ADS)

    Berger, Arjan

    In this work we solve the problem of the gauge dependence of molecular magnetic properties (magnetizabilities, circular dichroism) using time-dependent current-density functional theory [1]. We also present a new functional that accurately describes the optical absorption spectra of insulators, semiconductors and metals [2] N. Raimbault, P.L. de Boeij, P. Romaniello, and J.A. Berger Phys. Rev. Lett. 114, 066404 (2015) J.A. Berger, Phys. Rev. Lett. 115, 137402 (2015) This study has been partially supported through the Grant NEXT No. ANR-10-LABX-0037 in the framework of the Programme des Investissements d'Avenir.

  2. Thermal Corrections to Density Functional Simulations of Warm Dense Matter

    NASA Astrophysics Data System (ADS)

    Smith, Justin; Pribram-Jones, Aurora; Burke, Kieron

    Present density functional calculations of warm dense matter often use the Mermin-Kohn-Sham (MKS) scheme at finite temperature, but employ ground-state approximations to the exchange-correlation (XC) free energy. In the simplest solvable non-trivial model, an asymmetric Hubbard dimer, we calculate the exact many-body energies, the exact Mermin-Kohn-Sham functionals for this system, and extract the exact XC free energy. For moderate temperatures and weak correlation, we show this approximation is excellent, but fails for stronger correlations. Additionally, we use this system to test various conditions that must be satisfied.

  3. Atomistic force field for alumina fit to density functional theory

    SciTech Connect

    Sarsam, Joanne; Finnis, Michael W.; Tangney, Paul

    2013-11-28

    We present a force field for bulk alumina (Al{sub 2}O{sub 3}), which has been parametrized by fitting the energies, forces, and stresses of a large database of reference configurations to those calculated with density functional theory (DFT). We use a functional form that is simpler and computationally more efficient than some existing models of alumina parametrized by a similar technique. Nevertheless, we demonstrate an accuracy of our potential that is comparable to those existing models and to DFT. We present calculations of crystal structures and energies, elastic constants, phonon spectra, thermal expansion, and point defect formation energies.

  4. Linear Scaling Density Functional Calculations with Gaussian Orbitals

    NASA Technical Reports Server (NTRS)

    Scuseria, Gustavo E.

    1999-01-01

    Recent advances in linear scaling algorithms that circumvent the computational bottlenecks of large-scale electronic structure simulations make it possible to carry out density functional calculations with Gaussian orbitals on molecules containing more than 1000 atoms and 15000 basis functions using current workstations and personal computers. This paper discusses the recent theoretical developments that have led to these advances and demonstrates in a series of benchmark calculations the present capabilities of state-of-the-art computational quantum chemistry programs for the prediction of molecular structure and properties.

  5. Conformational Analysis of Thioether Musks Using Density Functional Theory

    PubMed Central

    Setzer, William N.

    2009-01-01

    A conformational analysis of nine macrocyclic thioether musks has been carried out using molecular mechanics (MMFF), density functional theory (DFT) using both B3LYP and M06 functionals, as well as Hartree-Fock and post-Hartree-Fock (MP2) ab initio methods. 6-Thia-, 10-thia- and 4-methyl-5-thia-14-tetradecananolide, 4-thia-, 7-thia-, 11-thia- and 12-thia-15-pentadecanolide and 6-thia- and 12-thia-16-hexadecanolide were modeled. Unfortunately, there was little agreement between the computational methods at the levels of theory used in this study. PMID:20111690

  6. A Wigner Monte Carlo approach to density functional theory

    NASA Astrophysics Data System (ADS)

    Sellier, J. M.; Dimov, I.

    2014-08-01

    In order to simulate quantum N-body systems, stationary and time-dependent density functional theories rely on the capacity of calculating the single-electron wave-functions of a system from which one obtains the total electron density (Kohn-Sham systems). In this paper, we introduce the use of the Wigner Monte Carlo method in ab-initio calculations. This approach allows time-dependent simulations of chemical systems in the presence of reflective and absorbing boundary conditions. It also enables an intuitive comprehension of chemical systems in terms of the Wigner formalism based on the concept of phase-space. Finally, being based on a Monte Carlo method, it scales very well on parallel machines paving the way towards the time-dependent simulation of very complex molecules. A validation is performed by studying the electron distribution of three different systems, a Lithium atom, a Boron atom and a hydrogenic molecule. For the sake of simplicity, we start from initial conditions not too far from equilibrium and show that the systems reach a stationary regime, as expected (despite no restriction is imposed in the choice of the initial conditions). We also show a good agreement with the standard density functional theory for the hydrogenic molecule. These results demonstrate that the combination of the Wigner Monte Carlo method and Kohn-Sham systems provides a reliable computational tool which could, eventually, be applied to more sophisticated problems.

  7. Density functional theory and chromium: Insights from the dimers

    SciTech Connect

    Würdemann, Rolf; Kristoffersen, Henrik H.; Moseler, Michael; Walter, Michael

    2015-03-28

    The binding in small Cr clusters is re-investigated, where the correct description of the dimer in three charge states is used as criterion to assign the most suitable density functional theory approximation. The difficulty in chromium arises from the subtle interplay between energy gain from hybridization and energetic cost due to exchange between s and d based molecular orbitals. Variations in published bond lengths and binding energies are shown to arise from insufficient numerical representation of electron density and Kohn-Sham wave-functions. The best functional performance is found for gradient corrected (GGA) functionals and meta-GGAs, where we find severe differences between functionals from the same family due to the importance of exchange. Only the “best fit” from Bayesian error estimation is able to predict the correct energetics for all three charge states unambiguously. With this knowledge, we predict small bond-lengths to be exclusively present in Cr{sub 2} and Cr{sub 2}{sup −}. Already for the dimer cation, solely long bond-lengths appear, similar to what is found in the trimer and in chromium bulk.

  8. Curvature and Frontier Orbital Energies in Density Functional Theory

    SciTech Connect

    Stein, Tamar; Autschbach, Jochen; Govind, Niranjan; Kronik, Leeor; Baer, Roi

    2012-12-20

    Perdew et al. [Phys. Rev. Lett 49, 1691 (1982)] discovered and proved two different properties that exact Kohn-Sham density functional theory (DFT) must obey: (i) The exact total energy versus particle number must be a series of linear segments between integer electron points; (ii) Across an integer number of electrons, the exchange-correlation potential may ``jump’’ by a constant, known as the derivative discontinuity (DD). Here, we show analytically that in both the original and the generalized Kohn-Sham formulation of DFT, the two are in fact two sides of the same coin. Absence of a derivative discontinuity necessitates deviation from piecewise linearity, and the latter can be used to correct for the former, thereby restoring the physical meaning of the orbital energies. Using selected small molecules, we show that this results in a simple correction scheme for any underlying functional, including semi-local and hybrid functionals as well as Hartree-Fock theory, suggesting a practical correction for the infamous gap problem of density functional theory. Moreover, we show that optimally-tuned range-separated hybrid functionals can inherently minimize both DD and curvature, thus requiring no correction, and show that this can be used as a sound theoretical basis for novel tuning strategies.

  9. Van der Waals Density Functional Theory with Applications

    NASA Astrophysics Data System (ADS)

    Langreth, David C.

    2004-03-01

    We discuss the development of electronic density functionals that are applicable for weakly bound systems where the van der Waals interaction and its ramifications become important. Our current functionals approach the correct asymptotic dependence at large distances and are seamless at small distances. The first form of the functional, appropriate for layered systems, has been recently applied to graphite, boron nitride, and molybdenum sulfide [H. Rydberg et al., Phys. Rev. Lett. 91, 126402 (2003) and D. C. Langreth, Int. J. Quant. Chem. (submitted), see http//:www.physics.rutgers.edu/ ˜langreth/preprints/dft2003.pdf]. The second form of the functional [M. Dion it et al. (to be published)] is appropriate for arbitrary geometries. Recent results on rare gas dimers and the benzene dimer suggest promise for this method as well.

  10. High-performance combinatorial algorithms

    SciTech Connect

    Pinar, Ali

    2003-10-31

    Combinatorial algorithms have long played an important role in many applications of scientific computing such as sparse matrix computations and parallel computing. The growing importance of combinatorial algorithms in emerging applications like computational biology and scientific data mining calls for development of a high performance library for combinatorial algorithms. Building such a library requires a new structure for combinatorial algorithms research that enables fast implementation of new algorithms. We propose a structure for combinatorial algorithms research that mimics the research structure of numerical algorithms. Numerical algorithms research is nicely complemented with high performance libraries, and this can be attributed to the fact that there are only a small number of fundamental problems that underlie numerical solvers. Furthermore there are only a handful of kernels that enable implementation of algorithms for these fundamental problems. Building a similar structure for combinatorial algorithms will enable efficient implementations for existing algorithms and fast implementation of new algorithms. Our results will promote utilization of combinatorial techniques and will impact research in many scientific computing applications, some of which are listed.

  11. Combinatorial optimization games

    SciTech Connect

    Deng, X.; Ibaraki, Toshihide; Nagamochi, Hiroshi

    1997-06-01

    We introduce a general integer programming formulation for a class of combinatorial optimization games, which immediately allows us to improve the algorithmic result for finding amputations in the core (an important solution concept in cooperative game theory) of the network flow game on simple networks by Kalai and Zemel. An interesting result is a general theorem that the core for this class of games is nonempty if and only if a related linear program has an integer optimal solution. We study the properties for this mathematical condition to hold for several interesting problems, and apply them to resolve algorithmic and complexity issues for their cores along the line as put forward in: decide whether the core is empty; if the core is empty, find an imputation in the core; given an imputation x, test whether x is in the core. We also explore the properties of totally balanced games in this succinct formulation of cooperative games.

  12. Generalized Pauli constraints in reduced density matrix functional theory

    SciTech Connect

    Theophilou, Iris; Helbig, Nicole; Lathiotakis, Nektarios N.; Marques, Miguel A. L.

    2015-04-21

    Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Coleman’s ensemble N-representability conditions. Recently, the topic of pure-state N-representability conditions, also known as generalized Pauli constraints, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble N-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone.

  13. Downlink Probability Density Functions for EOS-McMurdo Sound

    NASA Technical Reports Server (NTRS)

    Christopher, P.; Jackson, A. H.

    1996-01-01

    The visibility times and communication link dynamics for the Earth Observations Satellite (EOS)-McMurdo Sound direct downlinks have been studied. The 16 day EOS periodicity may be shown with the Goddard Trajectory Determination System (GTDS) and the entire 16 day period should be simulated for representative link statistics. We desire many attributes of the downlink, however, and a faster orbital determination method is desirable. We use the method of osculating elements for speed and accuracy in simulating the EOS orbit. The accuracy of the method of osculating elements is demonstrated by closely reproducing the observed 16 day Landsat periodicity. An autocorrelation function method is used to show the correlation spike at 16 days. The entire 16 day record of passes over McMurdo Sound is then used to generate statistics for innage time, outage time, elevation angle, antenna angle rates, and propagation loss. The levation angle probability density function is compared with 1967 analytic approximation which has been used for medium to high altitude satellites. One practical result of this comparison is seen to be the rare occurrence of zenith passes. The new result is functionally different than the earlier result, with a heavy emphasis on low elevation angles. EOS is one of a large class of sun synchronous satellites which may be downlinked to McMurdo Sound. We examine delay statistics for an entire group of sun synchronous satellites ranging from 400 km to 1000 km altitude. Outage probability density function results are presented three dimensionally.

  14. Revealing the density of encoded functions in a viral RNA

    PubMed Central

    Patel, Nikesh; Dykeman, Eric C.; Coutts, Robert H. A.; Lomonossoff, George P.; Rowlands, David J.; Phillips, Simon E. V.; Ranson, Neil; Twarock, Reidun; Tuma, Roman; Stockley, Peter G.

    2015-01-01

    We present direct experimental evidence that assembly of a single-stranded RNA virus occurs via a packaging signal-mediated mechanism. We show that the sequences of coat protein recognition motifs within multiple, dispersed, putative RNA packaging signals, as well as their relative spacing within a genomic fragment, act collectively to influence the fidelity and yield of capsid self-assembly in vitro. These experiments confirm that the selective advantages for viral yield and encapsidation specificity, predicted from previous modeling of packaging signal-mediated assembly, are found in Nature. Regions of the genome that act as packaging signals also function in translational and transcriptional enhancement, as well as directly coding for the coat protein, highlighting the density of encoded functions within the viral RNA. Assembly and gene expression are therefore direct molecular competitors for different functional folds of the same RNA sequence. The strongest packaging signal in the test fragment, encodes a region of the coat protein that undergoes a conformational change upon contact with packaging signals. A similar phenomenon occurs in other RNA viruses for which packaging signals are known. These contacts hint at an even deeper density of encoded functions in viral RNA, which if confirmed, would have profound consequences for the evolution of this class of pathogens. PMID:25646435

  15. Density-dependence of functional spiking networks in vitro

    SciTech Connect

    Ham, Michael I; Gintautuas, Vadas; Rodriguez, Marko A; Bettencourt, Luis M A; Bennett, Ryan; Santa Maria, Cara L

    2008-01-01

    During development, the mammalian brain differentiates into specialized regions with unique functional abilities. While many factors contribute to this functional specialization, we explore the effect neuronal density can have on neuronal interactions. Two types of networks, dense (50,000 neurons and glia support cells) and sparse (12,000 neurons and glia support cells), are studied. A competitive first response model is applied to construct activation graphs that represent pairwise neuronal interactions. By observing the evolution of these graphs during development in vitro we observe that dense networks form activation connections earlier than sparse networks, and that link-!llltropy analysis of the resulting dense activation graphs reveals that balanced directional connections dominate. Information theoretic measures reveal in addition that early functional information interactions (of order 3) are synergetic in both dense and sparse networks. However, during development in vitro, such interactions become redundant in dense, but not sparse networks. Large values of activation graph link-entropy correlate strongly with redundant ensembles observed in the dense networks. Results demonstrate differences between dense and sparse networks in terms of informational groups, pairwise relationships, and activation graphs. These differences suggest that variations in cell density may result in different functional specialization of nervous system tissue also in vivo.

  16. Current density partitioning in time-dependent current density functional theory

    SciTech Connect

    Mosquera, Martín A.; Wasserman, Adam; Department of Physics, Purdue University, West Lafayette, Indiana 47907

    2014-05-14

    We adapt time-dependent current density functional theory to allow for a fragment-based solution of the many-electron problem of molecules in the presence of time-dependent electric and magnetic fields. Regarding a molecule as a set of non-interacting subsystems that individually evolve under the influence of an auxiliary external electromagnetic vector-scalar potential pair, the partition 4-potential, we show that there are one-to-one mappings between this auxiliary potential, a sharply-defined set of fragment current densities, and the total current density of the system. The partition electromagnetic (EM) 4-potential is expressed in terms of the real EM 4-potential of the system and a gluing EM 4-potential that accounts for exchange-correlation effects and mutual interaction forces between fragments that are required to yield the correct electron dynamics. We prove the zero-force theorem for the fragmented system, establish a variational formulation in terms of action functionals, and provide a simple illustration for a charged particle in a ring.

  17. Subsystem real-time time dependent density functional theory.

    PubMed

    Krishtal, Alisa; Ceresoli, Davide; Pavanello, Michele

    2015-04-21

    We present the extension of Frozen Density Embedding (FDE) formulation of subsystem Density Functional Theory (DFT) to real-time Time Dependent Density Functional Theory (rt-TDDFT). FDE is a DFT-in-DFT embedding method that allows to partition a larger Kohn-Sham system into a set of smaller, coupled Kohn-Sham systems. Additional to the computational advantage, FDE provides physical insight into the properties of embedded systems and the coupling interactions between them. The extension to rt-TDDFT is done straightforwardly by evolving the Kohn-Sham subsystems in time simultaneously, while updating the embedding potential between the systems at every time step. Two main applications are presented: the explicit excitation energy transfer in real time between subsystems is demonstrated for the case of the Na4 cluster and the effect of the embedding on optical spectra of coupled chromophores. In particular, the importance of including the full dynamic response in the embedding potential is demonstrated. PMID:25903875

  18. Stochastic Optimally Tuned Range-Separated Hybrid Density Functional Theory.

    PubMed

    Neuhauser, Daniel; Rabani, Eran; Cytter, Yael; Baer, Roi

    2016-05-19

    We develop a stochastic formulation of the optimally tuned range-separated hybrid density functional theory that enables significant reduction of the computational effort and scaling of the nonlocal exchange operator at the price of introducing a controllable statistical error. Our method is based on stochastic representations of the Coulomb convolution integral and of the generalized Kohn-Sham density matrix. The computational cost of the approach is similar to that of usual Kohn-Sham density functional theory, yet it provides a much more accurate description of the quasiparticle energies for the frontier orbitals. This is illustrated for a series of silicon nanocrystals up to sizes exceeding 3000 electrons. Comparison with the stochastic GW many-body perturbation technique indicates excellent agreement for the fundamental band gap energies, good agreement for the band edge quasiparticle excitations, and very low statistical errors in the total energy for large systems. The present approach has a major advantage over one-shot GW by providing a self-consistent Hamiltonian that is central for additional postprocessing, for example, in the stochastic Bethe-Salpeter approach. PMID:26651840

  19. Density functional theory: Its origins, rise to prominence, and future

    NASA Astrophysics Data System (ADS)

    Jones, R. O.

    2015-07-01

    In little more than 20 years, the number of applications of the density functional (DF) formalism in chemistry and materials science has grown in an astonishing fashion. The number of publications alone shows that DF calculations make up a huge success story, and many younger colleagues are surprised to learn that the widespread application of density functional methods, particularly in chemistry, began only after 1990. This is indeed unexpected, because the origins are usually traced to the papers of Hohenberg, Kohn, and Sham more than a quarter of a century earlier. The DF formalism, its applications, and prospects were reviewed for this journal in 1989. About the same time, the combination of DF calculations with molecular dynamics promised to provide an efficient way to study structures and reactions in molecules and extended systems. This paper reviews the development of density-related methods back to the early years of quantum mechanics and follows the breakthrough in their application after 1990. The two examples from biochemistry and materials science are among the many current applications that were simply far beyond expectations in 1990. The reasons why—50 years after its modern formulation and after two decades of rapid expansion—some of the most cited practitioners in the field are concerned about its future are discussed.

  20. Applications of large-scale density functional theory in biology

    NASA Astrophysics Data System (ADS)

    Cole, Daniel J.; Hine, Nicholas D. M.

    2016-10-01

    Density functional theory (DFT) has become a routine tool for the computation of electronic structure in the physics, materials and chemistry fields. Yet the application of traditional DFT to problems in the biological sciences is hindered, to a large extent, by the unfavourable scaling of the computational effort with system size. Here, we review some of the major software and functionality advances that enable insightful electronic structure calculations to be performed on systems comprising many thousands of atoms. We describe some of the early applications of large-scale DFT to the computation of the electronic properties and structure of biomolecules, as well as to paradigmatic problems in enzymology, metalloproteins, photosynthesis and computer-aided drug design. With this review, we hope to demonstrate that first principles modelling of biological structure-function relationships are approaching a reality.

  1. Applications of large-scale density functional theory in biology.

    PubMed

    Cole, Daniel J; Hine, Nicholas D M

    2016-10-01

    Density functional theory (DFT) has become a routine tool for the computation of electronic structure in the physics, materials and chemistry fields. Yet the application of traditional DFT to problems in the biological sciences is hindered, to a large extent, by the unfavourable scaling of the computational effort with system size. Here, we review some of the major software and functionality advances that enable insightful electronic structure calculations to be performed on systems comprising many thousands of atoms. We describe some of the early applications of large-scale DFT to the computation of the electronic properties and structure of biomolecules, as well as to paradigmatic problems in enzymology, metalloproteins, photosynthesis and computer-aided drug design. With this review, we hope to demonstrate that first principles modelling of biological structure-function relationships are approaching a reality. PMID:27494095

  2. Van der Waals Interactions in Density Functional Theory: Intermolecular Complexes

    NASA Astrophysics Data System (ADS)

    Kannemann, Felix; Becke, Axel

    2010-03-01

    Conventional density functional theory (GGA and hybrid functionals) fails to account for dispersion interactions and is therefore not applicable to systems where van der Waals interactions play a dominant role, such as intermolecular complexes and biomolecules. The exchange-hole dipole moment (XDM) dispersion model of Becke and Johnson [A. D. Becke and E. R. Johnson, J. Chem. Phys. 127, 154108 (2007)] corrects for this deficiency. We have previously shown that the XDM dispersion model can be combined with standard GGA functionals (PW86 for exchange and PBE for correlation) to give accurate binding energy curves for rare-gas diatomics [F. O. Kannemann and A. D. Becke, J. Chem. Theory Comput. 5, 719 (2009)]. Here we present further tests of the GGA-XDM method using benchmark sets including hydrogen bonding, electrostatic, dispersion and stacking interactions, and systems ranging from rare-gas diatomics to biomolecular complexes.

  3. Nuclear chiral and magnetic rotation in covariant density functional theory

    NASA Astrophysics Data System (ADS)

    Meng, Jie; Zhao, Pengwei

    2016-05-01

    Excitations of chiral rotation observed in triaxial nuclei and magnetic and/or antimagnetic rotations (AMR) seen in near-spherical nuclei have attracted a lot of attention. Unlike conventional rotation in well-deformed or superdeformed nuclei, here the rotational axis is not necessary coinciding with any principal axis of the nuclear density distribution. Thus, tilted axis cranking (TAC) is mandatory to describe these excitations self-consistently in the framework of covariant density functional theory (CDFT). We will briefly introduce the formalism of TAC–CDFT and its application for magnetic and AMR phenomena. Configuration-fixed CDFT and its predictions for nuclear chiral configurations and for favorable triaxial deformation parameters are also presented, and the discoveries of the multiple chiral doublets in 133Ce and 103Rh are discussed.

  4. Hubbard operator density functional theory for Fermionic lattice models

    NASA Astrophysics Data System (ADS)

    Cheng, Zhengqian; Marianetti, Chris

    We formulate an effective action as a functional of Hubbard operator densities whose stationary point delivers all local static information of the interacting lattice model. Using the variational principle, we get a self-consistent equation for Hubbard operator densities. The computational cost of our approach is set by diagonalizing the local Fock space. We apply our method to the one and two band Hubbard model (including crystal field and on-site exchange) in infinite dimensions where the exact solution is known. Excellent agreement is obtained for the one-band model. In the two-band model, good agreement is obtained in the metallic region of the phase diagram in addition to the metal-insulator transition. While our approach does not address frequency dependent observables, it has a negligible computational cost as compared to dynamical mean field theory and could be highly applicable in the context total energies of strongly correlated materials and molecules.

  5. Density Functional Calculations for Atoms, Molecules and Clusters

    NASA Astrophysics Data System (ADS)

    Gunnarsson, O.; Jones, R. O.

    1980-01-01

    The density functional formalism provides a framework for including exchange and correlation effects in the calculation of ground state properties of many-electron systems. The reduction of the problem to the solution of single-particle equations leads to important numerical advantages over other ab initio methods of incorporating correlation effects. The essential features of the scheme are outlined and results obtained for atomic and molecular systems are surveyed. The local spin density (LSD) approximation gives generally good results for systems where the bonding involves s and p electrons, but results are less satisfactory for d-bonded systems. Non-local modifications to the LSD approximation have been tested on atomic systems yielding much improved total energies.

  6. Dipole polarizability of 120Sn and nuclear energy density functionals

    NASA Astrophysics Data System (ADS)

    Hashimoto, T.; Krumbholz, A. M.; Reinhard, P.-G.; Tamii, A.; von Neumann-Cosel, P.; Adachi, T.; Aoi, N.; Bertulani, C. A.; Fujita, H.; Fujita, Y.; Ganioǧlu, E.; Hatanaka, K.; Ideguchi, E.; Iwamoto, C.; Kawabata, T.; Khai, N. T.; Krugmann, A.; Martin, D.; Matsubara, H.; Miki, K.; Neveling, R.; Okamura, H.; Ong, H. J.; Poltoratska, I.; Ponomarev, V. Yu.; Richter, A.; Sakaguchi, H.; Shimbara, Y.; Shimizu, Y.; Simonis, J.; Smit, F. D.; Süsoy, G.; Suzuki, T.; Thies, J. H.; Yosoi, M.; Zenihiro, J.

    2015-09-01

    The electric dipole strength distribution in 120Sn between 5 and 22 MeV has been determined at the Research Center for Nuclear Physics, Osaka, from polarization transfer observables measured in proton inelastic scattering at E0=295 MeV and forward angles including 0∘. Combined with photoabsorption data, a highly precise electric dipole polarizability αD(120Sn) =8.93 (36 ) fm3 is extracted. The dipole polarizability as isovector observable par excellence carries direct information on nuclear symmetry energy and its density dependence. The correlation of the new value with the well-established αD(208Pb) serves as a test of its prediction by nuclear energy density functionals. Models based on modern Skyrme interactions describe the data fairly well while most calculations based on relativistic Hamiltonians cannot.

  7. Dynamic density functional theory of solid tumor growth: Preliminary models.

    PubMed

    Chauviere, Arnaud; Hatzikirou, Haralambos; Kevrekidis, Ioannis G; Lowengrub, John S; Cristini, Vittorio

    2012-03-01

    Cancer is a disease that can be seen as a complex system whose dynamics and growth result from nonlinear processes coupled across wide ranges of spatio-temporal scales. The current mathematical modeling literature addresses issues at various scales but the development of theoretical methodologies capable of bridging gaps across scales needs further study. We present a new theoretical framework based on Dynamic Density Functional Theory (DDFT) extended, for the first time, to the dynamics of living tissues by accounting for cell density correlations, different cell types, phenotypes and cell birth/death processes, in order to provide a biophysically consistent description of processes across the scales. We present an application of this approach to tumor growth. PMID:22489279

  8. Sublinear scaling for time-dependent stochastic density functional theory

    SciTech Connect

    Gao, Yi; Neuhauser, Daniel; Baer, Roi; Rabani, Eran

    2015-01-21

    A stochastic approach to time-dependent density functional theory is developed for computing the absorption cross section and the random phase approximation (RPA) correlation energy. The core idea of the approach involves time-propagation of a small set of stochastic orbitals which are first projected on the occupied space and then propagated in time according to the time-dependent Kohn-Sham equations. The evolving electron density is exactly represented when the number of random orbitals is infinite, but even a small number (≈16) of such orbitals is enough to obtain meaningful results for absorption spectrum and the RPA correlation energy per electron. We implement the approach for silicon nanocrystals using real-space grids and find that the overall scaling of the algorithm is sublinear with computational time and memory.

  9. Relativistic density functional theory for finite nuclei and neutron stars

    NASA Astrophysics Data System (ADS)

    Piekarewicz, Jorge

    In 1939 Oppenheimer and Volkoff demonstrated using Einstein's theory of general relativity that a neutron star supported exclusively by neutron degeneracy pressure will collapse into a black hole if its mass exceeds seven tenths of a solar mass. Seventy five years after such a pioneering prediction the existence of neutron stars with masses as large as two solar masses has been firmly established. This fact alone highlights the critical role that nuclear interactions play in explaining the structure of neutron stars. Indeed, a neutron star is a gold mine for the study of nuclear phenomena that span an enormous range of densities and neutron-proton asymmetries. Physical phenomena over such diverse scales are best described by a formalism based on Relativistic Density Functional Theory. In this contribution I focus on the synergy between theory, experiment, and observation that is needed to elucidate the myriad of exotic states of matter that are believed to exist in a neutron star.

  10. Progress at the interface of wave-function and density-functional theories

    SciTech Connect

    Gidopoulos, Nikitas I.

    2011-04-15

    The Kohn-Sham (KS) potential of density-functional theory (DFT) emerges as the minimizing effective potential in a variational scheme that does not involve fixing the unknown single-electron density. Using Rayleigh Schroedinger (RS) perturbation theory (PT), we construct ab initio approximations for the energy difference, the minimization of which determines the KS potential directly - thereby bypassing DFT's traditional algorithm to search for the density that minimizes the total energy. From second-order RS PT, we obtain variationally stable energy differences to be minimized, solving the severe problem of variational collapse of orbital-dependent exchange-correlation functionals based on second-order RS PT.