Use of the q-Gaussian mutation in evolutionary algorithms

This paper proposes the use of the q-Gaussian mutation with self-adaptation of the shape of the mutation distribution in evolutionary algorithms. The shape of the q-Gaussian mutation distribution is controlled by a real parameter q. In the proposed method, the real parameter q of the q-Gaussian mutation is encoded in the chromosome of individuals and hence is allowed to evolve during the evolutionary process. In order to test the new mutation operator, evolution strategy and evolutionary programming algorithms with self-adapted q-Gaussian mutation generated from anisotropic and isotropic distributions are presented. The theoretical analysis of the q-Gaussian mutation is also provided. In the experimental study, the q-Gaussian mutation is compared to Gaussian and Cauchy mutations in the optimization of a set of test functions. Experimental results show the efficiency of the proposed method of self-adapting the mutation distribution in evolutionary algorithms.

Analyzing the n→π* Electronic Transition of Formaldehyde in Water. A Sequential Monte Carlo/Time-Dependent Density Functional Theory

The n→π* absorption transition of formaldehyde in water is analyzed using combined and sequential classical Monte Carlo (MC) simulations and quantum mechanics (QM) calculations. MC simulations generate the liquid solute-solvent structures for subsequent QM calculations. Using time-dependent density functional theory in a localized set of gaussian basis functions (TD-DFT/6-311++G(d,p)) calculations are made on statistically relevant configurations to obtain the average solvatochromic shift. All results presented here use the electrostatic embedding of the solvent. The statistically converged average result obtained of 2300 cm-1 is compared to previous theoretical results available. Analysis is made of the effective dipole moment of the hydrogen-bonded shell and how it could be held responsible for the polarization of the solvent molecules in the outer solvation shells.

Data collapse, scaling functions, and analytical solutions of generalized growth models

We consider a nontrivial one-species population dynamics model with finite and infinite carrying capacities. Time-dependent intrinsic and extrinsic growth rates are considered in these models. Through the model per capita growth rate we obtain a heuristic general procedure to generate scaling functions to collapse data into a simple linear behavior even if an extrinsic growth rate is included. With this data collapse, all the models studied become independent from the parameters and initial condition. Analytical solutions are found when time-dependent coefficients are considered. These solutions allow us to perceive nontrivial transitions between species extinction and survival and to calculate the transition's critical exponents. Considering an extrinsic growth rate as a cancer treatment, we show that the relevant quantity depends not only on the intensity of the treatment, but also on when the cancerous cell growth is maximum.

Charm and longitudinal structure functions within the Kharzeev-Levin-Nardi model

We use the Kharzeev-Levin-Nardi (KLN) model of the low x gluon distributions to fit recent HERA data on F(L) and F(2)(c)(F(2)(b)). Having checked that this model gives a good description of the data, we use it to predict F(L) and F(2)(c) to be measured in a future electron-ion collider. The results are similar to those obtained with the de Florian-Sassot and Eskola-Paukkunen-Salgado nuclear gluon distributions. The conclusion of this exercise is that the KLN model, simple as it is, may still be used as an auxiliary tool to make estimates for both heavy-ion and electron-ion collisions.

High p(T) direct photon and pi(0) triggered azimuthal jet correlations and measurement of k(T) for isolated direct photons in p plus p collisions at root s=200 GeV

Correlations of charged hadrons of 1< p(T) < 10 Gev/c with high pT direct photons and pi(0) mesons in the range 5< p(T) < 15 Gev/c are used to study jet fragmentation in the gamma + jet and dijet channels, respectively. The magnitude of the partonic transverse momentum, k(T), is obtained by comparing to a model incorporating a Gaussian kT smearing. The sensitivity of the associated charged hadron spectra to the underlying fragmentation function is tested and the data are compared to calculations using recent global fit results. The shape of the direct photon-associated hadron spectrum as well as its charge asymmetry are found to be consistent with a sample dominated by quark-gluon Compton scattering. No significant evidence of fragmentation photon correlated production is observed within experimental uncertainties.

Use of correlated potential harmonic basis functions for the description of the (4)He trimer and small clusters

A correlated many-body basis function is used to describe the (4)He trimer and small helium clusters ((4)HeN) with N = 4-9. A realistic helium dimer potential is adopted. The ground state results of the (4)He dimer and trimer are in close agreement with earlier findings. But no evidence is found for the existence of Efimov state in the trimer for the actual (4)He-(4)He interaction. However, decreasing the potential strength we calculate several excited states of the trimer which exhibit Efimov character. We also solve for excited state energies of these clusters which are in good agreement with Monte Carlo hyperspherical description. (C) 2011 American Institute of Physics. [doi:10.1063/1.3583365]

Deformed Gaussian-orthogonal-ensemble description of small-world networks

