THE DETECTION OF HOT ETHANOL IN G34.3+0.15

Hot molecular cores in star-forming regions are known to have gas-phase chemical compositions determined by the evaporation of material from the icy mantles of interstellar grains, followed by subsequent reactions in the gas phase. Current models suggest that the evaporated material is rich in hydrogenated species, such as water, methane and methanol. In this paper, we report the detection of 14 rotational transitions of ethanol in the submillimetre spectrum of the molecular cloud associated with the ultra-compact H II region G34.3+0.15. We derive a rotation temperature of 125 K and a beam-averaged column density of 2.0x10(15) cm(-2), corresponding to a fractional abundance on the order of 4x10(-9). This large abundance, which is a lower limit due to the likelihood of beam dilution, cannot be made by purely gas-phase processes, and we conclude that the ethanol must be formed efficiently in the grain surface chemistry. Since it has been argued previously that methanol is formed via surface chemistry, it appears that alcohol formation may be a natural by-product of surface reactions.

THE CHEMISTRY OF PHOSPHORUS IN HOT MOLECULAR CORES

Hot molecular cores in star-forming regions are known to have gas-phase chemical compositions determined by the material evaporated from the icy mantles of interstellar grains, followed by subsequent reactions in the gas phase. Current models suggest that the evaporated material is rich in hydrogenated species. In this paper, we consider the chemistry induced in a hot core by the release of phosphine, PH3 from interstellar grains. We find that PH3 is rapidly destroyed by a series of reactions with atomic hydrogen and is converted, within 10(4) yr, into atomic P, and PO and PN, with P atoms being the most abundant species. Other P-bearing molecules can be formed in the hot gas, but on time-scales that are long compared to those of the hot cores.

ON THE ORIGIN OF NH IN DIFFUSE INTERSTELLAR CLOUDS

The observation by Meyer & Roch of NH in the interstellar clouds towards zeta Per and HD 27778 cannot be explained with conventional gas-phase chemistry models. A simple non-equilibrium model for the zeta Per cloud, which incorporates the grain-surface production of NH and OH or, alternatively, NH3 and H2O, is able to reproduce the abundances of all observed species (except CH+) quite accurately. Moreover, chemical models which include grain-surface reactions can reproduce the observed abundance not only of NH but also of CN, which is efficiently formed at low temperatures, initiated by the reaction of NH with C+. Pure gas-phase models and cloud interface models, in which NH and CH+ are formed in a warm and tenuous environment, fail to explain the observed high abundance of CN. Hence the observation of NH in zeta Per and HD 27778 provides evidence for the presence of grain-surface reactions leading to molecules other than H-2. It is predicted that NH2 and NH3 should have abundances not much below that of NH if NH3 instead of NH is formed on grains. With or without surface reactions, the column densities of H2O and C2H are expected to be about 10(13) cm-2, and these molecules may be detectable in the zeta Per cloud.

Multiferroic domain dynamics in strained strontium titanate

Multiferroicity can be induced in strontium titanate by applying biaxial strain. Using optical second harmonic generation, we report a transition from 4/mmm to the ferroelectric mm2 phase, followed by a transition to a ferroelastic-ferroelectric mm2 phase in a strontium titanate thin film. Piezoelectric force microscopy is used to study ferroelectric domain switching. Second harmonic generation, combined with phase-field modeling, is used to reveal the mechanism of coupled ferroelectric-ferroelastic domain wall motion. These studies have relevance to multiferroics with coupled polar and axial phenomena.

Nanoscale Origins of Nonlinear Behavior in Ferroic Thin Films

The nonlinear response of a ferroic to an applied field has been studied through the phenomenological Rayleigh Law for over a hundred years. Yet, despite this, the fundamental physical mechanisms at the nanoscale that lead to macroscopic Rayleigh behavior have remained largely elusive, and experimental evidence at small length scales is limited. Here, it is shown using a combination of scanning probe techniques and phase field modeling, that nanoscale piezoelectric response in prototypical Pb(Zr,Ti)O3 films appears to follow a distinctly non-Rayleigh regime. Through statistical analysis, it is found that an averaging of local responses can lead directly to Rayleigh-like behavior of the strain on a macroscale. Phase-field modeling confirms the twist of the ferroelastic interface is key in enhancing piezoelectric response. The studies shed light on the nanoscale origins of nonlinear behavior in disordered ferroics.

