Density functional energy decomposition into one- and two-atom contributions

The present work provides a generalization of Mayer's energy decomposition for the density-functional theory (DFT) case. It is shown that one- and two-atom Hartree-Fock energy components in Mayer's approach can be represented as an action of a one-atom potential VA on a one-atom density ρ A or ρ B. To treat the exchange-correlation term in the DFT energy expression in a similar way, the exchange-correlation energy density per electron is expanded into a linear combination of basis functions. Calculations carried out for a number of density functionals demonstrate that the DFT and Hartree-Fock two-atom energies agree to a reasonable extent with each other. The two-atom energies for strong covalent bonds are within the range of typical bond dissociation energies and are therefore a convenient computational tool for assessment of individual bond strength in polyatomic molecules. For nonspecific nonbonding interactions, the two-atom energies are low. They can be either repulsive or slightly attractive, but the DFT results more frequently yield small attractive values compared to the Hartree-Fock case. The hydrogen bond in the water dimer is calculated to be between the strong covalent and nonbonding interactions on the energy scale

One- and two-center physical space partitioning of the energy in the density functional theory

A conceptually new approach is introduced for the decomposition of the molecular energy calculated at the density functional theory level of theory into sum of one- and two-atomic energy components, and is realized in the "fuzzy atoms" framework. (Fuzzy atoms mean that the three-dimensional physical space is divided into atomic regions having no sharp boundaries but exhibiting a continuous transition from one to another.) The new scheme uses the new concept of "bond order density" to calculate the diatomic exchange energy components and gives them unexpectedly close to the values calculated by the exact (Hartree-Fock) exchange for the same Kohn-Sham orbitals

Electronic and vibrational contributions to first hyperpolarizability of donor-acceptor-substituted azobenzene

In this study we report on the electronic and vibrational (hyper)polarizabilities of donor–acceptorsubstituted azobenzene. It is observed that both electronic and vibrational contributions to the electric dipole first hyperpolarizability of investigated photoactive molecule substantially depend on the conformation. The contributions to the nuclear relaxation first hyperpolarizability are found to be quite important in the case of two considered isomers (cis and trans). Although the double-harmonic term is found to be the largest in terms of magnitude, it is shown that the total value of the nuclear relaxation contribution to vibrational first hyperpolarizability is a result of subtle interplay of higher-order contributions. As a part of the study, we also assess the performance of long-range-corrected density functional theory in determining vibrational contributions to electric dipole (hyper)polarizabilities. In most cases, the applied long-range-corrected exchange correlation potentials amend the drawbacks of their conventional counterparts

Sobre les relacions lineals d'energia lliure: mesures de semblança molecular quàntica sobre funcions de densitat electrònica modificad

En aquest treball es presenta l'ús de funcions de densitat electrònica de forat de Fermi per incrementar el paper que pren una regió molecular concreta, considerada com a responsable de la reactivitat molecular, tot i mantenir la mida de la funció de densitat original. Aquestes densitats s'utilitzen per fer mesures d'autosemblança molecular quàntica i es presenten com una alternativa a l'ús de fragments moleculars aillats en estudis de relació entre estructura i propietat. El treball es complementa amb un exemple pràctic, on es correlaciona l'autosemblanca molecular a partir de densitats modificades amb l'energia d'una reacció isodòsmica

Fuzzy set Theory and Quantum Chemistry

Es discuteixen breument algunes consideracions sobre l'aplicació de la Teoria dels Conjunts difusos a la Química quàntica. Es demostra aqui que molts conceptes químics associats a la teoria són adequats per ésser connectats amb l'estructura dels Conjunts difusos. També s'explica com algunes descripcions teoriques dels observables quàntics es potencien tractant-les amb les eines associades als esmentats Conjunts difusos. La funció densitat es pren com a exemple de l'ús de distribucions de possibilitat al mateix temps que les distribucions de probabilitat quàntiques

Sobre els mapes de semblança quàntica molecular

Partint de les definicions usuals de Mesures de Semblança Quàntica (MSQ), es considera la dependència d'aquestes mesures respecte de la superposició molecular. Pel cas particular en qnè els sistemes comparats siguin una molècula i un Àtom i que les mesures es calculin amb l'aproximació EASA, les MSQ esdevenen funcions de les tres coordenades de l'espai. Mantenint fixa una de les tres coordenades, es pot representar fàcilment la variació del valor de semblança en un pla determinat, i obtenir els anomenats mapes de semblança. En aquest article, es comparen els mapes de semblança obtinguts amb diferents MSQ per a sistemes senzills

