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.

Geographies of Communality, Colonialism, and Capitalism: Ecology and the World-System

Drawing upon recent reworkings of world systems theory and Marx’s concept of metabolic rift, this paper attempts to ground early nineteenth-century Ireland more clearly within these metanarratives, which take the historical-ecological dynamics of the development of capitalism as their point of departure. In order to unravel the socio-spatial complexities of Irish agricultural production throughout this time, attention must be given to the prevalence of customary legal tenure, institutions of communal governance, and their interaction with the colonial apparatus, as an essential feature of Ireland’s historical geography often neglected by famine scholars. This spatially differentiated legacy of communality, embedded within a country-wide system of colonial rent, and burgeoning capitalist system of global trade, gave rise to profound regional differentiations and ecological contradictions, which became central to the distribution of distress during the Great Famine (1845-1852). Contrary to accounts which depict it as a case of discrete transition from feudalism to capitalism, Ireland’s pre-famine ecology must be understood through an analysis which emphasises these socio-spatial complexities. Consequently, this structure must be conceptualised as one in which communality, colonialism, and capitalism interact dynamically, and in varying stages of development and devolution, according to space and time.

Dynamics of the entanglement spectrum in spin chains

We study the dynamics of the entanglement spectrum, that is the time evolution of the eigenvalues of the reduced density matrices after a bipartition of a one-dimensional spin chain. Starting from the ground state of an initial Hamiltonian, the state of the system is evolved in time with a new Hamiltonian. We consider both instantaneous and quasi adiabatic quenches of the system Hamiltonian across a quantum phase transition. We analyse the Ising model that can be exactly solved and the XXZ for which we employ the time-dependent density matrix renormalisation group algorithm. Our results show once more a connection between the Schmidt gap, i.e. the difference of the two largest eigenvalues of the reduced density matrix and order parameters, in this case the spontaneous magnetisation.

Appropriate complexity for the prediction of coastal and estuarine geomorphic behaviour at decadal to centennial scales

Coastal and estuarine landforms provide a physical template that not only accommodates diverse ecosystem functions and human activities, but also mediates flood and erosion risks that are expected to increase with climate change. In this paper, we explore some of the issues associated with the conceptualisation and modelling of coastal morphological change at time and space scales relevant to managers and policy makers. Firstly, we revisit the question of how to define the most appropriate scales at which to seek quantitative predictions of landform change within an age defined by human interference with natural sediment systems and by the prospect of significant changes in climate and ocean forcing. Secondly, we consider the theoretical bases and conceptual frameworks for determining which processes are most important at a given scale of interest and the related problem of how to translate this understanding into models that are computationally feasible, retain a sound physical basis and demonstrate useful predictive skill. In particular, we explore the limitations of a primary scale approach and the extent to which these can be resolved with reference to the concept of the coastal tract and application of systems theory. Thirdly, we consider the importance of different styles of landform change and the need to resolve not only incremental evolution of morphology but also changes in the qualitative dynamics of a system and/or its gross morphological configuration. The extreme complexity and spatially distributed nature of landform systems means that quantitative prediction of future changes must necessarily be approached through mechanistic modelling of some form or another. Geomorphology has increasingly embraced so-called ‘reduced complexity’ models as a means of moving from an essentially reductionist focus on the mechanics of sediment transport towards a more synthesist view of landform evolution. However, there is little consensus on exactly what constitutes a reduced complexity model and the term itself is both misleading and, arguably, unhelpful. Accordingly, we synthesise a set of requirements for what might be termed ‘appropriate complexity modelling’ of quantitative coastal morphological change at scales commensurate with contemporary management and policy-making requirements: 1) The system being studied must be bounded with reference to the time and space scales at which behaviours of interest emerge and/or scientific or management problems arise; 2) model complexity and comprehensiveness must be appropriate to the problem at hand; 3) modellers should seek a priori insights into what kind of behaviours are likely to be evident at the scale of interest and the extent to which the behavioural validity of a model may be constrained by its underlying assumptions and its comprehensiveness; 4) informed by qualitative insights into likely dynamic behaviour, models should then be formulated with a view to resolving critical state changes; and 5) meso-scale modelling of coastal morphological change should reflect critically on the role of modelling and its relation to the observable world.

General rules for predicting where a catalytic reaction should occur on metal surfaces: A density functional theory study of C-H and C-O bond breaking/making on flat, stepped, and kinked metal surfaces

