A variational framework for spectral discretization of the density matrix in Kohn Sham density functional theory

Kohn-Sham density functional theory (KSDFT) is currently the main work-horse of quantum mechanical calculations in physics, chemistry, and materials science. From a mechanical engineering perspective, we are interested in studying the role of defects in the mechanical properties in materials. In real materials, defects are typically found at very small concentrations e.g., vacancies occur at parts per million, dislocation density in metals ranges from $10^{10} m^{-2}$ to $10^{15} m^{-2}$, and grain sizes vary from nanometers to micrometers in polycrystalline materials, etc. In order to model materials at realistic defect concentrations using DFT, we would need to work with system sizes beyond millions of atoms. Due to the cubic-scaling computational cost with respect to the number of atoms in conventional DFT implementations, such system sizes are unreachable. Since the early 1990s, there has been a huge interest in developing DFT implementations that have linear-scaling computational cost. A promising approach to achieving linear-scaling cost is to approximate the density matrix in KSDFT. The focus of this thesis is to provide a firm mathematical framework to study the convergence of these approximations. We reformulate the Kohn-Sham density functional theory as a nested variational problem in the density matrix, the electrostatic potential, and a field dual to the electron density. The corresponding functional is linear in the density matrix and thus amenable to spectral representation. Based on this reformulation, we introduce a new approximation scheme, called spectral binning, which does not require smoothing of the occupancy function and thus applies at arbitrarily low temperatures. We proof convergence of the approximate solutions with respect to spectral binning and with respect to an additional spatial discretization of the domain. For a standard one-dimensional benchmark problem, we present numerical experiments for which spectral binning exhibits excellent convergence characteristics and outperforms other linear-scaling methods.

Infared observations of H II regions

Spectral data are presented, giving intensities of the Brackett ɤ (B₇) line at six positions in M 42 and of the Brackett ten through fourteen (B₁₀-B₁₄) lines plus the He 4d³D-3p³p⁰ line at three positions in M 42. Observations of the Brackett ɤ line are also given for the planetary nebulae NGC 7027 and IC 418. Brackett gamma is shown to exhibit an anomalous satellite line in NGC 7027. Broadband data are presented, giving intensities at effective wavelengths of 1.25 μ, 1.65 μ, 2.2 μ, 3.5 μ and 4.8 μ for three positions in M 42.

Comparisons with visual and radio data as well as 12 micron and 20 micron data are used to derive reddening, electron temperatures, and electron densities for M 42 and the two planetaries, as well as a helium abundance for M 42. A representative electron temperature of 8400°K ± 1000°K, an electron density of 1.5 ±0.1 x 10³ cm^-3 and a He/H number density ratio of 0.10 +0.10/-0.05 are derived for the central region of M 42. The electron temperature is found to increase slightly with distance from the Trapezium.

M 42 is shown to emit in excess of the predicted recombination radiation throughout the entire infrared spectrum. The variations in the excess with wavelength and with position are analyzed to determine which of several physical processes may be operating. The longer wavelength infrared excess is shown to be dominated by dust emission, while the shorter wavelength infrared excess is caused by dust scattering. The dust is shown to be larger than the average interstellar particle. A new feature of the Orion red star ORS-1 is found in that it appears to have a reflection nebula around it.

A theory of strong and weak scintillations with applications to astrophysics

The propagation of waves in an extended, irregular medium is studied under the "quasi-optics" and the "Markov random process" approximations. Under these assumptions, a Fokker-Planck equation satisfied by the characteristic functional of the random wave field is derived. A complete set of the moment equations with different transverse coordinates and different wavenumbers is then obtained from the characteristic functional. The derivation does not require Gaussian statistics of the random medium and the result can be applied to the time-dependent problem. We then solve the moment equations for the phase correlation function, angular broadening, temporal pulse smearing, intensity correlation function, and the probability distribution of the random waves. The necessary and sufficient conditions for strong scintillation are also given.

We also consider the problem of diffraction of waves by a random, phase-changing screen. The intensity correlation function is solved in the whole Fresnel diffraction region and the temporal pulse broadening function is derived rigorously from the wave equation.

The method of smooth perturbations is applied to interplanetary scintillations. We formulate and calculate the effects of the solar-wind velocity fluctuations on the observed intensity power spectrum and on the ratio of the observed "pattern" velocity and the true velocity of the solar wind in the three-dimensional spherical model. The r.m.s. solar-wind velocity fluctuations are found to be ~200 km/sec in the region about 20 solar radii from the Sun.

We then interpret the observed interstellar scintillation data using the theories derived under the Markov approximation, which are also valid for the strong scintillation. We find that the Kolmogorov power-law spectrum with an outer scale of 10 to 100 pc fits the scintillation data and that the ambient averaged electron density in the interstellar medium is about 0.025 cm^-3. It is also found that there exists a region of strong electron density fluctuation with thickness ~10 pc and mean electron density ~7 cm^-3 between the PSR 0833-45 pulsar and the earth.

Survival, growth and production of Penaeus monodon and P. indicus at different density combinations with milkfish

