Uso de la hoja de cálculo en la resolución de problemas sencillos de campo electrostático.

Resumen basado en el de la publicación

Electronic structure of point defects in controlled self-doping of the TiO2 (110) surface: Combined photoemission spectroscopy and density functional theory study

Point defects in metal oxides such as TiO2 are key to their applications in numerous technologies. The investigation of thermally induced nonstoichiometry in TiO2 is complicated by the difficulties in preparing and determining a desired degree of nonstoichiometry. We study controlled self-doping of TiO2 by adsorption of 1/8 and 1/16 monolayer Ti at the (110) surface using a combination of experimental and computational approaches to unravel the details of the adsorption process and the oxidation state of Ti. Upon adsorption of Ti, x-ray and ultraviolet photoemission spectroscopy (XPS and UPS) show formation of reduced Ti. Comparison of pure density functional theory (DFT) with experiment shows that pure DFT provides an inconsistent description of the electronic structure. To surmount this difficulty, we apply DFT corrected for on-site Coulomb interaction (DFT+U) to describe reduced Ti ions. The optimal value of U is 3 eV, determined from comparison of the computed Ti 3d electronic density of states with the UPS data. DFT+U and UPS show the appearance of a Ti 3d adsorbate-induced state at 1.3 eV above the valence band and 1.0 eV below the conduction band. The computations show that the adsorbed Ti atom is oxidized to Ti2+ and a fivefold coordinated surface Ti atom is reduced to Ti3+, while the remaining electron is distributed among other surface Ti atoms. The UPS data are best fitted with reduced Ti2+ and Ti3+ ions. These results demonstrate that the complexity of doped metal oxides is best understood with a combination of experiment and appropriate computations.

Surface and interstitial Ti diffusion at the rutile TiO2(110) surface

Diffusion of Ti through the TiO2 (110) rutile surface plays a key role in the growth and reactivity of TiO2. To understand the fundamental aspects of this important process, we present an analysis of the diffusion of Ti adspecies at the stoichiometric TiO2(110) surface using complementary computational methodologies of density functional theory corrected for on-site Coulomb interactions (DFT+U) and a charge equilibration (QEq) atomistic potential to identify minimum energy pathways. We find that diffusion of Ti from the surface to subsurface (and vice versa) follows an intersticialcy exchange mechanism, involving exchange of surface Ti with the 6-fold coordinated Ti below the bridging oxygen rows. Diffusion in the subsurface between layers also follows an interstitialcy mechanism. The diffusion of Ti is discussed in light of continued attempts to understand the re-oxidation of non-stoichiometric TiO2(110) surfaces.

Solar induced climate effects

Solar electromagnetic radiation powers Earth’s climate system and, consequently, it is often naively assumed that changes in this solar output must be responsible for changes in Earth’s climate. However, the Sun is close to a blackbody radiator and so emits according to its surface temperature and the huge thermal time constant of the outer part of the Sun limits the variability in surface temperature and hence output. As a result, on all timescales of interest, changes in total power output are limited to small changes in effective surface temperature (associated with magnetic fields) and potential, although as yet undetected, solar radius variations. Larger variations are seen in the UV part of the spectrum which is emitted from the lower solar atmosphere (the chromosphere) and which influences Earth’s stratosphere. There is interest in“top-down” mechanisms whereby solar UV irradiance modulates stratospheric temperatures and winds which, in turn, may influence the underlying troposphere where Earth’s climate and weather reside. This contrasts with “bottom-up” effects in which the small total solar irradiance (dominated by the visible and near-IR) variations cause surface temperature changes which drive atmospheric circulations. In addition to these electromagnetic outputs, the Sun modulates energetic particle fluxes incident on the Earth. Solar Energetic Particles (SEP) are emitted by solar flares and from the shock fronts ahead of supersonic (and super-Alfvenic) ejections of material from the solar atmosphere. These SEPs enhance the destruction of polar stratospheric ozone which could be an additional form of top-down climate forcing. Even more energetic are Galactic Cosmic Rays (GCRs). These particles are not generated by the Sun, rather they originate at the shock fronts emanating from violent galactic events such as supernovae explosions; however, the expansion of the solar magnetic field into interplanetary space means that the Sun modulates the number of GCRs reaching Earth. These play a key role in enabling Earth’s global electric (thunderstorm) circuit and it has been proposed that they also modulate the formation of clouds. Both electromagnetic and corpuscular solar effects are known to vary over the solar magnetic cycle which is typically between 10 and 14 yrs in length (with an average close to 11 yrs). The solar magnetic field polarity at any one phase of one of these activity cycles is opposite to that at the same phase of the next cycle and this influences some phenomena, for example GCRs, which therefore show a 22 yr (“Hale”) cycle on average. Other phenomena, such as irradiance modulation, do not depend on the polarity of the magnetic field and so show only the basic 11-yr activity cycle. However, any effects on climate are much more significant for solar drifts over centennial timescales. This chapter discusses and evaluates potential effects on Earth’s climate system of variations in these solar inputs. Because of the great variety of proposed mechanisms, the wide range of timescales studied (from days to millennia) and the many debates (often triggered by the application of inadequate statistical methods), the literature on this subject is vast, complex, divergent and rapidly changing: consequently the number of references cited in this review is very large (yet still only a small fraction of the total).

