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.

Design of antenna-coupled lumped-element titanium nitride KIDs for long-wavelength multi-band continuum imaging

Many applications in cosmology and astrophysics at millimeter wavelengths including CMB polarization, studies of galaxy clusters using the Sunyaev-Zeldovich effect (SZE), and studies of star formation at high redshift and in our local universe and our galaxy, require large-format arrays of millimeter-wave detectors. Feedhorn and phased-array antenna architectures for receiving mm-wave light present numerous advantages for control of systematics, for simultaneous coverage of both polarizations and/or multiple spectral bands, and for preserving the coherent nature of the incoming light. This enables the application of many traditional "RF" structures such as hybrids, switches, and lumped-element or microstrip band-defining filters.

Simultaneously, kinetic inductance detectors (KIDs) using high-resistivity materials like titanium nitride are an attractive sensor option for large-format arrays because they are highly multiplexable and because they can have sensitivities reaching the condition of background-limited detection. A KID is a LC resonator. Its inductance includes the geometric inductance and kinetic inductance of the inductor in the superconducting phase. A photon absorbed by the superconductor breaks a Cooper pair into normal-state electrons and perturbs its kinetic inductance, rendering it a detector of light. The responsivity of KID is given by the fractional frequency shift of the LC resonator per unit optical power.

However, coupling these types of optical reception elements to KIDs is a challenge because of the impedance mismatch between the microstrip transmission line exiting these architectures and the high resistivity of titanium nitride. Mitigating direct absorption of light through free space coupling to the inductor of KID is another challenge. We present a detailed titanium nitride KID design that addresses these challenges. The KID inductor is capacitively coupled to the microstrip in such a way as to form a lossy termination without creating an impedance mismatch. A parallel plate capacitor design mitigates direct absorption, uses hydrogenated amorphous silicon, and yields acceptable noise. We show that the optimized design can yield expected sensitivities very close to the fundamental limit for a long wavelength imager (LWCam) that covers six spectral bands from 90 to 400 GHz for SZE studies.

Excess phase (frequency) noise has been observed in KID and is very likely caused by two-level systems (TLS) in dielectric materials. The TLS hypothesis is supported by the measured dependence of the noise on resonator internal power and temperature. However, there is still a lack of a unified microscopic theory which can quantitatively model the properties of the TLS noise. In this thesis we derive the noise power spectral density due to the coupling of TLS with phonon bath based on an existing model and compare the theoretical predictions about power and temperature dependences with experimental data. We discuss the limitation of such a model and propose the direction for future study.

Magnetic field and pressure effects in a saturated gas laser amplifier

Theoretical and experimental studies of a gas laser amplifier are presented, assuming the amplifier is operating with a saturating optical frequency signal. The analysis is primarily concerned with the effects of the gas pressure and the presence of an axial magnetic field on the characteristics of the amplifying medium. Semiclassical radiation theory is used, along with a density matrix description of the atomic medium which relates the motion of single atoms to the macroscopic observables. A two-level description of the atom, using phenomenological source rates and decay rates, forms the basis of our analysis of the gas laser medium. Pressure effects are taken into account to a large extent through suitable choices of decay rate parameters.

Two methods for calculating the induced polarization of the atomic medium are used. The first method utilizes a perturbation expansion which is valid for signal intensities which barely reach saturation strength, and it is quite general in applicability. The second method is valid for arbitrarily strong signals, but it yields tractable solutions only for zero magnetic field or for axial magnetic fields large enough such that the Zeeman splitting is much larger than the power broadened homogeneous linewidth of the laser transition. The effects of pressure broadening of the homogeneous spectral linewidth are included in both the weak-signal and strong-signal theories; however the effects of Zeeman sublevel-mixing collisions are taken into account only in the weak-signal theory.

The behavior of a He-Ne gas laser amplifier in the presence of an axial magnetic field has been studied experimentally by measuring gain and Faraday rotation of linearly polarized resonant laser signals for various values of input signal intensity, and by measuring nonlinearity - induced anisotropy for elliptically polarized resonant laser signals of various input intensities. Two high-gain transitions in the 3.39-μ region were used for study: a J = 1 to J = 2 (3s₂ → 3p₄) transition and a J = 1 to J = 1 (3s₂ → 3p₂) transition. The input signals were tuned to the centers of their respective resonant gain lines.

The experimental results agree quite well with corresponding theoretical expressions which have been developed to include the nonlinear effects of saturation strength signals. The experimental results clearly show saturation of Faraday rotation, and for the J = 1 t o J = 1 transition a Faraday rotation reversal and a traveling wave gain dip are seen for small values of axial magnetic field. The nonlinearity induced anisotropy shows a marked dependence on the gas pressure in the amplifier tube for the J = 1 to J = 2 transition; this dependence agrees with the predictions of the general perturbational or weak signal theory when allowances are made for the effects of Zeeman sublevel-mixing collisions. The results provide a method for measuring the upper (neon 3s₂) level quadrupole moment decay rate, the dipole moment decay rates for the 3s₂ → 3p₄ and 3s₂ → 3p₂ transitions, and the effects of various types of collision processes on these decay rates.

