First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

I. Numerical solution of exact pair equations. II. Numerical solution of first-order pair equations

Part I

Numerical solutions to the S-limit equations for the helium ground state and excited triplet state and the hydride ion ground state are obtained with the second and fourth difference approximations. The results for the ground states are superior to previously reported values. The coupled equations resulting from the partial wave expansion of the exact helium atom wavefunction were solved giving accurate S-, P-, D-, F-, and G-limits. The G-limit is -2.90351 a.u. compared to the exact value of the energy of -2.90372 a.u.

Part II

The pair functions which determine the exact first-order wavefunction for the ground state of the three-electron atom are found with the matrix finite difference method. The second- and third-order energies for the (1s1s)¹S, (1s2s)³S, and (1s2s)¹S states of the two-electron atom are presented along with contour and perspective plots of the pair functions. The total energy for the three-electron atom with a nuclear charge Z is found to be E(Z) = -1.125•Z² +1.022805•Z-0.408138-0.025515•(1/Z)+O(1/Z²)a.u.

Spectral density of first order Piecewise linear systems excited by white noise

The Fokker-Planck (FP) equation is used to develop a general method for finding the spectral density for a class of randomly excited first order systems. This class consists of systems satisfying stochastic differential equations of form ẋ + f(x) = m/Ʃ/j = 1 h_j(x)n_j(t) where f and the h_j are piecewise linear functions (not necessarily continuous), and the n_j are stationary Gaussian white noise. For such systems, it is shown how the Laplace-transformed FP equation can be solved for the transformed transition probability density. By manipulation of the FP equation and its adjoint, a formula is derived for the transformed autocorrelation function in terms of the transformed transition density. From this, the spectral density is readily obtained. The method generalizes that of Caughey and Dienes, J. Appl. Phys., 32.11.

This method is applied to 4 subclasses: (1) m = 1, h₁ = const. (forcing function excitation); (2) m = 1, h₁ = f (parametric excitation); (3) m = 2, h₁ = const., h₂ = f, n₁ and n₂ correlated; (4) the same, uncorrelated. Many special cases, especially in subclass (1), are worked through to obtain explicit formulas for the spectral density, most of which have not been obtained before. Some results are graphed.

Dealing with parametrically excited first order systems leads to two complications. There is some controversy concerning the form of the FP equation involved (see Gray and Caughey, J. Math. Phys., 44.3); and the conditions which apply at irregular points, where the second order coefficient of the FP equation vanishes, are not obvious but require use of the mathematical theory of diffusion processes developed by Feller and others. These points are discussed in the first chapter, relevant results from various sources being summarized and applied. Also discussed is the steady-state density (the limit of the transition density as t → ∞).

The first fifty years of the Freshwater Biological Association

A brief history of the Freshwater Biological Association is given in this article. The association has been in existence for fifty years, having been founded in 1929.

Joint research project on the aquatic ecotoxicity of organo-aluminium: Chemical reactions. First interim report

The overall goal of the joint research project is to relate the chemical reactions involved in the formation of organo-aluminium complexes under acid conditions to their toxic effects on the physiology of aquatic organisms. Finally, this research is intended to predict toxic effects arising from acidity and aluminium under varying environmental conditions. This interim report examines the chemical modelling of ion-binding by humic substances where a computer model has been developed and is being tested using field data, and conditions required for the precipitation of aluminium in surface waters.

I. Site effects and exciton structure in molecular crystals - Benzene. II. The first and second triplet states of solid benzene

The effect of intermolecular coupling in molecular energy levels (electronic and vibrational) has been investigated in neat and isotopic mixed crystals of benzene. In the isotopic mixed crystals of C₆H₆, C₆H₅D, m-C₆H₄D₂, p-C₆H₄D₂, sym-C₆H₃D₃, C₆D₅H, and C₆D₆ in either a C₆H₆ or C₆D₆ host, the following phenomena have been observed and interpreted in terms of a refined Frenkel exciton theory: a) Site shifts; b) site group splittings of the degenerate ground state vibrations of C₆H₆, C₆D₆, and sym-C₆H₃D₃; c) the orientational effect for the isotopes without a trigonal axis in both the ¹B_2u electronic state and the ground state vibrations; d) intrasite Fermi resonance between molecular fundamentals due to the reduced symmetry of the crystal site; and e) intermolecular or intersite Fermi resonance between nearly degenerate states of the host and guest molecules. In the neat crystal experiments on the ground state vibrations it was possible to observe many of these phenomena in conjunction with and in addition to the exciton structure.

