Testing isotherm models and recovering empirical relationships for adsorption in microporous carbons using virtual carbon models and grand canonical Monte Carlo simulations

Using the plausible model of activated carbon proposed by Harris and co-workers and grand canonical Monte Carlo simulations, we study the applicability of standard methods for describing adsorption data on microporous carbons widely used in adsorption science. Two carbon structures are studied, one with a small distribution of micropores in the range up to 1 nm, and the other with micropores covering a wide range of porosity. For both structures, adsorption isotherms of noble gases (from Ne to Xe), carbon tetrachloride and benzene are simulated. The data obtained are considered in terms of Dubinin-Radushkevich plots. Moreover, for benzene and carbon tetrachloride the temperature invariance of the characteristic curve is also studied. We show that using simulated data some empirical relationships obtained from experiment can be successfully recovered. Next we test the applicability of Dubinin's related models including the Dubinin-Izotova, Dubinin-Radushkevich-Stoeckli, and Jaroniec-Choma equations. The results obtained demonstrate the limits and applications of the models studied in the field of carbon porosity characterization.

Imaging the atomic structure of activated carbon

The precise atomic structure of activated carbon is unknown, despite its huge commercial importance in the purification of air and water. Diffraction methods have been extensively applied to the study of microporous carbons, but cannot provide an unequivocal identification of their structure. Here we show that the structure of a commercial activated carbon can be imaged directly using aberration-corrected transmission electron microscopy. Images are presented both of the as-produced carbon and of the carbon following heat treatment at 2000 degrees C. In the 2000 degrees C carbon clear evidence is found for the presence of pentagonal rings, suggesting that the carbons have a fullerene-related structure. Such a structure would help to explain the properties of activated carbon, and would also have important implications for the modelling of adsorption on microporous carbons.

Hyper-parallel tempering Monte Carlo simulations of Ar adsorption in new models of microporous non-graphitizing activated carbon: effect of microporosity

The adsorption of gases on microporous carbons is still poorly understood, partly because the structure of these carbons is not well known. Here, a model of microporous carbons based on fullerene- like fragments is used as the basis for a theoretical study of Ar adsorption on carbon. First, a simulation box was constructed, containing a plausible arrangement of carbon fragments. Next, using a new Monte Carlo simulation algorithm, two types of carbon fragments were gradually placed into the initial structure to increase its microporosity. Thirty six different microporous carbon structures were generated in this way. Using the method proposed recently by Bhattacharya and Gubbins ( BG), the micropore size distributions of the obtained carbon models and the average micropore diameters were calculated. For ten chosen structures, Ar adsorption isotherms ( 87 K) were simulated via the hyper- parallel tempering Monte Carlo simulation method. The isotherms obtained in this way were described by widely applied methods of microporous carbon characterisation, i. e. Nguyen and Do, Horvath - Kawazoe, high- resolution alpha(a)s plots, adsorption potential distributions and the Dubinin - Astakhov ( DA) equation. From simulated isotherms described by the DA equation, the average micropore diameters were calculated using empirical relationships proposed by different authors and they were compared with those from the BG method.

Inferring convective responses to El Niño with atmospheric electricity measurements at Shetland

Pacific ocean temperature anomalies associated with the El Niño–Southern Oscillation (ENSO) modulate atmospheric convection and hence thunderstorm electrification. The generated current flows globally via the atmospheric electric circuit, which can be monitored anywhere on Earth. Atmospheric electricity measurements made at Shetland (in Scotland) display a mean global circuit response to ENSO that is characterized by strengthening during 'El Niño' conditions, and weakening during 'La Niña' conditions. Examining the hourly varying response indicates that a potential gradient (PG) increase around noon UT is likely to be associated with a change in atmospheric convection and resultant lightning activity over equatorial Africa and Eastern Asia. A secondary increase in PG just after midnight UT can be attributed to more shower clouds in the central Pacific ocean during an 'El Niño'.

The Dirichlet-to-Neumann map for the elliptic sine Gordon

We analyse the Dirichlet problem for the elliptic sine Gordon equation in the upper half plane. We express the solution $q(x,y)$ in terms of a Riemann-Hilbert problem whose jump matrix is uniquely defined by a certain function $b(\la)$, $\la\in\R$, explicitly expressed in terms of the given Dirichlet data $g_0(x)=q(x,0)$ and the unknown Neumann boundary value $g_1(x)=q_y(x,0)$, where $g_0(x)$ and $g_1(x)$ are related via the global relation $\{b(\la)=0$, $\la\geq 0\}$. Furthermore, we show that the latter relation can be used to characterise the Dirichlet to Neumann map, i.e. to express $g_1(x)$ in terms of $g_0(x)$. It appears that this provides the first case that such a map is explicitly characterised for a nonlinear integrable {\em elliptic} PDE, as opposed to an {\em evolution} PDE.