To predict where a catalytic reaction should occur is a fundamental issue scientifically. Technologically, it is also important because it can facilitate the catalyst's design. However, to date, the understanding of this issue is rather limited. In this work, two types of reactions, CH4 CH3 + H and CO C + 0 on two transition metal surfaces, were chosen as model systems aiming to address in general where a catalytic reaction should occur. The dissociations of CH4 - CH3 + H and CO --> C + O and their reverse reactions on flat, stepped, and kinked Rh and Pd surfaces were studied in detail. We find the following: First, for the CH4 Ch(3) + H reaction, the dissociation barrier is reduced by similar to0.3 eV on steps and kinks as compared to that on flat surfaces. On the other hand, there is essentially no difference in barrier for the association reaction of CH3 + H on the flat surfaces and the defects. Second, for the CO C + 0 reaction, the dissociation barrier decreases dramatically (more than 0.8 eV on Rh and Pd) on steps and kinks as compared to that on flat surfaces. In contrast to the CH3 + H reaction, the C + 0 association reaction also preferentially occurs on steps and kinks. We also present a detailed analysis of the reaction barriers in which each barrier is decomposed quantitatively into a local electronic effect and a geometrical effect. Our DFT calculations show that surface defects such as steps and kinks can largely facilitate bond breaking, while whether the surface defects could promote bond formation depends on the individual reaction as well as the particular metal. The physical origin of these trends is identified and discussed. On the basis of our results, we arrive at some simple rules with respect to where a reaction should occur: (i) defects such as steps are always favored for dissociation reactions as compared to flat surfaces; and (ii) the reaction site of the association reactions is largely related to the magnitude of the bonding competition effect, which is determined by the reactant and metal valency. Reactions with high valency reactants are more likely to occur on defects (more structure-sensitive), as compared to reactions with low valency reactants. Moreover, the reactions on late transition metals are more likely to proceed on defects than those on the early transition metals.

Insight into the solvent effect: A density functional theory study of cisplatin hydrolysis

The solvent effect on reactions in solutions is crucial for many systems. In this study, the reaction barrier with respect to the number of solvent molecules included in the system is systematically studied using density function theory calculations. Our results show that the barriers rapidly converge with respect to the number of solvent molecules. The solvent effect is investigated by calculating cisplatin hydrolysis in several types of solvents. The results are analyzed and a linear relationship between the reaction barrier and the interaction strength of solvent-reactants is found. Insight into the general solvent effect is obtained. (c) 2006 American Institute of Physics.

A density functional theory study on the water formation at high coverages and the water effect in the Fischer-Tropsch synthesis

The O removal through water formation is an important process in the Fischer-Tropsch synthesis. In this study, both steps in water formation (O + H --> OH, OH + H --> H2O) are studied on the stepped Co(0001) at high coverages using density functional theory. We find the following. (i) In both O-O and O-OH co-adsorption systems, two transition states (TSs) were located for the O hydrogenation: in one TS, both O and H are on the same terrace, and in the other they are at the interface between the step edge and the terrace below. (ii) In both the O-O and O-OH co-adsorption systems, the O hydrogenation at the interface is easier (E-a = 0.32 eV in the O-O system, E-a = 1.10 eV in the O-OH system) than that on the same terrace (E-a = 1.49 eV in the O-O system, E-a = 1.80 eV in the O-OH system). (iii) In both the O-O and O-OH co-adsorption systems, only one TS for the OH hydrogenation was located, in which both OH and H are on the same terrace. (iv) Compared to the OH hydrogenation in the O-OH system (E-a = 1.46 eV), the reaction in the OH-OH system (E-a = 0.64 eV) is much easier. The barrier differences and the water effect on the Fischer-Tropsch synthesis are discussed. A possible route with low barriers for water formation is proposed.

Macroscopic limit of time-dependent density-functional theory for adiabatic local approximations of the exchange-correlation kernel

Time-dependent density-functional theory is a rather accurate and efficient way to compute electronic excitations for finite systems. However, in the macroscopic limit (systems of increasing size), for the usual adiabatic random-phase, local-density, or generalized-gradient approximations, one recovers the Kohn-Sham independent-particle picture, and thus the incorrect band gap. To clarify this trend, we investigate the macroscopic limit of the exchange-correlation kernel in such approximations by means of an algebraical analysis complemented with numerical studies of a one-dimensional tight-binding model. We link the failure to shift the Kohn-Sham spectrum of these approximate kernels to the fact that the corresponding operators in the transition space act only on a finite subspace.

