Partial spectral projected gradient method with active-set strategy for linearly constrained optimization

A method for linearly constrained optimization which modifies and generalizes recent box-constraint optimization algorithms is introduced. The new algorithm is based on a relaxed form of Spectral Projected Gradient iterations. Intercalated with these projected steps, internal iterations restricted to faces of the polytope are performed, which enhance the efficiency of the algorithm. Convergence proofs are given and numerical experiments are included and commented. Software supporting this paper is available through the Tango Project web page: http://www.ime.usp.br/similar to egbirgin/tango/.

Upper limit on the cosmic-ray photon flux above 10(19) eV using the surface detector of the Pierre Auger Observatory

A method is developed to search for air showers initiated by photons using data recorded by the surface detector of the Auger Observatory. The approach is based on observables sensitive to the longitudinal shower development, the signal risetime and the curvature of the shower front. Applying this method to the data, tipper limits on the flux of photons of 3.8 x 10(-3), 2.5 x 10(-3), and 2.2 x 10(-3) km(-2) sr(-1) yr(-1) above 10(19) eV, 2 x 10(19) eV, and 4 x 10(19) eV are derived, with corresponding limits on the fraction of photons being 2.0%, 5.1%, and 31% (all limits at 95% c.l.). These photon limits disfavor certain exotic models of sources of cosmic rays. The results also show that the approach adopted by the Auger Observatory to calibrate the shower energy is not strongly biased by a contamination from photons. (C) 2008 Elsevier B.V. All rights reserved.

Correlation of the highest-energy cosmic rays with the positions of nearby active galactic nuclei

