Bidentates versus monodentates in asymmetric hydrogenation catalysis: synergic effects on rate and allosteric effects on enantioselectivity

C-1-Symmetric phosphino/phosphonite ligands are prepared by the reactions of Ph2P(CH2)(2)P(NMe2)(2) with (S)-1,11'-bi-2-naphthol (to give L-A) or (S)-10,10'-bi-9-phenanthrol (to give L-B). Racemic 10,10'-bi-9-phenanthrol is synthesized in three steps from phenanthrene in 44% overall yield. The complexes [PdCl2(L-A,L-B)] (1a,b), [PtCl2(L-A,L-B)] (2a,b), [Rh(cod)(L-A,L-B)]BF4 (3a,b) and [Rh(L-A,L-B)(2)]BF4 (4a,b) are reported and the crystal structure of la has been determined. A P-31 NMR study shows that M, a 1:1 mixture of the monodentates, PMePh2 and methyl monophosphonite L-1a (based on (S)-1,11'-bi-2-naphthol), reacts with 1 equiv of [Rh(cod)(2)]BF4 to give the heteroligand complex [Rh(cod)(PMePh2)(L-1a)]BF4 (5) and homoligand complexes [Rh(cod)(PMePh2)(2)]BF4 (6) and [Rh(cod)(L-1a)(2)]BF4 (7) in the ratio 2:1:1. The same mixture of 5-7 is obtained upon mixing the isolated homoligand complexes 6 and 7 although the equilibrium is only established rapidly in the presence of an excess of PMePh2. The predominant species 5 is a monodentate ligand complex analogue of the chelate 3a. When the mixture of 5-7 is exposed to 5 atm H-2 for 1 h (the conditions used for catalyst preactivation in the asymmetric hydrogenation studies), the products are identified as the solvento species [Rh(PMePh2)(L-1a)(S)(2)]BF4 (5'), [Rh(S)(2)(PMePh2)(2)]BF4 (6') and [Rh(S)(2)(L-1a)(2)]BF4 (7') and are formed in the same 2:1:1 ratio. The reaction of M with 0.5 equiv of [Rh(cod)(2)]BF4 gives exclusively the heteroligand complex cis-[Rh(PMePh2)(2)(L-1a)(2)]BF4 (8), an analogue of 4a. The asymmetric hydrogenation of dehydroamino acid derivatives catalyzed by 3a,b is reported, and the enantioselectivities are compared with those obtained with (a) chelate catalysts derived from analogous diphosphonite ligands L-2a and L-2b, (b) catalysts based on methyl monophosphonites L-1a and L-1b, and (c) catalysts derived from mixture M. For the cinnamate and acrylate substrates studied, the catalysts derived from the phosphino/phosphonite bidentates L-A,L-B generally give superior enantioselectivities to the analogous diphosphonites L-2a and L-2b; these results are rationalized in terms of delta/lambda-chelate conformations and allosteric effects of the substrates. The rate of hydrogenation of acrylate substrate A with heterochelate 3a is significantly faster than with the homochelate analogues [Rh(L-2a)(cod)]BF4 and [Rh(dppe)(cod)]BF4. A synergic effect on the rate is also observed with the monodentate analogues: the rate of hydrogenation with the mixture containing predominantly heteroligand complex 5 is faster than with the monophosphine complex 6 or monophosphonite complex 7. Thus the hydrogenation catalysis carried out with M and [Rh(cod)(2)]BF4 is controlled by the dominant and most efficient heteroligand complex 5. In this study, the heterodiphos chelate 3a is shown to be more efficient and gives the opposite sense of optical induction t the heteromonophos analogue

Nanostructure formation in poly(gamma-benzyl-L-glutamate)-poly(ethylene glycol)-poly(gamma-benzyl-L-glutamate) triblock copolymers in the solid state

The morphology in the solid state of a series of triblock copolymers comprising a poly(ethylene glycol) (PEG) midblock and symmetric poly(gamma-benzyl-L-glutamate) (PBLG) end blocks has been studied using X-ray scattering and microscopy techniques. Transmission electron microscopy (TEM) on samples selectively stained with uranyl acetate provided clear assignment of morphologies for as-cast and annealed samples. The thickness of both PEG and PBLG domains was in good agreement with calculations based on the conformations of the respective chains, allowing for the crystal or amorphous state of PEG and the a-helical or P-sheet structure of the PBLG. Atomic force microscopy provided complementary information on surface morphology for several samples that was in good agreement with the structure observed by TEM. A morphology diagram was constructed. Cylindrical structures were observed for ordered samples with low f(PBLG), whereas at higher f(PLBG) there was evidence for broken lamellar and "hockey puck" nanostructures. Regular lamellae were observed for intermediate compositions.

Grey-box radial basis function modelling: the art of incorporating prior knowledge

A basic principle in data modelling is to incorporate available a priori information regarding the underlying data generating mechanism into the modelling process. We adopt this principle and consider grey-box radial basis function (RBF) modelling capable of incorporating prior knowledge. Specifically, we show how to explicitly incorporate the two types of prior knowledge: the underlying data generating mechanism exhibits known symmetric property and the underlying process obeys a set of given boundary value constraints. The class of orthogonal least squares regression algorithms can readily be applied to construct parsimonious grey-box RBF models with enhanced generalisation capability.

A new solution of the rendering equation with stratified Monte Carlo approach

This paper is turned to the advanced Monte Carlo methods for realistic image creation. It offers a new stratified approach for solving the rendering equation. We consider the numerical solution of the rendering equation by separation of integration domain. The hemispherical integration domain is symmetrically separated into 16 parts. First 9 sub-domains are equal size of orthogonal spherical triangles. They are symmetric each to other and grouped with a common vertex around the normal vector to the surface. The hemispherical integration domain is completed with more 8 sub-domains of equal size spherical quadrangles, also symmetric each to other. All sub-domains have fixed vertices and computable parameters. The bijections of unit square into an orthogonal spherical triangle and into a spherical quadrangle are derived and used to generate sampling points. Then, the symmetric sampling scheme is applied to generate the sampling points distributed over the hemispherical integration domain. The necessary transformations are made and the stratified Monte Carlo estimator is presented. The rate of convergence is obtained and one can see that the algorithm is of super-convergent type.

Robustness and applicability of Markov chain Monte Carlo algorithms for eigenvalue problems

In this paper we analyse applicability and robustness of Markov chain Monte Carlo algorithms for eigenvalue problems. We restrict our consideration to real symmetric matrices. Almost Optimal Monte Carlo (MAO) algorithms for solving eigenvalue problems are formulated. Results for the structure of both - systematic and probability error are presented. It is shown that the values of both errors can be controlled independently by different algorithmic parameters. The results present how the systematic error depends on the matrix spectrum. The analysis of the probability error is presented. It shows that the close (in some sense) the matrix under consideration is to the stochastic matrix the smaller is this error. Sufficient conditions for constructing robust and interpolation Monte Carlo algorithms are obtained. For stochastic matrices an interpolation Monte Carlo algorithm is constructed. A number of numerical tests for large symmetric dense matrices are performed in order to study experimentally the dependence of the systematic error from the structure of matrix spectrum. We also study how the probability error depends on the balancing of the matrix. (c) 2007 Elsevier Inc. All rights reserved.

Ocean acoustic models for low frequency propagation in 2D and 3D environments

In this paper we are mainly concerned with the development of efficient computer models capable of accurately predicting the propagation of low-to-middle frequency sound in the sea, in axially symmetric (2D) and in fully 3D environments. The major physical features of the problem, i.e. a variable bottom topography, elastic properties of the subbottom structure, volume attenuation and other range inhomogeneities are efficiently treated. The computer models presented are based on normal mode solutions of the Helmholtz equation on the one hand, and on various types of numerical schemes for parabolic approximations of the Helmholtz equation on the other. A new coupled mode code is introduced to model sound propagation in range-dependent ocean environments with variable bottom topography, where the effects of an elastic bottom, of volume attenuation, surface and bottom roughness are taken into account. New computer models based on finite difference and finite element techniques for the numerical solution of parabolic approximations are also presented. They include an efficient modeling of the bottom influence via impedance boundary conditions, they cover wide angle propagation, elastic bottom effects, variable bottom topography and reverberation effects. All the models are validated on several benchmark problems and versus experimental data. Results thus obtained were compared with analogous results from standard codes in the literature.

On the spectra of certain integro-differential-delay problems with applications in neurodynamics

We investigate the spectrum of certain integro-differential-delay equations (IDDEs) which arise naturally within spatially distributed, nonlocal, pattern formation problems. Our approach is based on the reformulation of the relevant dispersion relations with the use of the Lambert function. As a particular application of this approach, we consider the case of the Amari delay neural field equation which describes the local activity of a population of neurons taking into consideration the finite propagation speed of the electric signal. We show that if the kernel appearing in this equation is symmetric around some point a= 0 or consists of a sum of such terms, then the relevant dispersion relation yields spectra with an infinite number of branches, as opposed to finite sets of eigenvalues considered in previous works. Also, in earlier works the focus has been on the most rightward part of the spectrum and the possibility of an instability driven pattern formation. Here, we numerically survey the structure of the entire spectra and argue that a detailed knowledge of this structure is important within neurodynamical applications. Indeed, the Amari IDDE acts as a filter with the ability to recognise and respond whenever it is excited in such a way so as to resonate with one of its rightward modes, thereby amplifying such inputs and dampening others. Finally, we discuss how these results can be generalised to the case of systems of IDDEs.

Droplets of structured fluid on a flat substrate

We study the equilibrium morphology of droplets of symmetric AB diblock copolymer on a flat substrate. Using self-consistent field theory (SCFT), we provide the first predictions for the equilibrium droplet shape and its internal structure. When the sustrate affinity for the A component, $\eta_A$, is small, the droplet adopts a nearly spherical shape much like that of simple fluids. Inside the spherical droplet, however, concentric circular lamellar layers stack on top of each other; hence the thickness of the droplet is effectively quantized by a half-integer or integer number of layers. At larger $\eta_A$ and smaller contact angle, the area of the upper-most layer becomes relatively large, resulting in a nearly flat, faceted top surface, followed by a semi-spherical slope. This geometry is remarkably reminiscent of the droplet shapes observed with smetic liquid crystals.

Acoustic trapped modes in a three-dimensional waveguide of slowly varying cross section

In this paper we develop an asymptotic scheme to approximate the trapped mode solutions to the time harmonic wave equation in a three-dimensional waveguide with a smooth but otherwise arbitrarily shaped cross section and a single, slowly varying `bulge', symmetric in the longitudinal direction. Extending the work in Biggs (2012), we first employ a WKBJ-type ansatz to identify the possible quasi-mode solutions which propagate only in the thicker region, and hence find a finite cut-on region of oscillatory behaviour and asymptotic decay elsewhere. The WKBJ expansions are used to identify a turning point between the cut-on and cut-on regions. We note that the expansions are nonuniform in an interior layer centred on this point, and we use the method of matched asymptotic expansions to connect the cut-on and cut-on regions within this layer. The behaviour of the expansions within the interior layer then motivates the construction of a uniformly valid asymptotic expansion. Finally, we use this expansion and the symmetry of the waveguide around the longitudinal centre, x = 0, to extract trapped mode wavenumbers, which are compared with those found using a numerical scheme and seen to be extremely accurate, even to relatively large values of the small parameter.

Equitability revisited: why the “equitable threat score” is not equitable

In the forecasting of binary events, verification measures that are “equitable” were defined by Gandin and Murphy to satisfy two requirements: 1) they award all random forecasting systems, including those that always issue the same forecast, the same expected score (typically zero), and 2) they are expressible as the linear weighted sum of the elements of the contingency table, where the weights are independent of the entries in the table, apart from the base rate. The authors demonstrate that the widely used “equitable threat score” (ETS), as well as numerous others, satisfies neither of these requirements and only satisfies the first requirement in the limit of an infinite sample size. Such measures are referred to as “asymptotically equitable.” In the case of ETS, the expected score of a random forecasting system is always positive and only falls below 0.01 when the number of samples is greater than around 30. Two other asymptotically equitable measures are the odds ratio skill score and the symmetric extreme dependency score, which are more strongly inequitable than ETS, particularly for rare events; for example, when the base rate is 2% and the sample size is 1000, random but unbiased forecasting systems yield an expected score of around −0.5, reducing in magnitude to −0.01 or smaller only for sample sizes exceeding 25 000. This presents a problem since these nonlinear measures have other desirable properties, in particular being reliable indicators of skill for rare events (provided that the sample size is large enough). A potential way to reconcile these properties with equitability is to recognize that Gandin and Murphy’s two requirements are independent, and the second can be safely discarded without losing the key advantages of equitability that are embodied in the first. This enables inequitable and asymptotically equitable measures to be scaled to make them equitable, while retaining their nonlinearity and other properties such as being reliable indicators of skill for rare events. It also opens up the possibility of designing new equitable verification measures.

Verification of cloud-fraction forecasts

Cloud radar and lidar can be used to evaluate the skill of numerical weather prediction models in forecasting the timing and placement of clouds, but care must be taken in choosing the appropriate metric of skill to use due to the non- Gaussian nature of cloud-fraction distributions. We compare the properties of a number of different verification measures and conclude that of existing measures the Log of Odds Ratio is the most suitable for cloud fraction. We also propose a new measure, the Symmetric Extreme Dependency Score, which has very attractive properties, being equitable (for large samples), difficult to hedge and independent of the frequency of occurrence of the quantity being verified. We then use data from five European ground-based sites and seven forecast models, processed using the ‘Cloudnet’ analysis system, to investigate the dependence of forecast skill on cloud fraction threshold (for binary skill scores), height, horizontal scale and (for the Met Office and German Weather Service models) forecast lead time. The models are found to be least skillful at predicting the timing and placement of boundary-layer clouds and most skilful at predicting mid-level clouds, although in the latter case they tend to underestimate mean cloud fraction when cloud is present. It is found that skill decreases approximately inverse-exponentially with forecast lead time, enabling a forecast ‘half-life’ to be estimated. When considering the skill of instantaneous model snapshots, we find typical values ranging between 2.5 and 4.5 days. Copyright c 2009 Royal Meteorological Society

Baroclinic stationary waves in aquaplanet models

An aquaplanet model is used to study the nature of the highly persistent low-frequency waves that have been observed in models forced by zonally symmetric boundary conditions. Using the Hayashi spectral analysis of the extratropical waves, the authors find that a quasi-stationary wave 5 belongs to a wave packet obeying a well-defined dispersion relation with eastward group velocity. The components of the dispersion relation with k ≥ 5 baroclinically convert eddy available potential energy into eddy kinetic energy, whereas those with k < 5 are baroclinically neutral. In agreement with Green’s model of baroclinic instability, wave 5 is weakly unstable, and the inverse energy cascade, which had been previously proposed as a main forcing for this type of wave, only acts as a positive feedback on its predominantly baroclinic energetics. The quasi-stationary wave is reinforced by a phase lock to an analogous pattern in the tropical convection, which provides further amplification to the wave. It is also found that the Pedlosky bounds on the phase speed of unstable waves provide guidance in explaining the latitudinal structure of the energy conversion, which is shown to be more enhanced where the zonal westerly surface wind is weaker. The wave’s energy is then trapped in the waveguide created by the upper tropospheric jet stream. In agreement with Green’s theory, as the equator-to-pole SST difference is reduced, the stationary marginally stable component shifts toward higher wavenumbers, while wave 5 becomes neutral and westward propagating. Some properties of the aquaplanet quasi-stationary waves are found to be in interesting agreement with a low frequency wave observed by Salby during December–February in the Southern Hemisphere so that this perspective on low frequency variability, apart from its value in terms of basic geophysical fluid dynamics, might be of specific interest for studying the earth’s atmosphere.

Vertex splitting and connectivity augmentation in hypergraphs

We consider problems of splitting and connectivity augmentation in hypergraphs. In a hypergraph G = (V +s, E), to split two edges su, sv, is to replace them with a single edge uv. We are interested in doing this in such a way as to preserve a defined level of connectivity in V . The splitting technique is often used as a way of adding new edges into a graph or hypergraph, so as to augment the connectivity to some prescribed level. We begin by providing a short history of work done in this area. Then several preliminary results are given in a general form so that they may be used to tackle several problems. We then analyse the hypergraphs G = (V + s, E) for which there is no split preserving the local-edge-connectivity present in V. We provide two structural theorems, one of which implies a slight extension to Mader’s classical splitting theorem. We also provide a characterisation of the hypergraphs for which there is no such “good” split and a splitting result concerned with a specialisation of the local-connectivity function. We then use our splitting results to provide an upper bound on the smallest number of size-two edges we must add to any given hypergraph to ensure that in the resulting hypergraph we have λ(x, y) ≥ r(x, y) for all x, y in V, where r is an integer valued, symmetric requirement function on V*V. This is the so called “local-edge-connectivity augmentation problem” for hypergraphs. We also provide an extension to a Theorem of Szigeti, about augmenting to satisfy a requirement r, but using hyperedges. Next, in a result born of collaborative work with Zoltán Király from Budapest, we show that the local-connectivity augmentation problem is NP-complete for hypergraphs. Lastly we concern ourselves with an augmentation problem that includes a locational constraint. The premise is that we are given a hypergraph H = (V,E) with a bipartition P = {P1, P2} of V and asked to augment it with size-two edges, so that the result is k-edge-connected, and has no new edge contained in some P(i). We consider the splitting technique and describe the obstacles that prevent us forming “good” splits. From this we deduce results about which hypergraphs have a complete Pk-split. This leads to a minimax result on the optimal number of edges required and a polynomial algorithm to provide an optimal augmentation.

On the existence of extreme waves and the Stokes conjecture with vorticity

This is a study of singular solutions of the problem of traveling gravity water waves on flows with vorticity. We show that, for a certain class of vorticity functions, a sequence of regular waves converges to an extreme wave with stagnation points at its crests. We also show that, for any vorticity function, the profile of an extreme wave must have either a corner of 120° or a horizontal tangent at any stagnation point about which it is supposed symmetric. Moreover, the profile necessarily has a corner of 120° if the vorticity is nonnegative near the free surface.

The impact of the state of the troposphere on the response to stratospheric heating in a simplified GCM

Previous studies have made use of simplified general circulation models (sGCMs) to investigate the atmospheric response to various forcings. In particular, several studies have investigated the tropospheric response to changes in stratospheric temperature. This is potentially relevant for many climate forcings. Here the impact of changing the tropospheric climatology on the modeled response to perturbations in stratospheric temperature is investigated by the introduction of topography into the model and altering the tropospheric jet structure. The results highlight the need for very long integrations so as to determine accurately the magnitude of response. It is found that introducing topography into the model and thus removing the zonally symmetric nature of the model’s boundary conditions reduces the magnitude of response to stratospheric heating. However, this reduction is of comparable size to the variability in the magnitude of response between different ensemble members of the same 5000-day experiment. Investigations into the impact of varying tropospheric jet structure reveal a trend with lower-latitude/narrower jets having a much larger magnitude response to stratospheric heating than higher-latitude/wider jets. The jet structures that respond more strongly to stratospheric heating also exhibit longer time scale variability in their control run simulations, consistent with the idea that a feedback between the eddies and the mean flow is both responsible for the persistence of the control run variability and important in producing the tropospheric response to stratospheric temperature perturbations.