Questioning Kenneth N. Waltz (1924-2013)

Two little-known academics examine the doctoral thesis of a young theorist, Kenneth Waltz. They conclude that his work is important, despite its ambiguities. Some years later, they catch up with Professor Waltz between classes and explore how his ideas have developed.

Structure determination from powder diffraction data

Advances made over the past decade in structure determination from powder diffraction data are reviewed with particular emphasis on algorithmic developments and the successes and limitations of the technique. While global optimization methods have been successful in the solution of molecular crystal structures, new methods are required to make the solution of inorganic crystal structures more routine. The use of complementary techniques such as NMR to assist structure solution is discussed and the potential for the combined use of X-ray and neutron diffraction data for structure verification is explored. Structures that have proved difficult to solve from powder diffraction data are reviewed and the limitations of structure determination from powder diffraction data are discussed. Furthermore, the prospects of solving small protein crystal structures over the next decade are assessed.

Rational modification of the hierarchy of intermolecular interactions in molecular crystal structures by using tunable halogen bonds

A family of 16 isomolecular salts (3-XpyH)(2)[MX'(4)] (3-XpyH=3-halopyridinium; M=Co, Zn; X=(F), Cl, Br, (I); X'=Cl, Br, I) each containing rigid organic cations and tetrahedral halometallate anions has been prepared and characterized by X-ray single crystal and/or powder diffraction. Their crystal structures reflect the competition and cooperation between non-covalent interactions: N-H center dot center dot center dot X'-M hydrogen bonds, C-X center dot center dot center dot X'-M halogen bonds and pi-pi stacking. The latter are essentially unchanged in strength across the series, but both halogen bonds and hydrogen bonds are modified in strength upon changing the halogens involved. Changing the organic halogen (X) from F to I strengthens the C-X center dot center dot center dot X'-M halogen bonds, whereas an analogous change of the inorganic halogen (X') weakens both halogen bonds and N-H center dot center dot center dot X'-M hydrogen bonds. By so tuning the strength of the putative halogen bonds from repulsive to weak to moderately strong attractive interactions, the hierarchy of the interactions has been modified rationally leading to systematic changes in crystal packing. Three classes of crystal structure are obtained. In type A (C-F center dot center dot center dot X'-M) halogen bonds are absent. The structure is directed by N-H center dot center dot center dot X'-M hydrogen bonds and pi-stacking interactions. In type B structures, involving small organic halogens (X) and large inorganic halogens (X'), long (weak) C-X center dot center dot center dot X'-M interactions are observed with type I halogen-halogen interaction geometries (C-X center dot center dot center dot X' approximate to X center dot center dot center dot X'-M approximate to 155 degrees), but hydrogen bonds still dominate. Thus, minor but quite significant perturbations from the type A structure arise. In type C, involving larger organic halogens (X) and smaller inorganic halogens (X'), stronger halogen bonds are formed with a type II halogen-halogen interaction geometry (C-X center dot center dot center dot X' approximate to 180 degrees; X center dot center dot center dot X'-M approximate to 110 degrees) that is electrostatically attractive. The halogen bonds play a major role alongside hydrogen bonds in directing the type C structures, which as a result are quite different from type A and B.

GDASH: a grid-enabled program for structure solution from powder diffraction data

The simulated annealing approach to structure solution from powder diffraction data, as implemented in the DASH program, is easily amenable to parallelization at the individual run level. Very large scale increases in speed of execution can therefore be achieved by distributing individual DASH runs over a network of computers. The GDASH program achieves this by packaging DASH in a form that enables it to run under the Univa UD Grid MP system, which harnesses networks of existing computing resources to perform calculations.

MDASH: a multi-core-enabled program for structure solution from powder diffraction data

The simulated annealing approach to structure solution from powder diffraction data, as implemented in the DASH program, is easily amenable to parallelization at the individual run level. Modest increases in speed of execution can therefore be achieved by executing individual DASH runs on the individual cores of CPUs.

Temperature- and Pressure-induced Proton Transfer in the 1:1 Adduct Formed between Squaric Acid and 4,4 '-Bipyridine

