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.

A study of the ^9 Be nucleus

The reaction ⁷Li(³He, p)⁹Be has been used to measure excitations and intrinsic widths of levels in ⁹Be below the ⁷Li + d threshold. Previously unreported levels have been found at excitations of (13.78 ± .03) MeV and (16.671 ± .008) MeV with widths of (590 ± 60) keV and (41 ± 4) keV respectively. Two overlapping levels have been found at (11.81 ± .02) MeV and (11.29 ± .03) MeV with widths of (400 ± 30) keV and (620 ± 70) keV respectively. Branching ratios from ⁹Be levels populated in this reaction to the ground and first excited states of ⁸Be have been measured by observing the associated protons in coincidence with the decay neutrons. Branching ratios were found to be:

Excitation in ⁹Be .... Branching Ratio.......... Final Nucleus.........

(MeV) .......................... (percent) .....................................

.. 2.43 ........................... 7.5 ± 1.5 .............. ⁸Be(g.s.)

.. 3.03 ........................... 87 ± 13......................................

.. 4.65 ........................... 13 ± 4.......................................

.. 6.76 .............................. ≤ 2 ......................................

.. 11.29 ...............................≤ 2 ......................................

.. 11.81 ...............................≤ 3 ......................................

.. 6.76 ........................... .41 ≤ B.R. ≤ .69 ....... ⁸Be(2⁺)

.. 11.29 ........................... 14 ± 4 .......................................

.. 11.81 ........................... 12 ± 4 .......................................

Corresponding reduced widths for neutron emission are calculated and a comparison of the results with the expectations of current nuclear models is made. In particular the measured branching ratio to ⁸Be(g.s.) from ⁹Be(2.43 MeV) corresponds to an f-wave reduced with θ²_f = 2.1 x 10^-2, in units of ħ²/mR², with R = 4.35 fm. A comparison of this value with that predicted by a Nilsson model calculation, in which ⁹Be is taken to be a deformed nucleus, is discussed. The measured value for θ²_f is found to be consistent with that expected on the basis of measured E2-transition rates between rotational levels in ⁹Be.

Cholinergic neurons within the nucleus basalis magnocellularis in bilaterally lesioned rats with 192IgG-saporin.

En el presente estudio se han estudiado las neuronas colinérgicas del núcleo basal magnocelular y sus áreas de proyección en ratas lesionadas bilateralmente con la inmunotoxina 192IgG-saporina.

The statistical bootstrap model

A review is presented of the statistical bootstrap model of Hagedorn and Frautschi. This model is an attempt to apply the methods of statistical mechanics in high-energy physics, while treating all hadron states (stable or unstable) on an equal footing. A statistical calculation of the resonance spectrum on this basis leads to an exponentially rising level density ρ(m) ~ cm^-3 e^βom at high masses.

In the present work, explicit formulae are given for the asymptotic dependence of the level density on quantum numbers, in various cases. Hamer and Frautschi's model for a realistic hadron spectrum is described.

A statistical model for hadron reactions is then put forward, analogous to the Bohr compound nucleus model in nuclear physics, which makes use of this level density. Some general features of resonance decay are predicted. The model is applied to the process of NN annihilation at rest with overall success, and explains the high final state pion multiplicity, together with the low individual branching ratios into two-body final states, which are characteristic of the process. For more general reactions, the model needs modification to take account of correlation effects. Nevertheless it is capable of explaining the phenomenon of limited transverse momenta, and the exponential decrease in the production frequency of heavy particles with their mass, as shown by Hagedorn. Frautschi's results on "Ericson fluctuations" in hadron physics are outlined briefly. The value of β_o required in all these applications is consistently around [120 MeV]^-1 corresponding to a "resonance volume" whose radius is very close to ƛ_π. The construction of a "multiperipheral cluster model" for high-energy collisions is advocated.

On the Lyapunov transformation for stable matrices

The matrices studied here are positive stable (or briefly stable). These are matrices, real or complex, whose eigenvalues have positive real parts. A theorem of Lyapunov states that A is stable if and only if there exists H ˃ 0 such that AH + HA* = I. Let A be a stable matrix. Three aspects of the Lyapunov transformation L_A :H → AH + HA* are discussed.

1. Let C₁ (A) = {AH + HA* :H ≥ 0} and C₂ (A) = {H: AH+HA* ≥ 0}. The problems of determining the cones C₁(A) and C₂(A) are still unsolved. Using solvability theory for linear equations over cones it is proved that C₁(A) is the polar of C₂(A*), and it is also shown that C₁ (A) = C₁(A^-1). The inertia assumed by matrices in C₁(A) is characterized.