What Monte Carlo models can do and cannot do efficiently?

The question "what Monte Carlo models can do and cannot do efficiently" is discussed for some functional spaces that define the regularity of the input data. Data classes important for practical computations are considered: classes of functions with bounded derivatives and Holder type conditions, as well as Korobov-like spaces. Theoretical performance analysis of some algorithms with unimprovable rate of convergence is given. Estimates of computational complexity of two classes of algorithms - deterministic and randomized for both problems - numerical multidimensional integration and calculation of linear functionals of the solution of a class of integral equations are presented. (c) 2007 Elsevier Inc. All rights reserved.

Correlation between regioselectivity and site charge in propene polymerisation catalysed by metallocene

The reactions of propene with [Zr(cyclopentadienyl)(2)Me](+) have been investigated using density functional theory in order to study the correlation between regioselectivity and site charge in propene polymerisation. The reaction paths of the 1,2 and 2,1 additions of the methyl group to propene have been established. The geometries and energies of the reactants, transition states and products have been obtained using both PBEPBE/LANL2DZ and B3LYP/LANL2DZ methodologies. The results with both density functionals show that the activation energy for 1,2-insertion is lower than that for 2,1-insertion (Fig. 5) and this is consistent with the experiment results. Also for both density functionals, the difference of the thermal dynamic driving forces between the 2,1 product named 2-21 and the 1,2 product named 2-12 is significantly lower than the difference between the energy barriers. It is noted that in the reactants, the Mulliken partial charge on the central carbon atom C2 is positive and it can be concluded that 1,2-insertion is favoured because it can proceed via a cationic reaction.

Smoothed histogram modification for image processing

The technique of constructing a transformation, or regrading, of a discrete data set such that the histogram of the transformed data matches a given reference histogram is commonly known as histogram modification. The technique is widely used for image enhancement and normalization. A method which has been previously derived for producing such a regrading is shown to be “best” in the sense that it minimizes the error between the cumulative histogram of the transformed data and that of the given reference function, over all single-valued, monotone, discrete transformations of the data. Techniques for smoothed regrading, which provide a means of balancing the error in matching a given reference histogram against the information lost with respect to a linear transformation are also examined. The smoothed regradings are shown to optimize certain cost functionals. Numerical algorithms for generating the smoothed regradings, which are simple and efficient to implement, are described, and practical applications to the processing of LANDSAT image data are discussed.

A Wiener-Hopf type factorization for the exponential functional of Levy processes

For a Lévy process ξ=(ξt)t≥0 drifting to −∞, we define the so-called exponential functional as follows: Formula Under mild conditions on ξ, we show that the following factorization of exponential functionals: Formula holds, where × stands for the product of independent random variables, H− is the descending ladder height process of ξ and Y is a spectrally positive Lévy process with a negative mean constructed from its ascending ladder height process. As a by-product, we generate an integral or power series representation for the law of Iξ for a large class of Lévy processes with two-sided jumps and also derive some new distributional properties. The proof of our main result relies on a fine Markovian study of a class of generalized Ornstein–Uhlenbeck processes, which is itself of independent interest. We use and refine an alternative approach of studying the stationary measure of a Markov process which avoids some technicalities and difficulties that appear in the classical method of employing the generator of the dual Markov process.

A general method for finding extremal states of Hamiltonian dynamical systems, with applications to perfect fluids

In addition to the Hamiltonian functional itself, non-canonical Hamiltonian dynamical systems generally possess integral invariants known as ‘Casimir functionals’. In the case of the Euler equations for a perfect fluid, the Casimir functionals correspond to the vortex topology, whose invariance derives from the particle-relabelling symmetry of the underlying Lagrangian equations of motion. In a recent paper, Vallis, Carnevale & Young (1989) have presented algorithms for finding steady states of the Euler equations that represent extrema of energy subject to given vortex topology, and are therefore stable. The purpose of this note is to point out a very general method for modifying any Hamiltonian dynamical system into an algorithm that is analogous to those of Vallis etal. in that it will systematically increase or decrease the energy of the system while preserving all of the Casimir invariants. By incorporating momentum into the extremization procedure, the algorithm is able to find steadily-translating as well as steady stable states. The method is applied to a variety of perfect-fluid systems, including Euler flow as well as compressible and incompressible stratified flow.

Adsorption of organic molecules at the TiO2(110) surface: the effect of van der Waals interactions