Photo-induced phase transitions in azobenzene-doped liquid crystals

The effect of irradiation (UV-visible light) on a nematic liquid crystal doped with a photoactive azobenzene derivative was investigated. The selective irradiation results in either an E implies Z or Z implies E isomerization of the azobenzene unit. The effect of the isomerization is to cause a reversible depression of the liquid crystal to isotropic (LC implies l) phase transition temperature of the doped mixture, which can be monitored optically as an isothermal phase transition. This depression also results in a biphasic liquid crystal+isotropic region which is discussed. The authors investigate the cause and magnitude of the phase depression as a function of the amount of doped 4-butyl-4'-methoxyazobenzene (photoactive unit) in 4-cyano-4'-n-pentylbiphenyl (liquid crystal unit), and as a function of the percentage conversion of E implies Z (caused by isomerization) in the azobenzene. The photostationary state of the doped mixtures achieved by Z implies E isomerization is considered and its effect upon the transition temperature of the mixture and response time of the system is discussed. They discuss the implications of the photostationary state with regards to the reversibility of the photo-induced phase transition and hence potential applications.

Local gratings due to angular hole burning in a photorefractive polymer

We investigate the processes involved in writing real-time holographic gratings in a photorefractive polymer (PRP) that incorporates an azo-dye. In such systems there may be gratings due to mechanisms associated with trans–cis isomerization (angular hole burning (AHB) and/or angular redistribution), which appear in addition to those arising from the photorefractive (PR) effect. The work presented here helps to understand the interactions which may occur between these different gratings. The formation of local gratings due to mechanisms associated with photoisomerization is studied, in a new PRP based on the photoconductor poly(N-vinylcarbazole):2, 4, 7-trinitro-9-fluorenone, plasticized with N-ethylcarbazole. The polymer includes the azo-dye 4-nitro-4'-pentyloxy-azobenzene and we observe both PR and photoisomerization gratings. The gratings are shown to be both polarization-sensitive and reversible. The presence of the photoisomerization gratings (which diffract almost as strongly as the PR gratings) significantly affects the field-dependent diffractive behaviour of the composite. A measurement of the lifetime of the cis state is made (τcis = 38 s) using photoinduced dichroism. This is close to the decay time constant of the local gratings (τdecay = 42 s), and it is suggested that the local grating mechanism is AHB of the azo-dye. This is the first time (to the knowledge of the authors) that a local grating due to AHB has been demonstrated in a PRP.

Modelling the structure of a polymer glass poly(methylmethacrylate) through neutron scattering experiments

Determination of the local structure of a polymer glass by scattering methods is complex due to the number of spatial and orientational correlations, both from within the polymer chain (intrachain) and between neighbouring chains (interchain), from which the scattering arises. Recently considerable advances have been made in the structural analysis of relatively simple polymers such as poly(ethylene) through the use of broad Q neutron scattering data tightly coupled to atomistic modelling procedures. This paper presents the results of an investigation into the use of these procedures for the analysis of the local structure of a-PMMA which is chemically more complex with a much greater number of intrachain structural parameters. We have utilised high quality neutron scattering data obtained using SANDALS at ISIS coupled with computer models representing both the single chain and bulk polymer system. Several different modelling approaches have been explored which encompass such techniques as Reverse Monte Carlo refinement and energy minimisation and their relative merits and successes are discussed. These different approaches highlight structural parameters which any realistic model of glassy atactic PMMA must replicate.

A robust controller design method for feedback substitution schemes using genetic algorithms

Controllers for feedback substitution schemes demonstrate a trade-off between noise power gain and normalized response time. Using as an example the design of a controller for a radiometric transduction process subjected to arbitrary noise power gain and robustness constraints, a Pareto-front of optimal controller solutions fulfilling a range of time-domain design objectives can be derived. In this work, we consider designs using a loop shaping design procedure (LSDP). The approach uses linear matrix inequalities to specify a range of objectives and a genetic algorithm (GA) to perform a multi-objective optimization for the controller weights (MOGA). A clonal selection algorithm is used to further provide a directed search of the GA towards the Pareto front. We demonstrate that with the proposed methodology, it is possible to design higher order controllers with superior performance in terms of response time, noise power gain and robustness.