To theoretically interpret these diverse experimental data, the concepts of interchange symmetry, the ideal mixed crystal, and site wave functions have been developed and are presented in detail. In the interpretation of the exciton data the relative signs of the intermolecular coupling constants have been emphasized, and in the limit of the ideal mixed crystal a technique is discussed for locating the exciton band center or unobserved exciton components. A differentiation between static and dynamic interactions is made in the Frenkel limit which enables the concepts of site effects and exciton coupling to be sharpened. It is thus possible to treat the crystal induced effects in such a fashion as to make their similarities and differences quite apparent.

A calculation of the ground state vibrational phenomena (site shifts and splittings, orientational effects, and exciton structure) and of the crystal lattice modes has been carried out for these systems. This calculation serves as a test of the approximations of first order Frenkel theory and the atom-atom, pair wise interaction model for the intermolecular potentials. The general form of the potential employed was V(r) = Be^-Cr - A/r⁶ ; the force constants were obtained from the potential by assuming the atoms were undergoing simple harmonic motion.

In part II the location and identification of the benzene first and second triplet states (³B_1u and ³E_1u) is given.

Decoding cosets of first-order Reed-Muller code

Proper encoding of transmitted information can improve the performance of a communication system. To recover the information at the receiver it is necessary to decode the received signal. For many codes the complexity and slowness of the decoder is so severe that the code is not feasible for practical use. This thesis considers the decoding problem for one such class of codes, the comma-free codes related to the first-order Reed-Muller codes.

A factorization of the code matrix is found which leads to a simple, fast, minimum memory, decoder. The decoder is modular and only n modules are needed to decode a code of length 2ⁿ. The relevant factorization is extended to any code defined by a sequence of Kronecker products.

The problem of monitoring the correct synchronization position is also considered. A general answer seems to depend upon more detailed knowledge of the structure of comma-free codes. However, a technique is presented which gives useful results in many specific cases.

Optimisation of the first principle code Octopus for massive parallel architectures: application to light harvesting complexes

150 p.

Proceedings of the first world meeting on bigeye tuna

The annual catches of big eye are exceeded by those of only two other species of tuna, skipjack, Katsuwonus pelamis, and yellowfin, Thunnus albacares. However, because most of the bigeye caught are consumed fresh, whereas most of the skipjack and yellowfin caught are canned, the economic value of big eye exceeds that of any other species of tuna. Despite its importance, less is known of the biology of bigeye than of the biology of any of the other principal market species of tunas. Historically, bigeye have been harvested mostly by longlines, which take only medium to large fish. During recent years, however, greater amounts of small bigeye have been caught by purse seines and other surface gear. This is a matter of concern for several reasons. First, long line fishermen are concerned that the harvesting of small bigeye will decrease the amounts of medium to large bigeye available to them. Second, since small bigeye are canned, rather than eaten fresh, consumers are concerned about the possible decrease in the supply of high-quality fresh fish. Third, economists are concerned about the possible economic loss associated with harvesting fish at less than their maximum economic value. Fourth, biologists are concerned about the possibility that harvesting of small bigeye could decrease the overall catches of that species. These concerns cannot be properly addressed until more knowledge of the biology of big eye is available. The purposes of the meeting were to review and discuss the information available and to make recommendations for further research.

Integrated treatment of first episode psychosis with online training (e-learning): study protocol for a randomised controlled trial

