Linear-scaling methods theory

Other Work on Water-Related Systems. Sonoda et al.61 have simulated a time-resolved optical Kerr effect experiment. In this model, which uses molecular dynamics to represent the behaviour of the extended medium, the principle intermolecular effects are generated by the dipole-induced-dipole (DID) mechanism, but the effect of the second order molecular response is also include through terms involving the static molecular / tensor, calculated by an MP2 method. Weber et al.6S have applied ab initio linear scaling response theory to water clusters. Skaf and Vechi69 have used MP2/6-311 ++ G(d,p) calculation of the a and y tensors of water and dimethylsulfoxide (DMSO) to carry out a molecular dynamics simulation of DMSO/Water mixtures. Frediani et al.70 have used a new development of the polarizable continuum model to study the polarizability of halides at the water/air interface. 

The present review has been very selective, stressing the rationale behind density-functional methods above their applications and excluding many important topics (both theoretical and computational). The interested reader may refer to anyone of the many books [91-93] or review articles [94-101] on density-functional theory for more details. Of special importance is the extension of density-functional theory to time-dependent external potentials [102-105], as this enables the dynamical behavior of molecules, including electronic excitation, to be addressed in the context of DFT [106-108]. As they are particularly relevant to the present discussion, we cite several articles related to the formal foundations of density-functional theory [85,100,109-111], linear-scaling methods [63,112-116], exchange-correlation energy functionals [25, 117-122], and qualitative tools for describing chemical reactions [123-126,126-132]. 

Traditional wavefunction and density functional theories are applicable to large systems, including nanomaterials however, their inqilementation often involves different algorithms. These include the various linear scaling methods, hybrid (often referred to as QM/MM) mediods, and sparse matrix mediods. 

Kussmann, J., 8c Ochsenfeld, C. (2007). Linear-scaling method for calculating nuclear magnetic resonance chemical shifts using gauge-including atomic orbitals within Hartree-Fock and density-functional theory. Journal of Chemical Physics, 127, 054103. 

Basis Sets Correlation Consistent Sets Configuration Interaction Coupled-cluster Theory Density Functional Applications Density Functional Theory Applications to Transition Metal Problems G2 Theory Integrals of Electron Repulsion Integrals Overlap Linear Scaling Methods for Electronic Structure Calculations Localized MO SCF Methods Mpller-Plesset Perturbation Theory Monte Carlo Quantum Methods for Electronic Structure Numerical Hartree-Fock Methods for Molecules Pseudospectral Methods in Ab Initio Quantum Chemistry Self-consistent Reaction Field Methods Symmetry in Hartree-Fock Theory. 

Configuration Interaction Density Functional Theory (DFT), Hartree-Fock (HF), and the Self-consistent Field Linear Scaling Methods for Electronic Structure Calculations Pseudospectral Methods in Ab Initio Quantum Chemistry. 

If affordable, there is a range of very accurate coupled-cluster and symmetry-adapted perturbation theories available which can approach spectroscopic accuracy [57, 200, 201]. However, these are only applicable to the smallest alcohol cluster systems using currently available computational resources. Near-linear scaling algorithms [192] and explicit correlation methods [57] promise to extend the applicability range considerably. Furthermore, benchmark results for small systems can guide both experimentalists and theoreticians in the characterization of larger molecular assemblies. 

The concept of purification is well known in the linear-scaling literature for one-particle theories like Hartree-Fock and density functional theory, where it denotes the iterative process by which an arbitrary one-particle density matrix is projected onto an idempotent 1-RDM [2,59-61]. An RDM is said to be pure A-representable if it arises from the integration of an Al-particle density matrix T T, where T (the preimage) is an Al-particle wavefiinction [3-5]. Any idempotent 1-RDM is N-representable with a unique Slater-determinant preimage. Within the linear-scaling literature the 1-RDM may be directly computed with unconstrained optimization, where iterative purification imposes the A-representabUity conditions [59-61]. Recently, we have shown that these methods for computing the 1 -RDM directly... 

Computational methods have accompanied the development of the Polarizable Continuum Model theory throughout its history. In the building of the molecular cavity and its sampling together with the resolution of the BEM equations we nowadays have a large choice of alternative algorithms, suitable for all kinds of molecular calculations. Linear scaling both in time and space is achieved in both fields. 

In this section, we briefly discuss some of the electronic structure methods which have been used in the calculations of the PE functions which are discussed in the following sections. There are variety of ab initio electronic structure methods which can be used for the calculation of the PE surface of the electronic ground state. Most widely used are Hartree-Fock (HF) based methods. In this approach, the electronic wavefunction of a closed-shell system is described by a determinant composed of restricted one-electron spin orbitals. The unrestricted HF (UHF) method can handle also open-shell electronic systems. The limitation of HF based methods is that they do not account for electron correlation effects. For the electronic ground state of closed-shell systems, electron correlation effects can be accounted for relatively easily by second-order Mpller-Plesset perturbation theory (MP2). In modern implementations of MP2, linear scaling with the size of the system has been achieved. It is thus possible to treat quite large molecules and clusters at this level of theory. 

Nam K, JL Gao, DM York (2005) An efficient linear-scaling Ewald method for long-range electrostatic interactions in combined QM/MM calculations. J. Chem. Theory Comput. 1 (1) 2-13... 

Kudin, K., Scuseria, G. (2000). Linear-scaling density-functional theory with Gaussian orbitals and periodic boundary conditions efficient evaluation of energy and forces via the fast multipole method. Phys. Rev. B61,16443. 

We believe J.W. Linnett, in his Methuen Monograph, Wave Mechanics and Valency Theory, 1956, was the first to use elliptical projections to display the phases of the spherical harmonics on the unit sphere. We adopt a projection in which both 0 and (f) coordinates are plotted on linear scales on the minor and major axes of a 30° ellipse of eccentricity a/3/2. This cartographic device is the one proposed by Apianus in 1524, and known as the Apianus II projection. In our early work on the Spherical Shell method we called this a modified Mollweide projection, reversing the historical sequence. 

---