Accurate quadratic-response approximation for the self-consistent pseudopotential of semiconductor nanostructures

Quadratic-response theory is shown to provide a conceptually simple but accurate approximation for the self-consistent one-electron potential of semiconductor nanostructures. Numerical examples are presented for GaAs/AlAs and InGaAs/InP (001) superlattices using the local-density approximation to density-functional theory and norm-conserving pseudopotentials without spin-orbit coupling. When the reference crystal is chosen to be the virtual-crystal average of the two bulk constituents, the absolute error in the quadratic-response potential for Gamma(15) valence electrons is about 2 meV for GaAs/AlAs and 5 meV for InGaAs/InP. Low-order multipole expansions of the electron density and potential response are shown to be accurate throughout a small neighborhood of each reciprocal lattice vector, thus providing a further simplification that is confirmed to be valid for slowly varying envelope functions. Although the linear response is about an order of magnitude larger than the quadratic response, the quadratic terms are important both quantitatively (if an accuracy of better than a few tens of meV is desired) and qualitatively (due to their different symmetry and long-range dipole effects).Comment: 16 pages, 20 figures; v2: new section on limitations of theor

Elimination of unoccupied state summations in it ab initio self-energy calculations for large supercells

We present a new method for the computation of self-energy corrections in large supercells. It eliminates the explicit summation over unoccupied states, and uses an iterative scheme based on an expansion of the Green's function around a set of reference energies. This improves the scaling of the computational time from the fourth to the third power of the number of atoms for both the inverse dielectric matrix and the self-energy, yielding improved efficiency for 8 or more silicon atoms per unit cell

Electron-phonon interaction in the solid form of the smallest fullerene C20_{20}

The electron-phonon coupling of a theoretically devised carbon phase made by assembling the smallest fullerenes C$_{20}$ is calculated from first principles. The structure consists of C$_{20}$ cages in an {\it fcc} lattice interlinked by two bridging carbon atoms in the interstitial tetrahedral sites ({\it fcc}-C$_{22}$). The crystal is insulating but can be made metallic by doping with interstitial alkali atoms. In the compound NaC$_{22}$ the calculated coupling constant $\lambda/N(0)$ is 0.28 eV, a value much larger than in C$_{60}$, as expected from the larger curvature of C$_{20}$. On the basis of the McMillan's formula, the calculated $\lambda$=1.12 and a $\mu^*$ assumed in the range 0.3-0.1 a superconducting T$_c$ in the range 15-55 K is predicted.Comment: 7 page

Integral quantification of seasonal soil moisture changes in farmland by cosmic-ray neutrons

Soil moisture at the plot or hill-slope scale is an important link between local vadose zone hydrology and catchment hydrology. However, so far only a few methods are on the way to close this gap between point measurements and remote sensing. One new measurement methodology that could determine integral soil moisture at this scale is the aboveground sensing of cosmic-ray neutrons, more precisely of ground albedo neutrons. The present study performed ground albedo neutron sensing (GANS) at an agricultural field in northern Germany. To test the method it was accompanied by other soil moisture measurements for a summer period with corn crops growing on the field and a later autumn-winter period without crops and a longer period of snow cover. Additionally, meteorological data and aboveground crop biomass were included in the evaluation. Hourly values of ground albedo neutron sensing showed a high statistical variability. Six-hourly values corresponded well with classical soil moisture measurements, after calibration based on one reference dry period and three wet periods of a few days each. Crop biomass seemed to influence the measurements only to minor degree, opposed to snow cover which has a more substantial impact on the measurements. The latter could be quantitatively related to a newly introduced field neutron ratio estimated from neutron counting rates of two energy ranges. Overall, our study outlines a procedure to apply the ground albedo neutron sensing method based on devices now commercially available, without the need for accompanying numerical simulations and suited for longer monitoring periods after initial calibration

Quadratic response theory for spin-orbit coupling in semiconductor heterostructures

This paper examines the properties of the self-energy operator in lattice-matched semiconductor heterostructures, focusing on nonanalytic behavior at small values of the crystal momentum, which gives rise to long-range Coulomb potentials. A nonlinear response theory is developed for nonlocal spin-dependent perturbing potentials. The ionic pseudopotential of the heterostructure is treated as a perturbation of a bulk reference crystal, and the self-energy is derived to second order in the perturbation. If spin-orbit coupling is neglected outside the atomic cores, the problem can be analyzed as if the perturbation were a local spin scalar, since the nonlocal spin-dependent part of the pseudopotential merely renormalizes the results obtained from a local perturbation. The spin-dependent terms in the self-energy therefore fall into two classes: short-range potentials that are analytic in momentum space, and long-range nonanalytic terms that arise from the screened Coulomb potential multiplied by a spin-dependent vertex function. For an insulator at zero temperature, it is shown that the electronic charge induced by a given perturbation is exactly linearly proportional to the charge of the perturbing potential. These results are used in a subsequent paper to develop a first-principles effective-mass theory with generalized Rashba spin-orbit coupling.Comment: 20 pages, no figures, RevTeX4; v2: final published versio