The study of spectral behavior of networks has gained enthusiasm over the last few years. In particular, random matrix theory (RMT) concepts have proven to be useful. In discussing transition from regular behavior to fully chaotic behavior it has been found that an extrapolation formula of the Brody type can be used. In the present paper we analyze the regular to chaotic behavior of small world (SW) networks using an extension of the Gaussian orthogonal ensemble. This RMT ensemble, coined the deformed Gaussian orthogonal ensemble (DGOE), supplies a natural foundation of the Brody formula. SW networks follow GOE statistics until a certain range of eigenvalue correlations depending upon the strength of random connections. We show that for these regimes of SW networks where spectral correlations do not follow GOE beyond a certain range, DGOE statistics models the correlations very well. The analysis performed in this paper proves the utility of the DGOE in network physics, as much as it has been useful in other physical systems.

Balance functions from Au+Au, d+Au, and p+p collisions at root s(NN)=200 GeV

Balance functions have been measured for charged-particle pairs, identified charged-pion pairs, and identified charged-kaon pairs in Au + Au, d + Au, and p + p collisions at root s(NN) = 200 GeV at the Relativistic Heavy Ion Collider using the STAR detector. These balance functions are presented in terms of relative pseudorapidity, Delta eta, relative rapidity, Delta y, relative azimuthal angle, Delta phi, and invariant relative momentum, q(inv). For charged-particle pairs, the width of the balance function in terms of Delta eta scales smoothly with the number of participating nucleons, while HIJING and UrQMD model calculations show no dependence on centrality or system size. For charged-particle and charged-pion pairs, the balance functions widths in terms of Delta eta and Delta y are narrower in central Au + Au collisions than in peripheral collisions. The width for central collisions is consistent with thermal blast-wave models where the balancing charges are highly correlated in coordinate space at breakup. This strong correlation might be explained by either delayed hadronization or limited diffusion during the reaction. Furthermore, the narrowing trend is consistent with the lower kinetic temperatures inherent to more central collisions. In contrast, the width of the balance function for charged-kaon pairs in terms of Delta y shows little centrality dependence, which may signal a different production mechanism for kaons. The widths of the balance functions for charged pions and kaons in terms of q(inv) narrow in central collisions compared to peripheral collisions, which may be driven by the change in the kinetic temperature.

Spin-polarized current and shot noise in the presence of spin flip in a quantum dot via nonequilibrium Green's functions

Using nonequilibrium Green's functions we calculate the spin-polarized current and shot noise in a ferromagnet-quantum-dot-ferromagnet system. Both parallel (P) and antiparallel (AP) magnetic configurations are considered. Coulomb interaction and coherent spin flip (similar to a transverse magnetic field) are taken into account within the dot. We find that the interplay between Coulomb interaction and spin accumulation in the dot can result in a bias-dependent current polarization p. In particular, p can be suppressed in the P alignment and enhanced in the AP case depending on the bias voltage. The coherent spin flip can also result in a switch of the current polarization from the emitter to the collector lead. Interestingly, for a particular set of parameters it is possible to have a polarized current in the collector and an unpolarized current in the emitter lead. We also found a suppression of the Fano factor to values well below 0.5.

INFERENCE FOR EIGENVALUES AND EIGENVECTORS OF GAUSSIAN SYMMETRIC MATRICES

This article presents maximum likelihood estimators (MLEs) and log-likelihood ratio (LLR) tests for the eigenvalues and eigenvectors of Gaussian random symmetric matrices of arbitrary dimension, where the observations are independent repeated samples from one or two populations. These inference problems are relevant in the analysis of diffusion tensor imaging data and polarized cosmic background radiation data, where the observations are, respectively, 3 x 3 and 2 x 2 symmetric positive definite matrices. The parameter sets involved in the inference problems for eigenvalues and eigenvectors are subsets of Euclidean space that are either affine subspaces, embedded submanifolds that are invariant under orthogonal transformations or polyhedral convex cones. We show that for a class of sets that includes the ones considered in this paper, the MLEs of the mean parameter do not depend on the covariance parameters if and only if the covariance structure is orthogonally invariant. Closed-form expressions for the MLEs and the associated LLRs are derived for this covariance structure.

LIPSCHITZ CLASSIFICATION OF FUNCTIONS ON A HOLDER TRIANGLE

The problem of semialgebraic Lipschitz classification of quasihomogeneous polynomials on a Holder triangle is studied. For this problem, the ""moduli"" are described completely in certain combinatorial terms.

Metastable BrO(2+) and NBr(2+) molecules in the gas phase