Manipulating Ferroelectric Domains in Nanostructures Under Electron Beams

Freestanding BaTiO3 nanodots exhibit domain structures characterized by distinct quadrants of ferroelastic 90 domains in transmission electron microscopy (TEM) observations. These differ significantly from flux-closure domain patterns in the same systems imaged by piezoresponse force microscopy. Based upon a series of phase field simulations of BaTiO3 nanodots, we suggest that the TEM patterns result from a radial electric field arising from electron beam charging of the nanodot. For sufficiently large charging, this converts flux-closure domain patterns to quadrant patterns with radial net polarizations. Not only does this explain the puzzling patterns that have been observed in TEM studies of ferroelectric nanodots, but also suggests how to manipulate ferroelectric domain patterns via electron beams.

Local probing of ferroelectric and ferroelastic switching through stress-mediated piezoelectric spectroscopy

Strain effects have a significant role in mediating classic ferroelectric behavior such as polarization switching and domain wall dynamics. These effects are of critical relevance if the ferroelectric order parameter is coupled to strain and is therefore, also ferroelastic. Here, switching spectroscopy piezoresponse force microscopy (SS-PFM) is combined with control of applied tip pressure to exert direct control over the ferroelastic and ferroelectric switching events, a modality otherwise unattainable in traditional PFM. As a proof of concept, stress-mediated SS-PFM is applied toward the study of polarization switching events in a lead zirconate titanate thin film, with a composition near the morphotropic phase boundary with co-existing rhombohedral and tetragonal phases. Under increasing applied pressure, shape modification of local hysteresis loops is observed, consistent with a reduction in the ferroelastic domain variants under increased pressure. These experimental results are further validated by phase field simulations. The technique can be expanded to explore more complex electromechanical responses under applied local pressure, such as probing ferroelectric and ferroelastic piezoelectric nonlinearity as a function of applied pressure, and electro-chemo-mechanical response through electrochemical strain microscopy.

Ferroelectric Vortices and Related Configurations

This chapter discusses that the theoretical studies, using both atomistic and phenomenological approaches, have made clear predictions about the existence and behaviour of ferroelectric (FE) vortices. Effective Hamiltonians can be implemented within both Monte Carlo (MC) and molecular dynamics (MD) simulations. In contrast to the effective Hamiltonian method, which is atomistic in nature, the phase field method employs a continuum approach, in which the polarization field is the order parameter. Properties of FE nanostructures are largely governed by the existence of a depolarization field, which is much stronger than the demagnetization field in magnetic nanosystems. The topological patterns seen in rare earth manganites are often referred to as vortices and yet this claim never seems to be explicitly justified. By inspection, the form of a vortex structure is such that there is a continuous rotation in the orientation of dipole vectors around the singularity at the centre of the vortex.

Parameters determination for the coupled electromagnetic transformer model, with particular reference to winding faults studies

In the past few years, a considerable research effort has been devoted to the development of transformer digital models in order to simulate its behaviour under transient and abnormal operating conditions. Although many three-phase transformer models have been presented in the literature, there is a surprisingly lack of studies regarding the incorporation of winding faults. This paper presents a coupled electromagnetic transformer model for the study of winding inter-turn short-circuits. Particular attention will be given to the model parameters determination, for both healthy and faulty operating conditions. Experimental and simulation test results are presented in the paper, demonstrating the adequacy of the model as well as the methodologies for the parameters determination.

Inverse problems in dynamic cognitive modeling

Inverse problems for dynamical system models of cognitive processes comprise the determination of synaptic weight matrices or kernel functions for neural networks or neural/dynamic field models, respectively. We introduce dynamic cognitive modeling as a three tier top-down approach where cognitive processes are first described as algorithms that operate on complex symbolic data structures. Second, symbolic expressions and operations are represented by states and transformations in abstract vector spaces. Third, prescribed trajectories through representation space are implemented in neurodynamical systems. We discuss the Amari equation for a neural/dynamic field theory as a special case and show that the kernel construction problem is particularly ill-posed. We suggest a Tikhonov-Hebbian learning method as regularization technique and demonstrate its validity and robustness for basic examples of cognitive computations.

Axonal velocity distributions in neural feld equations

