A generic multibody simulation

Described is a dynamic simulation package which can be configured for orbital test scenarios involving multiple bodies. The rotational and translational state integration methods are selectable for each individual body and may be changed during a run if necessary. Characteristics of the bodies are determined by assigning components consisting of mass properties, forces, and moments, which are the outputs of user-defined environmental models. Generic model implementation is facilitated by a transformation processor which performs coordinate frame inversions. Transformations are defined in the initialization file as part of the simulation configuration. The simulation package includes an initialization processor, which consists of a command line preprocessor, a general purpose grammar, and a syntax scanner. These permit specifications of the bodies, their interrelationships, and their initial states in a format that is not dependent on a particular test scenario

Theory of the Stark Effect for P donors in Si

We develop a multi-valley effective mass theory for substitutional donors in silicon in an inhomogeneous environment. Valley-orbit coupling is treated perturbatively. We apply the theory to the Stark effect in Si:P. The method becomes more accurate at high fields, and it is designed to give correct experimental binding energies at zero field. Unexpectedly, the ground state energy for the donor electron is found to increase with electric field as a consequence of spectrum narrowing of the 1s manifold. Our results are of particular importance for the Kane quantum computer.Comment: published versio

Nonuniqueness and derivative discontinuities in density-functional theories for current-carrying and superconducting systems

Current-carrying and superconducting systems can be treated within density-functional theory if suitable additional density variables (the current density and the superconducting order parameter, respectively) are included in the density-functional formalism. Here we show that the corresponding conjugate potentials (vector and pair potentials, respectively) are {\it not} uniquely determined by the densities. The Hohenberg-Kohn theorem of these generalized density-functional theories is thus weaker than the original one. We give explicit examples and explore some consequences.Comment: revised version (typos corrected, some discussion added) to appear in Phys. Rev.

Determination of Compton profiles at solid surfaces from first-principles calculations

Projected momentum distributions of electrons, i.e. Compton profiles above the topmost atomic layer have recently become experimentally accessible by kinetic electron emission in grazing-incidence scattering of atoms at atomically flat single crystal metal surfaces. Sub-threshold emission by slow projectiles was shown to be sensitive to high-momentum components of the local Compton profile near the surface. We present a method to extract momentum distribution, Compton profiles, and Wigner and Husimi phase space distributions from ab-initio density-functional calculations of electronic structure. An application for such distributions to scattering experiments is discussed.Comment: 13 pages, 5 figures, submitted to PR

Molecular effects in the ionization of N2_2, O2_2 and F2_2 by intense laser fields

In this paper we study the response in time of N$_2$, O$_2$ and F$_2$ to laser pulses having a wavelength of 390nm. We find single ionization suppression in O$_2$ and its absence in F$_2$, in accordance with experimental results at $\lambda = 800$nm. Within our framework of time-dependent density functional theory we are able to explain deviations from the predictions of Intense-Field Many-Body $S$-Matrix Theory (IMST). We confirm the connection of ionization suppression with destructive interference of outgoing electron waves from the ionized electron orbital. However, the prediction of ionization suppression, justified within the IMST approach through the symmetry of the highest occupied molecular orbital (HOMO), is not reliable since it turns out that, e.g. in the case of F$_2$, the electronic response to the laser pulse is rather complicated and does not lead to dominant depletion of the HOMO. Therefore, the symmetry of the HOMO is not sufficient to predict ionization suppression. However, at least for F$_2$, the symmetry of the dominantly ionized orbital is consistent with the non-suppression of ionization.Comment: 19 pages, 5 figure

Quantum-Dot Cellular Automata using Buried Dopants

The use of buried dopants to construct quantum-dot cellular automata is investigated as an alternative to conventional electronic devices for information transport and elementary computation. This provides a limit in terms of miniaturisation for this type of system as each potential well is formed by a single dopant atom. As an example, phosphorous donors in silicon are found to have good energy level separation with incoherent switching times of the order of microseconds. However, we also illustrate the possibility of ultra-fast quantum coherent switching via adiabatic evolution. The switching speeds are numerically calculated and found to be 10's of picoseconds or less for a single cell. The effect of decoherence is also simulated in the form of a dephasing process and limits are estimated for operation with finite dephasing. The advantages and limitations of this scheme over the more conventional quantum-dot based scheme are discussed. The use of a buried donor cellular automata system is also discussed as an architecture for testing several aspects of buried donor based quantum computing schemes.Comment: Minor changes in response to referees comments. Improved section on scaling and added plot of incoherent switching time