Selective Hydrogenation of Acetylene over Pd-Boron Catalysts: A Density Functional Theory Study

Boron-modified Pd catalysts have shown excellent performance for the selective hydrogenation of alkynes experimentally. In the current work, we investigated the hydrogenation of acetylene on boron-modified Pd(111) and Pd(211) surfaces, utilizing density functional theory calculations. The activity of acetylene hydrogenation has been studied by estimating the effective barrier of the whole process. The selectivity of ethylene formation is investigated from a comparison between the desorption and the hydrogenation of ethylene as well as comparison between the ethylene and the 1,3-butadiene formation. Formation of subsurface carbon and hydrogen on both boron-modified Pd(111) and Pd(211) surfaces has also been evaluated, since these have been reported to affect both the activity and the selectivity of acetylene hydrogenation to produce ethylene on Pd surfaces. Our results provide some important insights into the Pd B catalysts for selective hydrogenation of acetylene and also for more complex hydrogenation systems, such as stereoselective hydrogenation of longer chain alkynes and selective hydrogenation of vegetable oil.

Density functional theory study on the activation of molecular oxygen on a stepped gold surface in an aqueous environment: a new approach for simulating reactions in solution

The activation of oxygen molecules is an important issue in the gold-catalyzed partial oxidation of alcohols in aqueous solution. The complexity of the solution arising from a large number of solvent molecules makes it difficult to study the reaction in the system. In this work, O-2 activation on an Au catalyst is investigated using an effective approach to estimate the reaction barriers in the presence of solvent. Our calculations show that O-2 can be activated, undergoing OOH* in the presence of water molecules. The OOH* can readily be formed on Au(211) via four possible pathways with almost equivalent free energy barriers at the aqueous-solid interface: the direct or indirect activation of O-2 by surface hydrogen or the hydrolysis of O-2 following a Langmuir-Hinshelwood mechanism or an Eley-Rideal mechanism. Among them, the Eley-Rideal mechanism may be slightly more favorable due to the restriction of the low coverage of surface H on Au(211) in the other mechanisms. The results shed light on the importance of water molecules on the activation of oxygen in gold-catalyzed systems. Solvent is found to facilitate the oxygen activation process mainly by offering extra electrons and stabilizing the transition states. A correlation between the energy barrier and the negative charge of the reaction center is found. The activation barrier is substantially reduced by the aqueous environment, in which the first solvation shell plays the most important role in the barrier reduction. Our approach may be useful for estimating the reaction barriers in aqueous systems.

A density functional theory study of CH2 and H adsorption on Ni(111)

Ab initio total energy calculations within the density functional theory framework have been used to study the adsorption of CH2 and H as well as the coadsorption of CH2 and H on Ni(111). H binds strongly at threefold hollow sites with calculated adsorption energies of 2.60 and 2.54 eV at the face-centered-cubic (fcc) and hexagonal-close-packed (hcp) hollow sites, respectively. Adsorption energies and H-Ni distances are found to agree well with both experimental and theoretical results. CH2 adsorbs strongly at all high symmetry sites with calculated adsorption energies of 3.26, 3.22, 3.14 and 2.36 eV at the fcc, hcp, bridge and top sites, respectively. Optimized structures are reported at all sites, and, in the most stable hollow sites there is considerable internal reorganization of the CH2 fragment. The CH2 molecule is tilted, the hydrogens are inequivalent and the C-H bonds are lengthened relative to the gas phase. In the CH2-H coadsorption systems the adsorbates have a tendency to move toward bridge sites. The bonding of all adsorbates to the surface is analyzed in detail. (C) 2000 American Institute of Physics. [S0021-9606(00)71213-X].

Density functional theory study of the interaction between CO and on a Pt surface:CO/Pt(111), O/Pt(111), and CO/O/Pt(111)