Understanding the interaction of organic molecules with TiO2 surfaces is important for a wide range of technological applications. While density functional theory (DFT) calculations can provide valuable insight about these interactions, traditional DFT approaches with local exchange-correlation functionals suffer from a poor description of non-bonding van der Waals (vdW) interactions. We examine here the contribution of vdW forces to the interaction of small organic molecules (methane, methanol, formic acid and glycine) with the TiO2 (110) surface, based on DFT calculations with the optB88-vdW functional. The adsorption geometries and energies at different configurations were also obtained in the standard generalized gradient approximation (GGA-PBE) for comparison. We find that the optB88-vdW consistently gives shorter surface adsorbate-to-surface distances and slightly stronger interactions than PBE for the weak (physisorbed) modes of adsorption. In the case of strongly adsorbed (chemisorbed) molecules both functionals give similar results for the adsorption geometries, and also similar values of the relative energies between different chemisorption modes for each molecule. In particular both functionals predict that dissociative adsorption is more favourable than molecular adsorption for methanol, formic acid and glycine, in general agreement with experiment. The dissociation energies obtained from both functionals are also very similar, indicating that vdW interactions do not affect the thermodynamics of surface deprotonation. However, the optB88-vdW always predicts stronger adsorption than PBE. The comparison of the methanol adsorption energies with values obtained from a Redhead analysis of temperature programmed desorption data suggests that optB88-vdW significantly overestimates the adsorption strength, although we warn about the uncertainties involved in such comparisons.

Evolution PDEs and augmented eigenfunctions. half line.

The solution of an initial-boundary value problem for a linear evolution partial differential equation posed on the half-line can be represented in terms of an integral in the complex (spectral) plane. This representation is obtained by the unified transform introduced by Fokas in the 90's. On the other hand, it is known that many initial-boundary value problems can be solved via a classical transform pair, constructed via the spectral analysis of the associated spatial operator. For example, the Dirichlet problem for the heat equation can be solved by applying the Fourier sine transform pair. However, for many other initial-boundary value problems there is no suitable transform pair in the classical literature. Here we pose and answer two related questions: Given any well-posed initial-boundary value problem, does there exist a (non-classical) transform pair suitable for solving that problem? If so, can this transform pair be constructed via the spectral analysis of a differential operator? The answer to both of these questions is positive and given in terms of augmented eigenfunctions, a novel class of spectral functionals. These are eigenfunctions of a suitable differential operator in a certain generalised sense, they provide an effective spectral representation of the operator, and are associated with a transform pair suitable to solve the given initial-boundary value problem.

Morphological properties of superclusters of galaxies

We studied superclusters of galaxies in a volume-limited sample extracted from the Sloan Digital Sky Survey Data Release 7 and from mock catalogues based on a semi-analytical model of galaxy evolution in the Millennium Simulation. A density field method was applied to a sample of galaxies brighter than M(r) = -21+5 log h(100) to identify superclusters, taking into account selection and boundary effects. In order to evaluate the influence of the threshold density, we have chosen two thresholds: the first maximizes the number of objects (D1) and the second constrains the maximum supercluster size to similar to 120 h(-1) Mpc (D2). We have performed a morphological analysis, using Minkowski Functionals, based on a parameter, which increases monotonically from filaments to pancakes. An anticorrelation was found between supercluster richness (and total luminosity or size) and the morphological parameter, indicating that filamentary structures tend to be richer, larger and more luminous than pancakes in both observed and mock catalogues. We have also used the mock samples to compare supercluster morphologies identified in position and velocity spaces, concluding that our morphological classification is not biased by the peculiar velocities. Monte Carlo simulations designed to investigate the reliability of our results with respect to random fluctuations show that these results are robust. Our analysis indicates that filaments and pancakes present different luminosity and size distributions.

Soliton solutions for quasilinear Schrodinger equations with critical growth

In this paper we establish the existence of standing wave solutions for quasilinear Schrodinger equations involving critical growth. By using a change of variables, the quasilinear equations are reduced to semilinear one. whose associated functionals are well defined in the usual Sobolev space and satisfy the geometric conditions of the mountain pass theorem. Using this fact, we obtain a Cerami sequence converging weakly to a solution v. In the proof that v is nontrivial, the main tool is the concentration-compactness principle due to P.L. Lions together with some classical arguments used by H. Brezis and L. Nirenberg (1983) in [9]. (C) 2009 Elsevier Inc. All rights reserved.