The Lanczos algorithm with selective orthogonalization

A new stable and efficient implementation of the Lanczos algorithm is presented. The algorithm is a powerful method for finding a few eigenvalues and eigenvectors at one or both ends of the spectrum of a symmetric matrix A. The algorithm is particularly effective if A is large and sparse in that the only way in which A enters the calculation is through a subroutine which computes Av for any vector v. Thus the user is free to take advantage of any sparsity structure in A and A need not even be represented as a matrix et al

Multifunctional Magnetoelectric Materials for Device Applications

Mutiferroics are a novel class of next generation multifunctional materials, which display simultaneous magnetic spin, electric dipole, and ferroelastic ordering, and have drawn increasing interest due to their multi-functionality for a variety of device applications. Since single-phase materials exist rarely in nature with such cross-coupling properties, an intensive research activity is being pursued towards the discovery of new single-phase multiferroic materials and the design of new engineered materials with strong magneto-electric (ME) coupling. This review article summarizes the development of different kinds of multiferroic material: single-phase and composite ceramic, laminated composite, and nanostructured thin films. Thin-film nanostructures have higher magnitude direct ME coupling values and clear evidence of indirect ME coupling compared with bulk materials. Promising ME coupling coefficients have been reported in laminated composite materials in which signal to noise ratio is good for device fabrication. We describe the possible applications of these materials

Photoionization cross sections of O II, O III, O IV, and O V: benchmarking R-matrix theory and experiments

For crucial tests between theory and experiment, ab initio close coupling calculations are carried out for photoionization of O II, O III, O IV, O V. The relativistic fine structure and resonance effects are studied using the R-matrix and its relativistic variant the Breit Pauli R-matrix (BPRM) approximation. Detailed comparison is made with high resolution experimental measurements carried out in three different set-ups: Advanced Light Source at Berkeley, and synchrotron radiation experiments at University of Aarhus and University of Paris-Sud. The comparisons illustrate physical effects in photoionization such as (i) fine structure, (ii) resolution, and (iii) metastable components. Photoionization cross sections sigma{PI} of the ground and a few low lying excited states of these ions obtained in the experimental spectrum include combined features of these states. Theoretically calculated resonances need to be resolved with extremely fine energy mesh for precise comparison. In addition, prominent resonant features are observed in the measured spectra from transitions allowed with relativistic fine structure, but not in LS coupling. The sigma_{PI} are obtained for ground and metastable (i) 2s^22p^3(^4S^o, ^2D^o, ^2P^o) states of O II, (ii) 2s^22p^2(^3P,^1D,^1S) and 2s2p^3(^5S^o) states of O III, (iii) 2s^22p(^2P^o_J) and 2s2p^2(^4P_J) levels of O IV, and (iv) 2s^2(^1S) and 2s2p(^3P^o,^1P^o) states of O V. It is found that resonances in ground and metastable cross sections can be a diagnostic of experimental beam composition, with potential ap plications to astrophysical and laboratory plasma environments.Comment: 27 pages, 7 figs., submitted to Phys. Rev. A., text with high resolution figures at http://www.astronomy.ohio-state.edu/~pradhan/Oions.p

Dynamical Origin of Extrasolar Planet Eccentricity Distribution

We explore the possibility that the observed eccentricity distribution of extrasolar planets arose through planet-planet interactions, after the initial stage of planet formation was complete. Our results are based on ~3250 numerical integrations of ensembles of randomly constructed planetary systems, each lasting 100 Myr. We find that for a remarkably wide range of initial conditions the eccentricity distributions of dynamically active planetary systems relax towards a common final equilibrium distribution, well described by the fitting formula dn ~ e exp[-1/2 (e/0.3)^2] de. This distribution agrees well with the observed eccentricity distribution for e > 0.2, but predicts too few planets at lower eccentricities, even when we exclude planets subject to tidal circularization. These findings suggest that a period of large-scale dynamical instability has occurred in a significant fraction of newly formed planetary systems, lasting 1--2 orders of magnitude longer than the ~1 Myr interval in which gas-giant planets are assembled. This mechanism predicts no (or weak) correlations between semimajor axis, eccentricity, inclination, and mass in dynamically relaxed planetary systems. An additional observational consequence of dynamical relaxation is a significant population of planets (>10%) that are highly inclined (>25deg) with respect to the initial symmetry plane of the protoplanetary disk; this population may be detectable in transiting planets through the Rossiter-McLaughlin effect.Comment: Accepted to ApJ, conclusions updated to reflect the current observational constraint

A Migration Study of \u3ci\u3eStelidota Geminata\u3c/i\u3e (Coleoptera: Nitidulidae)

