A multiconfigurational time-dependent Hartree-Fock method for excited electronic states. I. General formalism and application to open-shell states

The solution of the time-dependent Schrodinger equation for systems of interacting electrons is generally a prohibitive task, for which approximate methods are necessary. Popular approaches, such as the time-dependent Hartree-Fock (TDHF) approximation and time-dependent density functional theory (TDDFT), are essentially single-configurational schemes. TDHF is by construction incapable of fully accounting for the excited character of the electronic states involved in many physical processes of interest; TDDFT, although exact in principle, is limited by the currently available exchange-correlation functionals. On the other hand, multiconfigurational methods, such as the multiconfigurational time-dependent Hartree-Fock (MCTDHF) approach, provide an accurate description of the excited states and can be systematically improved. However, the computational cost becomes prohibitive as the number of degrees of freedom increases, and thus, at present, the MCTDHF method is only practical for few-electron systems. In this work, we propose an alternative approach which effectively establishes a compromise between efficiency and accuracy, by retaining the smallest possible number of configurations that catches the essential features of the electronic wavefunction. Based on a time-dependent variational principle, we derive the MCTDHF working equation for a multiconfigurational expansion with fixed coefficients and specialise to the case of general open-shell states, which are relevant for many physical processes of interest. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3600397]

A multiconfigurational time-dependent Hartree-Fock method for excited electronic states. II. Coulomb interaction effects in single conjugated polymer chains

Conjugated polymers have attracted considerable attention in the last few decades due to their potential for optoelectronic applications. A key step that needs optimisation is charge carrier separation following photoexcitation. To understand better the dynamics of the exciton prior to charge separation, we have performed simulations of the formation and dynamics of localised excitations in single conjugated polymer strands. We use a nonadiabatic molecular dynamics method which allows for the coupled evolution of the nuclear degrees of freedom and of multiconfigurational electronic wavefunctions. We show the relaxation of electron-hole pairs to form excitons and oppositely charged polaron pairs and discuss the modifications to the relaxation process predicted by the inclusion of the Coulomb interaction between the carriers. The issue of charge photogeneration in conjugated polymers in dilute solution is also addressed. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3600404]

