Adaptive Finite Element Methods with Inexact Solvers for the Nonlinear Poisson-Boltzmann Equation

In this article we study adaptive finite element methods (AFEM) with inexact solvers for a class of semilinear elliptic interface problems. We are particularly interested in nonlinear problems with discontinuous diffusion coefficients, such as the nonlinear Poisson-Boltzmann equation and its regularizations. The algorithm we study consists of the standard SOLVE-ESTIMATE-MARK-REFINE procedure common to many adaptive finite element algorithms, but where the SOLVE step involves only a full solve on the coarsest level, and the remaining levels involve only single Newton updates to the previous approximate solution. We summarize a recently developed AFEM convergence theory for inexact solvers, and present a sequence of numerical experiments that give evidence that the theory does in fact predict the contraction properties of AFEM with inexact solvers. The various routines used are all designed to maintain a linear-time computational complexity.Comment: Submitted to DD20 Proceeding

Bound, virtual and resonance SS-matrix poles from the Schr\"odinger equation

A general method, which we call the potential $S$-matrix pole method, is developed for obtaining the $S$-matrix pole parameters for bound, virtual and resonant states based on numerical solutions of the Schr\"odinger equation. This method is well-known for bound states. In this work we generalize it for resonant and virtual states, although the corresponding solutions increase exponentially when $r\to\infty$. Concrete calculations are performed for the $1^+$ ground and the $0^+$ first excited states of $^{14}\rm{N}$, the resonance $^{15}\rm{F}$ states ($1/2^+$, $5/2^+$), low-lying states of $^{11}\rm{Be}$ and $^{11}\rm{N}$, and the subthreshold resonances in the proton-proton system. We also demonstrate that in the case the broad resonances their energy and width can be found from the fitting of the experimental phase shifts using the analytical expression for the elastic scattering $S$-matrix. We compare the $S$-matrix pole and the $R$-matrix for broad $s_{1/2}$ resonance in ${}^{15}{\rm F}$Comment: 14 pages, 5 figures (figures 3 and 4 consist of two figures each) and 4 table

Experimental study of the decay PHI(1020)--->ETA+GAMMA in multi-photon final state

In the SND experiment at VEPP-2M e^+e^- collider the phi(1020) ---> eta+gamma ---> 3pi^0+gamma decay was studied. The branching ratio B(phi--->eta+gamma) = (1.246 +- 0.025 +- 0.057)% was measured.Comment: 4 pages, 1 figure, 1 table, Submitted to JETP Let

Morphology and Orientation Selection of Non-Metallic Inclusions in Electrified Molten Metal

The effect of electric current on morphology and orientation selection of non-metallic inclusions in molten metal has been investigated using theoretical modelling and numerical calculation. Two geometric factors, namely the circularity (fc) and alignment ratio (fe) were introduced to describe the inclusions shape and configuration. Electric current free energy was calculated and the values were used to determine the thermodynamic preference between different microstructures. Electric current promotes the development of inclusion along the current direction by either expatiating directional growth or enhancing directional agglomeration. Reconfiguration of the inclusions to reduce the system electric resistance drives the phenomena. The morphology and orientation selection follows the routine to reduce electric free energy. The numerical results are in agreement with our experimental observations

Production of intermediate-mass dileptons in relativistic heavy ion collisions

The production of intermediate mass dileptons in ultrarelativistic nuclear collisions at SPS energies is studied. The acceptance and detector resolution inherent to measurements by the NA50 experimental collaboration are accurately modeled. The measured centrality dependence of the intermediate mass lepton pair excess is also addressed.Comment: 9 pages, 8 figures, ReVTe

Study of Δ(1232)\Delta(1232) isobar electroproduction at VEPP-2M e+e−e^+e^- collider

Results from the Spherical Nonmagnetic Detector (SND) on $\Delta (1232)$ isobar electroproduction in the collisions of beam electrons (positrons) and residual gas nuclei in the VEPP-2M $e^+e^-$ collider are presented. On the basis of the obtained data the expected counting rate of this process in future high luminosity $e^+e^-$ colliders (~$\phi$-, $c$-$\tau$- and $b$-factories) was estimated.Comment: 7 pages LATEX and 3 figure

Formalism for dilepton production via virtual photon bremsstrahlung in hadronic reactions

We derive a set of new formulas for various distributions in dilepton production via virtual photon bremsstrahlung from pseudoscalar mesons and unpolarized spin-one-half fermions. These formulas correspond to the leading and sub-leading terms in the Low-Burnett-Kroll expansion for real photon bremsstrahlung. The relation of our leading-term formulas to previous works is also shown. Existing formulas are examined in the light of Lorentz covariance and gauge invariance. Numerical comparison is made in a simple example, where an "exact" formula and real photon data exist. The results reveal large discrepancies among different bremsstrahlung formulas. Of all the leading-term bremsstrahlung formulas, the one derived in this work agrees best with the exact formula. The issues of M_T-scaling and event generators are also addressed.Comment: 37 pages, RevTeX, epsf.sty, 10 embedded figure

Evidence for Factorization in Three-body B --> D(*) K- K0 Decays

Motivated by recent experimental results, we use a factorization approach to study the three-body B --> D(*) K- K0 decay modes. Two mechanisms are proposed for kaon pair production: current-produced (from vacuum) and transition (from B meson). The Bbar0 --> D(*)+ K- K0 decay is governed solely by the current-produced mechanism. As the kaon pair can be produced only by the vector current, the matrix element can be extracted from e+ e- --> K Kbar processes via isospin relations. The decay rates obtained this way are in good agreement with experiment. Both current-produced and transition processes contribute to B- --> D(*)0 K- K0 decays. By using QCD counting rules and the measured B- --> D(*)0 K- K0 decay rates, the measured decay spectra can be understood.Comment: 17 pages, 6 figure

Study of the process e+e−→π+π−π0e^+e^- \to \pi^+\pi^-\pi^0 in the energy region s\sqrt[]{s} below 0.98 GeV

The cross section of the process $e^+e^-\to \pi^+\pi^-\pi^0$ was measured in the Spherical Neutral Detector (SND) experiment at the VEPP-2M collider in the energy region $\sqrt[]{s}$ below 980 MeV. This measurement was based on about $1.2 \times 10^6$ selected events. The obtained cross section was analyzed together with the SND and DM2 data in the energy region $\sqrt[]{s}$ up to 2 GeV. The $\omega$-meson parameters: $m_\omega=782.79\pm 0.08\pm 0.09$ MeV, $\Gamma_\omega=8.68\pm 0.04\pm 0.15$ MeV and $\sigma(\omega\to 3\pi)=1615\pm 9\pm 57$ nb were obtained. It was found that the experimental data cannot be described by a sum of only $\omega$, $\phi$, $\omega^\prime$ and $\omega^{\prime\prime}$ resonances contributions. This can be interpreted as a manifestation of $\rho\to 3\pi$ decay, suppressed by $G$-parity, with relative probability $B(\rho\to 3\pi) = (1.01\pm^{0.54}_{0.36}\pm 0.034) \times 10^{-4}$.Comment: 41 pages REVTEX and 34 figure