An efficient algorithm for calculation of the Luenberger canonical form

An algorithm is suggested to obtain the Luenberger canonical form for multivariable systems. The method computes the canonical form directly without having to compute the transformation matrix. In addition, there is a large reduction in the number of calculations. The reduced computations along with Gaussian techniques lend to greater inherent accuracy and the ability to refine the solution with additional computations

Phase-resolved HST/STIS spectroscopy of the exposed white dwarf in the high-field polar AR UMa

Phase-resolved HST/STIS ultraviolet spectroscopy of the high-field polar AR UMa confirms that the WD photospheric Ly alpha Zeeman features are formed in a magnetic field of ~200 MG. In addition to the Ly alpha pi and sigma+ components, we detect the forbidden hydrogen 1s0->2s0 transition, which becomes ``enabled'' in the presence of both strong magnetic and electric fields. Our attempt in fitting the overall optical+UV low state spectrum with single temperature magnetic WD models remains rather unsatisfactory, indicating either a shortcoming in the present models or a new physical process acting in AR UMa. As a result, our estimate of the WD temperature remains somewhat uncertain, Twd=20000+-5000K. We detect a broad emission bump centered at ~1445A and present throughout the entire binary orbit, and a second bump near ~1650A, which appears only near the inferior conjunction of the secondary star. These are suggestive of low harmonic cyclotron emission produced by low-level (M-dot~1e-13 Msun/yr) accretion onto both magnetic poles. However, there is no evidence in the power spectrum of light variations for accretion in gas blobs. The observed Ly alpha emission line shows a strong phase dependence with maximum flux and redshift near orbital phase phi~0.3, strongly indicating an origin on the trailing hemisphere of the secondary star. An additional Ly alpha absorption feature with similar phasing as the Ly alpha emission, but a \~700km/s blueshift could tentatively be ascribed to absorption of WD emission in a moderately fast wind. We derive a column density of neutral hydrogen of NH=(1.1+-1.0)1e18 cm**-2, the lowest of any known polar.Comment: 26 pages, 10 figures, AAS TeX 5.0, accepted for publication in the Astrophysical Journa

Electrostatics of ions inside the nanopores and trans-membrane channels

A model of a finite cylindrical ion channel through a phospholipid membrane of width $L$ separating two electrolyte reservoirs is studied. Analytical solution of the Poisson equation is obtained for an arbitrary distribution of ions inside the trans-membrane pore. The solution is asymptotically exact in the limit of large ionic strength of electrolyte on the two sides of membrane. However, even for physiological concentrations of electrolyte, the electrostatic barrier sizes found using the theory are in excellent agreement with the numerical solution of the Poisson equation. The analytical solution is used to calculate the electrostatic potential energy profiles for pores containing charged protein residues. Availability of a semi-exact interionic potential should greatly facilitate the study of ionic transport through nanopores and ion channels

Slope Instability of the Earthen Levee in Boston, UK: Numerical Simulation and Sensor Data Analysis

The paper presents a slope stability analysis for a heterogeneous earthen levee in Boston, UK, which is prone to occasional slope failures under tidal loads. Dynamic behavior of the levee under tidal fluctuations was simulated using a finite element model of variably saturated linear elastic perfectly plastic soil. Hydraulic conductivities of the soil strata have been calibrated according to piezometers readings, in order to obtain correct range of hydraulic loads in tidal mode. Finite element simulation was complemented with series of limit equilibrium analyses. Stability analyses have shown that slope failure occurs with the development of a circular slip surface located in the soft clay layer. Both models (FEM and LEM) confirm that the least stable hydraulic condition is the combination of the minimum river levels at low tide with the maximal saturation of soil layers. FEM results indicate that in winter time the levee is almost at its limit state, at the margin of safety (strength reduction factor values are 1.03 and 1.04 for the low-tide and high-tide phases, respectively); these results agree with real-life observations. The stability analyses have been implemented as real-time components integrated into the UrbanFlood early warning system for flood protection

Revisiting two-step Forbush decreases

