Application of the direct Liapunov method to the problem of symmetric stability in the atmosphere

The problem of symmetric stability is examined within the context of the direct Liapunov method. The sufficient conditions for stability derived by Fjørtoft are shown to imply finite-amplitude, normed stability. This finite-amplitude stability theorem is then used to obtain rigorous upper bounds on the saturation amplitude of disturbances to symmetrically unstable flows.By employing a virial functional, the necessary conditions for instability implied by the stability theorem are shown to be in fact sufficient for instability. The results of Ooyama are improved upon insofar as a tight two-sided (upper and lower) estimate is obtained of the growth rate of (modal or nonmodal) symmetric instabilities.The case of moist adiabatic systems is also considered.

The importance of surface pressure changes in the response of the atmosphere to zonally-symmetric thermal and mechanical forcing

he classical problem of the response of a balanced, axisymmetric vortex to thermal and mechanical forcing is re-examined, paying special attention to the lower boundary condition. The correct condition is DΦ/Dt = 0, where Φ is the geopotential and D/Dt the material derivative, which explicitly accounts for a mass redistribution as part of the mean-flow response. This redistribution is neglected when using the boundary condition Dp/Dt = 0, which has conventionally been applied in this problem. It is shown that applying the incorrect boundary condition, and thereby ignoring the surface pressure change, leads to a zonal wind acceleration δū/δt that is too strong, especially near the surface. The effect is significant for planetary-scale forcing even when applied at tropopause level. A comparison is made between the mean-flow evolution in a baroclinic life-cycle, as simulated in a fully nonlinear, primitive-equation model, and that predicted by using the simulated eddy fluxes in the zonally-symmetric response problem. Use of the correct lower boundary condition is shown to lead to improved agreement.

A modified PNN algorithm with optimal PD modeling using the orthogonal least squares method

In this paper a modified algorithm is suggested for developing polynomial neural network (PNN) models. Optimal partial description (PD) modeling is introduced at each layer of the PNN expansion, a task accomplished using the orthogonal least squares (OLS) method. Based on the initial PD models determined by the polynomial order and the number of PD inputs, OLS selects the most significant regressor terms reducing the output error variance. The method produces PNN models exhibiting a high level of accuracy and superior generalization capabilities. Additionally, parsimonious models are obtained comprising a considerably smaller number of parameters compared to the ones generated by means of the conventional PNN algorithm. Three benchmark examples are elaborated, including modeling of the gas furnace process as well as the iris and wine classification problems. Extensive simulation results and comparison with other methods in the literature, demonstrate the effectiveness of the suggested modeling approach.

Non-orthogonal version of the arbitrary polygonal C-grid and a new diamond grid

Quasi-uniform grids of the sphere have become popular recently since they avoid parallel scaling bottle- necks associated with the poles of latitude–longitude grids. However quasi-uniform grids of the sphere are often non- orthogonal. A version of the C-grid for arbitrary non- orthogonal grids is presented which gives some of the mimetic properties of the orthogonal C-grid. Exact energy conservation is sacrificed for improved accuracy and the re- sulting scheme numerically conserves energy and potential enstrophy well. The non-orthogonal nature means that the scheme can be used on a cubed sphere. The advantage of the cubed sphere is that it does not admit the computa- tional modes of the hexagonal or triangular C-grids. On var- ious shallow-water test cases, the non-orthogonal scheme on a cubed sphere has accuracy less than or equal to the orthog- onal scheme on an orthogonal hexagonal icosahedron. A new diamond grid is presented consisting of quasi- uniform quadrilaterals which is more nearly orthogonal than the equal-angle cubed sphere but with otherwise similar properties. It performs better than the cubed sphere in ev- ery way and should be used instead in codes which allow a flexible grid structure.

Morphology induced spinodal decomposition at the surface of symmetric diblock copolymer films

Atomic force microscopy is used to study the ordering dynamics of symmetric diblock copolymer films. The films order to form a lamellar structure which results in a frustration when the film thickness is incommensurate with the lamellae. By probing the morphology of incommensurate films in the early ordering stages, we discover an intermediate phase of lamellae arranged perpendicular to the film surface. This morphology is accompanied by a continuous growth in amplitude of the film surface topography with a characteristic wavelength, indicative of a spinodal process. Using selfconsistent field theory, we show that the observation of perpendicular lamellae suggests an intermediate state with parallel lamellae at the substrate and perpendicular lamellae at the free surface. The calculations confirm that the intermediate state is unstable to thickness fluctuations, thereby driving the spinodal growth of surface structures.

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̅.

Approximation by harmonic polynomials in star-shaped domains and exponential convergence of Trefftz hp-dGFEM

We study the approximation of harmonic functions by means of harmonic polynomials in two-dimensional, bounded, star-shaped domains. Assuming that the functions possess analytic extensions to a delta-neighbourhood of the domain, we prove exponential convergence of the approximation error with respect to the degree of the approximating harmonic polynomial. All the constants appearing in the bounds are explicit and depend only on the shape-regularity of the domain and on delta. We apply the obtained estimates to show exponential convergence with rate O(exp(−b square root N)), N being the number of degrees of freedom and b>0, of a hp-dGFEM discretisation of the Laplace equation based on piecewise harmonic polynomials. This result is an improvement over the classical rate O(exp(−b cubic root N )), and is due to the use of harmonic polynomial spaces, as opposed to complete polynomial spaces.

Elastic net orthogonal forward regression

