Lanczos subspace filter diagonalization: Homogeneous recursive filtering and a low-storage method for the calculation of matrix elements

We develop a new iterative filter diagonalization (FD) scheme based on Lanczos subspaces and demonstrate its application to the calculation of bound-state and resonance eigenvalues. The new scheme combines the Lanczos three-term vector recursion for the generation of a tridiagonal representation of the Hamiltonian with a three-term scalar recursion to generate filtered states within the Lanczos representation. Eigenstates in the energy windows of interest can then be obtained by solving a small generalized eigenvalue problem in the subspace spanned by the filtered states. The scalar filtering recursion is based on the homogeneous eigenvalue equation of the tridiagonal representation of the Hamiltonian, and is simpler and more efficient than our previous quasi-minimum-residual filter diagonalization (QMRFD) scheme (H. G. Yu and S. C. Smith, Chem. Phys. Lett., 1998, 283, 69), which was based on solving for the action of the Green operator via an inhomogeneous equation. A low-storage method for the construction of Hamiltonian and overlap matrix elements in the filtered-basis representation is devised, in which contributions to the matrix elements are computed simultaneously as the recursion proceeds, allowing coefficients of the filtered states to be discarded once their contribution has been evaluated. Application to the HO2 system shows that the new scheme is highly efficient and can generate eigenvalues with the same numerical accuracy as the basic Lanczos algorithm.

Subspace wavepacket evolution with Newton polynomials

Time-dependent wavepacket evolution techniques demand the action of the propagator, exp(-iHt/(h)over-bar), on a suitable initial wavepacket. When a complex absorbing potential is added to the Hamiltonian for combating unwanted reflection effects, polynomial expansions of the propagator are selected on their ability to cope with non-Hermiticity. An efficient subspace implementation of the Newton polynomial expansion scheme that requires fewer dense matrix-vector multiplications than its grid-based counterpart has been devised. Performance improvements are illustrated with some benchmark one and two-dimensional examples. (C) 2001 Elsevier Science B.V. All rights reserved.

Complex natural resonances of conducting planar objects buried in a dielectric half-space

A scheme is presented to incorporate a mixed potential integral equation (MPIE) using Michalski's formulation C with the method of moments (MoM) for analyzing the scattering of a plane wave from conducting planar objects buried in a dielectric half-space. The robust complex image method with a two-level approximation is used for the calculation of the Green's functions for the half-space. To further speed up the computation, an interpolation technique for filling the matrix is employed. While the induced current distributions on the object's surface are obtained in the frequency domain, the corresponding time domain responses are calculated via the inverse fast Fourier transform (FFT), The complex natural resonances of targets are then extracted from the late time response using the generalized pencil-of-function (GPOF) method. We investigate the pole trajectories as we vary the distance between strips and the depth and orientation of single, buried strips, The variation from the pole position of a single strip in a homogeneous dielectric medium was only a few percent for most of these parameter variations.

Fluorescence spectrum of a two-level atom driven by a multiple modulated field

We investigate the fluorescence spectrum of a two-level atom driven by a multiple amplitude-modulated field. The driving held is modeled as a polychromatic field composed of a strong central (resonant) component and a large number of symmetrically detuned sideband fields displaced from the central component by integer multiples of a constant detuning. Spectra obtained here differ qualitatively from those observed for a single pair of modulating fields [B. Blind, P.R. Fontana, and P. Thomann, J. Phys. B 13, 2717 (1980)]. In the case of a small number of the modulating fields, a multipeaked spectrum is obtained with the spectral features located at fixed frequencies that are independent of the number of modulating fields and their Rabi frequencies. As the number of the modulating fields increases, the spectrum ultimately evolves to the well-known Mellow triplet with the sidebands shifted from the central component by an effective Rabi frequency whose magnitude depends on the initial relative phases of the components of the driving held. For equal relative phases, the effective Rabi frequency of the driving field can be reduced to zero resulting in the disappearance of fluorescence spectrum, i.e., the atom can stop interacting with the field. When the central component and the modulating fields are 180 degrees out of phase, the spectrum retains its triplet structure with the sidebands located at frequencies equal to the sum of the Rabi frequencies of the component of the driving field. Moreover, we shaw that the frequency of spontaneous emission can be controlled and switched from one frequency to another when the Rabi frequency or initial phase of the modulating fields are varied.