We have applied a combination of spectroscopic and diffraction methods to study the adduct formed between squaric acid and bypridine, which has been postulated to exhibit proton transfer associated with a single-crystal to single-crystal phase transition at ca. 450 K. A combination of X-ray single-crystal and very-high flux powder neutron diffraction data confirmed that a proton does transfer from the acid to the base in the high-temperature form. Powder X-ray diffraction measurements demonstrated that the transition was reversible but that a significant kinetic energy barrier must be overcome to revert to the original structure. Computational modeling is consistent with these results. Modeling also revealed that, while the proton transfer event would be strongly discouraged in the gas phase, it occurs in the solid state due to the increase in charge state of the molecular ions and their arrangement inside the lattice. The color change is attributed to a narrowing of the squaric acid to bipyridine charge-transfer energy gap. Finally, evidence for the possible existence of two further phases at high pressure is also presented.

Vibrational spectroscopy of a compound with a CS7 ring

The product of the Asinger reaction between elemental sulfur, n-butylamine and acetophenone is 8-(n-butylaminophenylmethyliden)-1,2,3,4,5,6,7-heptathiocane which contains a CS7 ring. A combination of infrared, Raman and inelastic neutron scattering spectroscopies with periodic density functional theory calculations is used to provide a complete assignment of the vibrational spectra of this unusual species. The similarity between the Raman spectra of the compound and that of elemental sulfur is particularly striking. Copyright (C) 2009 John Wiley & Sons, Ltd.

Reactive Halogens in the Marine Boundary Layer (RHaMBLe): the tropical North Atlantic experiments

The NERC UK SOLAS-funded Reactive Halogens in the Marine Boundary Layer (RHaMBLe) programme comprised three field experiments. This manuscript presents an overview of the measurements made within the two simultaneous remote experiments conducted in the tropical North Atlantic in May and June 2007. Measurements were made from two mobile and one ground-based platforms. The heavily instrumented cruise D319 on the RRS Discovery from Lisbon, Portugal to São Vicente, Cape Verde and back to Falmouth, UK was used to characterise the spatial distribution of boundary layer components likely to play a role in reactive halogen chemistry. Measurements onboard the ARSF Dornier aircraft were used to allow the observations to be interpreted in the context of their vertical distribution and to confirm the interpretation of atmospheric structure in the vicinity of the Cape Verde islands. Long-term ground-based measurements at the Cape Verde Atmospheric Observatory (CVAO) on São Vicente were supplemented by long-term measurements of reactive halogen species and characterisation of additional trace gas and aerosol species during the intensive experimental period. This paper presents a summary of the measurements made within the RHaMBLe remote experiments and discusses them in their meteorological and chemical context as determined from these three platforms and from additional meteorological analyses. Air always arrived at the CVAO from the North East with a range of air mass origins (European, Atlantic and North American continental). Trace gases were present at stable and fairly low concentrations with the exception of a slight increase in some anthropogenic components in air of North American origin, though NOx mixing ratios during this period remained below 20 pptv. Consistency with these air mass classifications is observed in the time series of soluble gas and aerosol composition measurements, with additional identification of periods of slightly elevated dust concentrations consistent with the trajectories passing over the African continent. The CVAO is shown to be broadly representative of the wider North Atlantic marine boundary layer; measurements of NO, O3 and black carbon from the ship are consistent with a clean Northern Hemisphere marine background. Aerosol composition measurements do not indicate elevated organic material associated with clean marine air. Closer to the African coast, black carbon and NO levels start to increase, indicating greater anthropogenic influence. Lower ozone in this region is possibly associated with the increased levels of measured halocarbons, associated with the nutrient rich waters of the Mauritanian upwelling. Bromide and chloride deficits in coarse mode aerosol at both the CVAO and on D319 and the continuous abundance of inorganic gaseous halogen species at CVAO indicate significant reactive cycling of halogens. Aircraft measurements of O3 and CO show that surface measurements are representative of the entire boundary layer in the vicinity both in diurnal variability and absolute levels. Above the inversion layer similar diurnal behaviour in O3 and CO is observed at lower mixing ratios in the air that had originated from south of Cape Verde, possibly from within the ITCZ. ECMWF calculations on two days indicate very different boundary layer depths and aircraft flights over the ship replicate this, giving confidence in the calculated boundary layer depth.

