Further considerations of cosmic ray modulation of infra-red radiation in the atmosphere

Understanding effects of ionisation in the lower atmosphere is a new interdisciplinary area, crossing the traditionally distinct scientific boundaries between astro-particle and atmospheric physics and also requiring understanding of both heliospheric and magnetospheric influences on cosmic rays. Following the paper of Erlykin et al. (2014) we develop further the interpretation of our observed changes in long-wave (LW) radiation, Aplin and Lockwood (2013) by taking account of both cosmic ray ionisation yields and atmospheric radiative transfer. To demonstrate this, we show that the thermal structure of the whole atmosphere needs to be considered along with the vertical profile of ionisation. Allowing for, in particular, ionisation by all components of a cosmic ray shower and not just by the muons, reveals that the effect we have detected is certainly not inconsistent with laboratory observations of the LW absorption cross section. The analysis presented here, although very different from that of Erlykin et al., does come to the same conclusion that the events detected by AL were not caused by individual cosmic ray primaries – not because it is impossible on energetic grounds, but because events of the required energy are too infrequent for the 12 h_1 rate at which they were seen by the AL experiment. The present paper numerically models the effect of three different scenario changes to the primary GCR spectrum which all reproduce the required magnitude of the effect observed by AL. However, they cannot solely explain the observed delay in the peak effect which, if confirmed, would appear to open up a whole new and interesting area in the study of water oligomers and their effects on LW radiation. We argue that a technical artefact in the AL experiment is highly unlikely and that our initial observations merit both a wide-ranging follow-up experiment and more rigorous, self-consistent, three-dimensional radiative transfer modelling

DK-1D: a drift-kinetic simulation tool for modelling the shear Alfvén wave and its interaction with collisionless plasma

We present a highly accurate tool for the simulation of shear Alfven waves (SAW) in collisionless plasma. SAW are important in space plasma environments because for small perpendicular scale lengths they can support an electric field parallel to the ambient magnetic field. Electrons can be accelerated by the parallel electric field and these waves have been implicated as the source of vibrant auroral displays. However, the parallel electric field carried by SAW is small in comparison to the perpendicular electric field of the wave, making it difficult to measure directly in the laboratory, or by satellites in the near-Earth plasma environment. In this paper, we present a simulation code that provides a means to study in detail the SAW-particle interaction in both space and laboratory plasma. Using idealised, small-amplitude propagating waves with a single perpendicular wavenumber, the simulation code accurately reproduces the damping rates and parallel electric field amplitudes predicted by linear theory for varying temperatures and perpendicular scale lengths. We present a rigorous kinetic derivation of the parallel electric field strength for small-amplitude SAW and show that commonly-used inertial and kinetic approximations are valid except for where the ratio of thermal to Alfv\'{e}n speed is between 0.7 and 1.0. We also present nonlinear simulations of large-amplitude waves and show that in cases of strong damping, the damping rates and parallel electric field strength deviate from linear predictions when wave energies are greater than only a few percent of the plasma kinetic energy, a situation which is often observed in the magnetosphere. The drift-kinetic code provides reliable, testable predictions of the parallel electric field strength which can be investigated directly in the laboratory, and will help to bridge the gap between studies of SAW in man-made and naturally occuring plasma.

EISCAT observations of plasma convection and the high-latitude, winter F-region during substorm activity

A 24 h period of observations by the EISCAT radar and other ground-based instrumentation is used to study the role of plasma convection in determining the morphology of the high-latitude F-region during winter. It is suggested that, in the afternoon sector of the polar convection pattern, rapid zonal (westward) flows caused low F-region electron densities due to an extension of the mid-latitude trough far into the sunlit hemisphere. Low densities on the dawn side prior to 0600 UT may also have been due to a trough-like feature. Although the generation mechanism is unclear, the trough may be the fossil remnant of a substorm. Around midnight, high F-region densities were seen, probably due to plasma flow emerging from the cap through soft particle precipitation in the auroral oval. Two substorms occurred at times when the radar was south of the auroral oval. Both caused enhanced convection speeds, a swing to equatorward flow, enhanced E-region densities and a depleted F-region. The first was seen as a Westward Travelling Surge, and the swing to purely southward flow which followed the surge front did not return to westward flows until 80–110 min later. The Harang discontinuity was observed co-rotating eastwards between the substorms, 65 ± 30 min before the separatrix between the dawn and dusk convection cells.

Notations and conventions in molecular spectroscopy: part 3. Permutation and permutation-inversion symmetry notation