Determination of eutectic solidification mode in Sr-modified hypoeutectic Al-Si alloys by EBSD

The effect of eutectic modification by strontium on nucleation and growth of the eutectic in hypoeutectic Al-Si foundry alloys has been investigated by electron back-scattering diffraction (EBSD) mapping. Specimens were prepared from three hypoeutectic AlSi base alloys with 5, 7 and 10 mass%Si and with different strontium contents up to 740 ppm for modification of eutectic silicon. By comparing the orientation of the aluminium in the eutectic to that of the surrounding primary aluminium dendrites? the growth mode of the eutectic could be determined. The mapping results indicate that the eutectic grew from the primary phase in unmodified alloys. When the eutectic was modified by strontium, eutectic grains nucleated separately from the primary dendrites. However, in alloys with high strontium levels, the eutectic again grew from the primary phase. These observed effects of strontium additions on the eutectic solidification mode are independent of silicon content in the range between 5 and 10 mass%Si.

Eutectic growth mode in strontium, antimony and phosphorus modified hypoeutectic Al-Si foundry alloys

The effect of strontium (Sr), antimony (Sb) and phosphorus (P) on nucleation and growth mode of the eutectic in hypoeutectic Al-10 mass%Si alloys has been investigated by electron back-scattering diffraction (EBSD) mapping. Specimens were prepared from a hypoeutectic Al-10 mass%Si base alloy, adding different levels of strontium, antimony and phosphorus for modification of eutectic silicon. By comparing the orientation of the aluminium in the eutectic to that of the surrounding primary aluminium dendrites, the solidification mode of the eutectic could be determined. The results of these studies show that the eutectic nucleation mode, and subsequent growth mode, is strongly dependent on additive elements. The EBSD mapping results indicate that the eutectic grew from the primary phase in unmodified and phosphorus-containing alloys. When the eutectic was modified by strontium or antimony, eutectic grains nucleated and grew separately from the primary dendrites.

Integrability and exact spectrum of a pairing model for nucleons

A pairing model for nucleons, introduced by Richardson in 1966, which describes proton-neutron pairing as well as proton-proton and neutron-neutron pairing, is re-examined in the context of the quantum inverse scattering method. Specifically, this shows that the model is integrable by enabling the explicit construction of the conserved operators. We determine the eigenvalues of these operators in terms of the Bethe ansatz, which in turn leads to an expression for the energy eigenvalues of the Hamiltonian.

Integrable variant of the one-dimensional Hubbard model

A new integrable model which is a variant of the one-dimensional Hubbard model is proposed. The integrability of the model is verified by presenting the associated quantum R-matrix which satisfies the Yang-Baxter equation. We argue that the new model possesses the SO(4) algebra symmetry, which contains a representation of the eta-pairing SU(2) algebra and a spin SU(2) algebra. Additionally, the algebraic Bethe ansatz is studied by means of the quantum inverse scattering method. The spectrum of the Hamiltonian, eigenvectors, as well as the Bethe ansatz equations, are discussed. (C) 2002 American Institute of Physics.

Microstrip circular phased array design and development using microwave antenna CAD tools

The complex design and development of a planar multilayer phased array antenna in microstrip technology can be simplified using two commercially available design tools 1) Ansoft Ensemble and 2) HP-EEsof Touchstone. In the approach presented here, Touchstone is used to design RF switches and phase shifters whose scattering parameters are incorporated in Ensemble simulations using its black box tool. Using this approach, Ensemble is able to fully analyze the performance of radiating and beamforming layers of a phased array prior to its manufacturing. This strategy is demonstrated in a design example of a 12-element linearly-polarized circular phased array operating at L band. A comparison between theoretical and experimental results of the array is demonstrated.

Temperature dependence of polaronic transport through single molecules and quantum dots

Motivated by recent experiments on electric transport through single molecules and quantum dots, we investigate a model for transport that allows for significant coupling between the electrons and a boson mode isolated on the molecule or dot. We focus our attention on the temperature-dependent properties of the transport. In the Holstein picture for polaronic transport in molecular crystals the temperature dependence of the conductivity exhibits a crossover from coherent (band) to incoherent (hopping) transport. Here, the temperature dependence of the differential conductance on resonance does not show such a crossover, but is mostly determined by the lifetime of the resonant level on the molecule or dot.

