Non-Orthogonal Density Matrix Perturbation Theory

Density matrix perturbation theory [Phys. Rev. Lett. Vol. 92, 193001 (2004)] provides an efficient framework for the linear scaling computation of response properties [Phys. Rev. Lett. Vol. 92, 193002 (2004)]. In this article, we generalize density matrix perturbation theory to include properties computed with a perturbation dependent non-orthogonal basis. Such properties include analytic derivatives of the energy with respect to nuclear displacement, as well as magnetic response computed with a field dependent basis. The non-orthogonal density matrix perturbation theory is developed in the context of recursive purification methods, which are briefly reviewed.Comment: 8 pages, 2 figure

Expansion algorithm for the density matrix

A purification algorithm for expanding the single-particle density matrix in terms of the Hamiltonian operator is proposed. The scheme works with a predefined occupation and requires less than half the number of matrix-matrix multiplications compared to existing methods at low (90%) occupancy. The expansion can be used with a fixed chemical potential in which case it is an asymmetric generalization of and a substantial improvement over grand canonical McWeeny purification. It is shown that the computational complexity, measured as number of matrix multiplications, essentially is independent of system size even for metallic materials with a vanishing band gap.Comment: 5 pages, 4 figures, to appear in Phys. Rev.

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

Density Matrix Perturbation Theory

An expansion method for perturbation of the zero temperature grand canonical density matrix is introduced. The method achieves quadratically convergent recursions that yield the response of the zero temperature density matrix upon variation of the Hamiltonian. The technique allows treatment of embedded quantum subsystems with a computational cost scaling linearly with the size of the perturbed region, O(N_pert.), and as O(1) with the total system size. It also allows direct computation of the density matrix response functions to any order with linear scaling effort. Energy expressions to 4th order based on only first and second order density matrix response are given.Comment: 4 pages, 2 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

Carer's perception of and reaction to reassurance seeking in obsessive compulsive disorder

The purpose of the present study is to explore the experience of being asked for reassurance from the perspective of carers of OCD sufferers, and to examine its relationship to sufferers’ reassurance seeking by a direct comparison with data obtained from the person they normally offer reassurance to. Forty-two individuals with OCD and their carers completed alternate versions of the Reassurance Seeking Questionnaire. Result suggest that carers report most commonly providing reassurance when asked to do so, and the frequency of their reassurance provision is associated with how carefully sufferers seek reassurance, rather than their OCD symptom severity. The carer's perspectives on the impact of reassurance provision was accurate; both sufferers and carers perceive that reassurance works only temporarily, but even if the anxiety-relieving effect of reassurance decreases in the medium term, it is likely to be perceived as beneficial because carers accurately perceived that sufferers would feel much worse if they refuse to provide reassurance. The present study is the first to quantitatively investigate carer’s experiences of reassurance provision, and elucidate why carers feel the need to provide it

O(N) methods in electronic structure calculations

Linear scaling methods, or O(N) methods, have computational and memory requirements which scale linearly with the number of atoms in the system, N, in contrast to standard approaches which scale with the cube of the number of atoms. These methods, which rely on the short-ranged nature of electronic structure, will allow accurate, ab initio simulations of systems of unprecedented size. The theory behind the locality of electronic structure is described and related to physical properties of systems to be modelled, along with a survey of recent developments in real-space methods which are important for efficient use of high performance computers. The linear scaling methods proposed to date can be divided into seven different areas, and the applicability, efficiency and advantages of the methods proposed in these areas is then discussed. The applications of linear scaling methods, as well as the implementations available as computer programs, are considered. Finally, the prospects for and the challenges facing linear scaling methods are discussed.Comment: 85 pages, 15 figures, 488 references. Resubmitted to Rep. Prog. Phys (small changes

Evaluation of the efficacy of Alpron disinfectant for dental unit water lines

AIMS: To assess the efficacy of a disinfectant, Alpron, for controlling microbial contamination within dental unit water lines. METHODS: The microbiological quality of water emerging from the triple syringe, high speed handpiece, cup filler and surgery hand wash basin from six dental units was assessed for microbiological total viable counts at 22 degrees C and 37 degrees C before and after treatment with Alpron solutions. RESULTS: The study found that the use of Alpron disinfectant solutions could reduce microbial counts in dental unit water lines to similar levels for drinking water. This effect was maintained in all units for up to six weeks following one course of treatment. In four out of six units the low microbial counts were maintained for 13 weeks. CONCLUSIONS: Disinfectants may have a short term role to play in controlling microbial contamination of dental unit water lines to drinking water quality. However, in the longer term attention must be paid to redesigning dental units to discourage the build up of microbial biofilms

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

An attenuated strain of Bacillus anthracis (CDC 684) has a large chromosomal inversion and altered growth kinetics

BackgroundAn isolate originally labeled Bacillus megaterium CDC 684 was found to contain both pXO1 and pXO2, was non-hemolytic, sensitive to gamma-phage, and produced both the protective antigen and the poly-D-glutamic acid capsule. These phenotypes prompted Ezzell et al., (J. Clin. Microbiol. 28:223) to reclassify this isolate to Bacillus anthracis in 1990.ResultsWe demonstrate that despite these B. anthracis features, the isolate is severely attenuated in a guinea pig model. This prompted whole genome sequencing and closure. The comparative analysis of CDC 684 to other sequenced B. anthracis isolates and further analysis reveals: a) CDC 684 is a close relative of a virulent strain, Vollum A0488; b) CDC 684 defines a new B. anthracis lineage (at least 51 SNPs) that includes 15 other isolates; c) the genome of CDC 684 contains a large chromosomal inversion that spans 3.3 Mbp; d) this inversion has caused a displacement of the usual spatial orientation of the origin of replication (ori) to the termination of replication (ter) from 180° in wild-type B. anthracis to 120° in CDC 684 and e) this isolate also has altered growth kinetics in liquid media.ConclusionsWe propose two alternative hypotheses explaining the attenuated phenotype of this isolate. Hypothesis 1 suggests that the skewed ori/ter relationship in CDC 684 has altered its DNA replication and/or transcriptome processes resulting in altered growth kinetics and virulence capacity. Hypothesis 2 suggests that one or more of the single nucleotide polymorphisms in CDC 684 has altered the expression of a regulatory element or other genes necessary for virulence