Structured minimal-memory inexact quasi-Newton method and secant preconditioners for augmented Lagrangian optimization

Augmented Lagrangian methods for large-scale optimization usually require efficient algorithms for minimization with box constraints. On the other hand, active-set box-constraint methods employ unconstrained optimization algorithms for minimization inside the faces of the box. Several approaches may be employed for computing internal search directions in the large-scale case. In this paper a minimal-memory quasi-Newton approach with secant preconditioners is proposed, taking into account the structure of Augmented Lagrangians that come from the popular Powell-Hestenes-Rockafellar scheme. A combined algorithm, that uses the quasi-Newton formula or a truncated-Newton procedure, depending on the presence of active constraints in the penalty-Lagrangian function, is also suggested. Numerical experiments using the Cute collection are presented.

FROM TRISECTIONS IN MODULE CATEGORIES TO QUASI-DIRECTED COMPONENTS

In this paper, we define and study a special type of trisections in a module category, namely the compact trisections which characterize quasi-directed components. We apply this notion to the study of laura algebras and we use it to define a class of algebras with predictable Auslander-Reiten components.

Gravitational signature of Schwarzschild black holes in dynamical Chern-Simons gravity

Dynamical Chern-Simons gravity is an extension of general relativity in which the gravitational field is coupled to a scalar field through a parity-violating Chern-Simons term. In this framework, we study perturbations of spherically symmetric black hole spacetimes, assuming that the background scalar field vanishes. Our results suggest that these spacetimes are stable, and small perturbations die away as a ringdown. However, in contrast to standard general relativity, the gravitational waveforms are also driven by the scalar field. Thus, the gravitational oscillation modes of black holes carry imprints of the coupling to the scalar field. This is a smoking gun for Chern-Simons theory and could be tested with gravitational-wave detectors, such as LIGO or LISA. For negative values of the coupling constant, ghosts are known to arise, and we explicitly verify their appearance numerically. Our results are validated using both time evolution and frequency domain methods.

Electrical transport in single-crystalline Li(0.9)Mo(6)O(17): A two-band Luttinger liquid exhibiting Bose metal behavior

Temperature-dependent electrical resistance in quasi-one-dimensional Li(0.9)Mo(6)O(17) is described by two Luttinger liquid anomalous exponents alpha, each associated with a distinct one dimensional band. The band with alpha < 1 is argued to crossover to a higher dimension below the temperature T(M'), leading to superconductivity. Disorder and magnetic fields are shown to induce the Bose metal behavior in this bulk compound.

Observational evidences on the modulation of the South American Low Level Jet east of the Andes according the ENSO variability

The differences on the phase and wavelength of the quasi-stationary waves over the South America generated by El Nino (EN) and La Nina (LN) events seem to affect the daily evolution of the South American Low Level Jet east of the Andes (SALLJ). For the austral summer period of 1977 2004 the SALLJ episodes detected according to Bonner criterion 1 show normal to above-normal frequency in EN years, and in LN years the episodes show normal to below-normal frequency. During EN and LN years the SALLJ episodes were associated with positive rainfall anomalies over the La Plata Basin, but more intense during LN years. During EN years the increase in the SALLJ cases were associated to intensification of the Subtropical Jet (SJ) around 30 degrees S and positive Sea Level Pressure (SLP) anomalies over the western equatorial Atlantic and tropical South America, particularly over central Brazil. This favored the intensification of the northeasterly trade winds over the northern continent and it channeled by the Andes mountain to the La Plata Basin region where negative SLP are found. The SALLJ cases identified during the LN events were weaker and less frequent when compared to those for EN years. In this case the SJ was weaker than in EN years and the negative SLP anomalies over the tropical continent contributed to the inversion of the northeasterly trade winds. Also a southerly flow anomaly was generated by the geostrophic balance due to the anomalous blocking over southeast Pacific and the intense cyclonic transient over the southern tip of South America. As result the warm tropical air brought by the SALLJ encounters the cold extratropical air from the southerly winds over the La Plata basin. This configuration can increase the conditional instability over the La Plata basin and may explain the more intense positive rainfall anomalies in SALLJ cases during LN years than in EN years.

Continental passive microwave-based rainfall estimation algorithm: Application to the Amazon Basin

This paper presents a new statistical algorithm to estimate rainfall over the Amazon Basin region using the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI). The algorithm relies on empirical relationships derived for different raining-type systems between coincident measurements of surface rainfall rate and 85-GHz polarization-corrected brightness temperature as observed by the precipitation radar (PR) and TMI on board the TRMM satellite. The scheme includes rain/no-rain area delineation (screening) and system-type classification routines for rain retrieval. The algorithm is validated against independent measurements of the TRMM-PR and S-band dual-polarization Doppler radar (S-Pol) surface rainfall data for two different periods. Moreover, the performance of this rainfall estimation technique is evaluated against well-known methods, namely, the TRMM-2A12 [ the Goddard profiling algorithm (GPROF)], the Goddard scattering algorithm (GSCAT), and the National Environmental Satellite, Data, and Information Service (NESDIS) algorithms. The proposed algorithm shows a normalized bias of approximately 23% for both PR and S-Pol ground truth datasets and a mean error of 0.244 mm h(-1) ( PR) and -0.157 mm h(-1)(S-Pol). For rain volume estimates using PR as reference, a correlation coefficient of 0.939 and a normalized bias of 0.039 were found. With respect to rainfall distributions and rain area comparisons, the results showed that the formulation proposed is efficient and compatible with the physics and dynamics of the observed systems over the area of interest. The performance of the other algorithms showed that GSCAT presented low normalized bias for rain areas and rain volume [0.346 ( PR) and 0.361 (S-Pol)], and GPROF showed rainfall distribution similar to that of the PR and S-Pol but with a bimodal distribution. Last, the five algorithms were evaluated during the TRMM-Large-Scale Biosphere-Atmosphere Experiment in Amazonia (LBA) 1999 field campaign to verify the precipitation characteristics observed during the easterly and westerly Amazon wind flow regimes. The proposed algorithm presented a cumulative rainfall distribution similar to the observations during the easterly regime, but it underestimated for the westerly period for rainfall rates above 5 mm h(-1). NESDIS(1) overestimated for both wind regimes but presented the best westerly representation. NESDIS(2), GSCAT, and GPROF underestimated in both regimes, but GPROF was closer to the observations during the easterly flow.

