Perturbation theory for classical thermodynamic Green's functions

A systematic time-dependent perturbation scheme for classical canonical systems is developed based on a Wick's theorem for thermal averages of time-ordered products. The occurrence of the derivatives with respect to the canonical variables noted by Martin, Siggia, and Rose implies that two types of Green's functions have to be considered, the propagator and the response function. The diagrams resulting from Wick's theorem are "double graphs" analogous to those introduced by Dyson and also by Kawasaki, in which the response-function lines form a "tree structure" completed by propagator lines. The implication of a fluctuation-dissipation theorem on the self-energies is analyzed and compared with recent results by Deker and Haake.

New constraints on a light spinless particle coupled to photons

We obtain new stringent constraints on a light spinless particle f coupled only to photons at low energies, considering its effects on the extragalactic photon background, the black-body spectrum of the cosmic microwave background radiation and the cosmological abundance of deuterium.

Exact gravitational shockwaves and Planckian scattering on branes

We obtain a solution describing a gravitational shock wave propagating along a Randall-Sundrum brane. The interest of such a solution is twofold: on the one hand, it is the first exact solution for a localized source on a Randall-Sundrum three-brane. On the other hand, one can use it to study forward scattering at Planckian energies, including the effects of the continuum of Kaluza-Klein modes. We map out the different regimes for the scattering obtained by varying the center-of-mass energy and the impact parameter. We also discuss exact shock waves in ADD scenarios with compact extra dimensions.

A fully-automatic caudate nucleus segmentation of brain MRI: application in volumetric analysis of pediatric attention-deficit/hyperactivity disorder

Background Accurate automatic segmentation of the caudate nucleus in magnetic resonance images (MRI) of the brain is of great interest in the analysis of developmental disorders. Segmentation methods based on a single atlas or on multiple atlases have been shown to suitably localize caudate structure. However, the atlas prior information may not represent the structure of interest correctly. It may therefore be useful to introduce a more flexible technique for accurate segmentations. Method We present Cau-dateCut: a new fully-automatic method of segmenting the caudate nucleus in MRI. CaudateCut combines an atlas-based segmentation strategy with the Graph Cut energy-minimization framework. We adapt the Graph Cut model to make it suitable for segmenting small, low-contrast structures, such as the caudate nucleus, by defining new energy function data and boundary potentials. In particular, we exploit information concerning the intensity and geometry, and we add supervised energies based on contextual brain structures. Furthermore, we reinforce boundary detection using a new multi-scale edgeness measure. Results We apply the novel CaudateCut method to the segmentation of the caudate nucleus to a new set of 39 pediatric attention-deficit/hyperactivity disorder (ADHD) patients and 40 control children, as well as to a public database of 18 subjects. We evaluate the quality of the segmentation using several volumetric and voxel by voxel measures. Our results show improved performance in terms of segmentation compared to state-of-the-art approaches, obtaining a mean overlap of 80.75%. Moreover, we present a quantitative volumetric analysis of caudate abnormalities in pediatric ADHD, the results of which show strong correlation with expert manual analysis. Conclusion CaudateCut generates segmentation results that are comparable to gold-standard segmentations and which are reliable in the analysis of differentiating neuroanatomical abnormalities between healthy controls and pediatric ADHD.

The prediction of chemical reactivity from first principles: Collision and reaction dynamics

Aquest treball fa una revisió de mesures experimentals i càlculs teòrics sobre la dinàmica de col·lisions i reaccions moleculars. Els experiments se centren en col·lisions, a energies intermèdies, que involucren sistemes del tipus ió-àtom i iómolècula, per les quals es mesuren seccions eficaces totals, estat a estat, així com aquelles que discerneixen les diferents contribucions del moment angular d'espín. Els resultats obtinguts s'interpreten satisfactòriament en termes d'acoblaments no adiabàtics entre els diferents estats electrònics dels sistemes col·lisionants. Els càlculs teòrics utilitzen la metodologia quasiclàssica, així com metodologies mecanoquàntiques recentment desenvolupades, tant aproximades com exactes. S'han obtingut resultats totalment convergits per sistemes tipus, mentre que s'han analitzat, de manera detallada i extensiva, les característiques dinàmiques de sistemes triatòmic, tetraatòmic i pentaatòmic.

Chemisorption of atomic aluminum on Si(111): Evidence for an adsorbate-induced relaxation based on ab initio cluster-model calculations