By modelling the average activity of large neuronal populations, continuum mean field models (MFMs) have become an increasingly important theoretical tool for understanding the emergent activity of cortical tissue. In order to be computationally tractable, long-range propagation of activity in MFMs is often approximated with partial differential equations (PDEs). However, PDE approximations in current use correspond to underlying axonal velocity distributions incompatible with experimental measurements. In order to rectify this deficiency, we here introduce novel propagation PDEs that give rise to smooth unimodal distributions of axonal conduction velocities. We also argue that velocities estimated from fibre diameters in slice and from latency measurements, respectively, relate quite differently to such distributions, a significant point for any phenomenological description. Our PDEs are then successfully fit to fibre diameter data from human corpus callosum and rat subcortical white matter. This allows for the first time to simulate long-range conduction in the mammalian brain with realistic, convenient PDEs. Furthermore, the obtained results suggest that the propagation of activity in rat and human differs significantly beyond mere scaling. The dynamical consequences of our new formulation are investigated in the context of a well known neural field model. On the basis of Turing instability analyses, we conclude that pattern formation is more easily initiated using our more realistic propagator. By increasing characteristic conduction velocities, a smooth transition can occur from self-sustaining bulk oscillations to travelling waves of various wavelengths, which may influence axonal growth during development. Our analytic results are also corroborated numerically using simulations on a large spatial grid. Thus we provide here a comprehensive analysis of empirically constrained activity propagation in the context of MFMs, which will allow more realistic studies of mammalian brain activity in the future.

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.

Ice supersaturation and the potential for contrail formation in a changing climate

Ice supersaturation (ISS) in the upper troposphere and lower stratosphere is important for the formation of cirrus clouds and long-lived contrails. Cold ISS (CISS) regions (taken here to be ice-supersaturated regions with temperature below 233 K) are most relevant for contrail formation.We analyse projected changes to the 250 hPa distribution and frequency of CISS regions over the 21st century using data from the Representative Concentration Pathway 8.5 simulations for a selection of Coupled Model Intercomparison Project Phase 5 models. The models show a global-mean, annual-mean decrease in CISS frequency by about one-third, from 11 to 7% by the end of the 21st century, relative to the present-day period 1979–2005. Changes are analysed in further detail for three subregions where air traffic is already high and increasing (Northern Hemisphere mid-latitudes) or expected to increase (tropics and Northern Hemisphere polar regions). The largest change is seen in the tropics, where a reduction of around 9 percentage points in CISS frequency by the end of the century is driven by the strong warming of the upper troposphere. In the Northern Hemisphere mid-latitudes the multi-model-mean change is an increase in CISS frequency of 1 percentage point; however the sign of the change is dependent not only on the model but also on latitude and season. In the Northern Hemisphere polar regions there is an increase in CISS frequency of 5 percentage points in the annual mean. These results suggest that, over the 21st century, climate change may have large impacts on the potential for contrail formation; actual changes to contrail cover will also depend on changes to the volume of air traffic, aircraft technology and flight routing.

Bayesian model comparison with un-normalised likelihoods

Models for which the likelihood function can be evaluated only up to a parameter-dependent unknown normalizing constant, such as Markov random field models, are used widely in computer science, statistical physics, spatial statistics, and network analysis. However, Bayesian analysis of these models using standard Monte Carlo methods is not possible due to the intractability of their likelihood functions. Several methods that permit exact, or close to exact, simulation from the posterior distribution have recently been developed. However, estimating the evidence and Bayes’ factors for these models remains challenging in general. This paper describes new random weight importance sampling and sequential Monte Carlo methods for estimating BFs that use simulation to circumvent the evaluation of the intractable likelihood, and compares them to existing methods. In some cases we observe an advantage in the use of biased weight estimates. An initial investigation into the theoretical and empirical properties of this class of methods is presented. Some support for the use of biased estimates is presented, but we advocate caution in the use of such estimates.

Quantum current algebra for the AdS(5) x S(5) superstring

The sigma model describing the dynamics of the superstring in the AdS(5) x S(5) background can be constructed using the coset PSU(2, 2 vertical bar 4)/SO(4, 1) x SO(5). A basic set of operators in this two dimensional conformal field theory is composed by the left invariant currents. Since these currents are not (anti) holomorphic, their OPE`s is not determined by symmetry principles and its computation should be performed perturbatively. Using the pure spinor sigma model for this background, we compute the one-loop correction to these OPE`s. We also compute the OPE`s of the left invariant currents with the energy momentum tensor at tree level and one loop.