Parallel Shooting Sequential Quadratic Programming for Nonlinear MPC Problems

In this paper, we propose a parallel shooting algorithm for solving nonlinear model predictive control problems using sequential quadratic programming. This algorithm is built on a two-phase approach where we first test and assess sequential convergence over many initial trajectories in parallel. However, if none converge, the algorithm starts varying the Newton step size in parallel instead. Through this parallel shooting approach, it is expected that the number of iterations to converge to an optimal solution can be decreased. Furthermore, the algorithm can be further expanded and accelerated by implementing it on GPUs. We illustrate the effectiveness of the proposed Parallel Shooting Sequential Quadratic Programming (PS-SQP) method in some benchmark examples for nonlinear model predictive control. The developed PS-SQP parallel solver converges faster on average and especially when significant nonlinear behaviour is excited in the NMPC horizon.Comment: 7 pages, 6 figures, submitted and accepted for the 7th IEEE Conference on Control Technology and Applications (CCTA) 202

Why holes are not like electrons. II. The role of the electron-ion interaction

In recent work, we discussed the difference between electrons and holes in energy band in solids from a many-particle point of view, originating in the electron-electron interaction, and argued that it has fundamental consequences for superconductivity. Here we discuss the fact that there is also a fundamental difference between electrons and holes already at the single particle level, arising from the electron-ion interaction. The difference between electrons and holes due to this effect parallels the difference due to electron-electron interactions: {\it holes are more dressed than electrons}. We propose that superconductivity originates in 'undressing' of carriers from $both$ electron-electron and electron-ion interactions, and that both aspects of undressing have observable consequences.Comment: Continuation of Phys.Rev.B65, 184502 (2002) = cond-mat/0109385 (2001

Radio continuum and far-infrared emission from the galaxies in the Eridanus group

The Eridanus galaxies follow the well-known radio-FIR correlation. Majority (70%) of these galaxies have their star formation rates below that of the Milky Way. The galaxies having a significant excess of radio emission are identified as low luminosity AGNs based on their radio morphologies obtained from the GMRT observations. There are no powerful AGNs (L{20cm} > 10^{23} W Hz^{-1}) in the group. The two most far-infrared and radio luminous galaxies in the group have optical and HI morphologies suggestive of recent tidal interactions. The Eridanus group also has two far-infrared luminous but radio-deficient galaxies. It is believed that these galaxies are observed within a few Myr of the onset of an intense star formation episode after being quiescent for at least a 100 Myr. The upper end of the radio luminosity distribution of the Eridanus galaxies (L_{20cm} ~ 10^{22} W Hz^{-1}) is consistent with that of the field galaxies, other groups, and late-type galaxies in nearby clusters.Comment: 16 pages; Accepted for publication in Journal of Astroph. & Astron. March, 200

Low-temperature structural model of hcp solid C70_{70}

We report intermolecular potential-energy calculations for solid C_${70}$ and determine the optimum static orientations of the molecules at low temperature; we find them to be consistent with the monoclinic structural model proposed by us in an earlier report [Solid State Commun. {\bf 105), 247 (1998)]. This model indicates that the C_5 axis of the molecule is tilted by an angle $\approx$18^o from the monoclinic b axis in contrast with the molecular orientation proposed by Verheijen {\it et al.} [J. Chem. Phys. {\bf 166}, 287 (1992)] where the C_5 axis is parallel to the monoclinic b axis. In this calculation we have incorporated the effective bond charge Coulomb potential together with the Lennard-Jones potential between the molecule at the origin of the monoclinic unit cell and its six nearest neighbours, three above and three below. The minimum energy configuration for the molecular orientations turns out to be at $\theta$=18^o, $\phi$=8^o, and $\psi$=5^o, where $\theta$, $\phi$, and $\psi$ define the molecular orientations.Comment: ReVTeX (4 pages) + 2 PostScript figure

GMRT HI observations of the Eridanus group of galaxies

The GMRT HI 21cm-line observations of galaxies in the Eridanus group are presented. The Eridanus group, at a distance of ~23 Mpc, is a loose group of \~200 galaxies. The group extends more than 10 Mpc in projection. The velocity dispersion of the galaxies in the group is ~240 km/s. The galaxies are clustered into different sub-groups. The overall population mix of the group is 30% (E+S0) and 70% (Sp+Irr). The observations of 57 Eridanus galaxies were carried out with the GMRT for ~200 hour. HI emission was detected from 31 galaxies. The channel rms of ~1.0 mJy beam^{-1} was achieved for most of the image-cubes made with 4 hour of data. The corresponding HI column density sensitivity (3-sigma) is ~1x10^{20} cm^{-2} for a velocity-width of ~13.4 km/s. The 3-sigma detection limit of HI mass is ~1.2x10^{7} M_sun for a line-width of 50 km/s. Total HI images, HI velocity fields, global HI line profiles, HI mass surface densities, HI disk parameters and HI rotation curves are presented. The velocity fields are analysed separately for the approaching and the receding sides of the galaxies. This data will be used to study the HI and the radio continuum properties, the Tully-Fisher relations, the dark matter halos, and the kinematical and HI lopsidedness in galaxies.Comment: 75 pages including HI atlas; Accepted for publication in Journal of Astroph. & Astron. March, 200

Evidence for a nuclear compartment of transcription and splicing located at chromosome domain boundaries

The nuclear topography of splicing snRNPs, mRNA transcripts and chromosome domains in various mammalian cell types are described. The visualization of splicing snRNPs, defined by the Sm antigen, and coiled bodies, revealed distinctly different distribution patterns in these cell types. Heat shock experiments confirmed that the distribution patterns also depend on physiological parameters. Using a combination of fluorescencein situ hybridization and immunodetection protocols, individual chromosome domains were visualized simultaneously with the Sm antigen or the transcript of an integrated human papilloma virus genome. Three-dimensional analysis of fluorescence-stained target regions was performed by confocal laser scanning microscopy. RNA transcripts and components of the splicing machinery were found to be generally excluded from the interior of the territories occupied by the individual chromosomes. Based on these findings we present a model for the functional compartmentalization of the cell nucleus. According to this model the space between chromosome domains, including the surface areas of these domains, defines a three-dimensional network-like compartment, termed the interchromosome domain (ICD) compartment, in which transcription and splicing of mRNA occurs

IMAGES-III: The evolution of the Near-Infrared Tully-Fisher relation over the last 6 Gyr

Using the multi-integral field spectrograph GIRAFFE at VLT, we have derived the K-band Tully-Fisher relation (TFR) at z~0.6 for a representative sample of 65 galaxies with emission lines. We confirm that the scatter in the z~0.6 TFR is caused by galaxies with anomalous kinematics, and find a positive and strong correlation between the complexity of the kinematics and the scatter that they contribute to the TFR. Considering only relaxed-rotating disks, the scatter, and possibly also the slope of the TFR, do not appear to evolve with z. We detect an evolution of the K-band TFR zero point between z~0.6 and z=0, which, if interpreted as an evolution of the K-band luminosity of rotating disks, would imply that a brightening of 0.66+/-0.14 mag occurs between z~0.6 and z=0. Any disagreement with the results of Flores et al. (2006) are attributed to both an improvement of the local TFR and the more detailed accurate measurement of the rotation velocities in the distant sample. Most of the uncertainty can be explained by the relatively coarse spatial-resolution of the kinematical data. Because most rotating disks at z~0.6 are unlikely to experience further merging events, one may assume that their rotational velocity does not evolve dramatically. If true, our result implies that rotating disks observed at z~0.6 are rapidly transforming their gas into stars, to be able to double their stellar masses and be observed on the TFR at z=0. The rotating disks observed are indeed emission-line galaxies that are either starbursts or LIRGs, which implies that they are forming stars at a high rate. Thus, a significant fraction of the rotating disks are forming the bulk of their stars within 6 to 8 Gyr, in good agreement with former studies of the evolution of the M-Z relation.Comment: 17 pages, 11 figures, accepted for publication in A&A. v2 taking into account comments from language edito

CALIFA : a diameter-selected sample for an integral field spectroscopy galaxy survey

JMA acknowledges support from the European Research Council Starting Grant (SEDmorph; P.I. V. Wild).We describe and discuss the selection procedure and statistical properties of the galaxy sample used by the Calar Alto Legacy Integral Field Area (CALIFA) survey, a public legacy survey of 600 galaxies using integral field spectroscopy. The CALIFA "mother sample" was selected from the Sloan Digital Sky Survey (SDSS) DR7 photometric catalogue to include all galaxies with an r-band isophotal major axis between 45 '' and 79 : 2 '' and with a redshift 0 : 005 M-r > -23 : 1 and over a stellar mass range between 10(9.7) and 10(11.4) M-circle dot. In particular, within these ranges, the diameter selection does not lead to any significant bias against - or in favour of - intrinsically large or small galaxies. Only below luminosities of M-r = -19 (or stellar masses <10(9.7) M-circle dot) is there a prevalence of galaxies with larger isophotal sizes, especially of nearly edge-on late-type galaxies, but such galaxies form <10% of the full sample. We estimate volume-corrected distribution functions in luminosities and sizes and show that these are statistically fully compatible with estimates from the full SDSS when accounting for large-scale structure. For full characterization of the sample, we also present a number of value-added quantities determined for the galaxies in the CALIFA sample. These include consistent multi-band photometry based on growth curve analyses; stellar masses; distances and quantities derived from these; morphological classifications; and an overview of available multi-wavelength photometric measurements. We also explore different ways of characterizing the environments of CALIFA galaxies, finding that the sample covers environmental conditions from the field to genuine clusters. We finally consider the expected incidence of active galactic nuclei among CALIFA galaxies given the existing pre-CALIFA data, finding that the final observed CALIFA sample will contain approximately 30 Sey2 galaxies.Peer reviewe