Modelling monsoons: understanding and predicting current and future behaviour

The global monsoon system is so varied and complex that understanding and predicting its diverse behaviour remains a challenge that will occupy modellers for many years to come. Despite the difficult task ahead, an improved monsoon modelling capability has been realized through the inclusion of more detailed physics of the climate system and higher resolution in our numerical models. Perhaps the most crucial improvement to date has been the development of coupled ocean-atmosphere models. From subseasonal to interdecadal time scales, only through the inclusion of air-sea interaction can the proper phasing and teleconnections of convection be attained with respect to sea surface temperature variations. Even then, the response to slow variations in remote forcings (e.g., El Niño—Southern Oscillation) does not result in a robust solution, as there are a host of competing modes of variability that must be represented, including those that appear to be chaotic. Understanding the links between monsoons and land surface processes is not as mature as that explored regarding air-sea interactions. A land surface forcing signal appears to dominate the onset of wet season rainfall over the North American monsoon region, though the relative role of ocean versus land forcing remains a topic of investigation in all the monsoon systems. Also, improved forecasts have been made during periods in which additional sounding observations are available for data assimilation. Thus, there is untapped predictability that can only be attained through the development of a more comprehensive observing system for all monsoon regions. Additionally, improved parameterizations - for example, of convection, cloud, radiation, and boundary layer schemes as well as land surface processes - are essential to realize the full potential of monsoon predictability. A more comprehensive assessment is needed of the impact of black carbon aerosols, which may modulate that of other anthropogenic greenhouse gases. Dynamical considerations require ever increased horizontal resolution (probably to 0.5 degree or higher) in order to resolve many monsoon features including, but not limited to, the Mei-Yu/Baiu sudden onset and withdrawal, low-level jet orientation and variability, and orographic forced rainfall. Under anthropogenic climate change many competing factors complicate making robust projections of monsoon changes. Absent aerosol effects, increased land-sea temperature contrast suggests strengthened monsoon circulation due to climate change. However, increased aerosol emissions will reflect more solar radiation back to space, which may temper or even reduce the strength of monsoon circulations compared to the present day. Precipitation may behave independently from the circulation under warming conditions in which an increased atmospheric moisture loading, based purely on thermodynamic considerations, could result in increased monsoon rainfall under climate change. The challenge to improve model parameterizations and include more complex processes and feedbacks pushes computing resources to their limit, thus requiring continuous upgrades of computational infrastructure to ensure progress in understanding and predicting current and future behaviour of monsoons.

Characterisation of amorphous and nanocrystalline molecular materials by total scattering

The use of high-energy X-ray total scattering coupled with pair distribution function analysis produces unique structural fingerprints from amorphous and nanostructured phases of the pharmaceuticals carbamazepine and indomethacin. The advantages of such facility-based experiments over laboratory-based ones are discussed and the technique is illustrated with the characterisation of a melt-quenched sample of carbamazepine as a nanocrystalline (4.5 nm domain diameter) version of form III.

Skill Rating by Bayesian Inference

Systems Engineering often involves computer modelling the behaviour of proposed systems and their components. Where a component is human, fallibility must be modelled by a stochastic agent. The identification of a model of decision-making over quantifiable options is investigated using the game-domain of Chess. Bayesian methods are used to infer the distribution of players’ skill levels from the moves they play rather than from their competitive results. The approach is used on large sets of games by players across a broad FIDE Elo range, and is in principle applicable to any scenario where high-value decisions are being made under pressure.

Evidence for long-term regional changes in precipitation on the East Coast Mountains in Mauritius