Penaeus monodon and P. indicus juveniles were stocked and reared for about 3 months in earthen ponds at different density combinations with Chanos chanos. The presence of either Penaeus species at any density ratio did not affect significantly the C. chanos survival. Survival rates of the penaeids indicated that intraspecific and interspecific competition occurred and were reduced with the reduction in stocking rate. It is concluded that further studies on higher density ratios and feeding and economic consideration would be of help to the development of this kind of fishpond management system.

Electron microscopic examination of the viromatrix of Rana grylio virus in a fish cell line

Rana grylio virus (RGV), a Ranavirus belonging to the family Iridoviridae, assembles in the viromatrix which is a factory for viral genome replication and particle assembly. Ultrastructural studies of the viromatrix will clarify the pathway of assembly. The viromatrix and quantitative changes in RGV infected epithelipma papulosum cyprini (EPC) cells, one of fish cell lines, were studied by electron microscopy. It was shown that viromatrices were adjacent to the nucleus, and the electron density was lower than that of the surrounding cytoplasm. The viromatrix contained virus particles with different forms, electron-dense materials and amorphous structures which included tubules and membranous materials. Tubules were often observed in direct continuity with empty capsids. Several bundles of intermediate filaments were seen alongside the viromatrix and crystalline aggregates. Large clusters of mitochondria occurred in proximity to viromatrix. A total of 990 cells profiles were examined. The results showed that 394 cells contained viromatrix: 89.3% contained one, and 10.7% contained two to four viromatrices. The number of viromatrices increased gradually and reached a peak at 16 h p.i. The viromatrix area at 24 h p.i. increased up to 7.4 +/- 0.69 mu m(2) which was three-times lower than that at 6 h p.i. The number of empty capsids within viromatrix was generally more than that of "full" particles at different time points, and there was a strong positive correlation between them. (c) 2005 Elsevier B.V. All rights reserved.

Effects of polarization field on formation of two-dimensional electron gas in (0001) and (11(2)over-bar0) plane AlGaN/GaN heterostructures

Photoluminescence (PL) and temperature-dependent Hall effect measurements were carried out in (0001) and (11 (2) over bar0) AlGaN/GaN heterostructures grown on sapphire substrates by metalorganic chemical vapor deposition. There are strong spontaneous and piezoelectric electric fields (SPF) along the growth orientation of the (0001) AlGaN/GaN heterostructures. At the same time there are no corresponding SPF along that of the (1120) AlGaN/GaN. A strong PL peak related to the recombination between two-dimensional electron gas (2DEG) and photoexcited holes was observed at 3.258 eV at room temperature in (0001) AlGaN/GaN heterointerfaces while no corresponding PL peak was observed in (11 (2) over bar0). The existence of a 2DEG was observed in (0001) AlGaN/GaN multi-layers with a mobility saturated at 6000 cm(2)/V s below 80 K, whereas a much lower mobility was measured in (11 (2) over bar0). These results indicated that the SPF was the main element to cause the high mobility and high sheet-electron-density 2DEG in AlGaN/GaN heterostructures. (C) 2004 Elsevier B.V. All rights reserved.

Spin-polarized transport of two-dimensional electron gas embedded in a diluted magnetic semiconductor

The spin-polarized transport property of a diluted magnetic semiconductor two-dimensional electron gas is investigated theoretically at low temperature. A large current polarization can be found in this system even at small magnetic fields and oscillates with increasing magnetic field while the carrier polarization is vanishingly small. The magnitude as well as the sign of the current polarization can be tuned by varying magnetic field, the electron density and the Mn concentration. (c) 2005 American Institute of Physics.

Subband electron properties of InGaAs/InAlAs high-electron-mobility transistors with different channel chickness

Magnetotransport properties of In-0.53 GaAs/In-0.52 AlAs high electron mobility transistor (HEMT) structures with different channel thickness of 10-35 nm have been investigated in magnetic fields up to 13 T at 1.4 K. Fast Fourier transform has been employed to obtain the subband density and mobility of the two-dimensional electron gas in these HEMT structures. We found that the thickness of channel does not significantly enhance the electron density of the two-dimensional electron gas, however, it has strong effect on the proportion of electrons inhabited in different subbands. When the size of channel is 20 nm, the number of electrons occupying the excited subband, which have higher mobility, reaches the maximum. The experimental values obtained in this work are useful for the design and optimization of InGaAs/InAlAs HEMT devices.

Study of persistent photoconductivity and subband electronic properties of the two-dimensional electron gas in modulation doped GaAs/AlGaAs structure

We obtained the high mobility Of mu(2K) = 1.78 x 10(6) cm(2)/V . s in Si-doped GaAs/AlGaAs two-dimensional electron gas (2DEG) structures. After the sample was illuminated by a light-emitting diode in magnetic fields up to 6 T at T = 2K, we did observe the persistent photoconductivity effect and the electron density increased obviously. The electronic properties of 2DEG have been studied by Quantum-Hall-effect and Shubnikov-de Haas (SdH) oscillation measurements. We found that the electron concentrations of two subbands increase simultaneity with the increasing total electron concentration, and the electron mobility also increases obviously after being illuminated. At the same time, we also found that the electronic quantum lifetime becomes shorter, and a theoretical explunation is given through the widths of integral quantum Hall plateaus.

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.