First principles calculation of vibrational Raman spectra in large systems: signature of small rings in crystalline SiO2

We present an approach for the efficient calculation of vibrational Raman intensities in periodic systems within density functional theory. The Raman intensities are computed from the second order derivative of the electronic density matrix with respect to a uniform electric field. In contrast to previous approaches, the computational effort required by our method for the evaluation of the intensities is negligible compared to that required for the calculation of vibrational frequencies. As a first application, we study the signature of 3- and 4-membered rings in the the Raman spectra of several polymorphs of SiO2, including a zeolite having 102 atoms per unit cell.Comment: 4 pages, 2 figures, revtex4 Minor corrections; accepted in Phys. Rev. Let

Network inference using asynchronously updated kinetic Ising Model

Network structures are reconstructed from dynamical data by respectively naive mean field (nMF) and Thouless-Anderson-Palmer (TAP) approximations. For TAP approximation, we use two methods to reconstruct the network: a) iteration method; b) casting the inference formula to a set of cubic equations and solving it directly. We investigate inference of the asymmetric Sherrington- Kirkpatrick (S-K) model using asynchronous update. The solutions of the sets cubic equation depend of temperature T in the S-K model, and a critical temperature Tc is found around 2.1. For T < Tc, the solutions of the cubic equation sets are composed of 1 real root and two conjugate complex roots while for T > Tc there are three real roots. The iteration method is convergent only if the cubic equations have three real solutions. The two methods give same results when the iteration method is convergent. Compared to nMF, TAP is somewhat better at low temperatures, but approaches the same performance as temperature increase. Both methods behave better for longer data length, but for improvement arises, TAP is well pronounced.Comment: 6 pages, 4 figure

Band structure analysis of the conduction-band mass anisotropy in 6H and 4H SiC

The band structures of 6H and 4H SiC calculated by means of the FP-LMTO method are used to determine the effective mass tensors for their conduction-band minima. The results are shown to be consistent with recent optically detected cyclotron resonance measurements and predict an unusual band filling dependence for 6H-SiC.Comment: 5 pages including 4 postscript figures incorporated with epsfig figs. available as part 2: sicfig.uu self-extracting file to appear in Phys. Rev. B: Aug. 15 (Rapid Communications

Statistical mechanics of Floquet systems: the pervasive problem of near degeneracies

The statistical mechanics of periodically driven ("Floquet") systems in contact with a heat bath exhibits some radical differences from the traditional statistical mechanics of undriven systems. In Floquet systems all quasienergies can be placed in a finite frequency interval, and the number of near degeneracies in this interval grows without limit as the dimension N of the Hilbert space increases. This leads to pathologies, including drastic changes in the Floquet states, as N increases. In earlier work these difficulties were put aside by fixing N, while taking the coupling to the bath to be smaller than any quasienergy difference. This led to a simple explicit theory for the reduced density matrix, but with some major differences from the usual time independent statistical mechanics. We show that, for weak but finite coupling between system and heat bath, the accuracy of a calculation within the truncated Hilbert space spanned by the N lowest energy eigenstates of the undriven system is limited, as N increases indefinitely, only by the usual neglect of bath memory effects within the Born and Markov approximations. As we seek higher accuracy by increasing N, we inevitably encounter quasienergy differences smaller than the system-bath coupling. We therefore derive the steady state reduced density matrix without restriction on the size of quasienergy splittings. In general, it is no longer diagonal in the Floquet states. We analyze, in particular, the behavior near a weakly avoided crossing, where quasienergy near degeneracies routinely appear. The explicit form of our results for the denisty matrix gives a consistent prescription for the statistical mechanics for many periodically driven systems with N infinite, in spite of the Floquet state pathologies.Comment: 31 pages, 3 figure