Background: The integrated treatment of first episode psychosis has been shown to improve functionality and negative symptoms in previous studies. In this paper, we describe a study of integrated treatment (individual psychoeducation complementary to pharmacotherapy) versus treatment as usual, comparing results at baseline with those at 6-month re-assessment (at the end of the study) for these patients, and online training of professionals to provide this complementary treatment, with the following objectives: 1) to compare the efficacy of individual psychoeducation as add-on treatment versus treatment as usual in improving psychotic and mood symptoms; 2) to compare adherence to medication, functioning, insight, social response, quality of life, and brain-derived neurotrophic factor, between both groups; and 3) to analyse the efficacy of online training of psychotherapists. Methods/design: This is a single-blind randomised clinical trial including patients with first episode psychosis from hospitals across Spain, randomly assigned to either a control group with pharmacotherapy and regular sessions with their psychiatrist (treatment as usual) or an intervention group with integrated care including treatment as usual plus a psychoeducational intervention (14 sessions). Training for professionals involved at each participating centre was provided by the coordinating centre (University Hospital of Alava) through video conferences. Patients are evaluated with an extensive battery of tests assessing clinical and sociodemographic characteristics (Positive and Negative Syndrome Scale, State-Trait Anxiety Inventory, Liebowitz Social Anxiety Scale, Hamilton Rating Scale for Depression, Scale to Assess Unawareness of Mental Disorders, Strauss and Carpenter Prognostic Scale, Global Assessment of Functioning Scale, Morisky Green Adherence Scale, Functioning Assessment Short Test, World Health Organization Quality of Life instrument WHOQOL-BREF (an abbreviated version of the WHOQOL-100), and EuroQoL questionnaire), and brain-derived neurotrophic factor levels are measured in peripheral blood at baseline and at 6 months. The statistical analysis, including bivariate analysis, linear and logistic regression models, will be performed using SPSS. Discussion: This is an innovative study that includes the assessment of an integrated intervention for patients with first episode psychosis provided by professionals who are trained online, potentially making it possible to offer the intervention to more patients.

O Apostolado Positivista e o projeto da Ditadura Republicana no Brasil (1889-1891)

O objetivo da pesquisa foi buscar, por meio de um estudo aprofundado, informações relevantes da ação política do Apostolado Positivista do Brasil na transição entre a Monarquia e a República. O recorte temporal (1889-1891) remete o trabalho ao período do primeiro governo da república, chamado de governo provisório, que teve como Presidente, o Marechal Deodoro da Fonseca. O trabalho foi dividido em três capítulos: o primeiro ressalta o contexto do surgimento do pensamento positivista por intermédio da atuação de seu mentor Augusto Comte no debate intelectual da França. A partir disso, destacam-se aspectos que servem para esclarecer a cisão entre positivistas ortodoxos e heterodoxos e posteriormente, como as idéias positivistas se espalharam no contexto brasileiro. No mesmo capítulo, foi ressaltado ainda, o surgimento da Igreja Positivista do Brasil, as trajetórias de Miguel Lemos e Raimundo Teixeira Mendes até assumirem a direção e vice-direção do Apostolado. No segundo foi apresentada a ação política do Apostolado na tentativa de encaminhar o governo provisório para um regime político de ditadura definitiva. Tendo seu projeto político fracassado, o grupo passou a defender o cumprimento de artigos da constituição de 1891, que por diferentes motivos tivessem aspectos que se assemelhavam com a política tida pelo Apostolado como adequada. No terceiro foi mostrada a influência do positivismo em diferentes temas importantes no Brasil naquele momento, como, ensino, separação entre Estado e Igreja e movimento trabalhista. A persistência do lema positivista ordem e progresso na bandeira nacional até hoje, está entre os símbolos desta militância ortodoxa. Os membros do Apostolado Positivista não tiveram seu projeto de ditadura sacramentado apesar do forte influência da doutrina positivista no Brasil por dois aspectos, o primeiro por seguirem a recomendação de Comte de negar cargos políticos, o segundo pela doutrina positivista não ser orgânica a oligarquia cafeeira que iria nortear as bases de estruturação do país.