Cyclic variability of the circumstellar disk of the Be star zeta Tauri II. Testing the 2D global disk oscillation model

Context. About 2/3 of the Be stars present the so-called V/R variations, a phenomenon characterized by the quasi-cyclic variation in the ratio between the violet and red emission peaks of the HI emission lines. These variations are generally explained by global oscillations in the circumstellar disk forming a one-armed spiral density pattern that precesses around the star with a period of a few years. Aims. This paper presents self-consistent models of polarimetric, photometric, spectrophotometric, and interferometric observations of the classical Be star zeta Tauri. The primary goal is to conduct a critical quantitative test of the global oscillation scenario. Methods. Detailed three-dimensional, NLTE radiative transfer calculations were carried out using the radiative transfer code HDUST. The most up-to-date research on Be stars was used as input for the code in order to include a physically realistic description for the central star and the circumstellar disk. The model adopts a rotationally deformed, gravity darkened central star, surrounded by a disk whose unperturbed state is given by a steady-state viscous decretion disk model. It is further assumed that this disk is in vertical hydrostatic equilibrium. Results. By adopting a viscous decretion disk model for zeta Tauri and a rigorous solution of the radiative transfer, a very good fit of the time-average properties of the disk was obtained. This provides strong theoretical evidence that the viscous decretion disk model is the mechanism responsible for disk formation. The global oscillation model successfully fitted spatially resolved VLTI/AMBER observations and the temporal V/R variations in the H alpha and Br gamma lines. This result convincingly demonstrates that the oscillation pattern in the disk is a one-armed spiral. Possible model shortcomings, as well as suggestions for future improvements, are also discussed.

SOFT TISSUE INTEGRATION IN THE NECK AREA OF TITANIUM IMPLANTS - AN ANIMAL TRIAL

Dental implant materials are required to enable good apposition of bone and soft tissues. They must show sufficient resistance to chemical, physical and biological stress in the oral cavity to achieve good long-term outcomes. A critical issue is the apposition of the soft tissues, as they have provided a quasi-physiological closure of oral cavity. The present experiment was performed to study the peri-implant tissue response to non-submerged (1-stage) implant installation procedures. Two different implants types (NobelBiocare, NobelReplace (R) Tapered Groovy 4.3 x 10 mm and Replace (R) Select Tapered TiU RP 4.3 x 10 mm) were inserted into the right and left sides of 8 domestic pigs (Sus scrofa domestica) mandibles, between canines and premolars and immediately provided with a ceramic crown. Primary implant stability was determined using ressonance frequency analysis. Soft tissue parameters were assessed: sulcus depth (SDI) and junctional epithelium (JE). Following 70 days of healing, jaw sections were processed for histology and histomorphometric examination. Undecalcified histological sections demonstrated osseointegration with direct bone contact. The soft tissue parameters revealed no significant differences between the two implant types. The peri-implant soft tissues appear to behave similarly in both implant types.