Electronic states of a two-dimensional electron system in a lateral superlattice and a perpendicular magnetic field

The electronic state of a two-dimensional electron system (2DES) in the presence of a perpendicular uniform magnetic field and a lateral superlattice (LS) is investigated theoretically. A comparative study is made between a LS induced by a spatial electrostatic potential modulation (referred to as a PMLS) and that induced by a spatial magnetic-field modulation (referred ro asa MMLS). By utilizing a finite-temperature self-consistent Hartree-Fock approximation scheme; the dependence of the electronic state on different system parameters (e.g., the modulation period, the modulation strength, the effective electron-electron interaction strength, the averaged electron density, and the system temperature) is studied in detail. The inclusion of exchange effect is found to bring qualitative changes to the electronic state of a PMLS, leading generally to a nonuniform spin splitting, and consequently the behavior of the electronic state becomes similar to that of a MMLS. The Landau-level coupling is taken into account, and is found to introduce some interesting features not observed before. It is also found that, even in the regime of intermediate modulation strength, the density dependence of the spin splitting of energy levels, either for a PMLS or a MMLS, can be qualitatively understood within the picture of a 2DES in a perpendicular magnetic field with the modulation viewed as a perturbation. [S0163-1829(97)02248-0].

Self-consistent analysis of AlSb/InAs high electron mobility transistor structures

The influences of channel layer width, spacer layer width, and delta-doping density on the electron density and its distribution in the AlSb/InAs high electron mobility transistors (HEMTs) have been studied based on the self-consistent calculation of the Schrodinger and Poisson equations with both the strain and nonparabolicity effects being taken into account. The results show that, having little influence on the total two dimensional electron gas (2DEG) concentration in the channel, the HEMT's channel layer width has some influence on the electron mobility, with a channel as narrow as 100-130 angstrom being more beneficial. For the AlSb/InAs HEMT with a Te delta-doped layer, the 2DEG concentration as high as 9.1 X 10(12) cm(-2) can be achieved in the channel by enhancing the delta-doping concentration without the occurrence of the parallel conduction. When utilizing a Si delta-doped InAs layer as the electron-supplying layer of the AlSb/InAs HEMT, the effect of the InAs donor layer thickness is studied on the 2DEG concentration. To obtain a higher 2DEG concentration in the channel, it is necessary to use an InAs donor layer as thin as 4 monolayer. To test the validity of our calculation, we have compared our theoretical results (2DEG concentration and its distribution in different sub-bands of the channel) with the experimental ones done by other groups and show that our theoretical calculation is consistent with the experimental results.

Effects of polarization field on formation of two-dimensional electron gas in (0001) and (11(2)over-bar0) plane AlGaN/GaN heterostructures

Photoluminescence (PL) and temperature-dependent Hall effect measurements were carried out in (0001) and (11 (2) over bar0) AlGaN/GaN heterostructures grown on sapphire substrates by metalorganic chemical vapor deposition. There are strong spontaneous and piezoelectric electric fields (SPF) along the growth orientation of the (0001) AlGaN/GaN heterostructures. At the same time there are no corresponding SPF along that of the (1120) AlGaN/GaN. A strong PL peak related to the recombination between two-dimensional electron gas (2DEG) and photoexcited holes was observed at 3.258 eV at room temperature in (0001) AlGaN/GaN heterointerfaces while no corresponding PL peak was observed in (11 (2) over bar0). The existence of a 2DEG was observed in (0001) AlGaN/GaN multi-layers with a mobility saturated at 6000 cm(2)/V s below 80 K, whereas a much lower mobility was measured in (11 (2) over bar0). These results indicated that the SPF was the main element to cause the high mobility and high sheet-electron-density 2DEG in AlGaN/GaN heterostructures. (C) 2004 Elsevier B.V. All rights reserved.

Two-dimensional numerical simulation of the channel electron in an In0.52Al0.48As/In0.53Ga0.48As HEMT

A two-dimensional quantum model based on the solution of Schrodinger and Poisson equations is first presented for In0.52Al0.48As/In0.53Ga0.47As/InP HEMT. According to the model, the two-dimensional distributions of electron density and transverse electric field in the channel of InAlAs/InGaAs HEMT are discussed.

Low-compressibility and hard materials ReB2 and WB2: Prediction from first-principles study

First-principle calculations are performed to investigate the structural, elastic, and electronic properties of ReB2 and WB2. The calculated equilibrium structural parameters of ReB2 are consistent with the available experimental data. The calculations indicate that WB2 in the P6(3)/mmc space group is more energetically stable under the ambient condition than in the P6/mmm. Based on the calculated bulk modulus, shear modulus of polycrystalline aggregate, ReB2 and WB2 can be regarded as potential candidates of ultra-incompressible and hard materials. Furthermore, the elastic anisotropy is discussed by investigating the elastic stiffness constants. Density of states and electron density analysis unravel the covalent bonding between the transition metal atoms and the boron atoms as the driving force of the high bulk modulus and high shear modulus as well as small Poisson's ratio.