Canonical density matrix perturbation theory

Density matrix perturbation theory [Niklasson and Challacombe, Phys. Rev. Lett. 92, 193001 (2004)] is generalized to canonical (NVT) free energy ensembles in tight-binding, Hartree-Fock or Kohn-Sham density functional theory. The canonical density matrix perturbation theory can be used to calculate temperature dependent response properties from the coupled perturbed self-consistent field equations as in density functional perturbation theory. The method is well suited to take advantage of sparse matrix algebra to achieve linear scaling complexity in the computational cost as a function of system size for sufficiently large non-metallic materials and metals at high temperatures.Comment: 21 pages, 3 figure

Giant magnetic enhancement in Fe/Pd films and its influence on the magnetic interlayer coupling

The magnetic properties of thin Pd fcc(001) films with embedded monolayers of Fe are investigated by means of first principles density functional theory. The induced spin polarization in Pd is calculated and analyzed in terms of quantum interference within the Fe/Pd/Fe bilayer system. An investigation of the magnetic enhancement effects on the spin polarization is carried out and its consequences for the magnetic interlayer coupling are discussed. In contrast to {\it e.g.} the Co/Cu fcc(001) system we find a large effect on the magnetic interlayer coupling due to magnetic enhancement in the spacer material. In the case of a single embedded Fe monolayer we find aninduced Pd magnetization decaying with distance $n$ from the magnetic layer as ~$n^{-\alpha}$ with $\alpha \approx 2.4$. For the bilayer system we find a giant magnetic enhancement (GME) that oscillates strongly due to interference effects. This results in a strongly modified magnetic interlayer coupling, both in phase and magnitude, which may not be described in the pure Ruderman-Kittel-Kasuya-Yoshida (RKKY) picture. No anti-ferromagnetic coupling was found and by comparison with magnetically constrained calculations we show that the overall ferromagnetic coupling can be understood from the strong polarization of the Pd spacer

Time-reversible Born-Oppenheimer molecular dynamics

We present a time-reversible Born-Oppenheimer molecular dynamics scheme, based on self-consistent Hartree-Fock or density functional theory, where both the nuclear and the electronic degrees of freedom are propagated in time. We show how a time-reversible adiabatic propagation of the electronic degrees of freedom is possible despite the non-linearity and incompleteness of the self-consistent field procedure. Time-reversal symmetry excludes a systematic long-term energy drift for a microcanonical ensemble and the number of self-consistency cycles can be kept low (often only 2-4 cycles per nuclear time step) thanks to a good initial guess given by the adiabatic propagation of the electronic degrees of freedom. The time-reversible Born-Oppenheimer molecular dynamics scheme therefore combines a low computational cost with a physically correct time-reversible representation of the dynamics, which preserves a detailed balance between propagation forwards and backwards in time.Comment: 4 pages, 4 figure

Ab initio linear scaling response theory: Electric polarizability by perturbed projection

A linear scaling method for calculation of the static {\em ab inito} response within self-consistent field theory is developed and applied to calculation of the static electric polarizability. The method is based on density matrix perturbation theory [Niklasson and Challacombe, cond-mat/0311591], obtaining response functions directly via a perturbative approach to spectral projection. The accuracy and efficiency of the linear scaling method is demonstrated for a series of three-dimensional water clusters at the RHF/6-31G** level of theory. Locality of the response under a global electric field perturbation is numerically demonstrated by approximate exponential decay of derivative density matrix elements.Comment: 4.25 pages in PRL format, 2 figure

Extinction calculations of multi-sphere polycrystalline graphitic clusters - A comparison with the 2175 AA peak and between a rigorous solution and discrete-dipole approximations

Certain dust particles in space are expected to appear as clusters of individual grains. The morphology of these clusters could be fractal or compact. In this paper we study the light scattering by compact and fractal polycrystalline graphitic clusters consisting of touching identical spheres. We compare three general methods for computing the extinction of the clusters in the wavelength range 0.1 - 100 micron, namely, a rigorous solution (Gerardy & Ausloos 1982) and two different discrete-dipole approximation methods -- MarCODES (Markel 1998) and DDSCAT (Draine & Flatau 1994). We consider clusters of N = 4, 7, 8, 27,32, 49, 108 and 343 particles of radii either 10 nm or 50 nm, arranged in three different geometries: open fractal (dimension D = 1.77), simple cubic and face-centred cubic. The rigorous solution shows that the extinction of the fractal clusters, with N < 50 and particle radii 10 nm, displays a peak within 2% of the location of the observed interstellar extinction peak at ~4.6 inverse micron; the smaller the cluster, the closer its peak gets to this value. By contrast, the peak in the extinction of the more compact clusters lie more than 4% from 4.6 inverse micron. At short wavelengths (0.1 - 0.5 micron), all the methods show that fractal clusters have markedly different extinction from those of non-fractal clusters. At wavelengths > 5 micron, the rigorous solution indicates that the extinction from fractal and compact clusters are of the same order of magnitude. It was only possible to compute fully converged results of the rigorous solution for the smaller clusters, due to computational limitations, however, we find that both discrete-dipole approximation methods overestimate the computed extinction of the smaller fractal clusters.Comment: Corrections added in accordance with suggestions by the referee. 12 pages, 12 figures. Accepted for publication in Astronomy & Astrophysic

