Control of spiral waves and turbulent states in a cardiac model by travelling-wave perturbations

We propose a travelling-wave perturbation method to control the spatiotemporal dynamics in a cardiac model. It is numerically demonstrated that the method can successfully suppress the wave instability (alternans in action potential duration) in the one-dimensional case and convert spiral waves and turbulent states to the normal travelling wave states in the two-dimensional case. An experimental scheme is suggested which may provide a new design for a cardiac defibrillator.Comment: 9 pages, 5 figure

First-principles calculations of phase transition, elasticity, and thermodynamic properties for TiZr alloy

tructural transformation, pressure dependent elasticity behaviors, phonon, and thermodynamic properties of the equiatomic TiZr alloy are investigated by using first-principles density-functional theory. Our calculated lattice parameters and equation of state for $\alpha$ and $\omega$ phases as well as the phase transition sequence of $\alpha$$\mathtt{\rightarrow}$$\omega$$\mathtt{\rightarrow}$$\beta$ are consistent well with experiments. Elastic constants of $\alpha$ and $\omega$ phases indicate that they are mechanically stable. For cubic $\beta$ phase, however, it is mechanically unstable at zero pressure and the critical pressure for its mechanical stability is predicted to equal to 2.19 GPa. We find that the moduli, elastic sound velocities, and Debye temperature all increase with pressure for three phases of TiZr alloy. The relatively large $B/G$ values illustrate that the TiZr alloy is rather ductile and its ductility is more predominant than that of element Zr, especially in $\beta$ phase. Elastic wave velocities and Debye temperature have abrupt increase behaviors upon the $\alpha$$\mathtt{\rightarrow}$$\omega$ transition at around 10 GPa and exhibit abrupt decrease feature upon the $\omega$$\mathtt{\rightarrow}$$\beta$ transition at higher pressure. Through Mulliken population analysis, we illustrate that the increase of the \emph{d}-band occupancy will stabilize the cubic $\beta$ phase. Phonon dispersions for three phases of TiZr alloy are firstly presented and the $\beta$ phase phonons clearly indicate its dynamically unstable nature under ambient condition. Thermodynamics of Gibbs free energy, entropy, and heat capacity are obtained by quasiharmonic approximation and Debye model.Comment: 9 pages, 10 figure

Effects of a weakly interacting light U boson on the nuclear equation of state and properties of neutron stars in relativistic models

We investigate the effects of the light vector U-boson that couples weakly to nucleons in relativistic mean-field models on the equation of state and subsequently the consequence in neutron stars. It is analyzed that the U-boson can lead to a much clearer rise of the neutron star maximum mass in models with the much softer equation of state. The inclusion of the U-boson may thus allow the existence of the non-nucleonic degrees of freedom in the interior of large mass neutron stars initiated with the favorably soft EOS of normal nuclear matter. In addition, the sensitive role of the U-boson in the neutron star radius and its relation to the test of the non-Newtonian gravity that is herein addressed by the light U-boson are discussed.Comment: 8 pages,7 figure

An Integral Equation Method for the Cahn-Hilliard Equation in the Wetting Problem

We present an integral equation approach to solving the Cahn-Hilliard equation equipped with boundary conditions that model solid surfaces with prescribed Young's angles. The discretization of the system in time using convex splitting leads to a modified biharmonic equation at each time step. To solve it, we split the solution into a volume potential computed with free space kernels, plus the solution to a second kind integral equation (SKIE). The volume potential is evaluated with the help of a box-based volume-FMM method. For non-box domains, source density is extended by solving a biharmonic Dirichlet problem. The near-singular boundary integrals are computed using quadrature by expansion (QBX) with FMM acceleration. Our method has linear complexity in the number of surface/volume degrees of freedom and can achieve high order convergence with adaptive refinement to manage error from function extension