Ab initio total energy calculations within the Density Functional Theory framework were carried out for Pt(111), Pt(111)-p(2x2)-CO, Pt(111)-p(2x2)-O, and Pt(111)-p(2x2)-(CO+O) to provide an insight into the interaction between CO and O on metal surfaces, an important issue in CO oxidation, and also in promotion and poisoning effects of catalysis. The geometrical structures of these systems were optimized with respect to the total energy, the results of which agree with existing experimental values very well. It is found that (i) the local structures of Pt(111)-p(2x2)-(CO+O), such as the bond lengths of C-O, C-Pt, and O-Pt (chemisorbed O atom with Pt), are almost the same as that in Pt(111)-p(2x2)-CO and Pt(111)-p(2x2)-O, respectively, (ii) the total valence charge density distributions in Pt(111)-p(2x2)-(CO+O) are very similar to that in Pt(111)-p(2x2)-CO, except in the region of the chemisorbed oxygen atom, and also nearly identical to that in Pt(111)-p(2x2)-O, apart from in the region of the chemisorbed CO, and (iii) the chemisorption energy of CO on a precovered Pt(111)-p(2x2)-O and the chemisorption energy of O on a precovered Pt(111)-p(2x2)CO are almost equal to that in Pt(111)-p(2x2)-CO and Pt(111)-p(2x2)-O, respectively. These results indicate that the interaction between CO and chemisorbed oxygen on a metal surface is mainly shore range in nature. The discussions of Pt-CO and Pt-O bonding and the interaction between CO and the chemisorbed oxygen atom on Pt(111) are augmented by local densities of states and real space distributions of quantum states.

A density functional theory study on the interaction between chemisorbed CO and S on Rh(111)

Density functional theory calculations are carried out for Rh(111)-p(2 x 2)-CO, Rh(111)-p(2 x 2)-S, Rh(111)-p(2 x 2)-(S + CO), Rh(111)-p(3 x 3)-CO, Rh(111)-p(3 x 3)-S and Rh(111)-p(3 x 3)-(S + CO), aiming to shed some light on the S poisoning effect. Geometrical structures of these systems are optimized and chemisorption energies are determined. The presence of S does not significantly influence the geometrical structure and chemisorption energy of CO and vice versa, which strongly suggests that the interaction between CO and S on the Rh(111) surface is mainly short-range in nature. The long range electronic effect for the dramatic attenuation of the CO methanation activity by sulfur is likely to be incorrect. It is suggested that an ensemble effect may be dominant in the catalytic deactivation. (C) 1999 Elsevier Science B.V. All rights reserved.

A Density Functional Theory Study on the Active Center of Fe-Only Hydrogenase: Characterization and Electronic Structure of the Redox States

We have carried out extensive density functional theory (DFT) calculations for possible redox states of the active center in Fe-only hydrogenases. The active center is modeled by [(H(CH(3))S)(CO)(CN(-))Fe(p)(mu-DTN)(mu-CO)Fe(d)(CO)(CN(-))(L)](z) (z is the net charge in the complex; Fe(p)= the proximal Fe, Fe(d) = the distal Fe, DTN = (-SCH(2)NHCH(2)S-), L is the ligand that bonds with the Fed at the trans position to the bridging CO). Structures of possible redox states are optimized, and CO stretching frequencies are calculated. By a detailed comparison of all the calculated structures and the vibrational frequencies with the available experimental data, we find that (i) the fully oxidized, inactive state is an Fe(II)-Fe(II) state with a hydroxyl (OH(-)) group bonded at the Fe(d), (ii) the oxidized, active state is an Fe(II)-Fe(l) complex which is consistent with the assignment of Cao and Hall (J. Am. Chem. Soc. 2001, 123, 3734), and (iii) the fully reduced state is a mixture with the major component being a protonated Fe(l)-Fe(l) complex and the other component being its self-arranged form, Fe(II)-Fe(II) hydride, Our calculations also show that the exogenous CO can strongly bond with the Fe(II)-Fe(l) species, but cannot bond with the Fe(l)-Fe(l) complex. This result is consistent with experiments that CO tends to inhibit the oxidized, active state, but not the fully reduced state. The electronic structures of all the redox states have been analyzed. It is found that a frontier orbital which is a mixing state between the e(g) of Fe and the 2pi of the bridging CO plays a key role concerning the reactivity of Fe-only hydrogenases: (1) it is unoccupied in the fully oxidized, inactive state, half-occupied in the oxidized, active state, and fully occupied in the fully reduced state; (ii) the e(g)-2pi orbital is a bonding state, and this is the key reason for stability of the low oxidation states, such as Fe(l)-Fe(l) complexes; and (iii) in the e(g)-2pi orbital more charge accumulates between the bridging CO and the Fe(d) than between the bridging CO and the Fe(p), and the occupation increase in this orbital will enhance the bonding between the bridging CO and the Fe(d), leading to the bridging-CO shift toward the Fe(d).