An algorithm for arranging response surface designs in small blocks

We consider the problem of blocking response surface designs when the block sizes are prespecified to control variation efficiently and the treatment set is chosen independently of the block structure. We show how the loss of information due to blocking is related to scores defined by Mead and present an interchange algorithm based on scores to improve a given blocked design. Examples illustrating the performance of the algorithm are given and some comparisons with other designs are made. (C) 2000 Elsevier B.V. B.V. All rights reserved.

Some practical advice on polynomial regression analysis from blocked response surface designs

It is often necessary to run response surface designs in blocks. In this paper the analysis of data from such experiments, using polynomial regression models, is discussed. The definition and estimation of pure error in blocked designs are considered. It is recommended that pure error is estimated by assuming additive block and treatment effects, as this is more consistent with designs without blocking. The recovery of inter-block information using REML analysis is discussed, although it is shown that it has very little impact if thc design is nearly orthogonally blocked. Finally prediction from blocked designs is considered and it is shown that prediction of many quantities of interest is much simpler than prediction of the response itself.

A study of a numerical solution of the steady two dimensions Navier-Stokes equations in a constricted channel problem by a compact fourth order method

We present a numerical solution for the steady 2D Navier-Stokes equations using a fourth order compact-type method. The geometry of the problem is a constricted symmetric channel, where the boundary can be varied, via a parameter, from a smooth constriction to one possessing a very sharp but smooth corner allowing us to analyse the behaviour of the errors when the solution is smooth or near singular. The set of non-linear equations is solved by the Newton method. Results have been obtained for Reynolds number up to 500. Estimates of the errors incurred have shown that the results are accurate and better than those of the corresponding second order method. (C) 2002 Elsevier B.V. All rights reserved.

Statistics and kinetics of single-molecule electron transfer dynamics in complex environments: A simulation model study

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Critical temperature for first-order phase transitions in confined systems

We consider the Euclidean D-dimensional -lambda vertical bar phi vertical bar(4)+eta vertical bar rho vertical bar(6) (lambda,eta > 0) model with d (d <= D) compactified dimensions. Introducing temperature by means of the Ginzburg-Landau prescription in the mass term of the Hamiltonian, this model can be interpreted as describing a first-order phase transition for a system in a region of the D-dimensional space, limited by d pairs of parallel planes, orthogonal to the coordinates axis x(1), x(2),..., x(d). The planes in each pair are separated by distances L-1, L-2, ... , L-d. We obtain an expression for the transition temperature as a function of the size of the system, T-c({L-i}), i = 1, 2, ..., d. For D = 3 we particularize this formula, taking L-1 = L-2 = ... = L-d = L for the physically interesting cases d = 1 (a film), d = 2 (an infinitely long wire having a square cross-section), and for d = 3 (a cube). For completeness, the corresponding formulas for second-order transitions are also presented. Comparison with experimental data for superconducting films and wires shows qualitative agreement with our theoretical expressions.

Gauge formulation for higher order gravity

This work is an application of the second order gauge theory for the Lorentz group, where a description of the gravitational interaction is obtained that includes derivatives of the curvature. We analyze the form of the second field strength, G=partial derivative F+fAF, in terms of geometrical variables. All possible independent Lagrangians constructed with quadratic contractions of F and quadratic contractions of G are analyzed. The equations of motion for a particular Lagrangian, which is analogous to Podolsky's term of his generalized electrodynamics, are calculated. The static isotropic solution in the linear approximation was found, exhibiting the regular Newtonian behavior at short distances as well as a meso-large distance modification.

Difference variance dispersion graphs for comparing response surface designs with applications in food technology

Variance dispersion graphs have become a popular tool in aiding the choice of a response surface design. Often differences in response from some particular point, such as the expected position of the optimum or standard operating conditions, are more important than the response itself. We describe two examples from food technology. In the first, an experiment was conducted to find the levels of three factors which optimized the yield of valuable products enzymatically synthesized from sugars and to discover how the yield changed as the levels of the factors were changed from the optimum. In the second example, an experiment was conducted on a mixing process for pastry dough to discover how three factors affected a number of properties of the pastry, with a view to using these factors to control the process. We introduce the difference variance dispersion graph (DVDG) to help in the choice of a design in these circumstances. The DVDG for blocked designs is developed and the examples are used to show how the DVDG can be used in practice. In both examples a design was chosen by using the DVDG, as well as other properties, and the experiments were conducted and produced results that were useful to the experimenters. In both cases the conclusions were drawn partly by comparing responses at different points on the response surface.

Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.

Fractional polynomial response surface models

Second-order polynomial models have been used extensively to approximate the relationship between a response variable and several continuous factors. However, sometimes polynomial models do not adequately describe the important features of the response surface. This article describes the use of fractional polynomial models. It is shown how the models can be fitted, an appropriate model selected, and inference conducted. Polynomial and fractional polynomial models are fitted to two published datasets, illustrating that sometimes the fractional polynomial can give as good a fit to the data and much more plausible behavior between the design points than the polynomial model. © 2005 American Statistical Association and the International Biometric Society.

Braneworlds scenarios in a gravity model with higher order spatial three-curvature terms

