Accurate Charge-Dependent Nucleon-Nucleon Potential at Fourth Order of Chiral Perturbation Theory

We present the first nucleon-nucleon potential at next-to-next-to-next-to-leading order (fourth order) of chiral perturbation theory. Charge-dependence is included up to next-to-leading order of the isospin-violation scheme. The accuracy for the reproduction of the NN data below 290 MeV lab. energy is comparable to the one of phenomenological high-precision potentials. Since NN potentials of order three and less are known to be deficient in quantitative terms, the present work shows that the fourth order is necessary and sufficient for a reliable NN potential derived from chiral effective Lagrangians. The new potential provides a promising starting point for exact few-body calculations and microscopic nuclear structure theory (including chiral many-body forces derived on the same footing).Comment: 4 pages Revtex including one figur

Centers of Mass and Rotational Kinematics for the Relativistic N-Body Problem in the Rest-Frame Instant Form

In the Wigner-covariant rest-frame instant form of dynamics it is possible to develop a relativistic kinematics for the N-body problem. The Wigner hyperplanes define the intrinsic rest frame and realize the separation of the center-of-mass. Three notions of {\it external} relativistic center of mass can be defined only in terms of the {\it external} Poincar\'e group realization. Inside the Wigner hyperplane, an {\it internal} unfaithful realization of the Poincar\'e group is defined. The three concepts of {\it internal} center of mass weakly {\it coincide} and are eliminated by the rest-frame conditions. An adapted canonical basis of relative variables is found. The invariant mass is the Hamiltonian for the relative motions. In this framework we can introduce the same {\it dynamical body frames}, {\it orientation-shape} variables, {\it spin frame} and {\it canonical spin bases} for the rotational kinematics developed for the non-relativistic N-body problem.Comment: 78 pages, revtex fil

Updated tests of scaling and universality for the spin-spin correlations in the 2D and 3D spin-S Ising models using high-temperature expansions

We have extended, from order 12 through order 25, the high-temperature series expansions (in zero magnetic field) for the spin-spin correlations of the spin-S Ising models on the square, simple-cubic and body-centered-cubic lattices. On the basis of this large set of data, we confirm accurately the validity of the scaling and universality hypotheses by resuming several tests which involve the correlation function, its moments and the exponential or the second-moment correlation-lengths.Comment: 21 pages, 8 figure

Ground-State of Charged Bosons Confined in a Harmonic Trap

We study a system composed of N identical charged bosons confined in a harmonic trap. Upper and lower energy bounds are given. It is shown in the large N limit that the ground-state energy is determined within an accuracy of $\pm 8$ and that the mean field theory provides a reasonable result with relative error of less than 16% for the binding energy .Comment: 15 page

Shell Model Monte Carlo studies of neutron-rich nuclei in the 1s-0d-1p-0f shells

We demonstrate the feasibility of realistic Shell-Model Monte Carlo (SMMC) calculations spanning multiple major shells, using a realistic interaction whose bad saturation and shell properties have been corrected by a newly developed general prescription. Particular attention is paid to the approximate restoration of translational invariance. The model space consists of the full sd-pf shells. We include in the study some well-known T=0 nuclei and several unstable neutron-rich ones around N=20,28. The results indicate that SMMC can reproduce binding energies, B(E2) transitions, and other observables with an interaction that is practically parameter free. Some interesting insight is gained on the nature of deep correlations. The validity of previous studies is confirmed.Comment: 22 pages + 7 postscript figure

Extension to order β23\beta^{23} of the high-temperature expansions for the spin-1/2 Ising model on the simple-cubic and the body-centered-cubic lattices

Using a renormalized linked-cluster-expansion method, we have extended to order $\beta^{23}$ the high-temperature series for the susceptibility $\chi$ and the second-moment correlation length $\xi$ of the spin-1/2 Ising models on the sc and the bcc lattices. A study of these expansions yields updated direct estimates of universal parameters, such as exponents and amplitude ratios, which characterize the critical behavior of $\chi$ and $\xi$. Our best estimates for the inverse critical temperatures are $\beta^{sc}_c=0.221654(1)$ and $\beta^{bcc}_c=0.1573725(6)$. For the susceptibility exponent we get $\gamma=1.2375(6)$ and for the correlation length exponent we get $\nu=0.6302(4)$. The ratio of the critical amplitudes of $\chi$ above and below the critical temperature is estimated to be $C_+/C_-=4.762(8)$. The analogous ratio for $\xi$ is estimated to be $f_+/f_-=1.963(8)$. For the correction-to-scaling amplitude ratio we obtain $a^+_{\xi}/a^+_{\chi}=0.87(6)$.Comment: Misprints corrected, 8 pages, latex, no figure

Gene and pathway identification with Lp penalized Bayesian logistic regression