Kinetic effect on low frequency Alfven continuum in tokamaks

Using a quasi-toroidal set of coordinates in plasmas with coaxial circular magnetic surfaces, the Vlasov equation is solved, and dielectric tensor is found for large aspect ratio tokamaks in a low frequency band. Taking into account the q-profile and drift effects, Alfven wave continuum deformation by geodesic effects is analyzed. It is shown that the Alfven continuum has a minimum defined by the ion thermal velocity at the rational magnetic surfaces q(s)=-M/N, where M and N are the poloidal and toroidal mode numbers, respectively, and the parallel wave number is zero. Low frequency global Alfven waves are found below the continuum minimum. In hot ion plasmas, the geodesic term changes sign, provoking some deformation of Alfven velocity by a factor (1+q(2))(-1/2), and the continuum minimum disappears. (C) 2008 American Institute of Physics.

Localized modes of binary mixtures of Bose-Einstein condensates in nonlinear optical lattices

The properties of the localized states of a two-component Bose-Einstein condensate confined in a nonlinear periodic potential (nonlinear optical lattice) are investigated. We discuss the existence of different types of solitons and study their stability by means of analytical and numerical approaches. The symmetry properties of the localized states with respect to nonlinear optical lattices are also investigated. We show that nonlinear optical lattices allow the existence of bright soliton modes with equal symmetry in both components and bright localized modes of mixed symmetry type, as well as dark-bright bound states and bright modes on periodic backgrounds. In spite of the quasi-one-dimensional nature of the problem, the fundamental symmetric localized modes undergo a delocalizing transition when the strength of the nonlinear optical lattice is varied. This transition is associated with the existence of an unstable solution, which exhibits a shrinking (decaying) behavior for slightly overcritical (undercritical) variations in the number of atoms.

Phase transitions and regions of stability in Reissner-Nordstrom holographic superconductors

The phase transition of Reissner-Nordstrom AdS(4) interacting with a massive charged scalar field has been further revisited. We found exactly one stable and one unstable quasinormal mode region for the scalar field. The two of them are separated by the first marginally stable solution.

Passage of radiation through wormholes of arbitrary shape

We study quasinormal modes and scattering properties via calculation of the S matrix for scalar and electromagnetic fields propagating in the background of spherically symmetric and axially symmetric traversable Lorentzian wormholes of a generic shape. Such wormholes are described by the general Morris-Thorne ansatz. The properties of quasinormal ringing and scattering are shown to be determined by the behavior of the wormhole's shape function b(r) and shift factor Phi(r) near the throat. In particular, wormholes with the shape function b(r), such that b(dr) approximate to 1, have very long-lived quasinormal modes in the spectrum. We have proved that the axially symmetric traversable Lorentzian wormholes, unlike black holes and other compact rotating objects, do not allow for superradiance. As a by-product we have shown that the 6th order WKB formula used for scattering problems of black or wormholes gives quite high accuracy and thus can be used for quite accurate calculations of the Hawking radiation processes around various black holes.

Gravitational stability of simply rotating Myers-Perry black holes: Tensorial perturbations

We study the stability of D >= 7 asymptotically flat black holes rotating in a single two-plane against tensor-type gravitational perturbations. The extensive search of quasinormal modes for these black holes did not indicate any presence of growing modes, implying the stability of simply rotating Myers-Perry black holes against tensor-type perturbations.

Stability of higher dimensional Reissner-Nordstrom-anti-de Sitter black holes

We investigate stability of the D-dimensional Reissner-Nordstrom-anti-de Sitter metrics as solutions of the Einstein-Maxwell equations. We have shown that asymptotically anti-de Sitter (AdS) black holes are dynamically stable for all values of charge and anti-de Sitter radius in D=5,6...11 dimensional space-times. This does not contradict dynamical instability of RNAdS black holes found by Gubser in N=8 gauged supergravity, because the latter instability comes from the tachyon mode of the scalar field, coupled to the system. Asymptotically AdS black holes are known to be thermodynamically unstable for some region of parameters, yet, as we have shown here, they are stable against gravitational perturbations.

(In)stability of D-dimensional black holes in Gauss-Bonnet theory

We make an extensive study of evolution of gravitational perturbations of D-dimensional black holes in Gauss-Bonnet theory. There is an instability at higher multipoles l and large Gauss-Bonnet coupling alpha for D = 5, 6, which is stabilized at higher D. Although a small negative gap of the effective potential for the scalar type of gravitational perturbations exists for higher D and whatever alpha, it does not lead to any instability.