Strong-segregation limit of the self-consistent field theory for diblock copolymer melts

The self-consistent field theory (SCFT) introduced by Helfand for diblock copolymer melts is expected to converge to the strong-segregation theory (SST) of Semenov in the asymptotic limit, $\chi N \rightarrow \infty$. However, past extrapolations of the lamellar/cylinder and cylinder/sphere phase boundaries, within the standard unit-cell approximation, have cast some doubts on whether or not this is actually true. Here we push the comparison further by extending the SCFT calculations to $\chi N = 512,000$, by accounting for exclusion zones in the coronae of the cylindrical and spherical unit cells, and by examining finite-segregation corrections to SST. In doing so, we provide the first compelling evidence that SCFT does indeed reduce to SST.

Unit-cell approximation for diblock−copolymer brushes grafted to spherical particles

This paper examines the equilibrium phase behavior of thin diblock-copolymer films tethered to a spherical core, using numerical self-consistent field theory (SCFT). The computational cost of the calculation is greatly reduced by implementing the unit-cell approximation (UCA) routinely used in the study of bulk systems. This provides a tremendous reduction in computational time, permitting us to map out the phase behavior more extensively and allowing us to consider far larger particles. The main consequence of the UCA is that it omits packing frustration, but evidently the effect is minor for large particles. On the other hand, when the particles are small, the UCA calculation can be readily followed up with the full SCFT, the comparison to which conveniently allows one to quantitatively assess the effect of packing frustration.

A variational method to retrieve the extinction profile in liquid clouds using multiple field-of-view lidar

Liquid clouds play a profound role in the global radiation budget but it is difficult to remotely retrieve their vertical profile. Ordinary narrow field-of-view (FOV) lidars receive a strong return from such clouds but the information is limited to the first few optical depths. Wideangle multiple-FOV lidars can isolate radiation scattered multiple times before returning to the instrument, often penetrating much deeper into the cloud than the singly-scattered signal. These returns potentially contain information on the vertical profile of extinction coefficient, but are challenging to interpret due to the lack of a fast radiative transfer model for simulating them. This paper describes a variational algorithm that incorporates a fast forward model based on the time-dependent two-stream approximation, and its adjoint. Application of the algorithm to simulated data from a hypothetical airborne three-FOV lidar with a maximum footprint width of 600m suggests that this approach should be able to retrieve the extinction structure down to an optical depth of around 6, and total opticaldepth up to at least 35, depending on the maximum lidar FOV. The convergence behavior of Gauss-Newton and quasi-Newton optimization schemes are compared. We then present results from an application of the algorithm to observations of stratocumulus by the 8-FOV airborne “THOR” lidar. It is demonstrated how the averaging kernel can be used to diagnose the effective vertical resolution of the retrieved profile, and therefore the depth to which information on the vertical structure can be recovered. This work enables exploitation of returns from spaceborne lidar and radar subject to multiple scattering more rigorously than previously possible.

Inverse problems in neural field theory

We study inverse problems in neural field theory, i.e., the construction of synaptic weight kernels yielding a prescribed neural field dynamics. We address the issues of existence, uniqueness, and stability of solutions to the inverse problem for the Amari neural field equation as a special case, and prove that these problems are generally ill-posed. In order to construct solutions to the inverse problem, we first recast the Amari equation into a linear perceptron equation in an infinite-dimensional Banach or Hilbert space. In a second step, we construct sets of biorthogonal function systems allowing the approximation of synaptic weight kernels by a generalized Hebbian learning rule. Numerically, this construction is implemented by the Moore–Penrose pseudoinverse method. We demonstrate the instability of these solutions and use the Tikhonov regularization method for stabilization and to prevent numerical overfitting. We illustrate the stable construction of kernels by means of three instructive examples.

Detection of Chlamydia psittaci DNA in avian clinical samples by polymerase chain reaction