<p>Abstract</p> <p>Background</p> <p>Identifying genes and pathways associated with diseases such as cancer has been a subject of considerable research in recent years in the area of bioinformatics and computational biology. It has been demonstrated that the magnitude of differential expression does not necessarily indicate biological significance. Even a very small change in the expression of particular gene may have dramatic physiological consequences if the protein encoded by this gene plays a catalytic role in a specific cell function. Moreover, highly correlated genes may function together on the same pathway biologically. Finally, in sparse logistic regression with <it>L</it>_{<it>p </it>}(<it>p </it>< 1) penalty, the degree of the sparsity obtained is determined by the value of the regularization parameter. Usually this parameter must be carefully tuned through cross-validation, which is time consuming.</p> <p>Results</p> <p>In this paper, we proposed a simple Bayesian approach to integrate the regularization parameter out analytically using a new prior. Therefore, there is no longer a need for parameter selection, as it is eliminated entirely from the model. The proposed algorithm (BLpLog) is typically two or three orders of magnitude faster than the original algorithm and free from bias in performance estimation. We also define a novel similarity measure and develop an integrated algorithm to hunt the regulatory genes with low expression changes but having high correlation with the selected genes. Pathways of those correlated genes were identified with DAVID <url>http://david.abcc.ncifcrf.gov/</url>.</p> <p>Conclusion</p> <p>Experimental results with gene expression data demonstrate that the proposed methods can be utilized to identify important genes and pathways that are related to cancer and build a parsimonious model for future patient predictions.</p

Improved high-temperature expansion and critical equation of state of three-dimensional Ising-like systems

High-temperature series are computed for a generalized $3d$ Ising model with arbitrary potential. Two specific ``improved'' potentials (suppressing leading scaling corrections) are selected by Monte Carlo computation. Critical exponents are extracted from high-temperature series specialized to improved potentials, achieving high accuracy; our best estimates are: $\gamma=1.2371(4)$, $\nu=0.63002(23)$, $\alpha=0.1099(7)$, $\eta=0.0364(4)$, $\beta=0.32648(18)$. By the same technique, the coefficients of the small-field expansion for the effective potential (Helmholtz free energy) are computed. These results are applied to the construction of parametric representations of the critical equation of state. A systematic approximation scheme, based on a global stationarity condition, is introduced (the lowest-order approximation reproduces the linear parametric model). This scheme is used for an accurate determination of universal ratios of amplitudes. A comparison with other theoretical and experimental determinations of universal quantities is presented.Comment: 65 pages, 1 figure, revtex. New Monte Carlo data by Hasenbusch enabled us to improve the determination of the critical exponents and of the equation of state. The discussion of several topics was improved and the bibliography was update

Kernel based methods for accelerated failure time model with ultra-high dimensional data

<p>Abstract</p> <p>Background</p> <p>Most genomic data have ultra-high dimensions with more than 10,000 genes (probes). Regularization methods with <it>L</it>₁and <it>L_p</it>penalty have been extensively studied in survival analysis with high-dimensional genomic data. However, when the sample size <it>n </it>≪ <it>m </it>(the number of genes), directly identifying a small subset of genes from ultra-high (<it>m </it>> 10, 000) dimensional data is time-consuming and not computationally efficient. In current microarray analysis, what people really do is select a couple of thousands (or hundreds) of genes using univariate analysis or statistical tests, and then apply the LASSO-type penalty to further reduce the number of disease associated genes. This two-step procedure may introduce bias and inaccuracy and lead us to miss biologically important genes.</p> <p>Results</p> <p>The accelerated failure time (AFT) model is a linear regression model and a useful alternative to the Cox model for survival analysis. In this paper, we propose a nonlinear kernel based AFT model and an efficient variable selection method with adaptive kernel ridge regression. Our proposed variable selection method is based on the kernel matrix and dual problem with a much smaller <it>n </it>× <it>n </it>matrix. It is very efficient when the number of unknown variables (genes) is much larger than the number of samples. Moreover, the primal variables are explicitly updated and the sparsity in the solution is exploited.</p> <p>Conclusions</p> <p>Our proposed methods can simultaneously identify survival associated prognostic factors and predict survival outcomes with ultra-high dimensional genomic data. We have demonstrated the performance of our methods with both simulation and real data. The proposed method performs superbly with limited computational studies.</p

Nucleon-Nucleon Interaction: A Typical/Concise Review

Nearly a recent century of work is divided to Nucleon-Nucleon (NN) interaction issue. We review some overall perspectives of NN interaction with a brief discussion about deuteron, general structure and symmetries of NN Lagrangian as well as equations of motion and solutions. Meanwhile, the main NN interaction models, as frameworks to build NN potentials, are reviewed concisely. We try to include and study almost all well-known potentials in a similar way, discuss more on various commonly used plain forms for two-nucleon interaction with an emphasis on the phenomenological and meson-exchange potentials as well as the constituent-quark potentials and new ones based on chiral effective field theory and working in coordinate-space mostly. The potentials are constructed in a way that fit NN scattering data, phase shifts, and are also compared in this way usually. An extra goal of this study is to start comparing various potentials forms in a unified manner. So, we also comment on the advantages and disadvantages of the models and potentials partly with reference to some relevant works and probable future studies.Comment: 85 pages, 5 figures, than the previous v3 edition, minor changes, and typos fixe