Interplanetary coronal mass ejections (ICMEs) and their shocks can sweep out galactic cosmic rays (GCRs), thus creating Forbush decreases (FDs). The traditional model of FDs predicts that an ICME and its shock decrease the GCR intensity in a two-step profile. This model, however, has been the focus of little testing. Thus, our goal is to discover whether a passing ICME and its shock inevitably lead to a two-step FD, as predicted by the model. We use cosmic ray data from 14 neutron monitors and, when possible, high time resolution GCR data from the spacecraft International Gamma Ray Astrophysical Laboratory (INTEGRAL). We analyze 233 ICMEs that should have created two-step FDs. Of these, only 80 created FDs, and only 13 created two-step FDs. FDs are thus less common than predicted by the model. The majority of events indicates that profiles of FDs are more complicated, particularly within the ICME sheath, than predicted by the model. We conclude that the traditional model of FDs as having one or two steps should be discarded. We also conclude that generally ignored small-scale interplanetary magnetic field structure can contribute to the observed variety of FD profiles

Far Infrared Variability of Sagittarius A*: 25.5 Hours of Monitoring with HerschelHerschel

Variable emission from Sgr~A*, the luminous counterpart to the super-massive black hole at the center of our Galaxy, arises from the innermost portions of the accretion flow. Better characterization of the variability is important for constraining models of the low-luminosity accretion mode powering Sgr~A*, and could further our ability to use variable emission as a probe of the strong gravitational potential in the vicinity of the $4\times10^{6}\mathrm{M}_{\odot}$ black hole. We use the \textit{Herschel} Spectral and Photometric Imaging Receiver (SPIRE) to monitor Sgr~A* at wavelengths that are difficult or impossible to observe from the ground. We find highly significant variations at 0.25, 0.35, and 0.5 mm, with temporal structure that is highly correlated across these wavelengths. While the variations correspond to $<$1% changes in the total intensity in the \textit{Herschel} beam containing Sgr~A*, comparison to independent, simultaneous observations at 0.85 mm strongly supports the reality of the variations. The lowest point in the light curves, $\sim$0.5 Jy below the time-averaged flux density, places a lower bound on the emission of Sgr~A* at 0.25 mm, the first such constraint on the THz portion of the SED. The variability on few hour timescales in the SPIRE light curves is similar to that seen in historical 1.3 mm data, where the longest time series is available, but the distribution of variations in the sub-mm do not show a tail of large-amplitude variations seen at 1.3 mm. Simultaneous X-ray photometry from XMM-Newton shows no significant variation within our observing period, which may explain the lack of very large variations if X-ray and submillimeter flares are correlated.Comment: Accepted for publication in Ap

Affine maps of density matrices

For quantum systems described by finite matrices, linear and affine maps of matrices are shown to provide equivalent descriptions of evolution of density matrices for a subsystem caused by unitary Hamiltonian evolution in a larger system; an affine map can be replaced by a linear map, and a linear map can be replaced by an affine map. There may be significant advantage in using an affine map. The linear map is generally not completely positive, but the linear part of an equivalent affine map can be chosen to be completely positive and related in the simplest possible way to the unitary Hamiltonian evolution in the larger system.Comment: 4 pages, title changed, sentence added, reference update

Analysis of the Hydrogen-rich Magnetic White Dwarfs in the SDSS

We have calculated optical spectra of hydrogen-rich (DA) white dwarfs with magnetic field strengths between 1 MG and 1000 MG for temperatures between 7000 K and 50000 K. Through a least-squares minimization scheme with an evolutionary algorithm, we have analyzed the spectra of 114 magnetic DAs from the SDSS (95 previously published plus 14 newly discovered within SDSS, and five discovered by SEGUE). Since we were limited to a single spectrum for each object we used only centered magnetic dipoles or dipoles which were shifted along the magnetic dipole axis. We also statistically investigated the distribution of magnetic-field strengths and geometries of our sample.Comment: to appear in the proceedings of the 16th European Workshop on White Dwarfs, Barcelona, 200