Finite-size Effects in a Two-Dimensional Electron Gas with Rashba Spin-Orbit Interaction

Within the Kubo formalism, we estimate the spin-Hall conductivity in a two-dimensional electron gas with Rashba spin-orbit interaction and study its variation as a function of disorder strength and system size. The numerical algorithm employed in the calculation is based on the direct numerical integration of the time-dependent Schrodinger equation in a spin-dependent variant of the particle source method. We find that the spin-precession length, L_s controlled by the strength of the Rashba coupling, establishes the critical lengthscale that marks the significant reduction of the spin-Hall conductivity in bulk systems. In contrast, the electron mean free path, inversely proportional to the strength of disorder, appears to have only a minor effect.Comment: 5 pages, 3 figure

Pressure-induced isostructural phase transition of metal-doped silicon clathrates

We propose an atomistic model for the pressure-induced isostructural phase transition of metal-doped silicon clathrates, Ba8Si46 and K8Si46, that has been observed at 14 GPa and 23 GPa, respectively. The model explains successfully the equation of state, transition pressure, change of Raman spectra and dependence on the doped cations as well as the effects of substituting Si(6c) atoms with noble metals.Comment: 5 pages, two coumn, 5 figures. See http://www.iitaka.org/down.html for more informatio

Fast Algorithm for Finding the Eigenvalue Distribution of Very Large Matrices

A theoretical analysis is given of the equation of motion method, due to Alben et al., to compute the eigenvalue distribution (density of states) of very large matrices. The salient feature of this method is that for matrices of the kind encountered in quantum physics the memory and CPU requirements of this method scale linearly with the dimension of the matrix. We derive a rigorous estimate of the statistical error, supporting earlier observations that the computational efficiency of this approach increases with matrix size. We use this method and an imaginary-time version of it to compute the energy and the specific heat of three different, exactly solvable, spin-1/2 models and compare with the exact results to study the dependence of the statistical errors on sample and matrix size.Comment: 24 pages, 24 figure

Calculating response functions in time domain with non-orthonormal basis sets

We extend the recently proposed order-N algorithms (cond-mat/9703224) for calculating linear- and nonlinear-response functions in time domain to the systems described by nonorthonormal basis sets.Comment: 4 pages, no figure

Temperature dependence of ESR intensity for the nanoscale molecular magnet V15

The electron spin resonance (ESR) of nanoscale molecular magnet ${\rm V}_{15}$ is studied. Since the Hamiltonian of ${\rm V}_{15}$ has a large Hilbert space and numerical calculations of the ESR signal evaluating the Kubo formula with exact diagonalization method is difficult, we implement the formula with the help of the random vector technique and the Chebyshev polynominal expansion, which we name the double Chebyshev expansion method. We calculate the temperature dependence of the ESR intensity of ${\rm V}_{15}$ and compare it with the data obtained in experiment. As another complementary approach, we also implement the Kubo formula with the subspace iteration method taking only important low-lying states into account. We study the ESR absorption curve below $100{\rm K}$ by means of both methods. We find that side peaks appear due to the Dzyaloshinsky-Moriya interaction and these peaks grows as temperature decreases.Comment: 9 pages, 4 figures. To appear in J. Phys. Soc. Jpn. Supp

Algorithm for Linear Response Functions at Finite Temperatures: Application to ESR spectrum of s=1/2 Antiferromagnet Cu benzoate

We introduce an efficient and numerically stable method for calculating linear response functions $\chi(\vec{q},\omega)$ of quantum systems at finite temperatures. The method is a combination of numerical solution of the time-dependent Schroedinger equation, random vector representation of trace, and Chebyshev polynomial expansion of Boltzmann operator. This method should be very useful for a wide range of strongly correlated quantum systems at finite temperatures. We present an application to the ESR spectrum of s=1/2 antiferromagnet Cu benzoate.Comment: 4 pages, 4 figure

Quantum Dynamics of Spin Wave Propagation Through Domain Walls

Through numerical solution of the time-dependent Schrodinger equation, we demonstrate that magnetic chains with uniaxial anisotropy support stable structures, separating ferromagnetic domains of opposite magnetization. These structures, domain walls in a quantum system, are shown to remain stable if they interact with a spin wave. We find that a domain wall transmits the longitudinal component of the spin excitations only. Our results suggests that continuous, classical spin models described by LLG equation cannot be used to describe spin wave-domain wall interaction in microscopic magnetic systems

Fast algorithm for calculating two-photon absorption spectra

We report a numerical calculation of the two-photon absorption coefficient of electrons in a binding potential using the real-time real-space higher-order difference method. By introducing random vector averaging for the intermediate state, the task of evaluating the two-dimensional time integral is reduced to calculating two one-dimensional integrals. This allows the reduction of the computation load down to the same order as that for the linear response function. The relative advantage of the method compared to the straightforward multi-dimensional time integration is greater for the calculation of non-linear response functions of higher order at higher energy resolution.Comment: 4 pages, 2 figures. It will be published in Phys. Rev. E on 1, March, 199

First Principles Calculation of Elastic Properties of Solid Argon at High Pressures

The density and the elastic stiffness coefficients of fcc solid argon at high pressures from 1 GPa up to 80 GPa are computed by first-principles pseudopotential method with plane-wave basis set and the generalized gradient approximation (GGA). The result is in good agreement with the experimental result recently obtained with the Brillouin spectroscopy by Shimizu et al. [Phys. Rev. Lett. 86, 4568 (2001)]. The Cauchy condition was found to be strongly violated as in the experimental result, indicating large contribution from non-central many-body force. The present result has made it clear that the standard density functional method with periodic boundary conditions can be successfully applied for calculating elastic properties of rare gas solids at high pressures in contrast to those at low pressures where dispersion forces are important.Comment: 4 pages, 5 figures, submitted to PR

Time-dependent properties of proton decay from crossing single-particle metastable states in deformed nuclei

A dynamical study of the decay of a metastable state by quantum tunneling through an anisotropic, non separable, two-dimensional potential barrier is performed by the numerical solution of the time-dependent Schrodinger equation. Initial quasi- stationary proton states are chosen in the framework of a deformed Woods-Saxon single-particle model. The decay of two sets of states corresponding to true and quasi levels-crossing is studied and the evolution of their decay properties as a function of nuclear deformation is calculated around the crossing point. The results show that the investigation of the proton decay from metastable states in deformed nuclei can unambiguously distinguish between the two types of crossing and determine the structure of the nuclear states involved.Comment: 15 pages, 9 figures, submitted to Phys. Rev.