Influence of nano-structure on electrolytic properties in CeO2 based system

Doped ceria (CeO2) compounds are fluorite type oxides that show oxygen ionic conductivity higher than yttria stabilized zirconia, in oxidizing atmosphere. In order to improve the conductivity, the effective index was suggested to maximize the oxygen ionic conductivity in doped CeO2 based oxides. In addition, the true microstructure of doped CeO2 was observed at atomic scale for conclusion of conduction mechanism. Doped CeO2 had small domains (10-50 nm) with ordered structure in a grain. It is found that the electrolytic properties strongly depended on the nano-structural feature at atomic scale in doped CeO2 electrolyte.

Calculation of resonances and product state distributions for the unimolecular dissociation of H2S

Resonance phenomena associated with the unimolecular dissociation of H2S --> SH + H have been investigated quantum mechanically by the Lanczos homogeneous filter diagonalization method using a newly developed potential energy surface (J. Chem. Phys. 2001, 114, 320). Resonance energies, widths (rates), and product state distributions have been obtained. Both dissociation rates and product state distributions of SH show, strong fluctuations, indicating that the dissociation of H2S is essentially irregular. Statistical analysis of neighboring level spacing and width distributions also confirms this behavior. The dissociation rates and product state distributions are compared to the predictions of quantum phase space theory.

Secondary structure of the third extracellular loop responsible for ligand selectivity of a mammalian gonadotropin-releasing hormone receptor

The extracellular loop 3 (ECL3) of the mammalian gonadotropin-releasing hormone receptor (GnRH-R) contains an acidic amino acid (Glu(301) in the mouse GnRH-R,) that confers agonist selectivity for Are in mammalian GnRH. It is proposed that a specific conformation of ECL3 is necessary to orientate the carboxyl side chain of the acidic residue for interaction with Arg(8) of GnRH, which is supported by decreased affinity for Arg(8) GnRH but not Gln(8) GnRH when an adjacent Pro is mutated to Ala. To probe the structural contribution of the loop domain to the proposed presentation of the carboxyl side chain, we synthesized a model peptide (CGPEMLNRVSEPGC) representing residues 293-302 of mouse ECL3, where Cys and Gly residues are added symmetrically at the N and C termini, respectively, allowing the introduction of a disulfide bridge to simulate the distances at which the ECL3 is tethered to the transmembrane domains 6 and 7 of the receptor. The ability of the ECL3 peptide to bind GnRH with low affinity was demonstrated by its inhibition of GnRH stimulation of inositol phosphate production in cells expressing the GnRH-R. The CD bands of the ECL3 peptides exhibited a superposition of predominantly unordered structure and partial contributions from beta-sheet structure. Likewise, the analysis of the amide I and amide III bands from micro-Raman and FT Raman experiments revealed mainly unordered conformations of the cyclic and of the linear peptide. NMR data demonstrated the presence of a beta-hairpin among an ensemble of largely disordered structures in the cyclic peptide. The location of the turn linking the two strands of the hairpin was assigned to the three central residues L-296, N-297, and R-298. A small population of structured species among an ensemble of predominantly random coil conformation suggests that the unliganded receptor represents a variety of structural conformers, some of which have the potential to make contacts with the ligand. We propose a mechanism of receptor activation whereby binding of the agonist to the inactive receptor state induces and stabilizes a particular structural state of the loop domain, leading to further conformational rearrangements across the transmembrane domain and signal propagating interaction with G proteins. Interaction of the Glu(301) of the receptor with Arg(8) of GnRH induces a folded configuration of the ligand. Our proposal thus suggests that conformational changes of both ligand and receptor result from this interaction.

Quantum dynamical characterization of unimolecular resonances

We give a selective review of quantum mechanical methods for calculating and characterizing resonances in small molecular systems, with an emphasis on recent progress in Chebyshev and Lanczos iterative methods. Two archetypal molecular systems are discussed: isolated resonances in HCO, which exhibit regular mode and state specificity, and overlapping resonances in strongly bound HO2, which exhibit irregular and chaotic behavior. Future directions in this field are also discussed.