Parity-dependent State Engineering and Tomography in the ultrastrong coupling regime

Reaching the strong coupling regime of light-matter interaction has led to an impressive development in fundamental quantum physics and applications to quantum information processing. Latests advances in different quantum technologies, like superconducting circuits or semiconductor quantum wells, show that the ultrastrong coupling regime (USC) can also be achieved, where novel physical phenomena and potential computational benefits have been predicted. Nevertheless, the lack of effective decoupling mechanism in this regime has so far hindered control and measurement processes. Here, we propose a method based on parity symmetry conservation that allows for the generation and reconstruction of arbitrary states in the ultrastrong coupling regime of light-matter interactions. Our protocol requires minimal external resources by making use of the coupling between the USC system and an ancillary two-level quantum system.

Correlations between the performance characteristics of a liquid crystal laser and the macroscopic material properties.

In order to understand how the performance of a liquid-crystal laser depends on the physical properties of the low molar mass nematic host, we have studied the energy threshold and slope efficiency of ten optically pumped liquid-crystal lasers based on different hosts. Specifically, this leads to a variation in the birefringence, the orientational order parameter, and the order parameter of the transition dipole moment of the dye. It is found that low threshold energies and high slope efficiencies correlate with high order parameters and large birefringences. To a first approximation this can be understood by considering analytical expressions for the threshold and slope efficiency, which are derived from the space-independent rate equations for a two-level system, in terms of the macroscopic liquid crystal properties.

Process development on the monolithic fabrication of an ultra-compact 4×4 optical switch matrix on InP/InGaAsP material

We have fabricated an ultra-compact 4×4 optical matrix on InP/InGaAsP material. 1×4 MMI couplers and TIR mirrors are employed to produce a compact 1×2 mm2 device. A CH4/H2/O2 RIE dry etch process has been used to realize two-level dry etching: deep-etch for both the MMI couplers and the mirrors and shallow-etch for the rest of the routing waveguides. It was found that a metal/dielectric bilayer mask is essential for multi-dry-etch processes and high profile verticality. We have found a Ti intermediate mask for the deep-etch process which is removable by SF6 dry-etch before the following shallow process. Dry-etch removal of the intermediate mask is necessary to protect the deep-etched mirror sidewall.

Effects of control of C/N ratio by low-cost carbohydrate addition on water quality and pond ecology in freshwater prawn Macrobrachium rosenbergii post-larvae nursing system

An experiment was conducted to evaluate the effects of control of carbon/nitrogen ratio (C/N ratio) by addition of low cost carbohydrate to the water column on water quality and pond ecology in freshwater prawn Macrobrachium rosenbergii post-larvae nursing system. In this experiment, two level of dietary protein 20% and 35% without carbohydrate addition (‘P20' and ‘P35') and with carbohydrate addition (‘P20+CH' and ‘P35+CH') were compared in small ponds of 40 m² area stocked with 20 post-larvae (0.021 ± 0.001g) per m² . Maize flour was used as low cost carbohydrate and applied to the water column followed by the first feeding during the day. The addition of carbohydrate significantly reduced (p< 0.05) ammonia-nitrogen (NH sub(3)-N) and nitrite-nitrogen (NO sub(2) - N) of water in P20 + CH and P35 + CH treatments. It significantly increased (p< 0.05) the total heterotrophic bacteria (THB) population both in water and sediment. Fifty nine genera of plankton were identified belonging to the Bacillariophyceae (11), Chlorophyceae (21), Cyanophyceae (7), Dinophyceae (1), Rotifera (7) and Crustacea (9) without any significant difference (p>0.05) of total phytoplankton and zooplankton among the treatments. Survival rate of prawn was significantly lowest (p<0.05) in P20 and no significant difference (p>0.05) was observed between P20+CH and P35 treatments. Control of C/N ratio by the addition of low-cost carbohydrate to the pond water column benefited the freshwater prawn nursing practices in three ways (1) increased heterotrophic bacterial growth supplying bacterial protein augment the prawn post-larvae growth performances, (2) reduced demand for supplemental feed protein and subsequent reduction in feed cost and (3) reduced toxic NH sub(3)-N and NO sub(2)-N levels in pond nursing system.

Parallel implementation of Multilevel BDDC

In application of the Balancing Domain Decomposition by Constraints (BDDC) to a case with many substructures, solving the coarse problem exactly becomes the bottleneck which spoils scalability of the solver. However, it is straightforward for BDDC to substitute the exact solution of the coarse problem by another step of BDDC method with subdomains playing the role of elements. In this way, the algorithm of three-level BDDC method is obtained. If this approach is applied recursively, multilevel BDDC method is derived. We present a detailed description of a recently developed parallel implementation of this algorithm. The implementation is applied to an engineering problem of linear elasticity and a benchmark problem of Stokes flow in a cavity. Results by the multilevel approach are compared to those by the standard (two-level) BDDC method.

Generalized power method for sparse principal component analysis