Disentangling multi-level systems: averaging, correlations and memory

We consider two weakly coupled systems and adopt a perturbative approach based on the Ruelle response theory to study their interaction. We propose a systematic way of parameterizing the effect of the coupling as a function of only the variables of a system of interest. Our focus is on describing the impacts of the coupling on the long term statistics rather than on the finite-time behavior. By direct calculation, we find that, at first order, the coupling can be surrogated by adding a deterministic perturbation to the autonomous dynamics of the system of interest. At second order, there are additionally two separate and very different contributions. One is a term taking into account the second-order contributions of the fluctuations in the coupling, which can be parameterized as a stochastic forcing with given spectral properties. The other one is a memory term, coupling the system of interest to its previous history, through the correlations of the second system. If these correlations are known, this effect can be implemented as a perturbation with memory on the single system. In order to treat this case, we present an extension to Ruelle's response theory able to deal with integral operators. We discuss our results in the context of other methods previously proposed for disentangling the dynamics of two coupled systems. We emphasize that our results do not rely on assuming a time scale separation, and, if such a separation exists, can be used equally well to study the statistics of the slow variables and that of the fast variables. By recursively applying the technique proposed here, we can treat the general case of multi-level systems.

Relevance of sampling schemes in light of Ruelle's linear response theory

We reconsider the theory of the linear response of non-equilibrium steady states to perturbations. We �rst show that by using a general functional decomposition for space-time dependent forcings, we can de�ne elementary susceptibilities that allow to construct the response of the system to general perturbations. Starting from the de�nition of SRB measure, we then study the consequence of taking di�erent sampling schemes for analysing the response of the system. We show that only a speci�c choice of the time horizon for evaluating the response of the system to a general time-dependent perturbation allows to obtain the formula �rst presented by Ruelle. We also discuss the special case of periodic perturbations, showing that when they are taken into consideration the sampling can be �ne-tuned to make the de�nition of the correct time horizon immaterial. Finally, we discuss the implications of our results in terms of strategies for analyzing the outputs of numerical experiments by providing a critical review of a formula proposed by Reick.

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.

Sting jets in intense winter North-Atlantic windstorms

Extratropical cyclones dominate autumn and winter weather over western Europe. The strongest cyclones, often termed windstorms, have a large socio-economic impact due to the strong surface winds and associated storm surges in coastal areas. Here we show that sting jets are a common feature of windstorms; up to a third of the 100 most intense North Atlantic winter windstorms over the last two decades satisfy conditions for sting jets. The sting jet is a mesoscale descending airstream that can cause strong near-surface winds in the dry slot of the cyclone, a region not usually associated with strong winds. Despite their localized transient nature these sting jets can cause significant damage, a prominent example being the storm that devastated southeast England on 16 October 1987. We present the first regional climatology of windstorms with sting jets. Previously analysed sting jet cases appear to have been exceptional in their track over northwest Europe rather than in their strength.

Vocabulary is important for some, but not all reading skills

Although there is evidence for a close link between the development of oral vocabulary and reading comprehension, less clear is whether oral vocabulary skills relate to the development of word-level reading skills. This study investigated vocabulary and literacy in 81 children aged 8 to 10 years. In regression analyses, vocabulary accounted for unique variance in exception word reading and reading comprehension, but not text reading accuracy, decoding, or regular word reading. Consistent with these data, children with poor reading comprehension exhibited oral vocabulary weaknesses and read fewer exception words correctly. These findings demonstrate that oral vocabulary is associated with some, but not all, reading skills. Results are discussed in terms of current models of reading development.

The role of self-teaching in learning orthographic and semantic aspects of new words

This study explores how children learn the meaning (semantics) and spelling patterns (orthography) of novel words encountered in story context. English-speaking children (N = 88) aged 7 to 8 years read 8 stories and each story contained 1 novel word repeated 4 times. Semantic cues were provided by the story context such that children could infer the meaning of the word (specific context) or the category that the word belonged to (general context). Following story reading, posttests indicated that children showed reliable semantic and orthographic learning. Decoding was the strongest predictor of orthographic learning, indicating that self-teaching via phonological recoding was important for this aspect of word learning. In contrast, oral vocabulary emerged as the strongest predictor of semantic learning.