An efficient two-level model identification method aiming at maximising a model׳s generalisation capability is proposed for a large class of linear-in-the-parameters models from the observational data. A new elastic net orthogonal forward regression (ENOFR) algorithm is employed at the lower level to carry out simultaneous model selection and elastic net parameter estimation. The two regularisation parameters in the elastic net are optimised using a particle swarm optimisation (PSO) algorithm at the upper level by minimising the leave one out (LOO) mean square error (LOOMSE). There are two elements of original contributions. Firstly an elastic net cost function is defined and applied based on orthogonal decomposition, which facilitates the automatic model structure selection process with no need of using a predetermined error tolerance to terminate the forward selection process. Secondly it is shown that the LOOMSE based on the resultant ENOFR models can be analytically computed without actually splitting the data set, and the associate computation cost is small due to the ENOFR procedure. Consequently a fully automated procedure is achieved without resort to any other validation data set for iterative model evaluation. Illustrative examples are included to demonstrate the effectiveness of the new approaches.

State-resolved gas-surface reactivity of methane in the symmetric C-H stretch vibration on Ni(100)

The state-resolved reactivity of CH4 in its totally symmetric C-H stretch vibration (�1) has been measured on a Ni(100) surface. Methane molecules were accelerated to kinetic energies of 49 and 63:5 kJ=mol in a molecular beam and vibrationally excited to �1 by stimulated Raman pumping before surface impact at normal incidence. The reactivity of the symmetric-stretch excited CH4 is about an order of magnitude higher than that of methane excited to the antisymmetric stretch (�3) reported by Juurlink et al. [Phys. Rev. Lett. 83, 868 (1999)] and is similar to that we have previously observed for the excitation of the first overtone (2�3). The difference between the state-resolved reactivity for �1 and �3 is consistent with predictions of a vibrationally adiabatic model of the methane reaction dynamics and indicates that statistical models cannot correctly describe the chemisorption of CH4 on nickel.

Nonlinear identification using orthogonal forward regression with nested optimal regularization

An efficient data based-modeling algorithm for nonlinear system identification is introduced for radial basis function (RBF) neural networks with the aim of maximizing generalization capability based on the concept of leave-one-out (LOO) cross validation. Each of the RBF kernels has its own kernel width parameter and the basic idea is to optimize the multiple pairs of regularization parameters and kernel widths, each of which is associated with a kernel, one at a time within the orthogonal forward regression (OFR) procedure. Thus, each OFR step consists of one model term selection based on the LOO mean square error (LOOMSE), followed by the optimization of the associated kernel width and regularization parameter, also based on the LOOMSE. Since like our previous state-of-the-art local regularization assisted orthogonal least squares (LROLS) algorithm, the same LOOMSE is adopted for model selection, our proposed new OFR algorithm is also capable of producing a very sparse RBF model with excellent generalization performance. Unlike our previous LROLS algorithm which requires an additional iterative loop to optimize the regularization parameters as well as an additional procedure to optimize the kernel width, the proposed new OFR algorithm optimizes both the kernel widths and regularization parameters within the single OFR procedure, and consequently the required computational complexity is dramatically reduced. Nonlinear system identification examples are included to demonstrate the effectiveness of this new approach in comparison to the well-known approaches of support vector machine and least absolute shrinkage and selection operator as well as the LROLS algorithm.

A proof of van Aubel's theorem using orthogonal vectors

We show how two linearly independent vectors can be used to construct two orthogonal vectors of equal magnitude in a simple way. The proof that the constructed vectors are orthogonal and of equal magnitude is a good exercise for students studying properties of scalar and vector triple products. We then show how this result can be used to prove van Aubel's theorem that relates the two line segments joining the centres of squares on opposite sides of a plane quadrilateral.

Low rank representation on Riemannian manifold of symmetric positive definite matrices

Sparse coding aims to find a more compact representation based on a set of dictionary atoms. A well-known technique looking at 2D sparsity is the low rank representation (LRR). However, in many computer vision applications, data often originate from a manifold, which is equipped with some Riemannian geometry. In this case, the existing LRR becomes inappropriate for modeling and incorporating the intrinsic geometry of the manifold that is potentially important and critical to applications. In this paper, we generalize the LRR over the Euclidean space to the LRR model over a specific Rimannian manifold—the manifold of symmetric positive matrices (SPD). Experiments on several computer vision datasets showcase its noise robustness and superior performance on classification and segmentation compared with state-of-the-art approaches.

Equilibrium behavior of symmetric ABA triblock copolymer melts

Melts of ABA triblock copolymer molecules with identical end blocks are examined using self-consistent field theory (SCFT). Phase diagrams are calculated and compared with those of homologous AB diblock copolymers formed by snipping the triblocks in half. This creates additional end segments which decreases the degree of segregation. Consequently, triblock melts remain ordered to higher temperatures than their diblock counterparts. We also find that middle-block domains are easier to stretch than end-block domains. As a result, domain spacings are slightly larger, the complex phase regions are shifted towards smaller A-segment compositions, and the perforated-lamellar phase becomes more metastable in triblock melts as compared to diblock melts. Although triblock and diblock melts exhibit very similar phase behavior, their mechanical properties can differ substantially due to triblock copolymers that bridge between otherwise disconnected A domains. We evaluate the bridging fraction for lamellar, cylindrical, and spherical morphologies to be about 40%–45%, 60%–65%, and 75%–80%, respectively. These fractions only depend weakly on the degree of segregation and the copolymer composition.