Graph-based linear scaling electronic structure theory

We show how graph theory can be combined with quantum theory to calculate the electronic structure of large complex systems. The graph formalism is general and applicable to a broad range of electronic structure methods and materials, including challenging systems such as biomolecules. The methodology combines well-controlled accuracy, low computational cost, and natural low-communication parallelism. This combination addresses substantial shortcomings of linear scaling electronic structure theory, in particular with respect to quantum-based molecular dynamics simulations.Comment: 17 pages, 5 figure

Molecular-orbital-free algorithm for excited states in time-dependent perturbation theory

A non-linear conjugate gradient optimization scheme is used to obtain excitation energies within the Random Phase Approximation (RPA). The solutions to the RPA eigenvalue equation are located through a variational characterization using a modified Thouless functional, which is based upon an asymmetric Rayleigh quotient, in an orthogonalized atomic orbital representation. In this way, the computational bottleneck of calculating molecular orbitals is avoided. The variational space is reduced to the physically-relevant transitions by projections. The feasibility of an RPA implementation scaling linearly with system size, N, is investigated by monitoring convergence behavior with respect to the quality of initial guess and sensitivity to noise under thresholding, both for well- and ill-conditioned problems. The molecular- orbital-free algorithm is found to be robust and computationally efficient providing a first step toward a large-scale, reduced complexity calculation of time-dependent optical properties and linear response. The algorithm is extensible to other forms of time-dependent perturbation theory including, but not limited to, time-dependent Density Functional theory.Comment: 9 pages, 7 figure

Optomagnetic composite medium with conducting nanoelements

A new type of metal-dielectric composites has been proposed that is characterised by a resonance-like behaviour of the effective permeability in the infrared and visible spectral ranges. This material can be referred to as optomagnetic medium. The analytical formalism developed is based on solving the scattering problem for considered inclusions with impedance boundary condition, which yields the current and charge distributions within the inclusions. The presence of the effective magnetic permeability and its resonant properties lead to novel optical effects and open new possible applications.Comment: 48 pages, 13 figures. accepted to Phys. Rev. B; to appear vol. 66, 200

A Spatial Analysis of Rift Valley Fever Virus Seropositivity in Domestic Ruminants in Tanzania

Rift Valley fever (RVF) is an acute arthropod-borne viral zoonotic disease primarily occurring in Africa. Since RVF-like disease was reported in Tanzania in 1930, outbreaks of the disease have been reported mainly from the eastern ecosystem of the Great Rift Valley. This cross-sectional study was carried out to describe the variation in RVF virus (RVFV) seropositivity in domestic ruminants between selected villages in the eastern and western Rift Valley ecosystems in Tanzania, and identify potential risk factors. Three study villages were purposively selected from each of the two Rift Valley ecosystems. Serum samples from randomly selected domestic ruminants (n = 1,435) were tested for the presence of specific immunoglobulin G (IgG) and M (IgM), using RVF enzyme-linked immunosorbent assay methods. Mixed effects logistic regression modelling was used to investigate the association between potential risk factors and RVFV seropositivity. The overall RVFV seroprevalence (n = 1,435) in domestic ruminants was 25.8% and species specific seroprevalence was 29.7%, 27.7% and 22.0% in sheep (n = 148), cattle (n = 756) and goats (n = 531), respectively. The odds of seropositivity were significantly higher in animals sampled from the villages in the eastern than those in the western Rift Valley ecosystem (OR = 1.88, CI: 1.41, 2.51; p<0.001), in animals sampled from villages with soils of good than those with soils of poor water holding capacity (OR = 1.97; 95% CI: 1.58, 3.02; p< 0.001), and in animals which had been introduced than in animals born within the herd (OR = 5.08, CI: 2.74, 9.44; p< 0.001). Compared with animals aged 1-2 years, those aged 3 and 4-5 years had 3.40 (CI: 2.49, 4.64; p< 0.001) and 3.31 (CI: 2.27, 4.82, p< 0.001) times the odds of seropositivity. The findings confirm exposure to RVFV in all the study villages, but with a higher prevalence in the study villages from the eastern Rift Valley ecosystem

Применение тонкослойной хроматографии для индентификации дифенгидрамина гидрохлорида и прокаина гидрохлорида при их совместном присутствии

ТЕХНОЛОГИЯ ФАРМАЦЕВТИЧЕСКАЯДИМЕДРОЛНОВОКАИНХРОМАТОГРАФИЯ ТОНКОСЛОЙНАЯ /ИС