On the Nesterov-Todd Direction in Semidefinite Programming

On the Nesterov-Todd Direction in Semidefinite Programmin

Tidal signals in ocean-bottom magnetic measurements of the Northwestern Pacific: observation versus prediction

Motional induction in the ocean by tides has long been observed by both land and satellite measurements of magnetic fields. While these signals are weak (∼10 nT) when compared to the main magnetic field, their persistent nature makes them important for consideration during geomagnetic field modelling. Previous studies have reported several discrepancies between observations and numerical predictions of the tidal magnetic signals and those studies were inconclusive of the source of the error. We address this issue by (1) analysing magnetometer data from ocean-bottom stations, where the low-noise and high-signal environment is most suitable for detecting the weak tidal magnetic signals, (2) by numerically predicting the magnetic field with a spatial resolution that is 16times higher than the previous studies and (3) by using four different models of upper-mantle conductivity. We use vector magnetic data from six ocean-bottom electromagnetic (OBEM) stations located in the Northwestern Pacific Ocean. The OBEM tidal amplitudes were derived using an iteratively re-weighted least-squares (IRLS) method and by limiting the analysis of lunar semidiurnal (M2), lunar elliptic semidinurnal (N2) and diurnal (O1) tidal modes to the night-time. Using a 3-D electromagnetic induction solver and the TPX07.2 tidal model, we predict the tidal magnetic signal. We use earth models with non-uniform oceans and four 1-D mantle sections underneath taken from Kuvshinov and Olsen, Shimizu etal. and Baba etal. to compare the effect of upper-mantle conductivity. We find that in general, the predictions and observations match within 10-70 per cent across all the stations for each of the tidal modes. The median normalized percent difference (NPD) between observed and predicted amplitudes for the tidal modes M2, N2 and O1 were 15 per cent, 47 per cent and 98 per cent, respectively, for all the stations and models. At the majority of stations, and for each of the tidal modes, the higher resolution (0.25°×0.25°) modelling gave amplitudes consistently closer to the observations than the lower resolution (1°×1°) modelling. The difference in lithospheric resistance east and west of the Izu-Bonin trench system seems to be affecting the model response and observations in the O1 tidal mode. This response is not seen in the M2 and N2 modes, thereby indicating that the O1 mode is more sensitive to lithospheric resistanc

Stable periodic waves in coupled Kuramoto-Sivashinsky - Korteweg-de Vries equations

Periodic waves are investigated in a system composed of a Kuramoto-Sivashinsky - Korteweg-de Vries (KS-KdV) equation, which is linearly coupled to an extra linear dissipative equation. The model describes, e.g., a two-layer liquid film flowing down an inclined plane. It has been recently shown that the system supports stable solitary pulses. We demonstrate that a perturbation analysis, based on the balance equation for the field momentum, predicts the existence of stable cnoidal waves (CnWs) in the same system. It is found that the mean value U of the wave field u in the main subsystem, but not the mean value of the extra field, affects the stability of the periodic waves. Three different areas can be distinguished inside the stability region in the parameter plane (L,U), where L is the wave's period. In these areas, stable are, respectively, CnWs with positive velocity, constant solutions, and CnWs with negative velocity. Multistability, i.e., the coexistence of several attractors, including the waves with several maxima per period, appears at large value of L. The analytical predictions are completely confirmed by direct simulations. Stable waves are also found numerically in the limit of vanishing dispersion, when the KS-KdV equation goes over into the KS one.Comment: a latex text file and 16 eps files with figures. Journal of the Physical Society of Japan, in pres

Studies of Phase Turbulence in the One Dimensional Complex Ginzburg-Landau Equation

The phase-turbulent (PT) regime for the one dimensional complex Ginzburg-Landau equation (CGLE) is carefully studied, in the limit of large systems and long integration times, using an efficient new integration scheme. Particular attention is paid to solutions with a non-zero phase gradient. For fixed control parameters, solutions with conserved average phase gradient $\nu$ exist only for $|\nu|$ less than some upper limit. The transition from phase to defect-turbulence happens when this limit becomes zero. A Lyapunov analysis shows that the system becomes less and less chaotic for increasing values of the phase gradient. For high values of the phase gradient a family of non-chaotic solutions of the CGLE is found. These solutions consist of spatially periodic or aperiodic waves travelling with constant velocity. They typically have incommensurate velocities for phase and amplitude propagation, showing thereby a novel type of quasiperiodic behavior. The main features of these travelling wave solutions can be explained through a modified Kuramoto-Sivashinsky equation that rules the phase dynamics of the CGLE in the PT phase. The latter explains also the behavior of the maximal Lyapunov exponents of chaotic solutions.Comment: 16 pages, LaTeX (Version 2.09), 10 Postscript-figures included, submitted to Phys. Rev.

Stabilized Kuramoto-Sivashinsky system