Global climate change and its impacts are being increasingly studied and precipitation trends are one of the measures of quantifying climate change especially in the tropics. This study uses daily rainfall data to determine if there are changes in the long-term trends in rainfall variability in the East Coast Mountains of Mauritius during the last few decades, and to investigate the factors influencing the trends in the inter-annual to inter-decadal rainfall variability. Statistical modelling has been used to investigate the trends in total seasonal rainfall, the number of rain days and the mean amount of rain per rainy days and the local, regional and large-scale factors that affect them on inter-annual to inter-decadal time scales. The strongest inter-decadal trend was found in the number of rain days for both rainfall seasons, and the other variables were found to have weak or insignificant trends. Both local factors, such as the surrounding sea surface temperatures and large-scale phenomena such as Indian Monsoon and the El Niño Southern Oscillation were found to influence rainfall patterns.

Progress on “Changing coastlines: data assimilation for morphodynamic prediction and predictability”

The task of assessing the likelihood and extent of coastal flooding is hampered by the lack of detailed information on near-shore bathymetry. This is required as an input for coastal inundation models, and in some cases the variability in the bathymetry can impact the prediction of those areas likely to be affected by flooding in a storm. The constant monitoring and data collection that would be required to characterise the near-shore bathymetry over large coastal areas is impractical, leaving the option of running morphodynamic models to predict the likely bathymetry at any given time. However, if the models are inaccurate the errors may be significant if incorrect bathymetry is used to predict possible flood risks. This project is assessing the use of data assimilation techniques to improve the predictions from a simple model, by rigorously incorporating observations of the bathymetry into the model, to bring the model closer to the actual situation. Currently we are concentrating on Morecambe Bay as a primary study site, as it has a highly dynamic inter-tidal zone, with changes in the course of channels in this zone impacting the likely locations of flooding from storms. We are working with SAR images, LiDAR, and swath bathymetry to give us the observations over a 2.5 year period running from May 2003 – November 2005. We have a LiDAR image of the entire inter-tidal zone for November 2005 to use as validation data. We have implemented a 3D-Var data assimilation scheme, to investigate the improvements in performance of the data assimilation compared to the previous scheme which was based on the optimal interpolation method. We are currently evaluating these different data assimilation techniques, using 22 SAR data observations. We will also include the LiDAR data and swath bathymetry to improve the observational coverage, and investigate the impact of different types of observation on the predictive ability of the model. We are also assessing the ability of the data assimilation scheme to recover the correct bathymetry after storm events, which can dramatically change the bathymetry in a short period of time.

Pilot study of the short-term physico-chemical stability of atenolol tablets stored in a multi-compartment compliance aid

Study objectives: There is a possibility that lower air, moisture and light protection could impact on physico-chemical stability of medicines inside multi-compartment compliance aids (MCCAs), although this has not yet been proved. The objectives of the study were to examine the physico-chemical stability of atenolol tablets stored in a compliance aid at room temperature, and at elevated temperature and humidity to simulate practice conditions. Methods: Atenolol 100 mg tablets in 28-chamber, plastic compliance aids with transparent lids were stored for four weeks at room temperature and at 40°C with 75% relative humidity. Tablets were also stored at room temperature in original packaging and Petri dishes. Physical tests were conducted to standards as laid down in the British Pharmacopoeia 2005, and dissolution to those of the United States Pharmacopoeia volume 24. Chemical stability was assessed by a validated high-performance liquid chromatography (HPLC) method. Results: Tablets at room temperature in original packaging, in compliance aids and Petri dishes remained the same in appearance and passed physico-chemical tests. Tablets exposed to 40°C with 75% relative humidity in compliance aids passed tests for uniformity of weight, friability and chemical stability but became pale and moist, softer (82 newtons ± 4; p< 0.0001) than tablets in the original packaging (118 newtons ± 6), more friable (0.14% loss of mass) compared with other tablets (0.005%), and failed the tests for disintegration (>15 minutes) and dissolution (only 15% atenolol released at 30 minutes). Conclusion: Although chemical stability was unaffected, storage in compliance aids at 40°C with 75% relative humidity softened atenolol tablets, prolonged disintegration time and hindered dissolution which could significantly reduce bioavailability. This formulation could be suitable for storage in compliance aids at 25°C, but not in hotter, humid weather.