A polymerase chain reaction (PCR) assay was developed to detect Chlamydia psittaci DNA in faeces and tissue samples from avian species. Primers were designed to amplify a 264 bp product derived from part of the 5' non-translated region and part of the coding region of the ompA gene which encodes the major outer membrane protein. Amplified sequences were confirmed by Southern hybridization using an internal probe. The sensitivity of the combined assay was found to be between 60 to 600 fg of chlamydial DNA (approximately 6 to 60 genome copies). The specificity of the assay was confirmed since PCR product was not obtained from samples containing several serotypes of C. trachomatis, strains of C. pneumoniae, the type strain of C. pecorum, nor from samples containing microorganisms commonly found in the avian gut flora. In this study, 404 avian faeces and 141 avian tissue samples received by the Central Veterinary Laboratory over a 6 month period were analysed by PCR, antigen detection ELISA and where possible, cell culture isolation. PCR performed favourably compared with ELISA and cell culture, or with ELISA alone. The PCR assay was especially suited to the detection of C. psittaci DNA in avian faeces samples. The test was also useful when applied to tissue samples from small contact birds associated with a case of human psittacosis where ELISA results were negative and chlamydial isolation was a less favourable method due to the need for rapid diagnosis.

Metabifurcation analysis of a mean field model of the cortex

Mean field models (MFMs) of cortical tissue incorporate salient, average features of neural masses in order to model activity at the population level, thereby linking microscopic physiology to macroscopic observations, e.g., with the electroencephalogram (EEG). One of the common aspects of MFM descriptions is the presence of a high-dimensional parameter space capturing neurobiological attributes deemed relevant to the brain dynamics of interest. We study the physiological parameter space of a MFM of electrocortical activity and discover robust correlations between physiological attributes of the model cortex and its dynamical features. These correlations are revealed by the study of bifurcation plots, which show that the model responses to changes in inhibition belong to two archetypal categories or “families”. After investigating and characterizing them in depth, we discuss their essential differences in terms of four important aspects: power responses with respect to the modeled action of anesthetics, reaction to exogenous stimuli such as thalamic input, and distributions of model parameters and oscillatory repertoires when inhibition is enhanced. Furthermore, while the complexity of sustained periodic orbits differs significantly between families, we are able to show how metamorphoses between the families can be brought about by exogenous stimuli. We here unveil links between measurable physiological attributes of the brain and dynamical patterns that are not accessible by linear methods. They instead emerge when the nonlinear structure of parameter space is partitioned according to bifurcation responses. We call this general method “metabifurcation analysis”. The partitioning cannot be achieved by the investigation of only a small number of parameter sets and is instead the result of an automated bifurcation analysis of a representative sample of 73,454 physiologically admissible parameter sets. Our approach generalizes straightforwardly and is well suited to probing the dynamics of other models with large and complex parameter spaces.

Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular, we are able to treat “patchy” connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a “lattice-directed” traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs.

Monte Carlo field-theoretic simulations for melts of symmetric diblock copolymer

Monte Carlo field-theoretic simulations (MCFTS) are performed on melts of symmetric diblock copolymer for invariant polymerization indexes extending down to experimentally relevant values of N̅ ∼ 10^4. The simulations are performed with a fluctuating composition field, W_−(r), and a pressure field, W_+(r), that follows the saddle-point approximation. Our study focuses on the disordered-state structure function, S(k), and the order−disorder transition (ODT). Although shortwavelength fluctuations cause an ultraviolet (UV) divergence in three dimensions, this is readily compensated for with the use of an effective Flory−Huggins interaction parameter, χ_e. The resulting S(k) matches the predictions of renormalized one-loop (ROL) calculations over the full range of χ_eN and N̅ examined in our study, and agrees well with Fredrickson−Helfand (F−H) theory near the ODT. Consistent with the F−H theory, the ODT is discontinuous for finite N̅ and the shift in (χ_eN)_ODT follows the predicted N̅^−1/3 scaling over our range of N̅.

Reconnection at the high-latitude magnetopause during northward interplanetary magnetic field conditions

