A FORTRAN-90 low-energy electron diffraction program (LEED90 v1.1)

We describe a FORTRAN-90 program to compute low-energy electron diffraction I(V) curves. Plane-waves and layer doubling are used to compute the inter-layer multiple-scattering, while the intra-layer multiple-scattering is computed in the standard way expanding the wavefield on a basis of spherical waves. The program is kept as general as possible, in order to allow testing different parts of multiple-scattering calculations. In particular, it can handle non-diagonal t-matrices describing the scattering of non-spherical potentials, anisotropic vibrations, anharmonicity, etc. The program does not use old FORTRAN flavours, and has been written keeping in mind the advantage for parallelism brought forward by FORTRAN-90.

Surface chemistry of glycine on Pt{111} in different aqueous environments

Adsorption of glycine on Ptf111g under UHV conditions and in different aqueous environments was studied by XPS (UHV and ambient pressure) and NEXAFS. Under UHV conditions, glycine adsorbs in its neutral molecular state up to about 0.15 ML. Further deposition leads to the formation of an additional zwitterionic species, which is in direct contact with the substrate surface, followed by the growth of multilayers, which also consist of zwitterions. The neutral surface species is most stable and decomposes at 360 K through a multi-step process which includes the formation of methylamine and carbon monoxide. When glycine and water are co-adsorbed in UHV at low temperatures (< 170 K) inter-layer diffusion is inhibited and the surface composition depends on the adsorption sequence. Water adsorbed on top of a glycine layer does not lead to significant changes in its chemical state. When glycine is adsorbed on top of a pre-adsorbed chemisorbed water layer or thick ice layer, however, it is found in its zwitterionic state, even at low coverage. No difference is seen in the chemical state of glycine when the layers are exposed to ambient water vapor pressure up to 0.2 Torr at temperatures above 300 K. Also the decomposition temperature stays the same, 360 K, irrespective of the water vapor pressure. Only the reaction path of the decomposition products is affected by ambient water vapor.

The interplay of metal-atom ordering, Fermi leveltuning and thermoelectric properties in cobalt shandites Co3M2S2 (M = Sn, In)

A combination of structural, physical and computational techniques including powder X-ray and neutron diffraction, SQUID magnetometry, electrical and thermal transport measurements, DFT calculations and 119Sn Mössbauer and X-ray photoelec-tron spectroscopies has been applied to Co3Sn2-xInxS2 (0 ≤ x ≤ 2) in an effort to understand the relationship between metal-atom ordering and physical properties as the Fermi level is systematically varied. Whilst solid solution behavior is found throughout the composition region, powder neutron diffraction reveals that indium preferentially occupies an inter-layer site over an alternative kagome-like intra-layer site. DFT calculations indicate that this ordering, which leads to a lowering of energy, is related to the dif-fering bonding properties of tin and indium. Spectroscopic data suggest that throughout the composition range 0 ≤ x ≤ 2, all ele-ments adopt oxidation states that are significantly reduced from expectations based on formal charges. Chemical substitution ena-bles the electrical transport properties to be controlled through tuning of the Fermi level within a region of the density of states, which comprises narrow bands of predominantly Co d-character. This leads to a compositionally-induced double metal-to-semiconductor-to-metal transition. The marked increase in the Seebeck coefficient as the semiconducting region is approached leads to a substantial improvement in the thermoelectric figure of merit, ZT, which exhibits a maximum of ZT = 0.32 at 673 K. At 425 K, the figure of merit for phases in the region 0.8 ≤ x ≤ 0.85 is amongst the highest reported for sulphide phases, suggesting these materials may have applications in low-grade waste heat recovery.

The impact of charge transfer and structural disorder on the thermoelectric properties of cobalt intercalated TiS2

A family of phases, CoxTiS2 (0 ≤ x ≤ 0.75) has been prepared and characterised by powder X-ray and neutron diffraction, electrical and thermal transport property measurements, thermal analysis and SQUID magnetometry. With increasing cobalt content, the structure evolves from a disordered arrangement of cobalt ions in octahedral sites located in the van der Waals’ gap (x ≤ 0.2), through three different ordered vacancy phases, to a second disordered phase at x ≥ 0.67. Powder neutron diffraction reveals that both octahedral and tetrahedral inter-layer sites are occupied in Co0.67TiS2. Charge transfer from the cobalt guest to the TiS2 host affords a systematic tuning of the electrical and thermal transport properties. At low levels of cobalt intercalation (x < 0.1), the charge transfer increases the electrical conductivity sufficiently to offset the concomitant reduction in |S|. This, together with a reduction in the overall thermal conductivity leads to thermoelectric figures of merit that are 25 % higher than that of TiS2, ZT reaching 0.30 at 573 K for CoxTiS2 with 0.04 ≤ x ≤ 0.08. Whilst the electrical conductivity is further increased at higher cobalt contents, the reduction in |S| is more marked due to the higher charge carrier concentration. Furthermore both the charge carrier and lattice contributions to the thermal conductivity are increased in the electrically conductive ordered-vacancy phases, with the result that the thermoelectric performance is significantly degraded. These results illustrate the competition between the effects of charge transfer from guest to host and the disorder generated when cobalt cations are incorporated in the inter-layer space.

An enhancement of the ionospheric sporadic‐E layer in response to negative polarity cloud‐to‐ground lightning