A simple analysis of thermodynamic properties for classical plasmas: I. theory

By eliminating the short range negative divergence of the Debye–Hückel pair distribution function, but retaining the exponential charge screening known to operate at large interparticle separation, the thermodynamic properties of one-component plasmas of point ions or charged hard spheres can be well represented even in the strong coupling regime. Predicted electrostatic free energies agree within 5% of simulation data for typical Coulomb interactions up to a factor of 10 times the average kinetic energy. Here, this idea is extended to the general case of a uniform ionic mixture, comprising an arbitrary number of components, embedded in a rigid neutralizing background. The new theory is implemented in two ways: (i) by an unambiguous iterative algorithm that requires numerical methods and breaks the symmetry of cross correlation functions; and (ii) by invoking generalized matrix inverses that maintain symmetry and yield completely analytic solutions, but which are not uniquely determined. The extreme computational simplicity of the theory is attractive when considering applications to complex inhomogeneous fluids of charged particles.

Cyclic extrusion of a lava dome based on a stick-slip mechanism

Lava dome eruptions are sometimes characterised by large periodic fluctuations in extrusion rate over periods of hours that may be accompanied by Vulcanian explosions and pyroclastic flows. We consider a simple system of nonlinear equations describing a 1D flow of lava extrusion through a deep elastic dyke feeding a shallower cylindrical conduit in order to simulate this short-period cyclicity. Stick-slip conditions depending on a critical shear stress are assumed at the wall boundary of the cylindrical conduit. By analogy with the behaviour of industrial polymers in a plastic extruder, the elastic dyke acts like a barrel and the shallower cylindrical portion of the conduit as a die for the flow of magma acting as a polymer. When we applied the model to the Soufrière Hills Volcano, Montserrat, for which the key parameters have been evaluated from previous studies, cyclic extrusions with periods from 3 to 30 h were readily simulated, matching observations. The model also reproduces the reduced period of cycles observed when a major unloading event occurs due to lava dome collapse.

Pulsatory andesite lava flow at Bagana Volcano

Using a time series of TerraSAR-X spaceborne radar images we have measured the pulsatory motion of an andesite lava flow over a 14-month period at Bagana volcano, Papua New Guinea. Between October 2010 and December 2011, lava flowed continuously down the western flank of the volcano forming a 3 km-long blocky lava flow with a channel, levees, overflows and branches. We captured four successive pulses of lava advancing down the channel system, the first such behaviour of an andesite flow to be recorded using radar. Each pulse had a volume of the order of 107 m3 emplaced over many weeks. The average extrusion rate estimated from the radar data was 0.92 ± 0.35 m3 s-1 , and varied between 0.3 and 1.8 m3 s-1, with higher rates occurring earlier in each pulse. This, together with observations of sulphur dioxide emissions, explosions and incandescence suggest a variable supply rate of magma through Bagana’s conduit as the most likely source of the pulsatory behaviour.