On the origin of the forward peak and backward oscillations in the F + H2(v = 0) to HF(v' = 2) + H reaction

We present a semiclassical complex angular momentum (CAM) analysis of the forward scattering peak which occurs at a translational collision energy around 32 meV in the quantum mechanical calculations for the F + H2(v = 0, j = 0) ? HF(v' = 2, j' = 0) + H reaction on the Stark–Werner potential energy surface. The semiclassical CAM theory is modified to cover the forward and backward scattering angles. The peak is shown to result from constructive/destructive interference of the two Regge states associated with two resonances, one in the transition state region and the other in the exit channel van der Waals well. In addition, we demonstrate that the oscillations in the energy dependence of the backward differential cross section are caused by the interference between the direct backward scattering and the decay of the two resonance complexes returning to the backward direction after one full rotation.

Quasiparticle calculations of the electronic properties of ZrO2 and HfO2 polymorphs and their interface with Si

Quasiparticle calculations are performed to investigate the electronic band structures of various polymorphs of Hf and Zr oxides. The corrections with respect to density-functional-theory results are found to depend only weakly on the crystal structure. Based on these bulk calculations as well as those for bulk Si, the effect of quasiparticle corrections is also investigated for the band offsets at the interface between these oxides and Si assuming that the lineup of the potential at the interface is reproduced correctly within density-functional theory. On the one hand, the valence-band offsets are practically unchanged with a correction of a few tenths of electron volts. On the other hand, conduction-band offsets are raised by 1.3-1.5 eV. When applied to existing calculations for the offsets at the density-functional-theory level, our quasiparticle corrections provide results in good agreement with the experiment.

SrTiO3(001)(2x1) reconstructions: First-principles calculations of surface energy and atomic structure compared with scanning tunneling microscopy images

(1x1) and (2x1) reconstructions of the (001) SrTiO3 surface were studied using the first-principles full-potential linear muffin-tin orbital method. Surface energies were calculated as a function of TiO2 chemical potential, oxygen partial pressure and temperature. The (1x1) unreconstructed surfaces were found to be energetically stable for many of the conditions considered. Under conditions of very low oxygen partial pressure the (2x1) Ti2O3 reconstruction [Martin R. Castell, Surf. Sci. 505, 1 (2002)] is stable. The question as to why STM images of the (1x1) surfaces have not been obtained was addressed by calculating charge densities for each surface. These suggest that the (2x1) reconstructions would be easier to image than the (1x1) surfaces. The possibility that the presence of oxygen vacancies would destabilise the (1x1) surfaces was also investigated. If the (1x1) surfaces are unstable then there exists the further possibility that the (2x1) DL-TiO2 reconstruction [Natasha Erdman Nature (London) 419, 55 (2002)] is stable in a TiO2-rich environment and for p(O2)>10(-18) atm.

Free energy and molecular dynamics calculations for the cubic-tetragonal phase transition in zirconia

The high-temperature cubic-tetragonal phase transition of pure stoichiometric zirconia is studied by molecular dynamics (MD) simulations and within the framework of the Landau theory of phase transformations. The interatomic forces are calculated using an empirical, self-consistent, orthogonal tight-binding model, which includes atomic polarizabilities up to the quadrupolar level. A first set of standard MD calculations shows that, on increasing temperature, one particular vibrational frequency softens. The temperature evolution of the free-energy surfaces around the phase transition is then studied with a second set of calculations. These combine the thermodynamic integration technique with constrained MD simulations. The results seem to support the thesis of a second-order phase transition but with unusual, very anharmonic behavior above the transition temperature.

Density functional calculations for polymers and clusters - progress and limitations

The applicability of density functional (DF) methods has progressed greatly since the first workshop of this series ten years ago. Applications that show both the successes and the limitations can be found in the fields of: (a) the structures of the isomers of atomic clusters. and (b) the structure of organic molecules and polymers, and their reactions with additional molecules. We shall review some of the results and the lessons to be learned from them. (C) 2001 Elsevier Science B.V. All rights reserved.

Optimal basis set for electronic structure calculations in periodic systems

An efficient method for calculating the electronic structure of systems that need a very fine sampling of the Brillouin zone is presented. The method is based on the variational optimization of a single (i.e., common to all points in the Brillouin zone) basis set for the expansion of the electronic orbitals. Considerations from k.p-approximation theory help to understand the efficiency of the method. The accuracy and the convergence properties of the method as a function of the optimal basis set size are analyzed for a test calculation on a 16-atom Na supercell.

Structural and chemical embrittlement of grain boundaries by impurities: A general theory and first-principles calculations for copper

First-principles calculations of the Sigma 5(310)[001] symmetric tilt grain boundary in Cu with Bi, Na, and Ag substitutional impurities provide evidence that in the phenomenon of Bi embrittlement of Cu grain boundaries electronic effects do not play a major role; on the contrary, the embrittlement is mostly a structural or "size" effect. Na is predicted to be nearly as good an embrittler as Bi, whereas Ag does not embrittle the boundary in agreement with experiment. While we reject the prevailing view that "electronic" effects (i.e., charge transfer) are responsible for embrittlement, we do not exclude the role of chemistry. However, numerical results show a striking equivalence between the alkali metal Na and the semimetal Bi, small differences being accounted for by their contrasting "size" and "softness" (defined here). In order to separate structural and chemical effects unambiguously if not uniquely, we model the embrittlement process by taking the system of grain boundary and free surfaces through a sequence of precisely defined gedanken processes; each of these representing a putative mechanism. We thereby identify three mechanisms of embrittlement by substitutional impurities, two of which survive in the case of embrittlement or cohesion enhancement by interstitials. Two of the three are purely structural and the third contains both structural and chemical elements that by their very nature cannot be further unraveled. We are able to take the systems we study through each of these stages by explicit computer simulations and assess the contribution of each to the net reduction in intergranular cohesion. The conclusion we reach is that embrittlement by both Bi and Na is almost exclusively structural in origin; that is, the embrittlement is a size effect.