The strawberry sap beetle, Stelidota geminata (Say), is a major pest of strawberries in the northeastern United States. Further knowledge of the migratory habits of this insect pest can enhance the effectiveness of pest management strategies. This nitidulid was shown to migrate from its overwintering sites to one of its primary reproductive sites, strawberry fields, in late May. The beetle population peaked in the third week in July, 1993 in the strawberry field and then gradually declined. In 1994, the peak, as well as the total population, was much greater than in 1993. Furthermore, S. geminata was concentrated in the transition areas surrounding the strawberry fields prior to the ripening of the fruit

Classification of n-qubit states with minimum orbit dimension

The group of local unitary transformations acts on the space of n-qubit pure states, decomposing it into orbits. In a previous paper we proved that a product of singlet states (together with an unentangled qubit for a system with an odd number of qubits) achieves the smallest possible orbit dimension, equal to 3n/2 for n even and (3n + 1)/2 for n odd, where n is the number of qubits. In this paper we show that any state with minimum orbit dimension must be of this form, and furthermore, such states are classified up to local unitary equivalence by the sets of pairs of qubits entangled in singlets.Comment: 15 pages, latex, revision 2, conclusion added, some proofs shortene

Symmetric mixed states of nn qubits: local unitary stabilizers and entanglement classes

We classify, up to local unitary equivalence, local unitary stabilizer Lie algebras for symmetric mixed states into six classes. These include the stabilizer types of the Werner states, the GHZ state and its generalizations, and Dicke states. For all but the zero algebra, we classify entanglement types (local unitary equivalence classes) of symmetric mixed states that have those stabilizers. We make use of the identification of symmetric density matrices with polynomials in three variables with real coefficients and apply the representation theory of SO(3) on this space of polynomials.Comment: 10 pages, 1 table, title change and minor clarifications for published versio

Ultraviolet/X-ray variability and the extended X-ray emission of the radio-loud broad absorption line quasar PG 1004+130

We present the results of recent Chandra, XMM-Newton, and Hubble Space Telescope observations of the radio-loud (RL), broad absorption line (BAL) quasar PG 1004+130. We compare our new observations to archival X-ray and UV data, creating the most comprehensive, high signal-to-noise, multi-epoch, spectral monitoring campaign of a RL BAL quasar to date. We probe for variability of the X-ray absorption, the UV BAL, and the X-ray jet, on month-year timescales. The X-ray absorber has a low column density of $N_{H}=8\times10^{20}-4\times10^{21}$ cm$^{-2}$ when it is assumed to be fully covering the X-ray emitting region, and its properties do not vary significantly between the 4 observations. This suggests the observed absorption is not related to the typical "shielding gas" commonly invoked in BAL quasar models, but is likely due to material further from the central black hole. In contrast, the CIV BAL shows strong variability. The equivalent width (EW) in 2014 is EW=11.24$\pm$0.56 \AA, showing a fractional increase of $\Delta EW / \langle EW \rangle$=1.16$\pm$0.11 from the 2003 observation, 3183 days earlier in the rest-frame. This places PG 1004+130 among the most highly variable BAL quasars. By combining Chandra observations we create an exposure 2.5 times deeper than studied previously, with which to investigate the nature of the X-ray jet and extended diffuse X-ray emission. An X-ray knot, likely with a synchrotron origin, is detected in the radio jet ~8 arcsec (30 kpc) from the central X-ray source with a spatial extent of ~4 arcsec (15 kpc). No similar X-ray counterpart to the counterjet is detected. Asymmetric, non-thermal diffuse X-ray emission, likely due to inverse Compton scattering of Cosmic Microwave Background photons, is also detected.Comment: 15 pages, 7 figures, 3 tables. Accepted for publication in Ap

Reflectionless analytic difference operators II. Relations to soliton systems

This is the second part of a series of papers dealing with an extensive class of analytic difference operators admitting reflectionless eigenfunctions. In the first part, the pertinent difference operators and their reflectionless eigenfunctions are constructed from given ``spectral data'', in analogy with the IST for reflectionless Schr\"odinger and Jacobi operators. In the present paper, we introduce a suitable time dependence in the data, arriving at explicit solutions to a nonlocal evolution equation of Toda type, which may be viewed as an analog of the KdV and Toda lattice equations for the latter operators. As a corollary, we reobtain various known results concerning reflectionless Schr\"odinger and Jacobi operators. Exploiting a reparametrization in terms of relativistic Calogero--Moser systems, we also present a detailed study of $N$-soliton solutions to our nonlocal evolution equation