Lightning data, collected using a Boltek Storm Tracker system installed at Chilton, UK, were used to investigate the mean response of the ionospheric sporadic-E layer to lightning strokes in a superposed epoch study. The lightning detector can discriminate between positive and negative lightning strokes and between cloud-to-ground ( CG) and inter-cloud ( IC) lightning. Superposed epoch studies carried out separately using these subsets of lightning strokes as trigger events have revealed that the dominant cause of the observed ionospheric enhancement in the Es layer is negative cloud-to-ground lightning.

MetUM-GOML: a near-globally coupled atmosphere–ocean-mixed-layer model

Well-resolved air–sea interactions are simulated in a new ocean mixed-layer, coupled configuration of the Met Office Unified Model (MetUM-GOML), comprising the MetUM coupled to the Multi-Column K Profile Parameterization ocean (MC-KPP). This is the first globally coupled system which provides a vertically resolved, high near-surface resolution ocean at comparable computational cost to running in atmosphere-only mode. As well as being computationally inexpensive, this modelling framework is adaptable– the independent MC-KPP columns can be applied selectively in space and time – and controllable – by using temperature and salinity corrections the model can be constrained to any ocean state. The framework provides a powerful research tool for process-based studies of the impact of air–sea interactions in the global climate system. MetUM simulations have been performed which separate the impact of introducing inter- annual variability in sea surface temperatures (SSTs) from the impact of having atmosphere–ocean feedbacks. The representation of key aspects of tropical and extratropical variability are used to assess the performance of these simulations. Coupling the MetUM to MC-KPP is shown, for example, to reduce tropical precipitation biases, improve the propagation of, and spectral power associated with, the Madden–Julian Oscillation and produce closer-to-observed patterns of springtime blocking activity over the Euro-Atlantic region.

An assessment of upper ocean salinity content from the ocean reanalyses inter-comparison project (ORA-IP)

Many institutions worldwide have developed ocean reanalyses systems (ORAs) utilizing a variety of ocean models and assimilation techniques. However, the quality of salinity reanalyses arising from the various ORAs has not yet been comprehensively assessed. In this study, we assess the upper ocean salinity content (depth-averaged over 0–700 m) from 14 ORAs and 3 objective ocean analysis systems (OOAs) as part of the Ocean Reanalyses Intercomparison Project. Our results show that the best agreement between estimates of salinity from different ORAs is obtained in the tropical Pacific, likely due to relatively abundant atmospheric and oceanic observations in this region. The largest disagreement in salinity reanalyses is in the Southern Ocean along the Antarctic circumpolar current as a consequence of the sparseness of both atmospheric and oceanic observations in this region. The West Pacific warm pool is the largest region where the signal to noise ratio of reanalysed salinity anomalies is >1. Therefore, the current salinity reanalyses in the tropical Pacific Ocean may be more reliable than those in the Southern Ocean and regions along the western boundary currents. Moreover, we found that the assimilation of salinity in ocean regions with relatively strong ocean fronts is still a common problem as seen in most ORAs. The impact of the Argo data on the salinity reanalyses is visible, especially within the upper 500m, where the interannual variability is large. The increasing trend in global-averaged salinity anomalies can only be found within the top 0–300m layer, but with quite large diversity among different ORAs. Beneath the 300m depth, the global-averaged salinity anomalies from most ORAs switch their trends from a slightly growing trend before 2002 to a decreasing trend after 2002. The rapid switch in the trend is most likely an artefact of the dramatic change in the observing system due to the implementation of Argo.

What controls the vertical distribution of aerosol? Relationships between process sensitivity in HadGEM3–UKCA and inter-model variation from AeroCom Phase II

The vertical profile of aerosol is important for its radiative effects, but weakly constrained by observations on the global scale, and highly variable among different models. To investigate the controlling factors in one particular model, we investigate the effects of individual processes in HadGEM3–UKCA and compare the resulting diversity of aerosol vertical profiles with the inter-model diversity from the AeroCom Phase II control experiment. In this way we show that (in this model at least) the vertical profile is controlled by a relatively small number of processes, although these vary among aerosol components and particle sizes. We also show that sufficiently coarse variations in these processes can produce a similar diversity to that among different models in terms of the global-mean profile and, to a lesser extent, the zonal-mean vertical position. However, there are features of certain models' profiles that cannot be reproduced, suggesting the influence of further structural differences between models. In HadGEM3–UKCA, convective transport is found to be very important in controlling the vertical profile of all aerosol components by mass. In-cloud scavenging is very important for all except mineral dust. Growth by condensation is important for sulfate and carbonaceous aerosol (along with aqueous oxidation for the former and ageing by soluble material for the latter). The vertical extent of biomass-burning emissions into the free troposphere is also important for the profile of carbonaceous aerosol. Boundary-layer mixing plays a dominant role for sea salt and mineral dust, which are emitted only from the surface. Dry deposition and below-cloud scavenging are important for the profile of mineral dust only. In this model, the microphysical processes of nucleation, condensation and coagulation dominate the vertical profile of the smallest particles by number (e.g. total CN  >  3 nm), while the profiles of larger particles (e.g. CN  >  100 nm) are controlled by the same processes as the component mass profiles, plus the size distribution of primary emissions. We also show that the processes that affect the AOD-normalised radiative forcing in the model are predominantly those that affect the vertical mass distribution, in particular convective transport, in-cloud scavenging, aqueous oxidation, ageing and the vertical extent of biomass-burning emissions.