Implementation and testing of Lanczos-based algorithms for Random-Phase Approximation eigenproblems

The treatment of the Random-Phase Approximation Hamiltonians, encountered in different frameworks, like time-dependent density functional theory or Bethe-Salpeter equation, is complicated by their non-Hermicity. Compared to their Hermitian Hamiltonian counterparts, computational methods for the treatment of non-Hermitian Hamiltonians are often less efficient and less stable, sometimes leading to the breakdown of the method. Recently [Gruning et al. Nano Lett. 8 (2009) 28201, we have identified that such Hamiltonians are usually pseudo-Hermitian. Exploiting this property, we have implemented an algorithm of the Lanczos type for Random-Phase Approximation Hamiltonians that benefits from the same stability and computational load as its Hermitian counterpart, and applied it to the study of the optical response of carbon nanotubes. We present here the related theoretical grounds and technical details, and study the performance of the algorithm for the calculation of the optical absorption of a molecule within the Bethe-Salpeter equation framework. (C) 2011 Elsevier B.V. All rights reserved.

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.

Electronic properties of zircon and hafnon from many-body perturbation theory

The electronic properties of zircon and hafnon, two wide-gap high-kappa materials, are investigated using many-body perturbation theory (MBPT) combined with the Wannier interpolation technique. For both materials, the calculated band structures differ from those obtained within density-functional theory and MBPT by (i) a slight displacement of the highest valence-band maximum from the Gamma point and (ii) an opening of the indirect band gap to 7.6 and 8.0 eV for zircon and hafnon, respectively. The introduction of vertex corrections in the many-body self-energy does not modify the results except for a global rigid shift of the many-body corrections.

The failure of generalized gradient approximations (GGAs) and meta-GGAs for the two-center three-electron bonds in He-2(+), (H2O)(2)(+), and (NH3)(2)(+)

The radical cations He-2(+) (H2O)(2)(+), and (NH3)(2)(+) with two-center three-electron A-A bonds are investigated at the configuration interaction (CI), accurate Kohn-Sham (KS), generalized gradient approximation (GGA), and meta-GGA levels. Assessment of seven different GGA and six meta-GGA methods shows that the A(2)(+) systems remain a difficult case for density functional theory (DFT). All methods tested consistently overestimate the stability of A(2)(+): the corresponding D-e errors decrease for more diffuse valence densities in the series He-2(+) > (H2O)(2)(+) > (NH3)(2)(+). Upon comparison to the energy terms of the accurate Kohn-Sham solutions, the approximate exchange functionals are found to be responsible for the errors of GGA-type methods, which characteristically overestimate the exchange in A(2)(+). These so-called exchange functionals implicitly use localized holes. Such localized holes do occur if there is left-right correlation, i.e., the exchange functionals then also describe nondynamical correlation. However, in the hemibonded A(2)(+) systems the typical molecular (left-right, nondynamical) correlation of the two-electron pair bond is absent. The nondynamical correlation built into the exchange functionals is then spurious and yields too low energies.

Shape corrections to exchange-correlation potentials by gradient-regulated seamless connection of model potentials for inner and outer region

Shape corrections to the standard approximate Kohn-Sham exchange-correlation (xc) potentials are considered with the aim to improve the excitation energies (especially for higher excitations) calculated with time-dependent density functional perturbation theory. A scheme of gradient-regulated connection (GRAC) of inner to outer parts of a model potential is developed. Asymptotic corrections based either on the potential of Fermi and Amaldi or van Leeuwen and Baerends (LB) are seamlessly connected to the (shifted) xc potential of Becke and Perdew (BP) with the GRAC procedure, and are employed to calculate the vertical excitation energies of the prototype molecules N-2, CO, CH2O, C2H4, C5NH5, C6H6, Li-2, Na-2, K-2. The results are compared with those of the alternative interpolation scheme of Tozer and Handy as well as with the results of the potential obtained with the statistical averaging of (model) orbital potentials. Various asymptotically corrected potentials produce high quality excitation energies, which in quite a few cases approach the benchmark accuracy of 0.1 eV for the electronic spectra. Based on these results, the potential BP-GRAC-LB is proposed for molecular response calculations, which is a smooth potential and a genuine "local" density functional with an analytical representation. (C) 2001 American Institute of Physics.

Magnetic properties of copper hexadecaphthalocyanine (F16CuPc) thin films and powders

The structural and magnetic properties of F16CuPc thin films and powder, including x-ray diffraction (XRD), superconducting quantum interference device (SQUID) magnetometry, and theoretical modelling of exchange interactions are reported. Analysis of XRD from films, with thickness ranging between 100 and 160 nm, deposited onto Kapton and a perylene-3,4,9,10-tetracarboxylic-3,4,9,10-dianhydride (PTCDA) interlayer shows that the stacking angle (defined in the text) of the film is independent of the thickness, but that the texture is modified by both film thickness and substrate chemistry. The SQUID measurements suggest that all samples are paramagnetic, a result that is confirmed by our theoretical modelling including density functional theory calculations of one-dimensional molecular chains and Green's function perturbation theory calculations for a molecular dimer. By investigating theoretically a range of different geometries, we predict that the maximum possible exchange interaction between F16CuPc molecules is twice as large as that in unfluorinated copper-phthalocyanine (CuPc). This difference arises from the smaller intermolecular spacing in F16CuPc. Our density functional theory calculation for isolated F16CuPc molecule also shows that the energy levels of Kohn-Sham orbitals are rigidly shifted similar to 1 eV lower in F16CuPc compared to CuPc without a significant modification of the intramolecular spin physics, and that therefore the two molecules provide a suitable platform for independently varying magnetism and charge transport. 

Crystal and magnetic structure of Mn3IrSi

A new ternary Ir-Mn-Si phase with stoichiometry Mn₃IrSi has been synthesized and found to crystallize in the cubic AlAu₄-type structure, space group P2₁3 with Z=4, which is an ordered form of the beta-Mn structure. The unit cell dimension was determined by x-ray powder diffraction to a=6.4973(3) Angstrom. In addition to the crystal structure, we have determined the magnetic structure and properties using superconducting quantum interference device magnetometry and Rietveld refinements of neutron powder diffraction data. A complex noncollinear magnetic structure is found, with magnetic moments of 2.97(4)u(B) at 10 K only on the Mn atoms. The crystal structure consists of a triangular network built up by Mn atoms, on which the moments are rotated 120degrees around the triangle axes. The magnetic unit cell is the same as the crystallographic and carries no net magnetic moment. The Neel temperature was determined to be 210 K. A first-principles study, based on density functional theory in a general noncollinear formulation, reproduces the experimental results with good agreement. The observed magnetic structure is argued to be the result of frustration of antiferromagnetic couplings by the triangular geometry.

Combined studies of DFT atomistic modelling and in situ FTIR spectroscopy on surface oxidants and CO oxidation at Ru electrodes

We report the combined studies of density functional theory (DFT) calculations and electrochemical in situ FTIR spectroscopy on surface oxidants and mechanisms of CO oxidation at the Ru(0001) electrodes. It is shown that CO can co-adsorb with both O and OH species at lower potential region where a low coverage of the (2 x 2)-O/OH adlayer formed; the oxidation of CO adsorbates takes place at higher potentials where a high coverage of the (1 x 1)-O/OH adlayer formed. Surface O species are not the active oxidants under all coverages studied, due to the high reaction barriers between CO and O (>1 eV). However, surface OH species with higher coverage are identified as the active oxidants, and CO oxidation takes place via a two-steps' mechanism of CO + 3OH -> COOH + 2OH -> CO2 + H2O + OH, in which three nearby OH species are involved in the CO2 formation: CO reacts with OH, forming COOH; COOH then transfers the H to a nearby OH to form H2O and CO2, at the same time, another H in the H2O transfers to a nearby OH to form a weak adsorbed H2O and a new OH. The reaction barrier of these processes is reduced significantly to around 0.50 eV. These new results not only provide an insight into surface active oxidants on Ru, which is directly relevant to fuel cell catalysis, but also reveals the extra complexity of catalytic reactions taking place at solid/liquid electrochemical interface in comparison to the relatively simpler ones at solid/gas phase. 

Chiral thiouronium salts: synthesis, characterisation and application in NMR enantio-discrimination of chiral oxoanions

Chiral thioureas and functionalised chiral thiouronium salts were synthesised starting from the relatively cheap and easily available chiral amines: (S)-methylbenzylamine and rosin-derived (+)-dehydroabietylamine. The introduction of a delocalised positive charge to the thiourea functionality, by an alkylation reaction at the sulfur atom, enables dynamic rotameric processes: hindered rotations about the delocalised CN and CS bonds. Hence, four different rotamers/isomers may be recognised: syn-syn, syn-anti, anti-syn and anti-anti. Extensive H-1 and C-13 NMR studies have shown that in hydrogen-bond acceptor solvents, such as perdeuteriated dimethyl sulfoxide, the syn-syn conformation is preferable. On the other hand, when using non-polar solvents, such as CDCl3, the mixture of syn-syn and syn-anti isomers is detectable, with an excess of the latter. Apart from this, in the case of S-butyl-N,N'-bis(dehydroabietyl)thiouronium ethanoate in CDCl3, the H-1 NMR spectrum revealed that strong bifurcated hydrogen bonding between the anion and the cation causes global rigidity without signs of hindered rotamerism observable on the NMR time scale. This suggested that these new salts might be used as NMR discriminating agents for chiral oxoanions, and are indeed more effective than their archetypal guanidinium analogues or the neutral thioureas. The best results in recognition of a model substrate, mandelate, were obtained with S-butyl-N,N'-bis(dehydroabietyl) thiouronium bistriflamide. It was confirmed that the chiral recognition occurred not only for carboxylates but also for sulfonates and phosphonates. Further H-1 NMR studies confirmed a 1 : 1 recognition mode between the chiral agent (host) and the substrate (guest); binding constants were determined by H-1 NMR titrations in solutions of DMSO-d(6) in CDCl3. It was also found that the anion of the thiouronium salt had a significant influence on the recognition process: anions with poor hydrogen-bond acceptor abilities led to the best discrimination. The presence of host-guest hydrogen bonding was confirmed in the X-ray crystal structure of S-butyl-N,N'-bis(dehydroabietyl)thiouronium bromide and by computational studies (density functional theory).

Surface energy and surface proton order of the ice Ih basal and prism surfaces

Density-functional theory (DFT) is used to examine the basal and prism surfaces of ice Ih. Similar surface energies are obtained for the two surfaces; however, in each case a strong dependence of the surface energy on surface proton order is identified. This dependence, which can be as much as 50% of the absolute surface energy, is significantly larger than the bulk dependence (< 1%) on proton order, suggesting that the thermodynamic ground state of the ice surface will remain proton ordered well above the bulk order-disorder temperature of about 72 K. On the basal surface this suggestion is supported by Monte Carlo simulations with an empirical potential and solution of a 2D Ising model with nearest neighbor interactions taken from DFT. Order parameters that define the surface energy of each surface in terms of nearest neighbor interactions between dangling OH bonds (those which point out of the surface into vacuum) have been identified and are discussed. Overall, these results suggest that proton order-disorder effects have a profound impact on the stability of ice surfaces and will most likely have an effect on ice surface reactivity as well as ice crystal growth and morphology. S Supplementary data are available from stacks.iop.org/JPhysCM/22/074209/mmedia

The phase diagram of water at negative pressures: Virtual ices

The phase diagram of water at negative pressures as obtained from computer simulations for two models of water, TIP4P/2005 and TIP5P is presented. Several solid structures with lower densities than ice Ih, so-called virtual ices, were considered as possible candidates to occupy the negative pressure region of the phase diagram of water. In particular the empty hydrate structures sI, sII, and sH and another, recently proposed, low-density ice structure. The relative stabilities of these structures at 0 K was determined using empirical water potentials and density functional theory calculations. By performing free energy calculations and Gibbs-Duhem integration the phase diagram of TIP4P/2005 was determined at negative pressures. The empty hydrates sII and sH appear to be the stable solid phases of water at negative pressures. The phase boundary between ice Ih and sII clathrate occurs at moderate negative pressures, while at large negative pressures sH becomes the most stable phase. This behavior is in reasonable agreement with what is observed in density functional theory calculations.

Palladium-catalyzed liquid-phase hydrogenation/hydrogenolysis of disulfides

For the first time, the hydrogenation/hydrogenolysis of a range of disulfides has been achieved over a supported palladium catalyst using hydrogen under relatively benign conditions. These unexpected results demonstrate that it is possible to avoid the poisoning of the catalyst by either the nitrogen-containing groups or the sulfur species, allowing both efficient reaction and recycling of the catalyst under the proper conditions (e.g., at low temperatures). A slight loss in activity was found on recycling; however, the catalyst activity can be recovered using hydrogen pretreatment. The reaction mechanism for the hydrogenolysis and hydrogenation of ortho-, meta-, and para-dinitrodiphenyldisulfide to the corresponding aminothiophenol has been elucidated. Density functional theory calculations were used to investigate the adsorption mode of the dinitrodiphenyldisulfides; a clear dependence on adsorption geometry was found regarding whether the molecule is cleaved at the S-S bond before the reduction of the nitro group or vice versa. This study demonstrates the versatility of these catalysts for the hydrogenation/hydrogenolysis of sulfur-containing molecules, which normally are considered poisons, and will extend their use to a new family of substrates. (C) 2007 Elsevier Inc. All rights reserved.

The Role of Arg228 in the phosphorylation of galactokinase: The mechanism of GHMP kinases by QM/MM studies

GHMP kinases are a group of structurally-related small molecule kinases. They have been found in all kingdoms of life and are mostly responsible for catalysing the ATP-dependent phosphorylation of intermediary metabolites. Although the GHMP kinases are of clinical, pharmaceutical and biotechnological importance, the mechanism of GHMP-kinases is controversial. A catalytic base mechanism was suggested for mevalonate kinase that has a structural feature of the ?-phosphate of ATP close to an aspartate residue; however, for one GHMP member, homoserine kinase, where the residue acting as general base is absent, a direct phosphorylation mechanism was suggested. Furthermore, it has been proposed by some authors that all the GHMP kinases function via the direct phosphorylation mechanism. This controversy in mechanism has limited our ability to exploit these enzymes as drug targets and in biotechnology. Here the phosphorylation reaction mechanism of the human galactokinase, a member of GHMP kinase was investigated using molecular dynamics simulations and density functional theory-based QM/MM calculations (B3LYP-D/AMBER99). The reaction coordinates were localized by potential energy scan using adiabatic mapping method. Our results indicate that a highly conserved Glu174 captures Arg105 to the proximity of the a-phosphate of ATP forming a H-bond network, therefore the mobility of ATP in the large oxyanion hole is restricted. Arg228 functions to stabilize the negative charge developed at the ß,?-bridging oxygen of the ATP during bond cleavage. The reaction occurs via direct phosphorylation mechanism and the Asp186 in proximity of ATP does not directly participate in the reaction pathway. Since Arg228 is not conserved among GHMP kinases, reagents which form interactions with Arg228, and therefore can interrupt its function in phosphorylation may be developed into potential selective inhibitors for galactokinase.

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]

Performance of Nonlocal Optics When Applied to Plasmonic Nanostructures

Semiclassical nonlocal optics based on the hydrodynamic description of conduction electrons might be an adequate tool to study complex phenomena in the emerging field of nanoplasmonics. With the aim of confirming this idea, we obtain the local and nonlocal optical absorption spectra in a model nanoplasmonic device in which there are spatial gaps between the components at nanometric and subnanometric scales. After a comparison against time-dependent density functional calculations, we conclude that hydrodynamic nonlocal optics provides absorption spectra exhibiting qualitative agreement but not quantitative accuracy. This lack of accuracy, which is manifest even in the limit where induced electric currents are not established between the constituents of the device, is mainly due to the poor description of induced electron densities.