The field of Molecular Spectroscopy was surveyed in order to determine a set of conventions and symbols which are in common use in the spectroscopic literature. This document, which is Part 3 in a series, deals with symmetry notation referring to groups that involve nuclear permutations and the inversion operation. Further parts will follow, dealing inter alia with vibration-rotation spectroscopy and electronic spectroscopy.

Longitudinal selenium status in healthy British adults: assessment using biochemical and molecular biomarkers

Human selenium (Se) requirements are currently based on biochemical markers of Se status. In rats, tissue glutathione peroxidase-1 (Gpx1) mRNA levels can be used effectively to determine Se requirements; blood Gpx1 mRNA levels decrease in Se-deficient rats, so molecular biology-based markers have potential for human nutrition assessment. To study the efficacy of molecular biology markers for assessing Se status in humans, we conducted a longitudinal study on 39 subjects (age 45 +/- 11) in Reading, UK. Diet diaries (5 day) and blood were obtained from each subject at 2, 8, 17 and 23 weeks, and plasma Se, glutathione peroxidase (Gpx3) enzyme activity, and selenoprotein mRNA levels were determined. There were no significant longitudinal effects on Se biomarkers. Se intake averaged 48 +/- 14 mu g/d. Plasma Se concentrations averaged 1.13 +/- 0.16 mu mol/l. Plasma Se v. energy-corrected Se intake (ng Se/kJ/d) was significantly correlated, but neither Gpx3 activity v. Se intake (ng Se/kJ/d) nor Gpx3 activity v. plasma Se was significantly correlated. Collectively, this indicates that subjects were on the plateaus of the response curves. Selenoprotein mRNAs were quantitated in total RNA isolated from whole blood, but mRNA levels for Gpx1, selenoprotein H, and selenoprotein W (all highly regulated by Se in rodents), as well selenoprotein P, Gpx3, and phospholipid hydroperoxide glutathione peroxidase were also not significantly correlated with plasma Se. Thus selenoprotein molecular biomarkers, as well as traditional biochemical markers, are unable to further distinguish differences in Se status in these Se replete subjects. The efficacy of molecular biomarkers to detect Se deficiency needs to be tested in Se-deficient populations.

Insulin receptor substrate-2 gene variants in subjects with metabolic syndrome: association with plasma fatty acid levels and insulin resistance

Several insulin receptor substrate-2 (IRS-2) polymorphisms have been studied in relation to insulin resistance and type 2 diabetes. To examine whether the genetic variability at the IRS-2 gene locus was associated with the degree of insulin resistance and plasma fatty acid levels in metabolic syndrome (MetS) subjects. Methods and results: Insulin sensitivity, insulin secretion, glucose effectiveness, plasma fatty acid composition and three IRS-2 tag-single nucleotide polymorphisms (SNPs) were determined in 452 MetS subjects. Among subjects with the lowest level of monounsaturated (MUFA) (below the median), the rs2289046 A/A genotype was associated with lower glucose effectiveness (p<0.038), higher fasting insulin concentrations (p<0.028) and higher HOMA IR (p<0.038) as compared to subjects carrying the minor G-allele (A/G and G/G). In contrast, among subjects with the highest level of MUFA (above the median), the A/A genotype was associated with lower fasting insulin concentrations and HOMA-IR, whereas individuals carrying the G allele and with the highest level of ω-3 polyunsaturated fatty acids (above the median) showed lower fasting insulin (p<0.01) and HOMA-IR (p<0.02) as compared with A/A subjects. Conclusion: The rs2289046 polymorphism at the IRS2 gene locus may influence insulin sensitivity by interacting with certain plasma fatty acids in MetS subjects.

Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

EISCAT observations of ion composition and temperature anisotropy in the high-latitude F-region

The papers by Winser et al. [(1990) J. atmos. terr. Phys.52, 501] and Häggström and Collis [(1990) J. atmos. terr. Phys.52, 519] used plasma flows and ion temperatures, as measured by the EISCAT tristatic incoherent scatter radar, to investigate changes in the ion composition of the ionospheric F-layer at high latitudes, in response to increases in the speed of plasma convection. These studies reported that the ion composition rapidly changed from mainly O+ to almost completely (>90%) molecular ions, following rapid increases in ion drift speed by >1 km s−1. These changes appeared inconsisent with theoretical considerations of the ion chemistry, which could not account for the large fractions of molecular ions inferred from the obsevations. In this paper, we discuss two causes of this discrepancy. First, we reevaluate the theoretical calculations for chemical equilibrium and show that, if we correct the derived temperatures for the effect of the molecular ions, and if we employ more realistic dependences of the reaction rates on the ion temperature, the composition changes derived for the faster convection speeds can be explained. For the Winser et al. observations with the radar beam at an aspect angle of ϕ = 54.7° to the geomagnetic field, we now compute a change to 89% molecular ions in < 2 min, in response to the 3 km s−1 drift. This is broadly consistent with the observations. But for the two cases considered by Häggström and Collis, looking along the field line (ϕ = 0°), we compute the proportion of molecular ions to be only 4 and 16% for the observed plasma drifts of 1.2 and 1.6 km s−1, respectively. These computed proportions are much smaller than those derived experimentally (70 and 90%). We attribute the differences to the effects of non-Maxwellian, anisotropic ion velocity distribution functions. We also discuss the effect of ion composition changes on the various radar observations that report anisotropies of ion temperature.