Complex formation in solutions of oppositely charged polyelectrolytes at different polyion compositions and salt content

The formation of complexes in solutions containing positively charged polyions (polycations) and a variable amount of negatively charged polyions (polyanions) has been investigated by Monte Carlo simulations. The polyions were described as flexible chains of charged hard spheres interacting through a screened Coulomb potential. The systems were analyzed in terms of cluster compositions, structure factors, and radial distribution functions. At 50% charge equivalence or less, complexes involving two polycations and one polyanion were frequent, while closer to charge equivalence, larger clusters were formed. Small and neutral complexes dominated the solution at charge equivalence in a monodisperse system, while larger clusters again dominated the solution when the polyions were made polydisperse. The cluster composition and solution structure were also examined as functions of added salt by varying the electrostatic screening length. The observed formation of clusters could be rationalized by a few simple rules.

A Monte Carlo study of solutions of oppositely charged polyelectrolytes

The formation of complexes appearing in solutions containing oppositely charged polyelectrolytes has been investigated by Monte Carlo simulations using two different models. The polyions are described as flexible chains of 20 connected charged hard spheres immersed in a homogenous dielectric background representing water. The small ions are either explicitly included or their effect described by using a screened Coulomb potential. The simulated solutions contained 10 positively charged polyions with 0, 2, or 5 negatively charged polyions and the respective counterions. Two different linear charge densities were considered, and structure factors, radial distribution functions, and polyion extensions were determined. A redistribution of positively charged polyions involving strong complexes formed between the oppositely charged polyions appeared as the number of negatively charged polyions was increased. The nature of the complexes was found to depend on the linear charge density of the chains. The simplified model involving the screened Coulomb potential gave qualitatively similar results as the model with explicit small ions. Finally, owing to the complex formation, the sampling in configurational space is nontrivial, and the efficiency of different trial moves was examined.

Sensitivity of OMI SO2 measurements to variable eruptive behaviour at Soufrière Hills Volcano, Montserrat

Since 2004, the satellite-borne Ozone Mapping Instrument (OMI) has observed sulphur dioxide (SO2) plumes during both quiescence and effusive eruptive activity at Soufrière Hills Volcano, Montserrat. On average, OMI detected a SO2 plume 4-6 times more frequently during effusive periods than during quiescence in the 2008-2010 period. The increased ability of OMI to detect SO2 during eruptive periods is mainly due to an increase in plume altitude rather than a higher SO2 emission rate. Three styles of eruptive activity cause thermal lofting of gases (Vulcanian explosions; pyroclastic flows; a hot lava dome) and the resultant plume altitudes are estimated from observations and models. Most lofting plumes from Soufrière Hills are derived from hot domes and pyroclastic flows. Although Vulcanian explosions produced the largest plumes, some produced only negligible SO2 signals detected by OMI. OMI is most valuable for monitoring purposes at this volcano during periods of lava dome growth and during explosive activity.

Hydrodynamical simulations of Galactic fountains - I. Evolution of single fountains

The ejection of the gas out of the disc in late-type galaxies is related to star formation and is due mainly to Type II supernovae. In this paper, we studied in detail the development of the Galactic fountains in order to understand their dynamical evolution and their influence on the redistribution of the freshly delivered metals over the disc. To this aim, we performed a number of 3D hydrodynamical radiative cooling simulations of the gas in the Milky Way where the whole Galaxy structure, the Galactic differential rotation and the supernova explosions generated by a single OB association are considered. A typical fountain powered by 100 Type II supernovae may eject material up to similar to 2 kpc which than collapses back mostly in the form of dense, cold clouds and filaments. The majority of the gas lifted up by the fountains falls back on the disc remaining within a radial distance Delta R = 0.5 kpc from the place where the fountain originated. This localized circulation of disc gas does not influence the radial chemical gradients on large scale, as required by the chemical models of the Milky Way which reproduce the metallicity distribution without invoking large fluxes of metals. Simulations of multiple fountains fuelled by Type II supernovae of different OB associations will be presented in a companion paper.

