Monte Carlo Algorithm for Simulating Reversible Aggregation of Multisite Particles

We present an efficient and exact Monte Carlo algorithm to simulate reversible aggregation of particles with dedicated binding sites. This method introduces a novel data structure of dynamic bond tree to record clusters and sequences of bond formations. The algorithm achieves a constant time cost for processing cluster association and a cost between $\mathcal{O}(\log M)$ and $\mathcal{O}(M)$ for processing bond dissociation in clusters with $M$ bonds. The algorithm is statistically exact and can reproduce results obtained by the standard method. We applied the method to simulate a trivalent ligand and a bivalent receptor clustering system and obtained an average scaling of $\mathcal{O}(M^{0.45})$ for processing bond dissociation in acyclic aggregation, compared to a linear scaling with the cluster size in standard methods. The algorithm also demands substantially less memory than the conventional method.Comment: 8 pages, 3 figure

Semi-Inclusive B\to K(K^*) X Decays with Initial Bound State Effects

The effects of initial $b$ quark bound state for the semi-inclusive decays $B\to K(K^*) X$ are studied using light cone expansion and heavy quark effective theory methods. We find that the initial bound state effects on the branching ratios and CP asymmetries are small. In the light cone expansion approach, the CP-averaged branching ratios are increased by about 2% with respect to the free $b$-quark decay. For $\bar B^0 \to K^- (K^{*-}) X$, the CP-averaged branching ratios are sensitive to the phase $\gamma$ and the CP asymmetry can be as large as 7% (14%), whereas for $B^-\to \bar K^0 (\bar K^{*0})X$ the CP-averaged branching ratios are not sensitive to $\gamma$ and the CP asymmetries are small ($< 1$). The CP-averaged branching ratios are predicted to be in the ranges $(0.53 \sim 1.5)\times 10^{-4}$ [$(0.25 \sim 2.0)\times 10^{-4}$] for $\bar B^0 \to K^- (K^{*-})X$ and $(0.77 \sim 0.84)\times 10^{-4}$ [$(0.67 \sim 0.74)\times 10^{-4}$] for $B^-\to \bar K^0 (\bar K^{*0}) X$, depending on the value of the CP violating phase $\gamma$. In the heavy quark effective theory approach, we find that the branching ratios are decreased by about 10% and the CP asymmetries are not affected. These predictions can be tested in the near future.Comment: 29 pages, 12 ps figure

N_pN_n dependence of empirical formula for the lowest excitation energy of the 2^+ states in even-even nuclei

We examine the effects of the additional term of the type $\sim e^{- \lambda' N_pN_n}$ on the recently proposed empirical formula for the lowest excitation energy of the $2^+$ states in even-even nuclei. This study is motivated by the fact that this term carries the favorable dependence of the valence nucleon numbers dictated by the $N_pN_n$ scheme. We show explicitly that there is not any improvement in reproducing $E_x(2_1^+)$ by including the extra $N_pN_n$ term. However, our study also reveals that the excitation energies $E_x(2_1^+)$, when calculated by the $N_pN_n$ term alone (with the mass number $A$ dependent term), are quite comparable to those calculated by the original empirical formula.Comment: 14 pages, 5 figure

Evidence for Antiferromagnetic Order in La2−x_{2-x}Cex_{x}CuO4_{4} from Angular Magnetoresistance Measurements

We investigated the in-plane angular magnetoresistivity (AMR) of $T^{^{\prime}}$-phase La$_{2-x}$Ce$_{x}$CuO$_{4}$ (LCCO) thin films ($x=0.06-0.15$) fabricated by a pulsed laser deposition technique. The in-plane AMR with $\mathbf{H}\parallel ab$ shows a twofold symmetry instead of the fourfold behavior found in other electron-doped cuprates such as Pr$_{2-x}$Ce$_{x}$CuO$_{4}$ and Nd$_{2-x}$Ce$_{x}$CuO$_{4}$. The twofold AMR disappears above a certain temperature, $T_{D}$. The $T_{D}(x)$ is well above $T_{c}(x)$ for $x=0.06$ ($\sim 110$ K), and decreases with increasing doping, until it is no longer observed above $T_{c}(x)$ at $x=0.15$. This twofold AMR below $T_{D}(x)$ is suggested to originate from an antiferromagnetic or spin density wave order.Comment: to be published in Phys. Rev. B, Vol. 80 (2009

Analysis of hadronic invariant mass spectrum in inclusive charmless semileptonic B decays

We make an analysis of the hadronic invariant mass spectrum in inclusive charmless semileptonic B meson decays in a QCD-based approach. The decay width is studied as a function of the invariant mass cut. We examine their sensitivities to the parameters of the theory. The theoretical uncertainties in the determination of $|V_{ub}|$ from the hadronic invariant mass spectrum are investigated. A strategy for improving the theoretical accuracy in the value of $|V_{ub}|$ is described.Comment: 13 pages, 5 Postscript figure

Numerical approximation of the Euler-Poisson-Boltzmann model in the quasineutral limit

This paper analyzes various schemes for the Euler-Poisson-Boltzmann (EPB) model of plasma physics. This model consists of the pressureless gas dynamics equations coupled with the Poisson equation and where the Boltzmann relation relates the potential to the electron density. If the quasi-neutral assumption is made, the Poisson equation is replaced by the constraint of zero local charge and the model reduces to the Isothermal Compressible Euler (ICE) model. We compare a numerical strategy based on the EPB model to a strategy using a reformulation (called REPB formulation). The REPB scheme captures the quasi-neutral limit more accurately