2. The index of dissipation of A was defined to be the maximum number of equal eigenvalues of H, where H runs through all matrices in the interior of C₂(A). Upper and lower bounds, as well as some properties of this index, are given.

3. We consider the minimal eigenvalue of the Lyapunov transform AH+HA*, where H varies over the set of all positive semi-definite matrices whose largest eigenvalue is less than or equal to one. Denote it by ψ(A). It is proved that if A is Hermitian and has eigenvalues μ₁ ≥ μ₂…≥ μ_n ˃ 0, then ψ(A) = -(μ₁-μ_n)²/(4(μ₁ + μ_n)). The value of ψ(A) is also determined in case A is a normal, stable matrix. Then ψ(A) can be expressed in terms of at most three of the eigenvalues of A. If A is an arbitrary stable matrix, then upper and lower bounds for ψ(A) are obtained.

Design of a highly stable, high-conversion-efficiency pumping source for optical parametric amplifier by extending efficient crystal length

The objective of this study is to improve the stability of pumping source of optical parametric amplifier. Analysis by simulation leads to the conclusion that the stability of the second harmonic can be improved by using properly the intensity of fundamental light and corresponding length of the crystal. By the method of the noncollinear two-pass second harmonic or the tandem second harmonic, the efficient crystal length is extended to a proper value, and the stability of the second harmonic output has been improved two times more than that for the fundamental light, and the conversion-efficiency is about 70% in experiment. When the variation of the fundamental light is about 10%, the variation of the second harmonic intensity has been controlled within 5%. (c) 2006 Elsevier Ltd. All rights reserved.

Stable single frequency and single polarization DFB fiber lasers operated at 1053 nm

Stable single-frequency and single-polarization distributed-feedback (DFB) fiber laser was realized by giving a pressure on the phase shift region of the fiber grating. The output wavelength of the DFB fiber laser is 1053 nm. When the pump power of 980 nm laser diode is 100 and 254 mW, the output power can reach 8.3 and 37.1 mW and the polarization extinction ratio was 26 and 20 dB, respectively. After chopped by Acousto-optic modulator (0.3 Hz), the pulse peak value variance is 4.65%(peak to peak) and 1.098% (RMS) for 31 min. (C) 2006 Elsevier Ltd. All rights reserved.

Highly stable, diode-pumped Nd : YLF regenerative amplifier

We present the design and experimental results for a diode pumped Nd:YLF regenerative amplifier applied to amplify a nanosecond laser pulse. Numerical simulation shows that the maximum output energy and the best stability can be obtained when the regenerative amplifier operates in a saturated mode for all pulse duration and temporal profiles. Using extra post-pulse is a good method to decrease the square-pulse distortion caused by gain saturation effect. The amplifier shows output energy of 4.2mJ with a total energy gain of more than 107 and output energy stability of better than 1% rms. When extra post-pulse is added, square-pulse distortion is decreased from 1.33 to 1.17 for the amplifier that is seeded with an optical pulse of 3 ns.

Magnetic Cross-Linked Enzyme Aggregates (mCLEAs) of Candida antarctica Lipase: An Efficient and Stable Biocatalyst for Biodiesel Synthesis

Enzyme-catalyzed production of biodiesel is the object of extensive research due to the global shortage of fossil fuels and increased environmental concerns. Herein we report the preparation and main characteristics of a novel biocatalyst consisting of Cross-Linked Enzyme Aggregates (CLEAs) of Candida antarctica lipase B (CALB) which are covalently bound to magnetic nanoparticles, and tackle its use for the synthesis of biodiesel from non-edible vegetable and waste frying oils. For this purpose, insolubilized CALB was covalently cross-linked to magnetic nanoparticles of magnetite which the surface was functionalized with –NH2 groups. The resulting biocatalyst combines the relevant catalytic properties of CLEAs (as great stability and feasibility for their reutilization) and the magnetic character, and thus the final product (mCLEAs) are superparamagnetic particles of a robust catalyst which is more stable than the free enzyme, easily recoverable from the reaction medium and reusable for new catalytic cycles. We have studied the main properties of this biocatalyst and we have assessed its utility to catalyze transesterification reactions to obtain biodiesel from non-edible vegetable oils including unrefined soybean, jatropha and cameline, as well as waste frying oil. Using 1% mCLEAs (w/w of oil) conversions near 80% were routinely obtained at 30°C after 24 h of reaction, this value rising to 92% after 72 h. Moreover, the magnetic biocatalyst can be easily recovered from the reaction mixture and reused for at least ten consecutive cycles of 24 h without apparent loss of activity. The obtained results suggest that mCLEAs prepared from CALB can become a powerful biocatalyst for application at industrial scale with better performance than those currently available.