The Polar spacecraft had a prolonged encounter with the high-latitude dayside magnetopause on May 29, 1996. This encounter with the magnetopause occurred when the interplanetary magnetic field was directed northward. From the three-dimensional electron and ion distribution functions measured by the Hydra instrument, it has been possible to identify nearly all of the distinct boundary layer regions associated with high-latitude reconnection. The regions that have been identified are (1) the cusp; (2) the magnetopause current layer; (3) magnetosheath field lines that have interconnected in only the Northern Hemisphere; (4) magnetosheath field lines that have interconnected in only the Southern Hemisphere; (5) magnetosheath field lines that have interconnected in both the Northern and Southern Hemispheres; (6) magnetosheath that is disconnected from the terrestrial magnetic field; and (7) high-latitude plasma sheet field lines that are participating in magnetosheath reconnection. Reconnection over this time period was occurring at high latitudes over a broad local-time extent, interconnecting the magnetosheath and lobe and/or plasma sheet field lines in both the Northern and Southern Hemispheres. Newly closed boundary layer field lines were observed as reconnection occur-red first at high latitudes in one hemisphere and then later in the other. These observations establish the location of magnetopause reconnection during these northward interplanetary magnetic field conditions as being at high latitudes, poleward of the cusp, and further reinforce the general interpretation of electron and ion phase space density signatures as indicators of magnetic reconnection and boundary layer formation.

Cusp ion steps, field-aligned currents and poleward moving auroral forms

We predict the field-aligned currents around cusp ion steps produced by pulsed reconnection between the geomagnetic field and an interplanetary magnetic field (IMF) with a B-Y component that is large in magnitude. For B-Y > 0, patches of newly opened flux move westward and eastward in the Northern and Southern Hemispheres, respectively, under the influence of the magnetic curvature force. These flow directions are reversed for B-Y < 0. The speed of this longitudinal motion initially grows with elapsed time since reconnection, but then decays as the newly opened field lines straighten. We predict sheets of field-aligned current on the boundaries between the patches produced by successive reconnection pulses, associated with the difference in the speeds of their longitudinal motion. For low elapsed times since reconnection, near the equatorward edge of the cusp region where the field lines are accelerating, the field-aligned current sheets will be downward or upward in both hemispheres for positive or negative IMF B-Y, respectively. At larger elapsed times since reconnection, as events slow and evolve from the cusp into the mantle region, these field-aligned current directions will be reversed. Observations by the Polar spacecraft on August 26,1998, show the predicted upward current sheets at steps seen in the mantle region for IMF B-Y > 0. Mapped into the ionosphere, the steps coincide with poleward moving events seen by the CUTLASS HF radar. The mapped location of the largest step also coincides with a poleward moving arc seen by the UVI imager on Polar. We show that the arc is consistent with a region of upward field-aligned current that has become unstable, such that a potential drop of about 1 kV formed below the spacecraft. The importance of these observations is that they confirm that the poleward moving events, as seen by the HF radar and the UV imager, are due to pulsed magnetopause reconnection. Milan et al. [2000] noted that the great longitudinal extent of these events means that the required reconnection pulses would have contributed almost all the voltage placed across the magnetosphere at this time. The observations also show that auroral arcs can form on open field lines in response to the pulsed application of voltage at the magnetopause.

Spatial and temporal variability of the concentration field from localized releases in a regular building array

Spatial and temporal fluctuations in the concentration field from an ensemble of continuous point-source releases in a regular building array are analyzed from data generated by direct numerical simulations. The release is of a passive scalar under conditions of neutral stability. Results are related to the underlying flow structure by contrasting data for an imposed wind direction of 0 deg and 45 deg relative to the buildings. Furthermore, the effects of distance from the source and vicinity to the plume centreline on the spatial and temporal variability are documented. The general picture that emerges is that this particular geometry splits the flow domain into segments (e.g. “streets” and “intersections”) in each of which the air is, to a first approximation, well mixed. Notable exceptions to this general rule include regions close to the source, near the plume edge, and in unobstructed channels when the flow is aligned. In the oblique (45 deg) case the strongly three-dimensional nature of the flow enhances mixing of a scalar within the canopy leading to reduced temporal and spatial concentration fluctuations within the plume core. These fluctuations are in general larger for the parallel flow (0 deg) case, especially so in the long unobstructed channels. Due to the more complex flow structure in the canyon-type streets behind buildings, fluctuations are lower than in the open channels, though still substantially larger than for oblique flow. These results are relevant to the formulation of simple models for dispersion in urban areas and to the quantification of the uncertainties in their predictions.

