Tracing potential energy surfaces of electronic excitations via their transition origins: application to Oxirane

We show that the transition origins of electronic excitations identified by quantified natural transition orbital (QNTO) analysis can be employed to connect potential energy surfaces (PESs) according to their character across a widerange of molecular geometries. This is achieved by locating the switching of transition origins of adiabatic potential surfaces as the geometry changes. The transition vectors for analysing transition origins are provided by linear response time-dependent density functional theory (TDDFT) calculations under the Tamm-Dancoff approximation. We study the photochemical CO ring opening of oxirane as an example and show that the results corroborate the traditional Gomer-Noyes mechanism derived experimentally. The knowledge of specific states for the reaction also agrees well with that given by previous theoretical work using TDDFT surface-hopping dynamics that was validated by high-quality quantum Monte Carlo calculations. We also show that QNTO can be useful for considerably larger and more complex systems: by projecting the excitations to those of a reference oxirane molecule, the approach is able to identify and analyse specific excitations of a trans-2,3-diphenyloxirane molecule.

Ligand Discrimination in Myoglobin from Linear-Scaling DFT+U

Myoglobin modulates the binding of diatomic molecules to its heme group via hydrogen-bonding and steric interactions with neighboring residues, and is an important benchmark for computational studies of biomolecules. We have performed calculations on the heme binding site and a significant proportion of the protein environment (more than 1000 atoms) using linear-scaling density functional theory and the DFT+U method to correct for self-interaction errors associated with localized 3d states. We confirm both the hydrogen-bonding nature of the discrimination effect (3.6 kcal/mol) and assumptions that the relative strain energy stored in the protein is low (less than 1 kcal/mol). Our calculations significantly widen the scope for tackling problems in drug design and enzymology, especially in cases where electron localization, allostery or long-ranged polarization influence ligand binding and reaction.

Generalized Wannier functions: a comparison of molecular electric dipole polarizabilities

Localized Wannier functions provide an efficient and intuitive means by which to compute dielectric properties from first principles. They are most commonly constructed in a post-processing step, following total-energy minimization. Nonorthogonal generalized Wannier functions (NGWFs) [Skylaris et al., Phys. Rev. B 66, 035119 11 (2002); Skylaris et al., J. Chem. Phys. 122, 084119 (2005)] may also be optimized in situ, in the process of solving for the ground-state density. We explore the relationship between NGWFs and orthonormal, maximally localized Wannier functions (MLWFs) [Marzari and Vanderbilt, Phys. Rev. B 56, 12847 (1997); Souza, Marzari, and Vanderbilt, ibid. 65, 035109 (2001)], demonstrating that NGWFs may be used to compute electric dipole polarizabilities efficiently, with no necessity for post-processing optimization, and with an accuracy comparable to MLWFs.

Linear-scaling DFT+U with full local orbital optimization

We present an approach to the DFT+U method (Density Functional Theory + Hubbard model) within which the computational effort for calculation of ground state energies and forces scales linearly with system size. We employ a formulation of the Hubbard model using nonorthogonal projector functions to define the localized subspaces, and apply it to a local-orbital DFT method including in situ orbital optimization. The resulting approach thus combines linear-scaling and systematic variational convergence. We demonstrate the scaling of the method by applying it to nickel oxide nano-clusters with sizes exceeding 7,000 atoms.

Stress Development and Impurity Segregation during Oxidation of the Si(100) Surface

We have studied the segregation of P and B impurities during oxidation of the Si(100) surface by means of combined static and dynamical first-principles simulations based on density functional theory. In the bare surface, dopants segregate to chemically stable surface sites or to locally compressed subsurface sites. Surface oxidation is accompanied by development of tensile surface stress up to 2.9 N/m at a coverage of 1.5 monolayers of oxygen and by formation of oxidised Si species with charges increasing approximately linearly with the number of neighbouring oxygen atoms. Substitutional P and B defects are energetically unstable within the native oxide layer, and are preferentially located at or beneath the Si/SiOx interface. Consistently, first-principles molecular dynamics simulations of native oxide formation on doped surfaces reveal that dopants avoid the formation of P-O and B-O bonds, suggesting a surface oxidation mechanism whereby impurities remain trapped at the Si/SiOx interface. This seems to preclude a direct influence of impurities on the surface electrostatics and, hence, on the interactions with an external environment.

Development of a Classical Force Field for the Oxidised Si Surface: Application to Hydrophilic Wafer Bonding

We have developed a classical two- and three-body interaction potential to simulate the hydroxylated, natively oxidised Si surface in contact with water solutions, based on the combination and extension of the Stillinger-Weber potential and of a potential originally developed to simulate SiO2 polymorphs. The potential parameters are chosen to reproduce the structure, charge distribution, tensile surface stress and interactions with single water molecules of a natively oxidised Si surface model previously obtained by means of accurate density functional theory simulations. We have applied the potential to the case of hydrophilic silicon wafer bonding at room temperature, revealing maximum room temperature work of adhesion values for natively oxidised and amorphous silica surfaces of 97 mJ/m2 and 90mJ/m2, respectively, at a water adsorption coverage of approximately 1 monolayer. The difference arises from the stronger interaction of the natively oxidised surface with liquid water, resulting in a higher heat of immersion (203 mJ/m2 vs. 166 mJ/m2), and may be explained in terms of the more pronounced water structuring close to the surface in alternating layers of larger and smaller density with respect to the liquid bulk. The computed force-displacement bonding curves may be a useful input for cohesive zone models where both the topographic details of the surfaces and the dependence of the attractive force on the initial surface separation and wetting can be taken into account.