The doubly positively charged gas-phase molecules BrO(2+) and NBr(2+) have been produced by prolonged high-current energetic oxygen (17 keV (16)O(-)) ion surface bombardment (ion beam sputtering) of rubidium bromide (RbBr) and of ammonium bromide (NH(4)Br) powdered ionic salt samples, respectively, pressed into indium foil. These novel species were observed at half-integer m/z values in positive ion mass spectra for ion flight times of roughly similar to 12 mu s through a magnetic-sector secondary ion mass spectrometer. Here we present these experimental results and combine them with a detailed theoretical investigation using high level ab initio calculations of the ground states of BrO(2+) and NBr(2+), and a manifold of excited electronic states. NBr(2+) and BrO(2+), in their ground states, are long-lived metastable gas-phase molecules with well depths of 2.73 x 10(4) cm(-1) (3.38 eV) and 1.62 x 10(4) cm(-1) (2.01 eV); their fragmentation channels into two monocations lie 2.31 x 10(3) cm(-1) (0.29 eV) and 2.14 x 10(4) cm(-1) (2.65 eV) below the ground state minimum. The calculated lifetimes for NBr(2+) (v '' < 35) and BrO(2+) (v '' < 18) are large enough to be considered stable against tunneling. For NBr(2+), we predicted R(e) = 3.051 a(0) and omega(e) = 984 cm(-1); for BrO(2+), we obtained 3.033 a(0) and 916 cm(-1), respectively. The adiabatic double ionization energies of BrO and NBr to form metastable BrO(2+) and NBr(2+) are calculated to be 30.73 and 29.08 eV, respectively. The effect of spin-orbit interactions on the low-lying (Lambda + S) states is also discussed. (C) 2011 American Institute of Physics. [doi:10.1063/1.3562121]

The electronic states of SeF: A reinterpretation of the chemiluminescent emission of the reaction of selenium with fluorine

The low-lying doublet and quartet electronic states of the species SeF correlating with the first dissociation channel are investigated theoretically at a high-level of electronic correlation treatment, namely, the complete active space self-consistent field/multireference single and double excitations configuration interaction (CASSCF/MRSDCI) using a quintuple-zeta quality basis set including a relativistic effective core potential for the selenium atom. Potential energy curves for (Lambda+S) states and the corresponding spectroscopic properties are derived that allows for an unambiguous assignment of the only spectrum known experimentally as due to a spin-forbidden X (2)Pi-a (4)Sigma(-) transition, and not a A (2)Pi-X (2)Pi transition as assumed so far. For the bound excited doublets, yet unknown experimentally, this study is the first theoretical characterization of their spectroscopic properties. Also the spin-orbit coupling constant function for the X (2)Pi state is derived as well as the spin-orbit coupling matrix element between the X (2)Pi and a (4)Sigma(-) states. Dipole moment functions and vibrationally averaged dipole moments show SeF to be a very polar species. An overview of the lowest-lying spin-orbit (Omega) states completes this description. (C) 2010 American Institute of Physics. [doi: 10.1063/1.3426315]

Kernel machines for epilepsy diagnosis via EEG signal classification: A comparative study

Objective: We carry out a systematic assessment on a suite of kernel-based learning machines while coping with the task of epilepsy diagnosis through automatic electroencephalogram (EEG) signal classification. Methods and materials: The kernel machines investigated include the standard support vector machine (SVM), the least squares SVM, the Lagrangian SVM, the smooth SVM, the proximal SVM, and the relevance vector machine. An extensive series of experiments was conducted on publicly available data, whose clinical EEG recordings were obtained from five normal subjects and five epileptic patients. The performance levels delivered by the different kernel machines are contrasted in terms of the criteria of predictive accuracy, sensitivity to the kernel function/parameter value, and sensitivity to the type of features extracted from the signal. For this purpose, 26 values for the kernel parameter (radius) of two well-known kernel functions (namely. Gaussian and exponential radial basis functions) were considered as well as 21 types of features extracted from the EEG signal, including statistical values derived from the discrete wavelet transform, Lyapunov exponents, and combinations thereof. Results: We first quantitatively assess the impact of the choice of the wavelet basis on the quality of the features extracted. Four wavelet basis functions were considered in this study. Then, we provide the average accuracy (i.e., cross-validation error) values delivered by 252 kernel machine configurations; in particular, 40%/35% of the best-calibrated models of the standard and least squares SVMs reached 100% accuracy rate for the two kernel functions considered. Moreover, we show the sensitivity profiles exhibited by a large sample of the configurations whereby one can visually inspect their levels of sensitiveness to the type of feature and to the kernel function/parameter value. Conclusions: Overall, the results evidence that all kernel machines are competitive in terms of accuracy, with the standard and least squares SVMs prevailing more consistently. Moreover, the choice of the kernel function and parameter value as well as the choice of the feature extractor are critical decisions to be taken, albeit the choice of the wavelet family seems not to be so relevant. Also, the statistical values calculated over the Lyapunov exponents were good sources of signal representation, but not as informative as their wavelet counterparts. Finally, a typical sensitivity profile has emerged among all types of machines, involving some regions of stability separated by zones of sharp variation, with some kernel parameter values clearly associated with better accuracy rates (zones of optimality). (C) 2011 Elsevier B.V. All rights reserved.

Fuzzy pole placement based on piecewise Lyapunov functions

This paper presents a controller design method for fuzzy dynamic systems based on piecewise Lyapunov functions with constraints on the closed-loop pole location. The main idea is to use switched controllers to locate the poles of the system to obtain a satisfactory transient response. It is shown that the global fuzzy system satisfies the requirements for the design and that the control law can be obtained by solving a set of linear matrix inequalities, which can be efficiently solved with commercially available softwares. An example is given to illustrate the application of the proposed method. Copyright (C) 2009 John Wiley & Sons, Ltd.