Interaction models of atomic Al with Si4H9, Si4H7, and Si6H9 clusters have been studied to simulate Al chemisorption on the Si(111) surface in the atop, fourfold atop, and open sites. Calculations were carried out using nonempirical pseudopotentials in the framework of the ab initio Hartree-Fock procedure. Equilibrium bond distances, binding energies for adsorption, and vibrational frequencies of the adatoms are calculated. Several basis sets were used in order to show the importance of polarization effects, especially in the binding energies. Final results show the importance of considering adatom-induced relaxation effects to specify the order of energy stabilities for the three different sites, the fourfold atop site being the preferred one, in agreement with experimental findings.

Chemisorption of group-III metals on the Si(111) and Ge(111) surfaces: An ab initio study

Chemisorption of group-III metal adatoms on Si(111) and Ge(111) has been studied through the ab initio Hartree-Fock method including nonempirical pseudopotentials and using cluster models to simulate the surface. Three different high-symmetry sites (atop, eclipsed, and open) have been considered by using X4H9, X4H7, and X6H9 (X=Si,Ge) cluster models. In a first step, ideal surface geometries have been used. Metal-induced reconstruction upon chemisorption has also been taken into account. Equilibrium distances, binding energies, and vibrational frequencies have been obtained and compared with available experimental data. From binding-energy considerations, the atop and eclipsed sites seem to be the most favorable ones and thus a coadsorption picture may be suggested. Group-III metals exhibit a similar behavior and the same is true for Si(111) and Ge(111) surfaces when chemisorption is considered.

Evidence for two different bonding mechanisms of Al on Si(111)

By means of the ab initio cluster-model approach, we present theoretical evidence for two different mechanisms of bonding of atomic Al to Si(111). On the atop site (T1) the interaction of atomic Al to Si(111) is characteristic of an ionic bond whereas interaction above the threefold eclipsed site (T4) leads to the formation of a typical covalent bond. Moreover, both sites have a similar interaction energy if electronic correlation effects are included. While the conclusions regarding the nature of the chemisorption bond in the two sites do not depend either on the cluster-model size, the kind of embedding hydrogen atoms used, or the quality of the wave function (Hartree-Fock or configuration interaction), the chemisorption energy depends strongly on the wave function used. In fact, inclusion of correlation energy is necessary to properly describe the interaction energies.

Measures of ionicity of alkaline-earth oxides from the analysis of ab initio cluster wave functions

We present an analysis of the M-O chemical bonding in the binary oxides MgO, CaO, SrO, BaO, and Al2O3 based on ab initio wave functions. The model used to represent the local environment of a metal cation in the bulk oxide is an MO6 cluster which also includes the effect of the lattice Madelung potential. The analysis of the wave functions for these clusters leads to the conclusion that all the alkaline-earth oxides must be regarded as highly ionic oxides; however, the ionic character of the oxides decreases as one goes from MgO, almost perfectly ionic, to BaO. In Al2O3 the ionic character is further reduced; however, even in this case, the departure from the ideal, fully ionic, model of Al3+ is not exceptionally large. These conclusions are based on three measures, a decomposition of the Mq+-Oq- interaction energy, the number of electrons associated to the oxygen ions as obtained from a projection operator technique, and the analysis of the cation core-level binding energies. The increasing covalent character along the series MgO, CaO, SrO, and BaO is discussed in view of the existing theoretical models and experimental data.

Performance of correlation functionals in ab initio chemisorption cluster-model calculations: Alkali metals on Si(111)

The performance of different correlation functionals has been tested for alkali metals, Li to Cs, interacting with cluster models simulating different active sites of the Si(111) surface. In all cases, the ab initio Hartree-Fock density has been obtained and used as a starting point. The electronic correlation energy is then introduced as an a posteriori correction to the Hartree-Fock energy using different correlation functionals. By making use of the ionic nature of the interaction and of different dissociation limits we have been able to prove that all functionals tested introduce the right correlation energy, although to a different extent. Hence, correlation functionals appear as an effective and easy way to introduce electronic correlation in the ab initio Hartree-Fock description of the chemisorption bond in complex systems where conventional configuration interaction techniques cannot be used. However, the calculated energies may differ by some tens of eV. Therefore, these methods can be employed to get a qualitative idea of how important correlation effects are, but they have some limitations if accurate binding energies are to be obtained.

Neutral atoms in ionic lattices: Excited states of KCl:Ag(0)