In this paper we develop a new approach to sparse principal component analysis (sparse PCA). We propose two single-unit and two block optimization formulations of the sparse PCA problem, aimed at extracting a single sparse dominant principal component of a data matrix, or more components at once, respectively. While the initial formulations involve nonconvex functions, and are therefore computationally intractable, we rewrite them into the form of an optimization program involving maximization of a convex function on a compact set. The dimension of the search space is decreased enormously if the data matrix has many more columns (variables) than rows. We then propose and analyze a simple gradient method suited for the task. It appears that our algorithm has best convergence properties in the case when either the objective function or the feasible set are strongly convex, which is the case with our single-unit formulations and can be enforced in the block case. Finally, we demonstrate numerically on a set of random and gene expression test problems that our approach outperforms existing algorithms both in quality of the obtained solution and in computational speed. © 2010 Michel Journée, Yurii Nesterov, Peter Richtárik and Rodolphe Sepulchre.

Single-photon source in strong photon-atom interaction regime in quasiperiodic photonic crystals

The antibunching properties of the fluorescence from a two-level ideal system in a 12-fold quasiperiodic photonic crystal are investigated based on the calculated local density of states. We found that the antibunching phenomenon of the fluorescence from two-level ideal systems could be significantly changed by varying their positions, i.e., perfect antibunching and antibunching with damped Rabi oscillation phenomenon occurred in different positions and at different frequencies in photonic crystals as a result of the large differences in the local density of states. This study revealed that the multi-level coherence of fluorescence from a two-level ideal system could be manipulated by controlling the position of the two-level ideal system in photonic crystals and the emission frequency in the photonic band structure. Copyright (C) EPLA, 2008

Electromagnetic-field-induced cohesive force enhancement in a metallic nanowire

We theoretically study the conducting electronic contribution to the cohesive force in a metallic nanowire irradiated under a transversely polarized external electromagnetic field at low temperatures and in the ballistic regime. In the framework of the free-electron model, we have obtained a time-dependent two-level electronic wavefunction by means of a unitary transformation. Using a thermodynamic statistical approach with this wavefunction, we have calculated the cohesive force in the nanowire. We show that the cohesive force can be divided into two components, one of which is independent of the electromagnetic field (static component), which is consistent with the existing results in the literature. The magnitude of the other component is proportional to the electromagnetic field strength. This extra component of the cohesive force is originally from the coherent coupling between the two lateral energy levels of the wire and the electromagnetic field.

Intraband relaxation and its influences on quantum dot lasers

A comprehensive two-level numerical model is developed to describe carrier distribution in a quantum-dot laser. Light-emission spectra with different intraband relaxation rates (2ps, 7.5ps and 20ps) are calculated and analysed to investigate the influence of relaxation rates on performance of the quantum-dot laser. The results indicate that fast intraband relaxation favours not only the ground state single mode operation but also the higher injection efficiency.

Dephasing rate in an InAs/GaAs single-electron quantum dot qubit

We have obtained the parameter-phase diagram, which unambiguously defines the parameter region for the use of InAs/GaAs quantum dot as two-level quantum system in quantum computation in the framework of the effective-mass envelope function theory. Moreover, static electric field is found to efficiently prolong decoherence time. As a result, decoherence time may reach the order of magnitude of milli-seconds as external static electric field goes beyond 20 kV/cm if only vacuum fluctuation is taken as the main source for decoherence. Our calculated results are useful for guiding the solid-state implementation of quantum computing.

InAs/GaAs single-electron quantum dot qubit

The time evolution of the quantum mechanical state of an electron is calculated in the framework of the effective-mass envelope function theory for an InAs/GaAs quantum dot. The results indicate that the superposition state electron density oscillates in the quantum dot, with a period on the order of femtoseconds. The interaction energy E-ij between two electrons located in different quantum dots is calculated for one electron in the ith pure quantum state and another in the jth pure quantum state. We find that E-11]E-12]E-22, and E-ij decreases as the distance between the two quantum dots increases. We present a parameter-phase diagram which defines the parameter region for the use of an InAs/GaAs quantum dot as a two-level quantum system in quantum computation. A static electric field is found to efficiently prolong the decoherence time. Our results should be useful for designing the solid-state implementation of quantum computing. (C) 2001 American Institute of Physics.

Hydrostatic pressure effect on photoluminescence from a GaN0.015As0.985/GaAs quantum well

Photoluminescence from a GaN0.015As0.985/GaAs quantum well has been measured at 15 K under hydrostatic pressure up to 9 GPa. Both the emissions from the GaNAs well and GaAs barrier are observed. The GaNAs-related peak shows a much weaker pressure dependence compared to that of the GaAs band gap. A group of new peaks appear in the spectra when the pressure is beyond 2.5 GPa, which is attributed to the emissions from the N isoelectronic traps in GaAs. The pressure dependence of the GaNAs-related peaks was calculated using the two-level model with the measured pressure coefficients of the GaAs band gap and N level as fitting parameters. It is found that the calculated results deviate seriously from the experimental data. An increasing of the emission intensity and the linewidth of the GaNAs-related peaks was also observed and briefly discussed. (C) 2001 American Institute of Physics.