A numerical experimental study of inverse preconditioning for the parallel iterative solution to 3D finite element flow equations

Integration of the subsurface flow equation by finite elements (FE) in space and finite differences (FD) in time requires the repeated solution to sparse symmetric positive definite systems of linear equations. Iterative techniques based on preconditioned conjugate gradients (PCG) are one of the most attractive tool to solve the problem on sequential computers. A present challenge is to make PCG attractive in a parallel computing environment as well. To this aim a key factor is the development of an efficient parallel preconditioner. FSAI (factorized sparse approximate inverse) and enlarged FSAI relying on the approximate inverse of the coefficient matrix appears to be a most promising parallel preconditioner. In the present paper PCG using FSAI, diagonal and pARMS (parallel algebraic recursive multilevel solvers) preconditioners is implemented on the IBM SP4/512 and CLX/768 supercomputers with up to 32 processors to solve underground flow problems of a large size. The results show that FSAI may allow for a parallel relative efficiency larger than 50% on the largest problems with p=32 processors. Moreover, FSAI turns out to be significantly less expensive and more robust than pARMS. Finally, it is shown that for p in the upper range may be much improved if PCG–FSAI is implemented on CLX.     

Quantitative Comparison of REAPDOR and TRAPDOR Experiments by Numerical Simulations and Determination of H–Al Distances in Zeolites

Quantitative H–Al distances in acid sites of two zeolites with MFI and IFR framework topology were obtained by numerical simulation of ¹H{²⁷Al} rotational echo adiabatic passage double resonance (REAPDOR) experiments. A ²⁷Al offset-dependent data set yields for each resolved ¹H NMR line a corresponding nuclear electric quadrupole coupling constant of the neighboring ²⁷Al site. This information is used for analyzing a second data set for on-resonance irradiation, where the dipolar evolution time (number of rotor cycles) was varied, to yield the ¹H–²⁷Al dipolar coupling constant. Numerical simulations indicate that the REAPDOR method does not depend significantly on the polar angles, defining the orientation of the electric field gradient tensor of ²⁷Al with respect to the Al–H dipolar vector. In contrast, the transfer of populations in double resonance sequence is sensitive to these angles, and it can be thus used to measure them.     

A scheme for teleporting Schrdinger-cat states via the dispersive atom-cavity-field interaction

A proposal is presented for teleporting Schrding-cat states. The process of the teleportation is achieved through the dispersive atom-cavity-field interaction. In this proposal, only measurement on the cavity field and on the singlet atomic states are used.     

On diquark clustering in quark—gluon plasma

The possibility that pairs of quarks will form diquark clusters in the regime above deconfinement transition for hadron matter at finite density is revisited. Here we present the results on the diquark-diquark (dq-dq) interaction in the framework of constituent quark model taking account of spin, isospin and color degrees of freedom in the spirit of generalized Pauli principle. By constructing the appropriate spin and color states of the dq-dq clusters we compute the expectation values of the interaction Hamiltonian involving pairwise quark—quark interaction. We find that the effective interaction between two diquark clusters is quite sensitive to different configurations characterized by color and spin states, obtained after the coupling of two diquark states. The value of the coupling parameter for a particular color—spin state, i.e., -3, 1 is compared to the one obtained earlier by Donoghue and Sateesh,Phys. Rev. D38, 360 (1988) based on the effective Φ⁴-theory. This new value of λ derived for different color-spin dq-dq states, may lead to several important implications in the studies of diquark star and diquark gas.     

AN ITERATIVE EQUATION ON THE UNIT CIRCLE

A functional equation of nonlinear iterates is discussed on the circle S1 for its continuous solutions and differentiable solutions. By lifting to R, the existence, uniqueness and stability of those solutions are obtained. Techniques of continuation are used to guarantee the preservation of continuity and differentiability in lifting.     

The LPCTrap experiment: measurement of the β–ν angular correlation in 6He using a transparent Paul trap

We report on the status of the LPCTrap experiment, devoted to measure the β–ν angular correlation in the pure Gamow–Teller decay of ⁶He. This measurement is motivated by the search for the presence of tensor type contributions to the weak interaction. The ⁶He ions are confined in a novel transparent Paul trap. The β particles and the recoil ions are detected in coincidence to deduce the angular correlation parameter. The commissioning run performed in 2005 has given the proof of principle of this experiment. Up to 10⁵ coincidences were recorded during a second run in 2006.     

Intermolecular Interactions in 8-Aza-D-Homogonane Crystals and their Manifestation in the IR Spectra

Crystalline modifications of 8-aza-D-homogon-1,3,5(10),13-tetraen-12,17a-dione have been investigated by IR spectroscopy and x-ray structural analysis. It is shown that this compound crystallizes from solutions of chloroform with hexane in the form of a solvate comprised of chelate hydrogen bonds between the hydrogen atom of the chloroform molecule and the oxygen atoms of the carbonyl groups of the 8-azasteroid molecule. A relation between the changes in the characteristic absorption bands of the CH₂-, C=O-, and I=C groups and the structure of the crystalline modifications has been established.     

The elastic theory of a single DNA molecule

We study the elastic responses of double-(ds) and single-stranded (ss) DNA at external force fields. A double-strand-polymer elastic model is constructed and solved by path integral methods and Monte Carlo simulations to understand the entropic elasticity, cooperative extensibility, and supercoiling property of dsDNA. The good agreement with experiments indicates that short-ranged base-pair stacking interaction is crucial for the stability and the high deformability of dsDNA. Hairpin-coil transition in ssDNA is studied with generating function method. A threshold force is needed to pull the ssDNA hairpin patterns, stabilized by base pairing and base-pair stacking, into random coils. This phase transition is predicted to be of first order for stacking potential higher than some critical level, in accordance with experimental observations.     

Millimeter and submillimeter wave spectroscopic investigations into the rotation-tunneling spectrum of gGg ethylene glycol, HOCH2CH2OH

The rotation-tunneling spectrum of the second most stable gGg^′ conformer of ethylene glycol (1,2-ethanediol) in its ground vibrational state has been studied in selected regions between 77 and 579 GHz. Compared to the study of the more stable aGg^′ conformer, a much larger frequency range was studied, resulting in a much extended frequency list covering similar quantum numbers, J?55 and K_a?19. While the input data were reproduced within experimental uncertainties up to moderately high values of J and K_a larger residuals remain at higher quantum numbers. The severe mixing of the states caused by the Coriolis interaction between the two tunneling substates is suggested to provide a considerable part of the explanation. In addition, a Coriolis interaction of the gGg^′ ground vibrational state with an excited state of the aGg^′ conformer may also contribute. Relative intensities of closely spaced lines have been investigated to determine the signs of the Coriolis constants between the two tunneling substates relative to the dipole moment components and to estimate the magnitudes of the dipole moment components and the energy difference between the gGg^′ and the aGg^′ conformers. Results of ab initio calculations on the total dipole moment and the vibrational spectrum were needed for these estimates. The current analysis is limited to transitions with quantum numbers J?40 and K_a?6 plus those having J?22 and K_a?17 which could be reproduced within experimental uncertainties. The results are aimed at aiding radioastronomers to search for gGg^′ ethylene glycol in comets and in interstellar space.