Solvent Effects in Chemical Processes. Water-Assisted Proton Transfer Reaction of Pterin in Aqueous Environment

Pterins are members of a family of heterocyclic compounds present in a wide variety of biological systems and may exist in two forms, corresponding to an acid and a basic tautomer. In this work, the proton transfer reaction between these tautomeric forms was investigated in the gas phase and in aqueous solution. In gas phase, the intramolecular mechanism was carried out for die isolated pterin by quantum mechanical second-order Moller-Plesset Perturbation theory (MP2/aug-cc-pVDZ) calculations and it indicates that the acid form is more stable than the basic form by -1.4 kcal/mol with a barrier of 34.2 kcal/mol with respect to the basic form. In aqueous solution, the role of the water molecules in the proton transfer reaction was analyzed in two separated parts, the direct participation of one water molecule in the reaction path, called water-assisted mechanism, and the complementary participation of the aqueous solvation. The water-assisted mechanism was carried out for one pterin-water cluster by quantum mechanical calculations and it indicates that the acid form is still more stable by -3.3 kcal/mol with a drastic reduction of 70% of the barrier, The bulk solution effect on the intramolecular and water-assisted mechanisms was included by free energy perturbation implemented on Monte Carlo simulations. The bulk water effect is found to be substantial and decisive when the reaction path involves the water-assisted mechanism. In this case, the free energy barrier is only 6.7 kcal/mol and the calculated relative Gibbs free energy for the two tautomers is -11.2 kcal/mol. This value is used to calculate the pK(a) value of 8.2 +/- 0.6 that is in excellent agreement with the experimental result of 7.9.

Energy-momentum tensor in thermal strong-field QED with unstable vacuum

Themean value of the one-loop energy-momentum tensor in thermal QED with an electric-like background that creates particles from vacuum is calculated. The problem is essentially different from calculations of effective actions ( similar to the action of Heisenberg-Euler) in backgrounds that respect the stability of vacuum. The role of a constant electric background in the violation of both the stability of vacuum and the thermal character of particle distribution is investigated. Restrictions on the electric field and the duration over which one can neglect the back-reaction of created particles are established.

Gravitational field of a straight cosmic string in the framework of higher-derivative gravity

The gravitational properties of a straight cosmic string are studied in the linear approximation of higher-derivative gravity. These properties are shown to be very different from those found using linearized Einstein gravity: there exists a short range gravitational (anti-gravitational) force in the nonrelativistic limit; in addition, the derection angle of a light ray moving in a plane orthogonal to the string depends on the impact parameter.

Capacitive electrolyte-insulator-semiconductor structures functionalised with a polyelectrolyte/enzyme multilayer: New strategy for enhanced field-effect biosensing

A novel strategy for enhanced field-effect biosensing using capacitive electrolyte-insulator-semiconductor (EIS) structures functionalised with pH-responsive weak polyelectrolyte/enzyme or dendrimer/enzyme multilayers is presented. The feasibility of the proposed approach is exemplarily demonstrated by realising a penicillin biosensor based on a capacitive p-Si-SiO(2) EIS structure functionalised with a poly(allylamine hydrochloride) (PAH)/penicillinase and a poly(amidoamine) dendrimer/penicillinase multilayer. The developed sensors response to changes in both the local pH value near the gate surface and the charge of macromolecules induced via enzymatic reaction, resulting in a higher sensitivity. For comparison, an EIS penicillin biosensor with adsorptively immobilised penicillinase has been also studied. The highest penicillin sensitivity of 100 mV/dec has been observed for the EIS sensor functionalised with the PAH/penicillinase multilayer. The lower and upper detection limit was around 20 mu M and 10 mM, respectively. In addition, an incorporation of enzymes in a multilayer prepared by layer-by-layer technique provides a larger amount of immobilised enzymes per sensor area, reduces enzyme leaching effects and thus, enhances the biosensor lifetime (the loss of penicillin sensitivity after 2 months was 10-12%). (C) 2010 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim