Low-diffusivity scalar transport using a WENO scheme and dual meshing

Interfacial mass transfer of low-diffusive substances in an unsteady flow environment is marked by a very thin boundary layer at the interface and other regions with steep concentration gradients. A numerical scheme capable of resolving accurately most details of this process is presented. In this scheme, the fourth-order accurate WENO method developed by Liu et al. (1994) was implemented on a non-uniform staggered mesh to discretize the scalar convection while for the scalar diffusion a fourth-order accurate central discretization was employed. The discretization of the scalar convection-diffusion equation was combined with a fourth-order Navier-Stokes solver which solves the incompressible flow. A dual meshing strategy was employed, in which the scalar was solved on a finer mesh than the incompressible flow. The solver was tested by performing a number of two-dimensional simulations of an unstably stratified flow with low diffusivity scalar transport. The unstable stratification led to buoyant convection which was modelled using a Boussinesq approximation with a linear relationship between flow temperature and density. The order of accuracy for one-dimensional scalar transport on a stretched and uniform grid was also tested. The results show that for the method presented above a relatively coarse mesh is sufficient to accurately describe the fluid flow, while the use of a refined mesh for the low-diffusive scalars is found to be beneficial in order to obtain a highly accurate resolution with negligible numerical diffusion

Recent experimental results in sub- and near-barrier heavy ion fusion reactions

Recent advances obtained in the field of near and sub-barrier heavy-ion fusion reactions are reviewed. Emphasis is given to the results obtained in the last decade, and focus will be mainly on the experimental work performed concerning the influence of transfer channels on fusion cross sections and the hindrance phenomenon far below the barrier. Indeed, early data of sub-barrier fusion taught us that cross sections may strongly depend on the low-energy collective modes of the colliding nuclei, and, possibly, on couplings to transfer channels. The coupled-channels (CC) model has been quite successful in the interpretation of the experimental evidences. Fusion barrier distributions often yield the fingerprint of the relevant coupled channels. Recent results obtained by using radioactive beams are reported. At deep sub-barrier energies, the slope of the excitation function in a semi-logarithmic plot keeps increasing in many cases and standard CC calculations over-predict the cross sections. This was named a hindrance phenomenon, and its physical origin is still a matter of debate. Recent theoretical developments suggest that this effect, at least partially, may be a consequence of the Pauli exclusion principle. The hindrance may have far-reaching consequences in astrophysics where fusion of light systems determines stellar evolution during the carbon and oxygen burning stages, and yields important information for exotic reactions that take place in the inner crust of accreting neutron stars.Comment: 40 pages, 63 figures, review paper accepted for EPJ

Search for Invisible Decays of η\eta and η′\eta^\prime in J/ψ→ϕηJ/\psi \to \phi\eta and ϕη′\phi \eta^\prime

Using a data sample of $58\times 10^6$ $J/\psi$ decays collected with the BES II detector at the BEPC, searches for invisible decays of $\eta$ and $\eta^\prime$ in $J/\psi$ to $\phi\eta$ and $\phi\eta^\prime$ are performed. The $\phi$ signals, which are reconstructed in $K^+K^-$ final states, are used to tag the $\eta$ and $\eta^\prime$ decays. No signals are found for the invisible decays of either $\eta$ or $\eta^\prime$, and upper limits at the 90% confidence level are determined to be $1.65 \times 10^{-3}$ for the ratio $\frac{B(\eta\to \text{invisible})}{B(\eta\to\gamma\gamma)}$ and $6.69\times 10^{-2}$ for $\frac{B(\eta^\prime\to \text{invisible})}{B(\eta^\prime\to\gamma\gamma)}$. These are the first searches for $\eta$ and $\eta^\prime$ decays into invisible final states.Comment: 5 pages, 4 figures; Added references, Corrected typo

Semiparametric regression analysis for composite endpoints subject to componentwise censoring

Composite endpoints with censored data are commonly used as study outcomes in clinical trials. For example, progression-free survival is a widely used composite endpoint, with disease progression and death as the two components. Progression-free survival time is often defined as the time from randomization to the earlier occurrence of disease progression or death from any cause. The censoring times of the two components could be different for patients not experiencing the endpoint event. Conventional approaches, such as taking the minimum of the censoring times of the two components as the censoring time for progression-free survival time, may suffer from efficiency loss and could produce biased estimates of the treatment effect. We propose a new likelihood-based approach that decomposes the endpoints and models both the progression-free survival time and the time from disease progression to death. The censoring times for different components are distinguished. The approach makes full use of available information and provides a direct and improved estimate of the treatment effect on progression-free survival time. Simulations demonstrate that the proposed method outperforms several other approaches and is robust against various model misspecifications. An application to a prostate cancer clinical trial is provided

Negatively Charged Excitons and Photoluminescence in Asymmetric Quantum Well

We study photoluminescence (PL) of charged excitons ($X^-$) in narrow asymmetric quantum wells in high magnetic fields B. The binding of all $X^-$ states strongly depends on the separation $\delta$ of electron and hole layers. The most sensitive is the ``bright'' singlet, whose binding energy decreases quickly with increasing $\delta$ even at relatively small B. As a result, the value of B at which the singlet--triplet crossing occurs in the $X^-$ spectrum also depends on $\delta$ and decreases from 35 T in a symmetric 10 nm GaAs well to 16 T for $\delta=0.5$ nm. Since the critical values of $\delta$ at which different $X^-$ states unbind are surprisingly small compared to the well width, the observation of strongly bound $X^-$ states in an experimental PL spectrum implies virtually no layer displacement in the sample. This casts doubt on the interpretation of PL spectra of heterojunctions in terms of $X^-$ recombination

GLMMLasso: An Algorithm for High-Dimensional Generalized Linear Mixed Models Using L1-Penalization

We propose an L1-penalized algorithm for fitting high-dimensional generalized linear mixed models. Generalized linear mixed models (GLMMs) can be viewed as an extension of generalized linear models for clustered observations. This Lasso-type approach for GLMMs should be mainly used as variable screening method to reduce the number of variables below the sample size. We then suggest a refitting by maximum likelihood based on the selected variables only. This is an effective correction to overcome problems stemming from the variable screening procedure which are more severe with GLMMs. We illustrate the performance of our algorithm on simulated as well as on real data examples. Supplemental materials are available online and the algorithm is implemented in the R package glmmixedlasso