Measuring ion temperatures and studying the ion energy balance in the high-latitude ionosphere

Data are presented for a nighttime ion heating event observed by the EISCAT radar on 16 December 1988. In the experiment, the aspect angle between the radar beam and the geomagnetic field was fixed at 54.7°, which avoids any ambiguity in derived ion temperature caused by anisotropy in the ion velocity distribution function. The data were analyzed with an algorithm which takes account of the non-Maxwellian line-of-sight ion velocity distribution. During the heating event, the derived spectral distortion parameter (D∗) indicated that the distribution function was highly distorted from a Maxwellian form when the ion drift increased to 4 km s−1. The true three-dimensional ion temperature was used in the simplified ion balance equation to compute the ion mass during the heating event. The ion composition was found to change from predominantly O4 to mainly molecular ions. A theoretical analysis of the ion composition, using the MSIS86 model and published values of the chemical rate coefficients, accounts for the order-of-magnitude increase in the atomic/molecular ion ratio during the event, but does not successfully explain the very high proportion of molecular ions that was observed.

Comparisons between EISCAT observations and model calculations of the high latitude ionosphere

Calculations using a numerical model of the convection dominated high latitude ionosphere are compared with observations made by EISCAT as part of the UK-POLAR Special Programme. The data used were for 24–25 October 1984, which was characterized by an unusually steady IMF, with Bz < 0 and By > 0; in the calculations it was assumed that a steady IMF implies steady convection conditions. Using the electric field models of Heppner and Maynard (1983) appropriate to By > 0 and precipitation data taken from Spiroet al. (1982), we calculated the velocities and electron densities appropriate to the EISCAT observations. Many of the general features of the velocity data were reproduced by the model. In particular, the phasing of the change from eastward to westward flow in the vicinity of the Harang discontinuity, flows near the dayside throat and a region of slow flow at higher latitudes near dusk were well reproduced. In the afternoon sector modelled velocity values were significantly less than those observed. Electron density calculations showed good agreement with EISCAT observations near the F-peak, but compared poorly with observations near 211 km. In both cases, the greatest disagreement occurred in the early part of the observations, where the convection pattern was poorly known and showed some evidence of long term temporal change. Possible causes for the disagreement between observations and calculations are discussed and shown to raise interesting and, as yet, unresolved questions concerning the interpretation of the data. For the data set used, the late afternoon dip in electron density observed near the F-peak and interpreted as the signature of the mid-latitude trough is well reproduced by the calculations. Calculations indicate that it does not arise from long residence times of plasma on the nightside, but is the signature of a gap between two major ionization sources, viz. photoionization and particle precipitation.

Field-perpendicular and field-aligned plasma flows observed by EISCAT during a prolonged period of northward IMF

The effect of a prolonged period of strongly northward Interplanetary Magnetic Field (IMF) on the high-latitude F-region is studied using data from the EISCAT Common Programme Zero mode of operation on 11–12 August 1982. The analysis of the raw autocorrelation functions is kept to the directly derived parameters Ne, Te, Ti and velocity, and limits are defined for the errors introduced by assumptions about ion composition and by changes in the transmitted power and system constant. Simple data-cleaning criteria are employed to eliminate problems due to coherent signals and large background noise levels. The observed variations in plasma densities, temperatures and velocities are interpreted in terms of supporting data from ISEE-3 and local riometers and magnetometers. Both field-aligned and field-perpendicular plasma flows at Tromsø showed effects of the northward IMF: convection was slow and irregular and field-aligned flow profiles were characteristic of steady-state polar wind outflow with flux of order 1012 m−2 s−1. This period followed a strongly southward IMF which had triggered a substorm. The substorm gave enhanced convection, with a swing to equatorward flow and large (5 × 1012 m−2 s−1), steady-state field-aligned fluxes, leading to the possibility of O+ escape into the magnetosphere. The apparent influence of the IMF over both field-perpendicular and field-aligned flows is explained in terms of the cross-cap potential difference and the location of the auroral oval.