The optical-absorption spectrum of a cationic Ag0 atom in a KCl crystal has been studied theoretically by means of a series of cluster models of increasing size. Excitation energies have been determined by means of a multiconfigurational self-consistent field procedure followed by a second-order perturbation correlation treatment. Moreover results obtained within the density-functional framework are also reported. The calculations confirm the assignment of bands I and IV to transitions of the Ag-5s electron into delocalized states with mainly K-4s,4p character. Bands II and III have been assigned to internal transitions on the Ag atom, which correspond to the atomic Ag-4d to Ag-5s transition. We also determine the lowest charge transfer (CT) excitation energy and confirm the assignment of band VI to such a transition. The study of the variation of the CT excitation energy with the Ag-Cl distance R gives additional support to a large displacement of the Cl ions due to the presence of the Ag0 impurity. Moreover, from the present results, it is predicted that on passing to NaCl:Ag0 the CT onset would be out of the optical range while the 5s-5p transition would undergo a redshift of 0.3 eV. These conclusions, which underline the different character of involved orbitals, are consistent with experimental findings. The existence of a CT transition in the optical range for an atom inside an ionic host is explained by a simple model, which also accounts for the differences with the more common 3d systems. The present study sheds also some light on the R dependence of the s2-sp transitions due to s2 ions like Tl+.

Surface structure and stability of PdZn and PtZn alloys: Density-functional slab model studies

Geometric parameters of binary (1:1) PdZn and PtZn alloys with CuAu-L10 structure were calculated with a density functional method. Based on the total energies, the alloys are predicted to feature equal formation energies. Calculated surface energies of PdZn and PtZn alloys show that (111) and (100) surfaces exposing stoichiometric layers are more stable than (001) and (110) surfaces comprising alternating Pd (Pt) and Zn layers. The surface energy values of alloys lie between the surface energies of the individual components, but they differ from their composition weighted averages. Compared with the pure metals, the valence d-band widths and the Pd or Pt partial densities of states at the Fermi level are dramatically reduced in PdZn and PtZn alloys. The local valence d-band density of states of Pd and Pt in the alloys resemble that of metallic Cu, suggesting that a similar catalytic performance of these systems can be related to this similarity in the local electronic structures.

First principles LDA+U and GGA+U study of cerium oxides: Dependence on the effective U-parameter

The electronic structure and properties of cerium oxides (CeO2 and Ce2O3) have been studied in the framework of the LDA+U and GGA(PW91)+U implementations of density functional theory. The dependence of selected observables of these materials on the effective U parameter has been investigated in detail. The examined properties include lattice constants, bulk moduli, density of states, and formation energies of CeO2 and Ce2O3. For CeO2, the LDA+U results are in better agreement with experiment than the GGA+U results whereas for the computationally more demanding Ce2O3 both approaches give comparable accuracy. Furthermore, as expected, Ce2O3 is much more sensitive to the choice of the U value. Generally, the PW91 functional provides an optimal agreement with experiment at lower U energies than LDA does. In order to achieve a balanced description of both kinds of materials, and also of nonstoichiometric CeO2¿x phases, an appropriate choice of U is suggested for LDA+U and GGA+U schemes. Nevertheless, an optimum value appears to be property dependent, especially for Ce2O3. Optimum U values are found to be, in general, larger than values determined previously in a self-consistent way.

Wind turbine for developing countries

Aquest treball avarca des de reunions amb personal de qualitat d’una fàbrica decremalleres fins a la proposta d’un prototip d’una màquina que solucionés els seusproblemes concrets de producció. Pel camí s’ha fet recerca envers de solucionsprèvies que fossin factibles a realitzar en uns pocs mesos i amb els medis limitatsdels que es desponiaLes solucions, tant de mètodes d’assaig i observació com de disseny són un reflexclar de procés d’enginyeria en l’àmbit industrial espanyol. Els medis escassegen,però amb temps (tampoc gaire) i enginy es troben formes de satisfer lesnecessitats. Al cap i a la fi és el que s’espera de nosaltres com a enginyers.

Performance Analysis of Two Quantum Reaction Dynamics Codes: Time-Dependent and Time-Independent Strategies

The computer simulation of reaction dynamics has nowadays reached a remarkable degree of accuracy. Triatomic elementary reactions are rigorously studied with great detail on a straightforward basis using a considerable variety of Quantum Dynamics computational tools available to the scientific community. In our contribution we compare the performance of two quantum scattering codes in the computation of reaction cross sections of a triatomic benchmark reaction such as the gas phase reaction Ne + H2+ %12. NeH++ H. The computational codes are selected as representative of time-dependent (Real Wave Packet [ ]) and time-independent (ABC [ ]) methodologies. The main conclusion to be drawn from our study is that both strategies are, to a great extent, not competing but rather complementary. While time-dependent calculations advantages with respect to the energy range that can be covered in a single simulation, time-independent approaches offer much more detailed information from each single energy calculation. Further details such as the calculation of reactivity at very low collision energies or the computational effort related to account for the Coriolis couplings are analyzed in this paper.