Galactic outflows and the pollution of the galactic environment by supernovae

We here explore the effects of the SN explosions into the halo of star-forming galaxies like the Milky Way. Successive randomly distributed and clustered SNe explosions cause the formation of hot superbubbles that drive either fountains or galactic winds above the galactic disk, depending on the amount and concentration of energy that is injected by the SNe. In a galactic fountain, the ejected gas is re-captured by the gravitational potential and falls back onto the disk. From 3D non-equilibrium radiative cooling hydrodynamical simulations of these fountains, we find that they may reach altitudes up to about 5 kpc in the halo and thus allow for the formation of the so called intermediate-velocity-clouds (IVCs) which are often observed in the halos of disk galaxies. The high-velocity-clouds that are also observed but at higher altitudes (of up to 12 kpc) require another mechanism to explain their production. We argue that they could be formed either by the capture of gas from the intergalactic medium and/or by the action of magnetic fields that are carried to the halo with the gas in the fountains. Due to angular momentum losses to the halo, we find that the fountain material falls back to smaller radii and is not largely spread over the galactic disk. Instead, the SNe ejecta fall nearby the region where the fountain was produced, a result which is consistent with recent chemical models of the galaxy. The fall back material leads to the formation of new generations of molecular clouds and to supersonic turbulence feedback in the disk. (C) 2009 COSPAR. Published by Elsevier Ltd. All rights reserved.

NU-PROCESS IN EXOTIC MODELS

The exact physical conditions generating the abundances of r-elements in environments such as supernovae explosions are still under debate. We evaluated the characteristics expected for the neutrino wind in the proposed model of type-II supernova driven by conversion of nuclear matter to strange matter. Neutrinos will change the final abundance of elements after freeze out of r-process nucleosynthesis, specially those close to mass peaks.

The intraplate Porto dos Gauchos seismic zone in the Amazon craton - Brazil

The largest earthquake observed in the stable continental interior of the South American plate occurred in Serra do Tombador, Mato Grosso state - Brazil, on January 31,1955 with a magnitude of 6.2 m(b). Since then no other earthquake has been located near the 1955 epicentre. However, in Porto dos Gauchos, 100 km northeast of Serra do Tombador, a recurrent seismicity has been observed since 1959. Both Serra do Tombador and Porto dos Gauchos are located in the Phanerozoic Parecis basin. Two magnitude 5 earthquakes occurred in Porto dos Gauchos, in 1998 and 2005, with intensities up to VI and V, respectively. These two main shocks were followed by aftershock sequences lasting more than three years each. Local seismic stations have been deployed by the Seismological Observatory of the University of Brasilia since 1998 to study the ""Porto dos Gauchos"" seismic zone (PGSZ). A local seismic refraction survey was carried out with two explosions to help define the seismic velocity model. Both the 1998 and 2005 earthquake sequences occurred in the same WSW-ENE oriented fault zone with right-lateral strike-slip mechanisms. The epicentral zone is in the Parecis basin, near its northern border where there are buried grabens, generally trending WNW-ESE, such as the deep Mesoproterozoic Caiabis graben which lies partly beneath the Parecis basin. However, the epicentral distribution indicates that the 1998 and 2005 sequences are related to a N60 degrees E fault which probably crosses the entire Caiabis graben. The 1955 earthquake, despite the uncertainty in its epicentre, does not seem to be directly related to any buried graben either. The seismicity in the Porto dos Gauchos seismic zone, therefore, is not directly related to rifted crust. The probable direction of the maximum horizontal stress near Porto dos Gauchos is roughly E-W, consistent with other focal mechanisms further south in the Pantanal basin and Paraguay. but seems to be different from the NW-SE direction observed further north in the Amazon basin. The recurrent seismicity observed in Porto dos Gauchos, and the large 1955 earthquake nearby, make this area of the Parecis basin one of the most important seismic zones of Brazil. (C) 2009 Elsevier B.V. All rights reserved.