Biomedical subjects

Anders M N Niklasson

Publications and source records attributed to Anders M N Niklasson.

5 recordsLinked to original sources

PubMed2006

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 nonlinearity and incompleteness of the self-consistent field procedure. With a time-reversible lossless propagation the simulated dynamics is stabilized with respect to a systematic long-term energy drift and the number of self-consistency cycles can be kept low thanks to a good initial guess given from the electronic propagation. The proposed molecular dynamics scheme therefore combines a low computational cost with a physically correct time-reversible representation, which preserves a detailed balance between propagation forwards and backwards in time.

Journal Article↗

PubMed2005

Higher-order response in O(N) by perturbed projection.

Perturbed projection for linear scaling solution of the coupled-perturbed self-consistent-field equations [V. Weber, A.M.N. Niklasson, and M. Challacombe, Phys. Rev. Lett. 92, 193002 (2004)] is extended to the computation of higher-order static response properties. Although generally applicable, perturbed projection is further developed here in the context of the self-consistent first and second electric hyperpolarizabilities at the Hartree-Fock level of theory. Nonorthogonal, density-matrix analogs of Wigner's 2n+1 rule valid for linear one-electron perturbations are given up to fourth order. Linear scaling and locality of the higher-order response densities under perturbation by a global electric field are demonstrated for three-dimensional water clusters.

Journal Article↗

PubMed2005

Nonorthogonal density-matrix perturbation theory.

Recursive density-matrix perturbation theory [A.M.N. Niklasson and M. Challacombe, Phys. Rev. Lett. 92, 193001 (2004)] provides an efficient framework for the linear scaling computation of materials response properties [V. Weber, A.M.N. Niklasson, and M. Challacombe, Phys. Rev. Lett. 92, 193002 (2004)]. In this article, we generalize the density-matrix perturbation theory to include properties computed with a perturbation-dependent nonorthogonal 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 theory is developed in the context of linear scaling purification methods, which are briefly reviewed.

Journal Article↗

PubMed2004

Density matrix perturbation theory.

An orbital-free quantum perturbation theory is proposed. It gives the response of the density matrix upon variation of the Hamiltonian by quadratically convergent recursions based on perturbed projections. 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. The method allows efficient high order perturbation expansions, as demonstrated with an example involving a 10th order expansion. Density matrix analogs of Wigner's 2n+1 rule are also presented.

Journal Article↗

PubMed2004

Ab initio linear scaling response theory: electric polarizability by perturbed projection.

A linear scaling method for calculation of the static ab initio response within self-consistent field theory is developed and applied to the calculation of the static electric polarizability. The method is based on the density matrix perturbation theory [Phys. Rev. Lett. 92, 193001 (2004)]], 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. The locality of the response under a global electric field perturbation is numerically demonstrated by the approximate exponential decay of derivative density matrix elements.

Journal Article↗