In this work we study a Hořava-like 5-dimensional model in the context of braneworld theory. The equations of motion of such model are obtained and, within the realm of warped geometry, we show that the model is consistent if and only if λ takes its relativistic value 1. Furthermore, we show that the elimination of problematic terms involving the warp factor second order derivatives are eliminated by imposing detailed balance condition in the bulk. Afterwards, Israel's junction conditions are computed, allowing the attainment of an effective Lagrangian in the visible brane. In particular, we show that the resultant effective Lagrangian in the brane corresponds to a (3 + 1)-dimensional Hořava-like model with an emergent positive cosmological constant but without detailed balance condition. Now, restoration of detailed balance condition, at this time imposed over the brane, plays an interesting role by fitting accordingly the sign of the arbitrary constant β, insuring a positive brane tension and a real energy for the graviton within its dispersion relation. Also, the brane consistency equations are obtained and, as a result, the model admits positive brane tensions in the compactification scheme if, and only if, β is negative and the detailed balance condition is imposed. © 2013 Springer-Verlag Berlin Heidelberg and Società Italiana di Fisica.

Coherence characteristics of random lasing in a dye doped hybrid powder

The photon statistics of the random laser emission of a Rhodamine B doped di-ureasil hybrid powder is investigated to evaluate its degree of coherence above threshold. Although the random laser emission is a weighted average of spatially uncorrelated radiation emitted at different positions in the sample, a spatial coherence control was achieved due to an improved detection configuration based on spatial filtering. By using this experimental approach, which also allows for fine mode discrimination and timeresolved analysis of uncoupled modes from mode competition, an area not larger than the expected coherence size of the random laser is probed. Once the spectral and temporal behavior of nonoverlapping modes is characterized, an assessment of the photon-number probability distribution and the resulting second-order correlation coefficient as a function of time delay and wavelength was performed. The outcome of our single photon counting measurements revealed a high degree of temporal coherence at the time of maximum pump intensity and at wavelengths around the Rhodamine B gain maximum.

Avaliação psicométrica da escala de atitudes em relação à estatística

RESUMO: O objetivo do estudo foi estimar a validade e confiabilidade da Escala de atitudes em relação à Estatística (EAE) quando aplicada a estudantes de Ciências Farmacêuticas. A amostra de 253 estudantes foi subdividida em duas partes. Sessenta por cento da amostra foi utilizada para explorar a estrutura fatorial e 40% para confirmá-la. Para verificar a reprodutibilidade da escala a mesma foi aplicada em duplicata a 40 estudantes. Aplicou-se o Teste de esfericidade de Bartlett e o índice Kaiser-Meyer-Olkin (KMO). A extração dos fatores foi realizada pela Análise de Componentes Principais. Realizou-se rotação ortogonal Varimax. Calculou-se o Coeficiente alfa de Cronbach (α) e o Coeficiente de Correlação Intraclasse (ρ). Realizou-se análise fatorial confirmatória. Elaborou-se um modelo hierárquico de segunda ordem (MHSO). O teste de esfericidade de Bartlett e o índice KMO foram excelentes (χ 2 =1835,815, p<0,001; KMO=0,935). Verificou-se dois fatores com valores próprios acima de 1 (λ=9,748; λ=2,086) explicando 59,2% da variância total. A questão 2 foi removida. Observou-se excelente consistência interna e reprodutibilidade. O modelo fatorial apresentou índices de qualidade de ajustamento bons (λ=0,59-0,86, χ 2 /gl=1,691, CFI=0,919, GFI=0,804, RMSEA=0,079). A validade discriminante dos fatores foi adequada. A EAE apresentou estrutura bifatorial na amostra com níveis de validade e confiabilidade adequados.

Stripe-tetragonal phase transition in the two-dimensional Ising model with dipole interactions: Partition function zeros approach

We have performed multicanonical simulations to study the critical behavior of the two-dimensional Ising model with dipole interactions. This study concerns the thermodynamic phase transitions in the range of the interaction delta where the phase characterized by striped configurations of width h = 1 is observed. Controversial results obtained from local update algorithms have been reported for this region, including the claimed existence of a second-order phase transition line that becomes first order above a tricritical point located somewhere between delta = 0.85 and 1. Our analysis relies on the complex partition function zeros obtained with high statistics from multicanonical simulations. Finite size scaling relations for the leading partition function zeros yield critical exponents. that are clearly consistent with a single second-order phase transition line, thus excluding such a tricritical point in that region of the phase diagram. This conclusion is further supported by analysis of the specific heat and susceptibility of the orientational order parameter.

Higher-order perturbative relativistic corrections to energies and properties

Relativistic effects need to be considered in quantum-chemical calculations on systems including heavy elements or when aiming at high accuracy for molecules containing only lighter elements. In the latter case, consideration of relativistic effects via perturbation theory is an attractive option. Among the available techniques, Direct Perturbation Theory (DPT) in its lowest order (DPT2) has become a standard tool for the calculation of relativistic corrections to energies and properties.In this work, the DPT treatment is extended to the next order (DPT4). It is demonstrated that the DPT4 correction can be obtained as a second derivative of the energy with respect to the relativistic perturbation parameter. Accordingly, differentiation of a suitable Lagrangian, thereby taking into account all constraints on the wave function, provides analytic expressions for the fourth-order energy corrections. The latter have been implemented at the Hartree-Fock level and within second-order Møller-Plesset perturbaton theory using standard analytic second-derivative techniques into the CFOUR program package. For closed-shell systems, the DPT4 corrections consist of higher-order scalar-relativistic effects as well as spin-orbit corrections with the latter appearing here for the first time in the DPT series.Relativistic corrections are reported for energies as well as for first-order electrical properties and compared to results from rigorous four-component benchmark calculations in order to judge the accuracy and convergence of the DPT expansion for both the scalar-relativistic as well as the spin-orbit contributions. Additionally, the importance of relativistic effects to the bromine and iodine quadrupole-coupling tensors is investigated in a joint experimental and theoretical study concerning the rotational spectra of CH2BrF, CHBrF2, and CH2FI.