Disproportionation Phenomena on Free and Strained Sn/Ge(111) and Sn/Si(111) Surfaces

Distortions of the $\sqrt3\times\sqrt3$ Sn/Ge(111) and Sn/Si(111) surfaces are shown to reflect a disproportionation of an integer pseudocharge, $Q$, related to the surface band occupancy. A novel understanding of the $(3\times3)$-1U (``1 up, 2 down'') and 2U (``2 up, 1 down'') distortions of Sn/Ge(111) is obtained by a theoretical study of the phase diagram under strain. Positive strain keeps the unstrained value Q=3 but removes distorsions. Negative strain attracts pseudocharge from the valence band causing first a $(3\times3)$-2U distortion (Q=4) on both Sn/Ge and Sn/Si, and eventually a $(\sqrt3\times\sqrt3)$-3U (``all up'') state with Q=6. The possibility of a fluctuating phase in unstrained Sn/Si(111) is discussed.Comment: Revtex, 5 pages, 3 figure

Relation Between EEG Measures and Upper Limb Motor Recovery in Stroke Patients: A Scoping Review

Current clinical practice does not leverage electroencephalography (EEG) measurements in stroke patients, despite its potential to contribute to post-stroke recovery predictions. We review the literature on the effectiveness of various quantitative and qualitative EEG-based measures after stroke as a tool to predict upper limb motor outcome, in relation to stroke timeframe and applied experimental tasks. Moreover, we aim to provide guidance on the use of EEG in the assessment of upper limb motor recovery after stroke, suggesting a high potential for some metrics in the appropriate context. We identified relevant papers (N = 16) from databases ScienceDirect, Web of Science and MEDLINE, and assessed their methodological quality with the Joanna Briggs Institute (JBI) Critical Appraisal. We applied the Preferred Reporting Systems for Systematic Reviews and Meta-Analyses Extension for Scoping Reviews (PRISMA-ScR) Framework. Identified works used EEG to identify properties including event-related activation, spectral power in physiologically relevant bands, symmetry in brain dynamics, functional connectivity, cortico-muscular coherence and rhythmic coordination. EEG was acquired in resting state or in relation to behavioural conditions. Motor outcome was mainly evaluated with the Upper Limb Fugl-Meyer Assessment. Despite great variability in the literature, data suggests that the most promising EEG quantifiers for predicting post-stroke motor outcome are event-related measures. Measures of spectral power in physiologically relevant bands and measures of brain symmetry also show promise. We suggest that EEG measures may improve our understanding of stroke brain dynamics during recovery, and contribute to establishing a functional prognosis and choosing the rehabilitation approach

First-principles envelope-function theory for lattice-matched semiconductor heterostructures

In this paper a multi-band envelope-function Hamiltonian for lattice-matched semiconductor heterostructures is derived from first-principles norm-conserving pseudopotentials. The theory is applicable to isovalent or heterovalent heterostructures with macroscopically neutral interfaces and no spontaneous bulk polarization. The key assumption -- proved in earlier numerical studies -- is that the heterostructure can be treated as a weak perturbation with respect to some periodic reference crystal, with the nonlinear response small in comparison to the linear response. Quadratic response theory is then used in conjunction with k.p perturbation theory to develop a multi-band effective-mass Hamiltonian (for slowly varying envelope functions) in which all interface band-mixing effects are determined by the linear response. To within terms of the same order as the position dependence of the effective mass, the quadratic response contributes only a bulk band offset term and an interface dipole term, both of which are diagonal in the effective-mass Hamiltonian. Long-range multipole Coulomb fields arise in quantum wires or dots, but have no qualitative effect in two-dimensional systems beyond a dipole contribution to the band offsets.Comment: 25 pages, no figures, RevTeX4; v3: final published versio

Order-N Density-Matrix Electronic-Structure Method for General Potentials

A new order-N method for calculating the electronic structure of general (non-tight-binding) potentials is presented. The method uses a combination of the ``purification''-based approaches used by Li, Nunes and Vanderbilt, and Daw, and a representation of the density matrix based on ``travelling basis orbitals''. The method is applied to several one-dimensional examples, including the free electron gas, the ``Morse'' bound-state potential, a discontinuous potential that mimics an interface, and an oscillatory potential that mimics a semiconductor. The method is found to contain Friedel oscillations, quantization of charge in bound states, and band gap formation. Quantitatively accurate agreement with exact results is found in most cases. Possible advantages with regard to treating electron-electron interactions and arbitrary boundary conditions are discussed.Comment: 13 pages, REVTEX, 7 postscript figures (not quite perfect