Data collected by the Pierre Auger Observatory provide evidence for anisotropy in the arrival directions of the cosmic rays with the highest-energies, which are correlated with the positions of relatively nearby active galactic nuclei (AGN) [Pierre Auger Collaboration, Science 318 (2007) 938]. The correlation has maximum significance for cosmic rays with energy greater than similar to 6 x 10(19) eV and AGN at a distance less than similar to 75 Mpc. We have confirmed the anisotropy at a confidence level of more than 99% through a test with parameters specified a priori, using an independent data set. The observed correlation is compatible with the hypothesis that cosmic rays with the highest-energies originate from extra-galactic sources close enough so that their flux is not significantly attenuated by interaction with the cosmic background radiation (the Greisen-Zatsepin-Kuz`min effect). The angular scale of the correlation observed is a few degrees, which suggests a predominantly light composition unless the magnetic fields are very weak outside the thin disk of our galaxy. Our present data do not identify AGN as the sources of cosmic rays unambiguously, and other candidate sources which are distributed as nearby AGN are not ruled out. We discuss the prospect of unequivocal identification of individual sources of the highest-energy cosmic rays within a few years of continued operation of the Pierre Auger Observatory. (C) 2008 Elsevier B.V. All rights reserved.

The Stern-Gerlach Phenomenon According to Classical Electrodynamics

We present a description of the Stem-Gerlach type experiments using only the concepts of classical electrodynamics and the Newton`s equations of motion. The quantization of the projections of the spin (or the projections of the magnetic dipole) is not introduced in our calculations. The main characteristic of our approach is a quantitative analysis of the motion of the magnetic atoms at the entrance of the magnetic field region. This study reveals a mechanism which modifies continuously the orientation of the magnetic dipole of the atom in a very short time interval, at the entrance of the magnetic field region. The mechanism is based on the conservation of the total energy associated with a magnetic dipole which moves in a non uniform magnetic field generated by an electromagnet. A detailed quantitative comparison with the (1922) Stem-Gerlach experiment and the didactical (1967) experiment by J.R. Zacharias is presented. We conclude, contrary to the original Stern-Gerlach statement, that the classical explanations are not ruled out by the experimental data.

Study of the spectroscopic properties and first hyperpolarizabilities of disperse azo dyes derived from 2-amino-5-nitrothiazole

The solvatochromism and other spectroscopic and photophysical characteristics of four azo disperse dyes, derived from 2-amino-5-nitrothiazole, were evaluated and interpreted with the aid of experimental data and quantum mechanical calculations. For the non-substituted compound two conformers, E and Z, were proposed for the isolated molecules, being the second one considerably less stable. The optimization of these structures in combination with a SCRF methodology (IEFPCM, Simulating the molecules in a continuum dielectric with characteristics of methanol), suggests that the Z form is not stable in solution. This same behaviour is expected for the substituted compounds, which is corroborated by experimental data presented in previous investigations [A.E.H. Machado, L.M. Rodrigues, S. Gupta, A.M.F. Oliveira-Campos, A.M.S. Silva, J. Mol. Struct. 738 (2005) 239-245]. For the substituted compounds, two forms derived from E conformer (A and R) are possible. Quantum mechanical data suggest for the isolated molecules, that the low energy absorption hand of the E conformers involve at least two close electronic states. having the low-lying excited state a (1)(n,pi*) nature, and being the S-2 state attributed to a (1)(pi,pi*) transition. The data also suggest a small energy gap between the absorption peaks of A and B, related to the easy conversion between these forms. For the structures optimized in combination with the applied SCRF methodology, an states inversion is observed for the Substituted compounds, with a considerable diminish of the energy gap between A and B absorption peaks. The electronic spectra of these compounds are quite sensitive to changes in the solvent polarity. The positive solvatochromism is more evident in aprotic solvents, probably due to the polarization induced by the solute. These compounds do not fluoresce at 298 K, but present a small but perceptible fluorescence at 77 K, which seems to be favoured by the nature of the group in the 2 `-position of the phenyl ring. Moreover, such compounds present expressive values for first hyperpolarizability, which implies in good non-linear optics (NLO) responses and photoswitching capability. (C) 2008 Elsevier B.V. All rights reserved.

Atmospheric effects on extensive air showers observed with the surface detector of the Pierre Auger observatory

Atmospheric parameters, Such as pressure (P), temperature (T) and density (rho proportional to P/T), affect the development of extensive air showers initiated by energetic cosmic rays. We have Studied the impact of atmospheric variations on extensive air showers by means of the surface detector of the Pierre Auger Observatory. The rate of events shows a similar to 10% seasonal modulation and similar to 2% diurnal one. We find that the observed behaviour is explained by a model including the effects associated with the variations of P and rho. The former affects the longitudinal development of air showers while the latter influences the Moliere radius and hence the lateral distribution of the shower particles. The model is validated with full simulations of extensive air showers using atmospheric profiles measured at the site of the Pierre Auger Observatory. (C) 2009 Elsevier B.V. All rights reserved.

Upper limit on the cosmic-ray photon fraction at EeV energies from the Pierre Auger Observatory

From direct observations of the longitudinal development of ultra-high energy air showers performed with the Pierre Auger Observatory, upper limits of 3.8%, 2.4%, 3.5% and 11.7% (at 95% c.l.) are obtained on the fraction of cosmic-ray photons above 2, 3, 5 and 10 EeV (1 EeV equivalent to 10(18) eV), respectively. These are the first experimental limits on ultra-high energy photons at energies below 10 EeV. The results complement previous constraints on top-down models from array data and they reduce systematic uncertainties in the interpretation of shower data in terms of primary flux, nuclear composition and proton-air cross-section. (C) 2009 Elsevier B.V. All rights reserved.

Novel 6-methanesulfonamide-3,4-methylenedioxyphenyl-N-acylhydrazones: Orally effective anti-inflammatory drug candidates

We described herein the molecular design of novel in vivo anti-inflammatory 6-methanesulfonamide-3,4-methylenedioxyphenyl-N-acylhydrazone derivatives (1) planned by applying the molecular hybridization approach. This work also points out to the discovery of LASSBio-930 (1c) as a novel anti-inflammatory and anti-hyperalgesic prototype, which was able to reduce carrageenan-induced rat paw edema with an ED(50) of 97.8 mu mol/kg, acting mainly as a non-selective COX inhibitor. (C) 2009 Elsevier Ltd. All rights reserved.

The fluorescence detector of the Pierre Auger Observatory

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions. (C) 2010 Elsevier B.V All rights reserved.

Exact penalties for variational inequalities with applications to nonlinear complementarity problems

In this paper, we present a new reformulation of the KKT system associated to a variational inequality as a semismooth equation. The reformulation is derived from the concept of differentiable exact penalties for nonlinear programming. The best theoretical results are presented for nonlinear complementarity problems, where simple, verifiable, conditions ensure that the penalty is exact. We close the paper with some preliminary computational tests on the use of a semismooth Newton method to solve the equation derived from the new reformulation. We also compare its performance with the Newton method applied to classical reformulations based on the Fischer-Burmeister function and on the minimum. The new reformulation combines the best features of the classical ones, being as easy to solve as the reformulation that uses the Fischer-Burmeister function while requiring as few Newton steps as the one that is based on the minimum.

New and improved results for packing identical unitary radius circles within triangles, rectangles and strips

The focus of study in this paper is the class of packing problems. More specifically, it deals with the placement of a set of N circular items of unitary radius inside an object with the aim of minimizing its dimensions. Differently shaped containers are considered, namely circles, squares, rectangles, strips and triangles. By means of the resolution of non-linear equations systems through the Newton-Raphson method, the herein presented algorithm succeeds in improving the accuracy of previous results attained by continuous optimization approaches up to numerical machine precision. The computer implementation and the data sets are available at http://www.ime.usp.br/similar to egbirgin/packing/. (C) 2009 Elsevier Ltd, All rights reserved.

Second-order negative-curvature methods for box-constrained and general constrained optimization

A Nonlinear Programming algorithm that converges to second-order stationary points is introduced in this paper. The main tool is a second-order negative-curvature method for box-constrained minimization of a certain class of functions that do not possess continuous second derivatives. This method is used to define an Augmented Lagrangian algorithm of PHR (Powell-Hestenes-Rockafellar) type. Convergence proofs under weak constraint qualifications are given. Numerical examples showing that the new method converges to second-order stationary points in situations in which first-order methods fail are exhibited.

Improving ultimate convergence of an augmented Lagrangian method

Optimization methods that employ the classical Powell-Hestenes-Rockafellar augmented Lagrangian are useful tools for solving nonlinear programming problems. Their reputation decreased in the last 10 years due to the comparative success of interior-point Newtonian algorithms, which are asymptotically faster. In this research, a combination of both approaches is evaluated. The idea is to produce a competitive method, being more robust and efficient than its `pure` counterparts for critical problems. Moreover, an additional hybrid algorithm is defined, in which the interior-point method is replaced by the Newtonian resolution of a Karush-Kuhn-Tucker (KKT) system identified by the augmented Lagrangian algorithm. The software used in this work is freely available through the Tango Project web page:http://www.ime.usp.br/similar to egbirgin/tango/.

Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.

A novel nanostructured composite formed by interaction of copper octa(3-aminopropyl)octasilsesquioxane with azide ligands: Preparation, characterization and a voltammetric application

This study presents the preparation, characterization and application of copper octa(3-aminopropyl)octasilsesquioxane following its subsequent reaction with azide ions (ASCA). The precursor (AC) and the novel compound (ASCA) were characterized by Fourier transform infrared spectra (FTIR), nuclear magnetic resonance (NMR), electron paramagnetic resonance (EPR), scanning electronic microscopy (SEM), X-ray diffraction (XRD), Thermogravimetric analyses and voltammetric technique. The cyclic voltammogram of the modified graphite paste electrode with ASCA (GPE-ASCA), showed one redox couple with formal potential (E(1/2)(ox)) = 0.30 V and an irreversible process at 1.1 V (vs. Ag/AgCl; NaCl 1.0 M; v = 20 mV s(-1)). The material is very sensitive to nitrite concentrations. The modified graphite paste electrode (GPE-ASCA) gives a linear range from 1.0 x 10(-4) to 4.0 x 10(-3) mol L(-1) for the determination of nitrite, with a detection limit of 2.1 x 10(-4) mol L(-1) and the amperometric sensitivity of 8.04 mA/mol L(-1). (C) 2010 Elsevier Ltd. All rights reserved.