A model consisting of a mixed Kuramoto - Sivashinsky - KdV equation, linearly coupled to an extra linear dissipative equation, is proposed. The model applies to the description of surface waves on multilayered liquid films. The extra equation makes its possible to stabilize the zero solution in the model, opening way to the existence of stable solitary pulses (SPs). Treating the dissipation and instability-generating gain in the model as small perturbations, we demonstrate that balance between them selects two steady-state solitons from their continuous family existing in the absence of the dissipation and gain. The may be stable, provided that the zero solution is stable. The prediction is completely confirmed by direct simulations. If the integration domain is not very large, some pulses are stable even when the zero background is unstable. Stable bound states of two and three pulses are found too. The work was supported, in a part, by a joint grant from the Israeli Minsitry of Science and Technology and Japan Society for Promotion of Science.Comment: A text file in the latex format and 20 eps files with figures. Physical Review E, in pres

Stable two-dimensional solitary pulses in linearly coupled dissipative Kadomtsev-Petviashvili equations

A two-dimensional (2D) generalization of the stabilized Kuramoto - Sivashinsky (KS) system is presented. It is based on the Kadomtsev-Petviashvili (KP) equation including dissipation of the generic (Newell -- Whitehead -- Segel, NWS) type and gain. The system directly applies to the description of gravity-capillary waves on the surface of a liquid layer flowing down an inclined plane, with a surfactant diffusing along the layer's surface. Actually, the model is quite general, offering a simple way to stabilize nonlinear waves in media combining the weakly-2D dispersion of the KP type with gain and NWS dissipation. Parallel to this, another model is introduced, whose dissipative terms are isotropic, rather than of the NWS type. Both models include an additional linear equation of the advection-diffusion type, linearly coupled to the main KP-NWS equation. The extra equation provides for stability of the zero background in the system, opening a way to the existence of stable localized pulses. The consideration is focused on the case when the dispersive part of the system of the KP-I type, admitting the existence of 2D localized pulses. Treating the dissipation and gain as small perturbations and making use of the balance equation for the field momentum, we find that the equilibrium between the gain and losses may select two 2D solitons, from their continuous family existing in the conservative counterpart of the model (the latter family is found in an exact analytical form). The selected soliton with the larger amplitude is expected to be stable. Direct simulations completely corroborate the analytical predictions.Comment: a latex text file and 16 eps files with figures; Physical Review E, in pres

Molecular dynamics simulations of lead clusters

Molecular dynamics simulations of nanometer-sized lead clusters have been performed using the Lim, Ong and Ercolessi glue potential (Surf. Sci. {\bf 269/270}, 1109 (1992)). The binding energies of clusters forming crystalline (fcc), decahedron and icosahedron structures are compared, showing that fcc cuboctahedra are the most energetically favoured of these polyhedral model structures. However, simulations of the freezing of liquid droplets produced a characteristic form of ``shaved'' icosahedron, in which atoms are absent at the edges and apexes of the polyhedron. This arrangement is energetically favoured for 600-4000 atom clusters. Larger clusters favour crystalline structures. Indeed, simulated freezing of a 6525-atom liquid droplet produced an imperfect fcc Wulff particle, containing a number of parallel stacking faults. The effects of temperature on the preferred structure of crystalline clusters below the melting point have been considered. The implications of these results for the interpretation of experimental data is discussed.Comment: 11 pages, 18 figues, new section added and one figure added, other minor changes for publicatio

A Simple Model for Anisotropic Step Growth

We consider a simple model for the growth of isolated steps on a vicinal crystal surface. It incorporates diffusion and drift of adatoms on the terrace, and strong step and kink edge barriers. Using a combination of analytic methods and Monte Carlo simulations, we study the morphology of growing steps in detail. In particular, under typical Molecular Beam Epitaxy conditions the step morphology is linearly unstable in the model and develops fingers separated by deep cracks. The vertical roughness of the step grows linearly in time, while horizontally the fingers coarsen proportional to $t^{0.33}$. We develop scaling arguments to study the saturation of the ledge morphology for a finite width and length of the terrace.Comment: 20 pages, 12 figures; [email protected]

Inferring Multiple Graphical Structures

Gaussian Graphical Models provide a convenient framework for representing dependencies between variables. Recently, this tool has received a high interest for the discovery of biological networks. The literature focuses on the case where a single network is inferred from a set of measurements, but, as wetlab data is typically scarce, several assays, where the experimental conditions affect interactions, are usually merged to infer a single network. In this paper, we propose two approaches for estimating multiple related graphs, by rendering the closeness assumption into an empirical prior or group penalties. We provide quantitative results demonstrating the benefits of the proposed approaches. The methods presented in this paper are embeded in the R package 'simone' from version 1.0-0 and later

Computational Aspects of Protein Functionality

The purpose of this short article is to examine certain aspects of protein functionality with relation to some key organizing ideas. This is important from a computational viewpoint in order to take account of modelling both biological systems and knowledge of these systems. We look at some of the lexical dimensions of the function and how certain constructs can be related to underlying ideas. The pervasive computational metaphor is then discussed in relation to protein multifunctionality, and the specific case of von Willebrand factor as a ‘smart’ multifunctional protein is briefly considered. Some diagrammatic techniques are then introduced to better articulate protein function