Partitioning scheme for density functional calculations of extended systems

Frozen density embedding with non-integer subsystems' particle numbers

We extend the frozen density embedding theory to non-integer subsystems' particles numbers. Different features of this formulation are discussed, with special concern for approximate embedding calculations. In particular, we highlight the relation between the non-integer particle-number partition scheme and the resulting embedding errors. Finally, we provide a discussion of the implications of the present theory for the derivative discontinuity issue and the calculation of chemical reactivity descriptors.

Density functional calculations of extended, periodic systems using Coulomb corrected molecular fractionation with conjugated caps method (CC-MFCC)

A consistent DFT based formulation of the order-N molecular fractionation with conjugated caps method in which a molecular system is calculated considering a set of finite fragments, is proposed. Here we apply the method and test its performance on a periodic metal–organic framework system.

On the accuracy of frozen density embedding calculations with hybrid and orbital-dependent functionals for non-bonded interaction energies

We analyze the accuracy of the frozen density embedding (FDE) method, with hybrid and orbital-dependent exchange-correlation functionals, for the calculation of the total interaction energies of weakly interacting systems. Our investigation is motivated by the fact that these approaches require, in addition to the non-additive kinetic energy approximation, also approximate non-additive exact-exchange energies. Despite this further approximation, we find that the hybrid/orbital-dependent FDE approaches can reproduce the total energies with the same accuracy (about 1 mHa) as the one of conventional semi-local functionals. In many cases, thanks to error cancellation effects, hybrid/orbital-dependent approaches yield even the smallest error. A detailed energy-decomposition investigation is presented. Finally, the Becke-exchange functional is found to reproduce accurately the non-additive exact-exchange energies also for non-equilibrium geometries. These performances are rationalized in terms of a reduced-gradient decomposition of the non-additive exchange energy.

Characterization of hydrogen dissociation over aluminium-doped zinc oxide using an efficient massively parallel framework for QM/MM calculations

A task-farm parallelization framework has been implemented in the ChemShell computational chemistry environment to provide a facility for parallelizing common chemical calculations, including finite-difference Hessian evaluation, the nudged elastic band method for reaction path optimization, and population-based methods for global optimization. The optimization methods are provided by a parallel interface to the DL-FIND optimization library. As ChemShell can already exploit parallel external programs for energy and gradient evaluations, the new methods result in a two-level approach to parallelization that gives significantly improved performance for massively parallel calculations. For typical systems, speed-up factors of five to eight times have been observed compared with non-task-farmed calculations. The task-farming version of ChemShell has been used to study the heterolytic dissociation of a hydrogen molecule over a polar oxygen-terminated surface of aluminium-doped zinc oxide using an embedded cluster hybrid QM/MM approach. We calculate a 42 kcal mol −1 heat of reaction and a 30 kcal mol −1 activation energy, which is equivalent to a high backward reaction barrier of 72 kcal mol −1 per H 2 molecule, in close agreement with temperature programmed desorption experiments. The dissociation path includes a stable intermediate comprising a hydride ion in an oxygen vacancy and physisorbed atomic hydrogen.

Modelling components of future molecular devices

We discuss challenges involved in modelling different components of molecular devices and give several examples that demonstrate how computer modelling evolved over the last few years to become a comprehensive tool for designing molecules, predicting their adsorption and diffusion at surfaces, simulating atomic force microscopy imaging and manipulation of atoms and molecules at insulating surfaces and studying electron conduction in prototype molecular devices. We describe some of the computational techniques used for modelling adsorption, diffusion, imaging and manipulation of organic molecules at surfaces and challenges pertaining to these studies, give several examples of applications and discuss further prospects for theoretical modelling of complex organic molecules at surfaces.

The role of the exchange in the embedding electrostatic potential for the fragment molecular orbital method

We have examined the role of the exchange in describing the electrostatic potential in the fragment molecular orbital method and showed that it should be included in the total Fock matrix to obtain an accurate one-electron spectrum; however, adding it to the Fock matrices of individual fragments and pairs leads to very large errors. For the error analysis we have used the virial theorem; numerical tests have been performed for solvated phenol at the Hartree-Fock level with the 6-31G( *) and 6-311G( * *) basis sets.