title: The geometric origin of Lie point symmetries of the Schr\"{o}dinger and the Klein Gordon equations

abstract: We determine the Lie point symmetries of the Schr\"{o}dinger and the Klein Gordon equations in a general Riemannian space. It is shown that these symmetries are related with the homothetic and the conformal algebra of the metric of the space respectively. We consider the kinematic metric defined by the classical Lagrangian and show how the Lie point symmetries of the Schr\"{o}dinger equation and the Klein Gordon equation are related with the Noether point symmetries of this Lagrangian. The general results are applied to two practical problems a. The classification of all two and three dimensional potentials in a Euclidian space for which the Schr\"{o}dinger equation and the Klein Gordon equation admit Lie point symmetries and b. The application of Lie point symmetries of the Klein Gordon equation in the exterior Schwarzschild spacetime and the determination of the metric by means of conformally related Lagrangians.


title: Kepler-91b: a planet at the end of its life. Planet and giant host star properties via light-curve variations

abstract: The evolution of planetary systems is intimately linked to the evolution of their host star. Our understanding of the whole planetary evolution process is based on the large planet diversity observed so far. To date, only few tens of planets have been discovered orbiting stars ascending the Red Giant Branch. Although several theories have been proposed, the question of how planets die remains open due to the small number statistics. In this work we study the giant star Kepler-91 (KOI-2133) in order to determine the nature of a transiting companion. This system was detected by the Kepler Space Telescope. However, its planetary confirmation is needed. We confirm the planetary nature of the object transiting the star Kepler-91 by deriving a mass of $ M_p=0.88^{+0.17}_{-0.33} ~M_{\rm Jup}$ and a planetary radius of $R_p=1.384^{+0.011}_{-0.054} ~R_{\rm Jup}$. Asteroseismic analysis produces a stellar radius of $R_{\star}=6.30\pm 0.16 ~R_{\odot}$ and a mass of $M_{\star}=1.31\pm 0.10 ~ M_{\odot} $. We find that its eccentric orbit ($e=0.066^{+0.013}_{-0.017}$) is just $1.32^{+0.07}_{-0.22} ~ R_{\star}$ away from the stellar atmosphere at the pericenter. Kepler-91b could be the previous stage of the planet engulfment, recently detected for BD+48 740. Our estimations show that Kepler-91b will be swallowed by its host star in less than 55 Myr. Among the confirmed planets around giant stars, this is the planetary-mass body closest to its host star. At pericenter passage, the star subtends an angle of $48^{\circ}$, covering around 10% of the sky as seen from the planet. The planetary atmosphere seems to be inflated probably due to the high stellar irradiation.


title: Discovery of a transparent sightline at {\rho} < 20 kpc from an interacting pair of galaxies

abstract: We report the discovery of a transparent sightline at projected distances of {\rho} < 20 kpc to an interacting pair of mature galaxies at z = 0.12. The sightline of the UV-bright quasar PG1522+101 at z = 1.328 passes at {\rho} = 11.5 kpc from the higher-mass galaxy (M_* = 10^10.6 M_Sun) and {\rho} = 20.4 kpc from the lower-mass one (M_* = 10^10.0 M_Sun). The two galaxies are separated by 9 kpc in projected distance and 30 km/s in line-of-sight velocity. Deep optical images reveal tidal features indicative of close interactions. Despite the small projected distances, the quasar sightline shows little absorption associated with the galaxy pair with a total HI column density it no greater than log N(HI) = 13.65. This limiting HI column density is already two orders-of-magnitude less than what is expected from previous halo gas studies. In addition, we detect no heavy-element absorption features associated with the galaxy pair with 3-{\sigma} limits of log N(MgII) < 12.2 and log N(OVI) < 13.7. The probability of seeing such little absorption in a sightline passing at a small projected distance from two non-interacting galaxies is 0.2%. The absence of strong absorbers near the close galaxy pair suggests that the cool gas reservoirs of the galaxies have been significantly depleted by the galaxy interaction. These observations therefore underscore the potential impact of galaxy interactions on the gaseous halos around galaxies.


title: Combined Flux and Anisotropy Searches Improve Sensitivity to Gamma Rays from Dark Matter

abstract: Searches for dark matter annihilation signals in the diffuse gamma-ray background are promising. We present the first comprehensive study utilizing both flux and anisotropy, using the example of a spectral line signal. Besides improving sensitivity, a combined analysis will help separately determine the particle properties of dark matter and the cosmological aspects of its clumping into halo substructure. The significance of a signal in a shot-noise-dominated anisotropy analysis increases linearly with time $t$, as opposed to $\sqrt{t}$ for a flux analysis, so a flux hint might be confirmed with an anisotropy signal. A first combined line search with Fermi-LAT would provide powerful new sensitivity; one with the proposed GAMMA-400 would be dramatically better.

categories: astro-ph.HE astro-ph.CO astro-ph.GA hep-ex hep-ph


title: Evaluating Gyrochronology on the Zero-Age-Main-Sequence: Rotation Periods in the Southern Open Cluster Blanco 1 from the KELT-South Survey

abstract: We report periods for 33 members of Blanco 1 as measured from KELT-South light curves, the first reported rotation periods for this benchmark zero-age-main-sequence open cluster. The distribution of these stars spans from late-A or early-F dwarfs to mid-K with periods ranging from less than a day to ~8 days. The rotation period distribution has a morphology similar to the coeval Pleiades cluster, suggesting the universal nature of stellar rotation distributions. Employing two different gyrochronology methods, we find an age of 146+13-14 Myr for the cluster. Using the same techniques, we infer an age of 134+9-10 Myr for the Pleiades measured from existing literature rotation periods. These rotation-derived ages agree with independently determined cluster ages based on the lithium depletion boundary technique. Additionally, we evaluate different gyrochronology models, and quantify levels of agreement between the models and the Blanco 1/Pleiades rotation period distributions, including incorporating the rotation distributions of clusters at ages up to 1.1 Gyr. We find the Skumanich-like spin-down rate sufficiently describes the rotation evolution of stars hotter than the Sun; however, we find cooler stars rotating faster than predicted by a Skumanich-law, suggesting a mass dependence in the efficiency of stellar angular momentum loss rate. Finally, we compare the Blanco 1 and Pleiades rotation period distributions to available non-linear angular momentum evolution models. We find they require a significant mass dependence on the initial rotation rate of solar-type stars to reproduce the observed range of rotation periods at a given stellar mass, and are furthermore unable to predict the observed over-density of stars along the upper-envelope of the clusters' rotation distributions.


title: Soft X-ray Excess in the Coma Cluster from a Cosmic Axion Background

abstract: We show that the soft X-ray excess in the Coma cluster can be explained by a cosmic background of relativistic axions converting into photons in the cluster magnetic field. We provide a detailed self-contained review of the cluster soft X-ray excess, the proposed astrophysical explanations and the problems they face, and explain how a 0.1-1 keV axion background naturally arises at reheating in many string theory models of the early universe. We study the morphology of the soft excess by numerically propagating axions through stochastic, multi-scale magnetic field models that are consistent with observations of Faraday rotation measures from Coma. By comparing to ROSAT observations of the 0.2-0.4 keV soft excess, we find that the overall excess luminosity is easily reproduced for $g_{a\gamma\gamma} \sim 2 \times 10^{-13}$ GeV$^{-1}$. The resulting morphology is highly sensitive to the magnetic field power spectrum. For Gaussian magnetic field models, the observed soft excess morphology prefers magnetic field spectra with most power in coherence lengths on ${\cal O}(3 {\rm ~kpc})$ scales over those with most power on ${\cal O}(12 {\rm ~kpc})$ scales. Within this scenario, we bound the mean energy of the axion background to $50\, {\rm eV}\lesssim \langle E_a \rangle \lesssim 250\, {\rm eV}$, the axion mass to $m_a \lesssim 10^{-12}\,\hbox{eV}$, and derive a lower bound on the axion-photon coupling $g_{a\gamma\gamma} \gtrsim \sqrt{0.5/\Delta N_{\rm eff}}\, 1.4 \times 10^{-13}$ GeV$^{-1}$.


title: Compressed convolution

abstract: We introduce the concept of compressed convolution, a technique to convolve a given data set with a large number of non-orthogonal kernels. In typical applications our technique drastically reduces the effective number of computations. The new method is applicable to convolutions with symmetric and asymmetric kernels and can be easily controlled for an optimal trade-off between speed and accuracy. It is based on linear compression of the collection of kernels into a small number of coefficients in an optimal eigenbasis. The final result can then be decompressed in constant time for each desired convolved output. The method is fully general and suitable for a wide variety of problems. We give explicit examples in the context of simulation challenges for upcoming multi-kilo-detector cosmic microwave background (CMB) missions. For a CMB experiment with O(10,000) detectors with similar beam properties, we demonstrate that the algorithm can decrease the costs of beam convolution by two to three orders of magnitude with negligible loss of accuracy. Likewise, it has the potential to allow the reduction of disk space required to store signal simulations by a similar amount. Applications in other areas of astrophysics and beyond are optimal searches for a large number of templates in noisy data, e.g. from a parametrized family of gravitational wave templates; or calculating convolutions with highly overcomplete wavelet dictionaries, e.g. in methods designed to uncover sparse signal representations.


title: The Next-to-Minimal Two Higgs Doublet Model

abstract: The simplest extension of the Two Higgs Doublet Model is the addition of a real scalar singlet, S. The effects of mixing between the singlet and the doublets can be manifested in two ways. It can modify the couplings of the 126 GeV Higgs boson, h, and it can lead to direct detection of the heavy Higgs at the LHC. In this paper, we show that in the type-I model, for heavy Higgs masses in the 200-600 GeV range, the latter effect will be detected earlier than the former for most of parameter space. Should no such Higgs be discovered in this mass range, then the upper limit on the mixing will be sufficiently strong such that there will be no significant effects on the couplings of the h for most of parameter space. The reverse is true in the type-II model, the limits from measurements of the couplings of the h will dominate over the limits from non-observation of the heavy Higgs.


title: X-ray properties of K-selected galaxies at 0.5<z<2.0: Investigating trends with stellar mass, redshift and spectral type

abstract: We examine how the total X-ray luminosity correlates with stellar mass, stellar population, and redshift for a K-band limited sample of ~3500 galaxies at 0.5<z<2.0 from the NEWFIRM Medium Band Survey in the COSMOS field. The galaxy sample is divided into 32 different galaxy types, based on similarities between the spectral energy distributions. For each galaxy type, we further divide the sample into bins of redshift and stellar mass, and perform an X-ray stacking analysis using the Chandra COSMOS (C-COSMOS) data. We find that full band X-ray luminosity is primarily increasing with stellar mass, and at similar mass and spectral type is higher at larger redshifts. When comparing at the same stellar mass, we find that the X-ray luminosity is slightly higher for younger galaxies (i.e., weaker 4000\AA breaks), but the scatter in this relation is large. We compare the observed X-ray luminosities to those expected from low and high mass X-ray binaries (XRBs). For blue galaxies, XRBs can almost fully account for the observed emission, while for older galaxies with larger 4000\AA breaks, active galactic nuclei (AGN) or hot gas dominate the measured X-ray flux. After correcting for XRBs, the X-ray luminosity is still slightly higher in younger galaxies, although this correlation is not significant. AGN appear to be a larger component of galaxy X-ray luminosity at earlier times, as the hardness ratio increases with redshift. Together with the slight increase in X-ray luminosity this may indicate more obscured AGNs or higher accretion rates at earlier times.


title: A Sub-Earth-Mass Moon Orbiting a Gas Giant Primary or a High Velocity Planetary System in the Galactic Bulge

abstract: We present the first microlensing candidate for a free-floating exoplanet-exomoon system, MOA-2011-BLG-262, with a primary lens mass of M_host ~ 4 Jupiter masses hosting a sub-Earth mass moon. The data are well fit by this exomoon model, but an alternate star+planet model fits the data almost as well. Nevertheless, these results indicate the potential of microlensing to detect exomoons, albeit ones that are different from the giant planet moons in our solar system. The argument for an exomoon hinges on the system being relatively close to the Sun. The data constrain the product M pi_rel, where M is the lens system mass and pi_rel is the lens-source relative parallax. If the lens system is nearby (large pi_rel), then M is small (a few Jupiter masses) and the companion is a sub-Earth-mass exomoon. The best-fit solution has a large lens-source relative proper motion, mu_rel = 19.6 +- 1.6 mas/yr, which would rule out a distant lens system unless the source star has an unusually high proper motion. However, data from the OGLE collaboration nearly rule out a high source proper motion, so the exoplanet+exomoon model is the favored interpretation for the best fit model. However, the alternate solution has a lower proper motion, which is compatible with a distant (so stellar) host. A Bayesian analysis does not favor the exoplanet+exomoon interpretation, so Occam's razor favors a lens system in the bulge with host and companion masses of M_host = 0.12 (+0.19 -0.06) M_solar and m_comp = 18 (+28 -100 M_earth, at a projected separation of a_perp ~ 0.84 AU. The existence of this degeneracy is an unlucky accident, so current microlensing experiments are in principle sensitive to exomoons. In some circumstances, it will be possible to definitively establish the low mass of such lens systems through the microlensing parallax effect. Future experiments will be sensitive to less extreme exomoons.


title: On a shadow system of the S-K-T competition system

abstract: We study a boundary value problem with an integral constraint that arises from the modelings of species competition proposed by Lou and Ni in \cite{LN2}. Through bifurcation theories, we obtain the existence of non-constant positive solutions over one-dimensional domain, which are small perturbations from the positive constant solution of the system. Moreover, we determine the stability of the bifurcating solutions. Finally, for the diffusion rate being sufficiently small, we construct infinitely many positive solutions with single transition layer, which is represented as an approximation of a step function that. The transition-layer solution can be used to model the phenomenon of species segregation through competitions.


title: Green Bank Telescope observations of low column density HI around NGC 2997 and NGC 6946

abstract: Observations of ongoing HI accretion in nearby galaxies have only identified about 10% of the needed fuel to sustain star formation in these galaxies. Most of these observations have been conducted using interferometers and may have missed lower column density, diffuse, HI gas that may trace the missing 90% of gas. Such gas may represent the so-called "cold flows" predicted by current theories of galaxy formation to have never been heated above the virial temperature of the dark matter halo. As a first attempt to identify such cold flows around nearby galaxies and complete the census of HI down to N(HI)~10^18 cm^-2, I used the Robert C. Byrd Green Bank Telescope (GBT) to map the circumgalactic (r < 100-200 kpc) HI environment around NGC 2997 and NGC 6946. The resulting GBT observations cover a four square degree area around each galaxy with a 5-sigma detection limit of N(HI)~10^18 cm^-2 over a 20 km/s linewidth. This project complements absorption line studies, which are well-suited to the regime of lower N(HI). Around NGC 2997, the GBT HI data reveal an extended HI disk and all of its surrounding gas-rich satellite galaxies, but no filamentary features. Furthermore, the HI mass as measured with the GBT is only 7% higher than past interferometric measurements. After correcting for resolution differences, the HI extent of the galaxy is 23% larger at the N(HI)~1.2x10^18 cm^-2 level as measured by the GBT. On the other hand, the HI observations of NGC 6946 reveal a filamentary feature apparently connecting NGC 6946 with its nearest companions. This HI filament has N(HI)~10^18 cm^-2 and a FWHM of 55+-5 km/s and was invisible in past interferometer observations. The properties of this filament are broadly consistent with being a cold flow or debris from a past tidal interaction between NGC 6946 and its satellites.


title: Geometric mutual information at classical critical points

abstract: A practical use of the entanglement entropy in a 1d quantum system is to identify the conformal field theory describing its critical behavior. It is exactly $(c/3)\ln \ell$ for an interval of length $\ell$ in an infinite system, where $c$ is the central charge of the conformal field theory. Here we define the geometric mutual information, an analogous quantity for classical critical points. We compute this for 2d conformal field theories in an arbitrary geometry, and show in particular that for a rectangle cut into two rectangles, it is proportional to $c$. This makes it possible to extract $c$ in classical simulations, which we demonstrate for the critical Ising and 3-state Potts models.


title: On the Asymptotic Statistics of the Number of Occurrences of Multiple Permutation Patterns

abstract: We study statistical properties of the random variables $X_{\sigma}(\pi)$, the number of occurrences of the pattern $\sigma$ in the permutation $\pi$. We present two contrasting approaches to this problem: traditional probability theory and the ``less traditional'' computational approach. Through the perspective of the first one, we prove that for any pair of patterns $\sigma$ and $\tau$, the random variables $X_{\sigma}$ and $X_{\tau}$ are jointly asymptotically normal (when the permutation is chosen from $S_{n}$). From the other perspective, we develop algorithms that can show asymptotic normality and joint asymptotic normality (up to a point) and derive explicit formulas for quite a few moments and mixed moments empirically, yet rigorously. The computational approach can also be extended to the case where permutations are drawn from a set of pattern avoiders to produce many empirical moments and mixed moments. This data suggests that some random variables are not asymptotically normal in this setting.


title: Cosmological singularity theorems and black holes

abstract: An extension of Penrose's singularity theorem is proved for spacetimes where black holes are allowed to form from non-singular initial data. With standard assumptions about the spacetime, and assuming the existence of a trapped surface which lies outside of black hole horizons and is not completely surrounded by horizons, we show that the spacetime region outside (or on) the horizons must contain singularities. If the trapped surface is surrounded by horizons, we show that the horizons divide spacetime into causally disconnected pieces. Unlike the original Penrose's theorem, our theorems provide some information about the location of singularities. We illustrate how they can be used to rule out some cosmological scenarios.


title: Unit cell orientation of tetragonal-like BiFeO$_3$ thin films grown on highly miscut LaAlO$_3$ substrates

abstract: Synchrotron and lab-scale x-ray diffraction shows that tetragonal-like T'-BiFeO3 films on miscut LaAlO3 substrates (miscut < 5 deg) exhibit (00l)-planes tilted away from those of the substrate as predicted by the "Nagai model" (except for miscut <0.2 deg). Tilts as large as 1 deg are achieved even in 100 nm thick films, strikingly larger than those observed in other perovskites. We attribute this to the large c/a ratio and the high crystalline coherency of the T'-BiFeO3/LaAlO3 interface. This coherency is possible through an observed "diagonal-on-diagonal" alignment. Interestingly, the substrate miscut does not influence the relative population of monoclinic domains.


title: Evidence synthesis for count distributions based on heterogeneous and incomplete aggregated data

abstract: The analysis of count data is commonly done using Poisson models. Negative binomial models are a straightforward and readily motivated generalization for the case of overdispersed data, i.e., when the observed variance is greater than expected under a Poissonian model. Rate and overdispersion parameters then need to be considered jointly, which in general is not trivial. Here we are concerned with evidence synthesis in the case where the reporting of data is rather heterogeneous, i.e., events are reported either in terms of mean event counts, the proportion of event-free patients, or rate estimates and standard errors. Either figure carries some information about the relevant parameters, and it is the joint modeling that allows for coherent inference on the parameters of interest. The methods are motivated and illustrated by a systematic review in chronic obstructive pulmonary disease.


title: Terrestrial Planet Formation in a protoplanetary disk with a local mass depletion: A successful scenario for the formation of Mars

abstract: Models of terrestrial planet formation for our solar system have been successful in producing planets with masses and orbits similar to those of Venus and Earth. However, these models have generally failed to produce Mars-sized objects around 1.5 AU. The body that is usually formed around Mars' semimajor axis is, in general, much more massive than Mars. Only when Jupiter and Saturn are assumed to have initially very eccentric orbits (e $\sim$ 0.1), which seems fairly unlikely for the solar system, or alternately, if the protoplanetary disk is truncated at 1.0 AU, simulations have been able to produce Mars-like bodies in the correct location. In this paper, we examine an alternative scenario for the formation of Mars in which a local depletion in the density of the protosolar nebula results in a non-uniform formation of planetary embryos and ultimately the formation of Mars-sized planets around 1.5 AU. We have carried out extensive numerical simulations of the formation of terrestrial planets in such a disk for different scales of the local density depletion, and for different orbital configurations of the giant planets. Our simulations point to the possibility of the formation of Mars-sized bodies around 1.5 AU, specifically when the scale of the disk local mass-depletion is moderately high (50-75%) and Jupiter and Saturn are initially in their current orbits. In these systems, Mars-analogs are formed from the protoplanetary materials that originate in the regions of disk interior or exterior to the local mass-depletion. Results also indicate that Earth-sized planets can form around 1 AU with a substantial amount of water accreted via primitive water-rich planetesimals and planetary embryos. We present the results of our study and discuss their implications for the formation of terrestrial planets in our solar system.


title: Radiative effects on the thermoelectric problems

abstract: There are two main directions in this paper. One is to find sufficient conditions to ensure the existence of weak solutions to thermoelectric problems. At the steady-state, these problems consist by a coupled system of elliptic equations of the divergence form, commonly accomplished with nonlinear radiation-type conditions on at least on a nonempty part of the boundary of a $C^1$ domain. The model under study takes the thermoelectric Peltier and Seebeck effects into account, whose describe the Joule-Thomson effect. The proof method makes recourse of a fixed point argument. To this end, well-determined estimates are our main concern. The paper is in the second direction for the derivation of explicit $W^{1,p}$-estimates $(p>2)$ for solutions of nonlinear radiation-type problems, where the leading coefficient is assumed to be a discontinuous function on the space variable. In particular, the behavior of the leading coefficient is conveniently explicit on the estimate of any solution. This regularity result is sufficiently general to contribute to other problems, in which the dependence on the values of the involved constants is essential, instead of the problem under study only.


title: Fundamental Limits of Caching with Secure Delivery

abstract: Caching is emerging as a vital tool for alleviating the severe capacity crunch in modern content-centric wireless networks. The main idea behind caching is to store parts of popular content in end-users' memory and leverage the locally stored content to reduce peak data rates. By jointly designing content placement and delivery mechanisms, recent works have shown order-wise reduction in transmission rates in contrast to traditional methods. In this work, we consider the secure caching problem with the additional goal of minimizing information leakage to an external wiretapper. The fundamental cache memory vs. transmission rate trade-off for the secure caching problem is characterized. Rather surprisingly, these results show that security can be introduced at a negligible cost, particularly for large number of files and users. It is also shown that the rate achieved by the proposed caching scheme with secure delivery is within a constant multiplicative factor from the information-theoretic optimal rate for almost all parameter values of practical interest.


title: A class of transient acceleration models consistent with Big Bang Cosmology

abstract: Is it possible that the current cosmic accelerating expansion will turn into a decelerating one? Can this transition be realized by some viable theoretical model that is consistent with the standard Big Bang cosmology? We study a class of phenomenological models of a transient acceleration, based on a dynamical dark energy with a very general form of equation of state $p_{de}=\alpha\rho_{de}-\beta\rho_{de}^m$. It mimics the cosmological constant $\rho_{de}\rightarrow$ const for small scale factor $a$, and behaves as a barotropic gas with $\rho_{de}\rightarrow a^{-3(\alpha+1)}$ with $\alpha\ge 0$ for large $a$. The cosmic evolution of four models in the class has been examined in details, and all yields a smooth transient acceleration. Depending on the specific model, the future universe may be dominated either by the dark energy or by the matter. In two models, the dynamical dark energy can be explicitly realized by a scalar field with an analytical potential $V(\phi)$. Moreover, the statistical analysis shows that the models can be as robust as $\Lambda$CDM in confronting the observational data of SN Ia, CMB, and BAO. As improvements over the previous studies, our models overcome the over-abundance problem of dark energy during early eras, and satisfy the constraints on the dark energy from WMAP observations of CMB.


title: Optical selection of quasars: SDSS and LSST

abstract: Over the last decade, quasar sample sizes have increased from several thousand to several hundred thousand, thanks mostly to SDSS imaging and spectroscopic surveys. LSST, the next-generation optical imaging survey, will provide hundreds of detections per object for a sample of more than ten million quasars with redshifts of up to about seven. We briefly review optical quasar selection techniques, with emphasis on methods based on colors, variability properties and astrometric behavior.


title: The Marangoni flow of soluble amphiphiles

abstract: Surfactant distribution heterogeneities at a fluid/fluid interface trigger the Marangoni effect, i.e. a bulk flow due to a surface tension gradient. The influence of surfactant solubility in the bulk on these flows remains incompletely characterized. Here we study Marangoni flows sustained by injection of hydrosoluble surfactants at the air/water interface. We show that the flow extent increases with a decrease of the critical micelle concentration, i.e. the concentration at which these surfactants self-assemble in water. We document the universality of the surface velocity field and predict scaling laws based on hydrodynamics and surfactant physicochemistry that capture the flow features.


title: Comparison of quenched and annealed invariance principles for random conductance model: Part II

abstract: We show that there exists an ergodic conductance environment such that the weak (annealed) invariance principle holds for the corresponding continuous time random walk but the quenched invariance principle does not hold. In the present paper we give a proof of the full scaling limit for the weak invariance principle, improving the result in an earlier paper where we obtained a subsequential limit.


title: Optical variability of quasars: a damped random walk

abstract: A damped random walk is a stochastic process, defined by an exponential covariance matrix that behaves as a random walk for short time scales and asymptotically achieves a finite variability amplitude at long time scales. Over the last few years, it has been demonstrated, mostly but not exclusively using SDSS data, that a damped random walk model provides a satisfactory statistical description of observed quasar variability in the optical wavelength range, for rest-frame timescales from 5 days to 2000 days. The best-fit characteristic timescale and asymptotic variability amplitude scale with the luminosity, black hole mass, and rest wavelength, and appear independent of redshift. In addition to providing insights into the physics of quasar variability, the best-fit model parameters can be used to efficiently separate quasars from stars in imaging surveys with adequate long-term multi-epoch data, such as expected from LSST.


title: Identifying high-redshift GRBs with RATIR

abstract: We present a template fitting algorithm for determining photometric redshifts, $z_{\rm phot}$, of candidate high-redshift gamma-ray bursts (GRBs). Using afterglow photometry, obtained by the Reionization And Transients InfraRed (RATIR) camera, this algorithm accounts for the intrinsic GRB afterglow spectral energy distribution (SED), host dust extinction and the effect of neutral hydrogen (local and cosmological) along the line of sight. We present the results obtained by this algorithm and RATIR photometry of GRB 130606A, finding a range of best fit solutions $5.6 < z_{\rm phot} < 6.0$ for models of several host dust extinction laws (none, MW, LMC and SMC), consistent with spectroscopic measurements of the redshift of this GRB. Using simulated RATIR photometry, we find our algorithm provides precise measures of $z_{\rm phot}$ in the ranges $4 < z_{\rm phot} \lesssim 8$ and $9 < z_{\rm phot} < 10$ and can robustly determine when $z_{\rm phot}>4$. Further testing highlights the required caution in cases of highly dust extincted host galaxies. These tests also show that our algorithm does not erroneously find $z_{\rm phot} < 4$ when $z_{\rm sim}>4$, thereby minimizing false negatives and allowing us to rapidly identify all potential high-redshift events.


title: Generalized Approximate Message Passing for Cosparse Analysis Compressive Sensing

abstract: In cosparse analysis compressive sensing (CS), one seeks to estimate a non-sparse signal vector from noisy sub-Nyquist linear measurements by exploiting the knowledge that a given linear transform of the signal is cosparse, i.e., has sufficiently many zeros. We propose a novel approach to cosparse analysis CS based on the generalized approximate message passing (GAMP) algorithm. Unlike other AMP-based approaches to this problem, ours works with a wide range of analysis operators and regularizers. In addition, we propose a novel $\ell_0$-like soft-thresholder based on MMSE denoising for a spike-and-slab distribution with an infinite-variance slab. Numerical demonstrations on synthetic and practical datasets demonstrate advantages over existing AMP-based, greedy, and reweighted-$\ell_1$ approaches.


title: Background Electromagnetic Fields and NRQED Matching: Scalar Case

abstract: The low-energy structure of hadrons can be described systematically using effective field theory, and the parameters of the effective theory can be determined from lattice QCD computations. Recent work, however, points to inconsistencies between the background field method in lattice QCD and effective field theory matching conditions. We show that the background field problem necessitates inclusion of operators related by equations of motion. In the presence of time-dependent electromagnetic fields, for example, such operators modify Green's functions, thereby complicating the isolation of hadronic parameters which enter on-shell scattering amplitudes. The particularly simple case of a scalar hadron coupled to uniform electromagnetic fields is investigated in detail. At the level of the relativistic effective theory, operators related by equations of motion are demonstrated to be innocuous. The same result does not hold in the non-relativistic effective theory, and inconsistencies in matching are resolved by carefully treating operators related by equations of motion. As uniform external fields potentially allow for surface terms, the problem is additionally analyzed on a torus where such terms are absent. Finite-size corrections are derived for charged scalar correlation functions in uniform electric fields as a useful byproduct.


title: An Extensive Evaluation of Filtering Misclassified Instances in Supervised Classification Tasks

abstract: Removing or filtering outliers and mislabeled instances prior to training a learning algorithm has been shown to increase classification accuracy. A popular approach for handling outliers and mislabeled instances is to remove any instance that is misclassified by a learning algorithm. However, an examination of which learning algorithms to use for filtering as well as their effects on multiple learning algorithms over a large set of data sets has not been done. Previous work has generally been limited due to the large computational requirements to run such an experiment, and, thus, the examination has generally been limited to learning algorithms that are computationally inexpensive and using a small number of data sets. In this paper, we examine 9 learning algorithms as filtering algorithms as well as examining the effects of filtering in the 9 chosen learning algorithms on a set of 54 data sets. In addition to using each learning algorithm individually as a filter, we also use the set of learning algorithms as an ensemble filter and use an adaptive algorithm that selects a subset of the learning algorithms for filtering for a specific task and learning algorithm. We find that for most cases, using an ensemble of learning algorithms for filtering produces the greatest increase in classification accuracy. We also compare filtering with a majority voting ensemble. The voting ensemble significantly outperforms filtering unless there are high amounts of noise present in the data set. Additionally, we find that a majority voting ensemble is robust to noise as filtering with a voting ensemble does not increase the classification accuracy of the voting ensemble.


title: Balancing bike sharing systems (BBSS): instance generation from the CitiBike NYC data

abstract: Bike sharing systems are a very popular means to provide bikes to citizens in a simple and cheap way. The idea is to install bike stations at various points in the city, from which a registered user can easily loan a bike by removing it from a specialized rack. After the ride, the user may return the bike at any station (if there is a free rack). Services of this kind are mainly public or semi-public, often aimed at increasing the attractiveness of non-motorized means of transportation, and are usually free, or almost free, of charge for the users. Depending on their location, bike stations have specific patterns regarding when they are empty or full. For instance, in cities where most jobs are located near the city centre, the commuters cause certain peaks in the morning: the central bike stations are filled, while the stations in the outskirts are emptied. Furthermore, stations located on top of a hill are more likely to be empty, since users are less keen on cycling uphill to return the bike, and often leave their bike at a more reachable station. These issues result in substantial user dissatisfaction which may eventually cause the users to abandon the service. This is why nowadays most bike sharing system providers take measures to rebalance them. Over the last few years, balancing bike sharing systems (BBSS) has become increasingly studied in optimization. As such, generating meaningful instance to serve as a benchmark for the proposed approaches is an important task. In this technical report we describe the procedure we used to generate BBSS problem instances from data of the CitiBike NYC bike sharing system.


title: Disentangling radiative corrections using high-mass Drell-Yan at the LHC

abstract: We present a detailed numerical study of lepton-pair production via the Drell-Yan process above the Z-peak at the LHC. Our results consistently combine next-to-next-to-leading order QCD corrections and next-to-leading order electroweak effects, and include the leading photon-initiated processes using a recent extraction of the photon distribution function. We focus on the effects of electroweak corrections and of photon-photon scattering contributions, and demonstrate which kinematic distributions exhibit sensitivity to these corrections. We show that a combination of measurements allows them to be disentangled and separately determined.


title: Automation of The Guiding Center Expansion

abstract: We report on the use of the recently-developed Mathematica package \emph{VEST} (Vector Einstein Summation Tools) to automatically derive the guiding center transformation. Our Mathematica code employs a recursive procedure to derive the transformation order-by-order. This procedure has several novel features. (1) It is designed to allow the user to easily explore the guiding center transformation's numerous non-unique forms or representations. (2) The procedure proceeds entirely in cartesian position and velocity coordinates, thereby producing manifestly gyrogauge invariant results; the commonly-used perpendicular unit vector fields $e_1,e_2$ are never even introduced. (3) It is easy to apply in the derivation of higher-order contributions to the guiding center transformation without fear of human error. Our code therefore stands as a useful tool for exploring subtle issues related to the physics of toroidal momentum conservation in tokamaks.


title: The Hamiltonian Mechanics of Stochastic Acceleration

abstract: We show how to find the physical Langevin equation describing the trajectories of particles undergoing collisionless stochastic acceleration. These stochastic differential equations retain not only one-, but two-particle statistics, and inherit the Hamiltonian nature of the underlying microscopic equations. This opens the door to using stochastic variational integrators to perform simulations of stochastic interactions such as Fermi acceleration. We illustrate the theory by applying it to two example problems.


title: Gyrosymmetry: Global Considerations

abstract: In the guiding center theory, smooth unit vectors perpendicular to the magnetic field are required to define the gyrophase. The question of global existence of these vectors is addressed using a general result from the theory of characteristic classes. It is found that there is, in certain cases, an obstruction to global existence. In these cases, the gyrophase cannot be defined globally. The implications of this fact on the basic structure of the guiding center theory are discussed. In particular it is demonstrated that the guiding center asymptotic expansion of the equations of motion can still be performed in a globally consistent manner when a single global convention for measuring gyrophase is unavailable. The latter fact is demonstrated directly by deriving a new expression for the guiding-center Poincar\'e-Cartan form exhibiting no dependence on the choice of perpendicular unit vectors.


title: A Randomized Proper Orthogonal Decomposition Technique

abstract: In this paper, we consider the problem of model reduction of large scale systems, such as those obtained through the discretization of PDEs. We propose a randomized proper orthogonal decomposition (RPOD) technique to obtain the reduced order models by randomly choosing a subset of the inputs/outputs of the system to construct a suitable small sized Hankel matrix from the full Hankel matrix. It is shown that the RPOD technique is computationally orders of magnitude cheaper when compared to techniques such as the Eigensystem Realization algorithm (ERA)/Balanced POD (BPOD) while obtaining the same information in terms of the number and accuracy of the dominant modes. The method is tested on several different advection-diffusion equations.


title: Development of a proton Computed Tomography (pCT) scanner at NIU

abstract: We describe the development of a proton Computed Tomography (pCT) scanner at Northern Illinois University (NIU) in collaboration with Fermilab and Delhi University. This paper provides an overview of major components of the scanner and a detailed description of the data acquisition system (DAQ).


title: Application of chiral resonance Lagrangian theories to the muon g-2

abstract: We think that phenomenological resonance Lagrangian models, constrained by global fits from low energy hadron reaction data, can help to improve muon g-2 predictions. The main issue are those contributions which cannot be calculated by perturbative means: the hadronic vacuum polarization (HVP) effects and the hadronic light--by--light (HLbL) scattering contribution. I review recent progress in the evaluation of the HVP contribution within the broken Hidden Local Symmetry (HLS) framework, worked out in collaboration with M. Benayoun, P. David and L. DelBuono. Our HLS driven estimate reads a_mu (LO had)=(688.60+\-4.24) 10^{-10} and we find a_mu (theory)= (11 659 177.65 +/- 5.76) 10^{-10}.


title: Accelerating Parametric Probabilistic Verification

abstract: We present a novel method for computing reachability probabilities of parametric discrete-time Markov chains whose transition probabilities are fractions of polynomials over a set of parameters. Our algorithm is based on two key ingredients: a graph decomposition into strongly connected subgraphs combined with a novel factorization strategy for polynomials. Experimental evaluations show that these approaches can lead to a speed-up of up to several orders of magnitude in comparison to existing approaches


title: $\sigma$-Biderivations and $\sigma$-commuting maps of triangular algebras

abstract: Let $\A$ be an algebra and $\sigma$ an automorphism of $\A$. A linear map $d$ of $\A$ is called a $\sigma$-derivation of $\A$ if $d(xy) = d(x)y + \sigma(x)d(y)$, for all $x, y \in \A$. A bilinear map $D: \A \times \A \to \A$ is said to be a $\sigma$-biderivation of $\A$ if it is a $\sigma$-derivation in each component. An additive map $\Theta$ of $\A$ is $\sigma$-commuting if it satisfies $\Theta(x)x - \sigma(x)\Theta(x) = 0$, for all $x \in \A$. In this paper, we introduce the notions of inner and extremal $\sigma$-biderivations and of proper $\sigma$-commuting maps. One of our main results states that, under certain assumptions, every $\sigma$-biderivation of a triangular algebras is the sum of an extremal $\sigma$-biderivation and an inner $\sigma$-biderivation. Sufficient conditions are provided on a triangular algebra for all of its $\sigma$-biderivations (respectively, $\sigma$-commuting maps) to be inner (respectively, proper). A precise description of $\sigma$-commuting maps of triangular algebras is also given. A new class of automorphisms of triangular algebras is introduced and precisely described. We provide many classes of triangular algebras whose automorphisms can be precisely described.


title: Joint multi-mode dispersion extraction in Fourier and space time domains

abstract: In this paper we present a novel broadband approach for the extraction of dispersion curves of multiple time frequency overlapped dispersive modes such as in borehole acoustic data. The new approach works jointly in the Fourier and space time domains and, in contrast to existing space time approaches that mainly work for time frequency separated signals, efficiently handles multiple signals with significant time frequency overlap. The proposed method begins by exploiting the slowness (phase and group) and time location estimates based on frequency-wavenumber (f-k) domain sparsity penalized broadband dispersion extraction method as presented in \cite{AeronTSP2011}. In this context we first present a Cramer Rao Bound (CRB) analysis for slowness estimation in the (f-k) domain and show that for the f-k domain broadband processing, group slowness estimates have more variance than the phase slowness estimates and time location estimates. In order to improve the group slowness estimates we exploit the time compactness property of the modes to effectively represent the data as a linear superposition of time compact space time propagators parameterized by the phase and group slowness. A linear least squares estimation algorithm in the space time domain is then used to obtain improved group slowness estimates. The performance of the method is demonstrated on real borehole acoustic data sets.


title: Full Counting Statistics of Generic Spin Entangler with Quantum Dot-Ferromagnet detectors

abstract: Entanglement between spatially separated electrons in nanoscale transport is a fundamental property, yet to be demonstrated experimentally. Here we propose and analyse theoretically the transport statistics of a generic spin entangler coupled to a hybrid quantum dot-ferromagnet detector system. We show that the full distribution of charges arriving at the ferromagnetic terminals provide complete information on the spin state of the particles emitted by the entangler. This provides means for spin entanglement detection via electrical current correlations, with optimal measurement strategies depending on the a priori knowledge of ferromagnet polarization and spin-flip rates in the detector dots. The scheme is exemplified by applying it to Andreev and triple dot entanglers.


title: Nodal Areas and Structured Darkness

abstract: Generalized beams of light that are made to turn inside out, creating black nodal areas of total destructive interference are described. As an example we experimentally created an elliptical node from a uniformly illuminated elliptical aperture by use of a lossless phase mask called a q-plate. We demonstrate how a modified phase retrieval algorithm may be used to design phase masks that achieve this transformation for an arbitrary aperture shape when analytical methods are not available. This generic wave phenomenon may find uses in both optical and non-optical systems.


title: The Type IIb Supernova 2013df and Its Cool Supergiant Progenitor

abstract: We have obtained early-time photometry and spectroscopy of Supernova (SN) 2013df in NGC 4414. The SN is clearly of Type IIb, with notable similarities to SN 1993J. From its luminosity at secondary maximum light, it appears that less $^{56}$Ni ($\lesssim 0.06\ M_{\odot}$) was synthesized in the SN 2013df explosion than was the case for the SNe IIb 1993J, 2008ax, and 2011dh. Based on a comparison of the light curves, the SN 2013df progenitor must have been more extended in radius prior to explosion than the progenitor of SN 1993J. The total extinction for SN 2013df is estimated to be $A_V=0.30$ mag. The metallicity at the SN location is likely to be solar. We have conducted Hubble Space Telescope (HST) Target of Opportunity observations of the SN with the Wide Field Camera 3, and from a precise comparison of these new observations to archival HST observations of the host galaxy obtained 14 years prior to explosion, we have identified the progenitor of SN 2013df to be a yellow supergiant, somewhat hotter than a red supergiant progenitor for a normal Type II-Plateau SN. From its observed spectral energy distribution, assuming that the light is dominated by one star, the progenitor had effective temperature $T_{\rm eff} = 4250 \pm 100$ K and a bolometric luminosity $L_{\rm bol}=10^{4.94 \pm 0.06}\ L_{\odot}$. This leads to an effective radius $R_{\rm eff} = 545 \pm 65\ R_{\odot}$. The star likely had an initial mass in the range of 13 to 17 $M_{\odot}$; however, if it was a member of an interacting binary system, detailed modeling of the system is required to estimate this mass more accurately. The progenitor star of SN 2013df appears to have been relatively similar to the progenitor of SN 1993J.


title: Reexamining the half-lives of 195Os and 195Ir

abstract: Currently the half-life of 195Os is listed as unknown in most databases because the value of the only available measurement had been reassigned. We argue that the original assignment is correct and re-evaluate the half-life of 195Os to be 6.5(11)min, consistent with the original measurement. We also suggest to reassign the half-life of 195Ir to 2.29(17)h.


title: Correlations between user voting data, budget, and box office for films in the Internet Movie Database

abstract: The Internet Movie Database (IMDb) is one of the most-visited websites in the world and the premier source for information on films. Like Wikipedia, much of IMDb's information is user contributed. IMDb also allows users to voice their opinion on the quality of films through voting. We investigate whether there is a connection between this user voting data and certain economic film characteristics. To this end, we perform distribution and correlation analysis on a set of films chosen to mitigate effects of bias due to the language and country of origin of films. We show that production budget, box office gross, and total number of user votes for films are consistent with double-log normal distributions for certain time periods. Both total gross and user votes are consistent with a double-log normal distribution from the late 1980s onward, while for budget, it extends from 1935 to 1979. In addition, we find a strong correlation between number of user votes and the economic statistics, particularly budget. Remarkably, we find no evidence for a correlation between number of votes and average user rating. As previous studies have found a strong correlation between production budget and marketing expenses, our results suggest that total user votes is an indicator of a film's prominence or notability, which can be quantified by its promotional costs.


title: Classical Light Beams and Geometric Phases

abstract: We present a study of geometric phases in classical wave and polarisation optics using the basic mathematical framework of quantum mechanics. Important physical situations taken from scalar wave optics, pure polarisation optics, and the behaviour of polarisation in the eikonal or ray limit of Maxwell's equations in a transparent medium are considered. The case of a beam of light whose propagation direction and polarisation state are both subject to change is dealt with, attention being paid to the validity of Maxwell's equations at all stages. Global topological aspects of the space of all propagation directions are discussed using elementary group theoretical ideas, and the effects on geometric phases are elucidated.


title: Disentangling degenerate solutions from primary transit and secondary eclipse spectroscopy of exoplanets

abstract: Infrared transmission and emission spectroscopy of exoplanets, recorded from primary transit and secondary eclipse measurements, indicate the presence of the most abundant carbon and oxygen molecular species (H2O, CH4, CO, and CO2) in a few exoplanets. However, efforts to constrain the molecular abundances to within several orders of magnitude are thwarted by the broad range of degenerate solutions that fit the data. Here we explore, with radiative transfer models and analytical approximations, the nature of the degenerate solution sets resulting from the sparse measurements of "hot Jupiter" exoplanets. As demonstrated with simple analytical expressions, primary transit measurements probe roughly 4 atmospheric scale heights at each wavelength band. Derived mixing ratios from these data are highly sensitive to errors in the radius in planet (at a reference pressure), which are approximately a few percent. For example, an uncertainty of 1% in the radius of a 1000 K and H2-based exoplanet with Jupiter's radius and mass, causes an uncertainty of a factor of ~100-10000 in the derived gas mixing ratios. The degree of sensitivity depends on how the line strength increases with the optical depth (i.e. the curve of growth) and the atmospheric scale height. Temperature degeneracies in the solutions of the primary transit data, which manifest their effects through scale height and absorption coefficients, are smaller. We argue that these challenges can be partially surmounted by a combination of selected wavelength sampling of optical and IR measurements and, when possible, the joint analysis of transit and secondary eclipse data of exoplanets.


title: Classifiers With a Reject Option for Early Time-Series Classification

abstract: Early classification of time-series data in a dynamic environment is a challenging problem of great importance in signal processing. This paper proposes a classifier architecture with a reject option capable of online decision making without the need to wait for the entire time series signal to be present. The main idea is to classify an odor/gas signal with an acceptable accuracy as early as possible. Instead of using posterior probability of a classifier, the proposed method uses the "agreement" of an ensemble to decide whether to accept or reject the candidate label. The introduced algorithm is applied to the bio-chemistry problem of odor classification to build a novel Electronic-Nose called Forefront-Nose. Experimental results on wind tunnel test-bed facility confirms the robustness of the forefront-nose compared to the standard classifiers from both earliness and recognition perspectives.


title: ECOC-Based Training of Neural Networks for Face Recognition

abstract: Error Correcting Output Codes, ECOC, is an output representation method capable of discovering some of the errors produced in classification tasks. This paper describes the application of ECOC to the training of feed forward neural networks, FFNN, for improving the overall accuracy of classification systems. Indeed, to improve the generalization of FFNN classifiers, this paper proposes an ECOC-Based training method for Neural Networks that use ECOC as the output representation, and adopts the traditional Back-Propagation algorithm, BP, to adjust weights of the network. Experimental results for face recognition problem on Yale database demonstrate the effectiveness of our method. With a rejection scheme defined by a simple robustness rate, high reliability is achieved in this application.


title: Biological mechanism and identifiability of a class of stationary conductance model for Voltage-gated Ion channels

abstract: The physiology of voltage gated ion channels is complex and insights into their gating mechanism is incomplete. Their function is best represented by Markov models with relatively large number of distinct states that are connected by thermodynamically feasible transitions. On the other hand, popular models such as the one of Hodgkin and Huxley have empirical assumptions that are generally unrealistic. Experimental protocols often dictate the number of states in proposed Markov models, thus creating disagreements between various observations on the same channel. Here we aim to propose a limit to the minimum number of states required to model ion channels by employing a paradigm to define stationary conductance in a class of ion-channels. A simple expression is generated using concepts in elementary thermodynamics applied to protein conformational transitions. Further, it matches well many published channel current-voltage characteristics and parameters of the model are found to be identifiable and easily determined from usual experimental protocols.


title: Strict self-adjointness and shallow water models

abstract: We consider a class of third order equations from the point of view of strict self-adjointness. Necessary and sufficient conditions to the investigated class be strictly self-adjoint are obtained. Then, from a strictly self-adjoint subclass we consider those who admits a suitable scaling transformation. Consequently it is derived a family of equations including the Benjamin-Bona-Mahony, Camassa-Holm and Novikov equation. By a suitable choice of the parameters, we deduce an one-parameter family of equations unifying the last two mentioned equations. Then, using some recent techniques for constructing conserved vectors, we show that from the scale invariance it is obtained as a conserved density the same quantity employed to construct one of the well known Hamiltonians for the cited integrable equations.


title: Identities of symmetry for (h;q)-extension of higher-order Euler polynomials

abstract: In this paper, we study some symmetric properties of the multiple q-Euler zeta function. From these properties, we derive several identities of symmetry for the (h;q)-extension of higher-order Euler polynomials.


title: Plasmon resonance and heat generation model in nanostructures

abstract: In this paper, we investigate the photothermal effects of the plasmon resonance. Metal nanoparticles efficiently generate heat in the presence of electromagnetic radiation. The process is strongly enhanced when a fixed frequency of the incident wave illuminate on nanoprticles such that plasmon resonance happen. We shall introduce the electromagnetic radiation model and show exactly how and when the plasmon resonance happen. We then construct the heat generation and transfer model and derive the heat effect induced by plasmon resonance. Finally, we consider the heat generation under plasmon resonance in a concentric nanoshell structure.


title: Non-Hermitian quantum dynamics and entanglement of coupled nonlinear resonators

abstract: We consider a generalization of recently proposed non-Hermitian model for resonant cavities coupled by a chiral mirror by taking into account number non-conservation and nonlinear interactions. We analyze non-Hermitian quantum dynamics of populations and entanglement of the cavity modes. We find that an interplay of initial coherence and non-Hermitian coupling leads to a counterintuitive population transfer. While an initially coherent cavity mode is depleted, the other empty cavity can be populated more or less than the initially filled one. Moreover, presence of nonlinearity yields population collapse and revival as well as bipartite entanglement of the cavity modes. In addition to coupled cavities, we point out that similar models can be found in $\mathcal{PT}$ symmetric Bose-Hubbard dimers of Bose-Einstein condensates or in coupled soliton-plasmon waveguides. We specifically illustrate quantum dynamics of populations and entanglement in a heuristic model that we propose for a soliton-plasmon system with soliton amplitude dependent asymmetric interaction. Degree of asymmetry, nonlinearity and coherence are examined to control plasmon excitations and soliton-plasmon entanglement. Relations to $\mathcal{PT}$ symmetric lasers and Jahn-Teller systems are pointed out.


title: CANFAR+Skytree: A Cloud Computing and Data Mining System for Astronomy

abstract: At the Canadian Astronomy Data Centre, we have combined our cloud computing system, CANFAR, with the world's most advanced machine learning software, Skytree, to create the world's first cloud computing system for data mining in astronomy. CANFAR provides a generic environment for the storage and processing of large datasets, removing the requirement to set up and maintain a computing system when implementing an extensive undertaking such as a survey pipeline. 500 processor cores and several hundred terabytes of persistent storage are currently available to users. The storage is implemented via the International Virtual Observatory Alliance's VOSpace protocol, and is accessible both interactively, and to all processing jobs. The user interacts with CANFAR by utilizing virtual machines, which appear to them as equivalent to a desktop. Each machine is replicated as desired to perform large-scale parallel processing. Such an arrangement enables the user to immediately install and run the same astronomy code that they already utilize, in the same way as on a desktop. In addition, unlike many cloud systems, batch job scheduling is handled for the user on multiple virtual machines by the Condor job queueing system. Skytree is installed and run just as any other software on the system, and thus acts as a library of command line data mining functions that can be integrated into one's wider analysis. Thus we have created a generic environment for large-scale analysis by data mining, in the same way that CANFAR itself has done for storage and processing. Because Skytree scales to large data in linear runtime, this allows the full sophistication of the huge fields of data mining and machine learning to be applied to the hundreds of millions of objects that make up current large datasets. We demonstrate the utility of the CANFAR+Skytree system by showing science results obtained. [Abridged]


title: Focus Demo: CANFAR+Skytree: A Cloud Computing and Data Mining System for Astronomy

abstract: This is a companion Focus Demonstration article to the CANFAR+Skytree poster (Ball 2012), demonstrating the usage of the Skytree machine learning software on the Canadian Advanced Network for Astronomical Research (CANFAR) cloud computing system. CANFAR+Skytree is the world's first cloud computing system for data mining in astronomy.


title: Exploring the blazar zone in High Energy flares of FSRQs

abstract: The gamma-ray emission offers a powerful diagnostic tool to probe jets and their surroundings in flat spectrum radio quasars (FSRQ). In particular, sources emitting at high energies (>10 GeV) give us the strongest constraints. This motivates us to start a systematic study of flares with bright emission above 10 GeV, examining archival data of Fermi-LAT gamma-ray telescope. At the same time, we began to trigger Target of Opportunity observations to the Swift observatory at the occurrence of high-energy flares, obtaining a wide coverage of the spectral energy distributions for several FSRQs during flares. Among the others we investigate the SED of a peculiar flare of 3C 454.3, showing a remarkable hard gamma-ray spectrum, quite different from the brightest flares of this source, and a bright flare of CTA 102. We modeled the SED in the framework of the one--zone leptonic model, using also archival optical spectroscopic data to derive the luminosity of the broad lines and thus estimate the disk luminosity, from which the structural parameters of the FSRQ nucleus can be inferred. The model allowed us to evaluate the magnetic field intensity in the blazar zone, and to locate the emitting region of gamma rays in the particular case in which gamma-ray spectra show neither absorption from the BLR, nor the Klein-Nishina curvature expected in leptonic models assuming the BLR as source of seed photons for the External Compton. For FSRQs bright above 10 GeV, we where able to identify short periods lasting less than 1 day characterized by high rate of high energy gamma rays, and hard gamma-ray spectra. We discussed the observed spectra and variability timescales in terms of injection and cooling of energetic particles, arguing that these flares could be triggered by magnetic reconnections events or turbulence in the flow.


title: Electron Delocalization Determined Anomalous Stability in Small Water Rings

abstract: Water clusters are known to form through hydrogen bonding. However, this study shows that the formation of very small water clusters significantly deviates from this mechanism and instead involves both hydrogen bonding and electron delocalization. Our density functional theory calculations show that small water rings (H2O)n of n=3 or 4 show strong electron delocalization originating from both the hydrogen and oxygen atoms and extending to the ring center. This is very different from larger rings. Further energy decomposition of rings with n=3-6 demonstrates an upward trend in the polarization component but an decrease in the electrostatic and exchange repulsion components, presenting a minimum and accounting for 33% of interaction energy at n=3. This significantly promotes stability of the small water rings. Our findings provide a comprehensive analysis and improve our understanding of the stability characteristics of water clusters.


title: Average clock times for scattering through asymmetric barriers

abstract: The reflection and transmission Salecker-Wigner-Peres clock times averaged over the post-selected reflected and transmitted sub-ensembles, respectively, are investigated for the one dimensional scattering of a localized wave packet through an asymmetric barrier. The dwell time averaged over the same post-selected sub-ensembles is also considered. The emergence of negative average reflection times is examined and we show that while the average over the reflected sub-ensemble eliminates the negative peaks at resonance for the clock time, it still allows negative values for transparent barriers. The saturation of the average times with the barrier width (Hartman effect) is also addressed.


title: A Generalization of the Marshal-Olkin Scheme and Related Processes

abstract: A generalization of the Marshal-Olkin parametrization scheme is developed and stochastic models related to it are discussed here.


title: A K-theory approach to the tangent invariants of blow-ups

abstract: We extend the formula for the Chern classes of blow-ups of algebraic varieties due to Porteous and Lascu-Scott, and of symplectic and complex manifolds due to Geiges and Pasquotto, to the blow-ups of almost complex manifolds. Our approach is based on a concrete partition for the tangent bundle of a blow-up. The use of topological K-theory of vector bundles simplifies the classical approaches.


title: Asynchronous Physical-Layer Network Coding with Quasi-Cyclic Codes

abstract: Communication in the presence of bounded timing asynchronism which is known to the receiver but cannot be easily compensated is studied. Examples of such situations include point-to-point communication over inter-symbol interference (ISI) channels and asynchronous wireless networks. In these scenarios, although the receiver may know all the delays, it is often not be an easy task for the receiver to compensate the delays as the signals are mixed together. A novel framework called interleave/deinterleave transform (IDT) is proposed to deal with this problem. It is shown that the IDT allows one to design the delays so that quasi-cyclic (QC) codes with a proper shifting constraint can be used accordingly. When used in conjunction with QC codes, IDT provides significantly better performance than existing schemes relying solely on cyclic codes. Two instances of asynchronous physical-layer network coding, namely the integer-forcing equalization for ISI channels and asynchronous compute-and-forward, are then studied. For integer-forcing equalization, the proposed scheme provides improved performance over using cyclic codes. For asynchronous compute-and-forward, the proposed scheme shows that there is no loss in the achievable information due to delays which are integer multiples of the symbol duration. Further, the proposed approach shows that delays introduced by the channel can sometimes be exploited to obtain higher information rates than those obtainable in the synchronous case. The proposed IDT can be thought of as a generalization of the interleaving/deinterleaving idea proposed by Wang et al. which allows the use of QC codes thereby substantially increasing the design space.


title: Production and decays of a light $\phi^0$ in the LRTH model under the LHC Higgs data

abstract: In this paper we study the production and decays of a light pseudoscalar boson $\phi^0$ with $m_{\phi^0} \leq m_h/2$ appeared in the left-right twin Higgs (LRTH) model, and explore its phenomenological consequences when the latest LHC Higgs data are taken into account. We found that (a) the decay rate $Br(h\to \phi^0\phi^0)$ can be as large as $80\%$ and can suppress significantly the visible $\gamma\gamma$ signal rate, but the latest LHC Higgs data put a strong constraint on it: $Br(h \to \phi^0\phi^0) \leq 30\%$ at $3\sigma$ level; (b) the $p$-value of the LRTH model is around $0.6$, smaller than that of the SM in most of the parameter space and approaches the SM value $0.8$ for a sufficiently large $f$ parameter; (c) the neutral pseudoscalar $\phi^0$ dominantly decay into $b\bar{b}$ and the decay rate $Br(\phi^0\to b\bar{b})$ can be larger than $80\%$ for $m_{\phi^0}\leq 60$ GeV, and the second main decay mode is $\phi^0\to \tau^{+}\tau^{-}$ with a branching ratio about $14\%$; and (d) at the future $e^-e^+$ collider with $\sqrt{s}=250$ GeV, the processes $e^{+}e^{-}\to Zh\to Z(\phi^0\phi^0)\to Z(4b,2b2\tau)$ are promising for discovering such a light pseudoscalar $\phi^0$.


title: Tuning Magnetic Coupling in Sr$_2$IrO$_4$ Thin Films with Epitaxial Strain

abstract: We report x-ray resonant magnetic scattering (XRMS) and resonant inelastic x-ray scattering (RIXS) studies of epitaxially-strained $\mathrm{Sr_2IrO_4}$ thin films. The films were grown on $\mathrm{SrTiO_3}$ and $\mathrm{(LaAlO_3)_{0.3}(Sr_2AlTaO_6)_{0.7}}$ substrates, under slight tensile and compressive strains, respectively. Although the films develop a magnetic structure reminiscent of bulk $\mathrm{Sr_2IrO_4}$, the magnetic correlations are extremely anisotropic, with in-plane correlation lengths significantly longer than the out-of-plane correlation lengths. In addition, the compressive (tensile) strain serves to suppress (enhance) the magnetic ordering temperature $\mathrm{T_N}$, while raising (lowering) the energy of the zone boundary magnon. Quantum chemical calculations show that the tuning of magnetic energy scales can be understood in terms of strain-induced change in bond lengths.


title: Three-dimensional Mid-air Acoustic Manipulation by Ultrasonic Phased Arrays

abstract: The essence of levitation technology is the countervailing of gravity. It is known that an ultrasound standing wave is capable of suspending small particles at its sound pressure nodes. The acoustic axis of the ultrasound beam in conventional studies was parallel to the gravitational force, and the levitated objects were manipulated along the fixed axis (i.e. one-dimensionally) by controlling the phases or frequencies of bolted Langevin-type transducers. In the present study, we considered extended acoustic manipulation whereby millimetre-sized particles were levitated and moved three-dimensionally by localised ultrasonic standing waves, which were generated by ultrasonic phased arrays. Our manipulation system has two original features. One is the direction of the ultrasound beam, which is arbitrary because the force acting toward its centre is also utilised. The other is the manipulation principle by which a localised standing wave is generated at an arbitrary position and moved three-dimensionally by opposed and ultrasonic phased arrays. We experimentally confirmed that expanded-polystyrene particles of 0.6 mm and 2 mm in diameter could be manipulated by our proposed method.


title: On Real and Virtual Photons in the Davies Theory of Time-Symmetric Quantum Electrodynamics

abstract: This paper explores the distinction between virtual and real photons in the context of the Davies quantum relativistic extension of the Wheeler-Feynman classical electromagnetic theory. An alternative way of understanding this distinction is proposed, based on the transactional model.


title: Gauge Equivalence and the Inverse Spectral Problem for the Magnetic Schr\"odinger Operator on the Torus

abstract: We study the inverse spectral problem for the Schr\"odinger operator $H$ on the two-dimensional torus with even magnetic field $B(x)$ and even electric potential $V(x)$. V.Guillemin [11] proved that the spectrum of $H$ determines $B(x)$ and $V(x)$. A simple proof of Guillemin's results was given by the authors in [3]. In the present paper we consider gauge equivalent classes of magnetic potentials and give conditions which imply that the gauge equivalence class and the spectrum of $H$ determine the magnetic field and the electric potential. We also show that generically the spectrum and the magnetic field determine the "extended" gauge equivalence class of the magnetic potential. The proof is a modification of the proof in [3] with some corrections and clarifications.


title: Small knot mosaics and partition matrices

abstract: Lomonaco and Kauffman introduced knot mosaic system to give a definition of quantum knot system. This definition is intended to represent an actual physical quantum system. A knot $(m,n)$-mosaic is an $m \times n$ matrix of mosaic tiles which are $T_0$ through $T_{10}$ depicted as below, representing a knot or a link by adjoining properly that is called suitably connected. An interesting question in studying mosaic theory is how many knot $(m,n)$-mosaics are there. $D_{m,n}$ denotes the total number of all knot $(m,n)$-mosaics. This counting is very important because the total number of knot mosaics is indeed the dimension of the Hilbert space of these quantum knot mosaics. In this paper, we find a table of the precise values of $D_{m,n}$ for $4 \leq m \leq n \leq 6$ as below. Mainly we use a partition matrix argument which turns out to be remarkably efficient to count small knot mosaics. \begin{center} \begin{tabular}{|c|r|r|r|} \hline $D_{m,n}$ & $n=4$ & $n=5$ & $n=6$ \\ \hline $m=4$ & $2594$ & $54,226$ & $1,144,526$ \\ \hline $m=5$ & & $4,183,954$ & $331,745,962$ \\ \hline $m=6$ & & & $101,393,411,126$ \\ \hline \end{tabular} \end{center}


title: Higgs Precision Measurements and Flavor Physics: A Supersymmetric Example

abstract: [Review for "Chinese Science Bulletin" and CEPC+SPPC Proposal] Rare decays in flavor physics often suffer from Helicity suppress and Loop suppress. Helicity flip is a direct consequence of chiral $U(3)$ symmetry breaking and electroweak symmetry breaking. The identical feature is also shared by the mass generation of SM fermions. In this review, we use MSSM as an example to illustrate an explicit connection between bottom Yukawa coupling and rare decay process of $b\to s\gamma$. We take a symmetry approach to study the common symmetry breaking in supersymmetric correction to bottom quark mass generation and $b\to s\gamma$. We show that Large Peccei-Quinn symmetry breaking effect and $R$-symmetry breaking effect required by $b\to s\gamma$ inevitably lead to significant reduction of bottom Yukawa ${y}_{b}$. To compromise the reduction in $b\bar{b}$, a new decay is also needed to keep the Higgs total width as the SM value.


title: Spin-Orbit Driven Transitions Between Mott Insulators and Finite Momentum Superfluids of Bosons in Optical Lattices

abstract: Synthetic spin-orbit coupling in ultracold atomic gases can be taken to extremes rarely found in solids. We study a two dimensional Hubbard model of bosons in an optical lattice in the presence of spin-orbit coupling strong enough to drive direct transitions from Mott insulators to superfluids. Here we find phase-modulated superfluids with finite momentum that are generated entirely by spin-orbit coupling. We investigate the rich phase patterns of the superfluids, which may be directly probed using time-of-flight imaging of the spin-dependent momentum distribution.


title: Oblivious Query Processing

abstract: Motivated by cloud security concerns, there is an increasing interest in database systems that can store and support queries over encrypted data. A common architecture for such systems is to use a trusted component such as a cryptographic co-processor for query processing that is used to securely decrypt data and perform computations in plaintext. The trusted component has limited memory, so most of the (input and intermediate) data is kept encrypted in an untrusted storage and moved to the trusted component on ``demand.'' In this setting, even with strong encryption, the data access pattern from untrusted storage has the potential to reveal sensitive information; indeed, all existing systems that use a trusted component for query processing over encrypted data have this vulnerability. In this paper, we undertake the first formal study of secure query processing, where an adversary having full knowledge of the query (text) and observing the query execution learns nothing about the underlying database other than the result size of the query on the database. We introduce a simpler notion, oblivious query processing, and show formally that a query admits secure query processing iff it admits oblivious query processing. We present oblivious query processing algorithms for a rich class of database queries involving selections, joins, grouping and aggregation. For queries not handled by our algorithms, we provide some initial evidence that designing oblivious (and therefore secure) algorithms would be hard via reductions from two simple, well-studied problems that are generally believed to be hard. Our study of oblivious query processing also reveals interesting connections to database join theory.


title: Integrable hydrodynamic chains for WZNW model

abstract: The new integrable hydrodynamic equations obtained for WZNW model with $SU(2)$, $SO(3)$, $SP(2)$ and $SU(\infty )$ constant torsions


title: A Simple Method to Produce Algorithmic MIDI Music based on Randomness, Simple Probabilities and Multi-Threading

abstract: This paper introduces a simple method for producing multichannel MIDI music that is based on randomness and simple probabilities. One distinctive feature of the method is that it produces and sends in parallel to the sound card more than one unsynchronized channels by exploiting the multi-threading capabilities of general purpose programming languages. As consequence the derived sound offers a quite ``full" and ``unpredictable" acoustic experience to the listener. Subsequently the paper reports the results of an evaluation with users. The results were very surprising: the majority of users responded that they could tolerate this music in various occasions.


title: The Garnir Relations for Weyl Groups

abstract: The so-called Garnir relations are very important in the representation theory of the symmetric groups. In this paper, we give a possible extension of Garnir relations for Weyl groups in general.


title: Hybrid Radio-map for Noise Tolerant Wireless Indoor Localization

abstract: In wireless networks, radio-map based locating techniques are commonly used to cope the complex fading feature of radio signal, in which a radio-map is built by calibrating received signal strength (RSS) signatures at training locations in the offline phase. However, in severe hostile environments, such as in ship cabins where severe shadowing, blocking and multi-path fading effects are posed by ubiquitous metallic architecture, even radio-map cannot capture the dynamics of RSS. In this paper, we introduced multiple feature radio-map location method for severely noisy environments. We proposed to add low variance signature into radio map. Since the low variance signatures are generally expensive to obtain, we focus on the scenario when the low variance signatures are sparse. We studied efficient construction of multi-feature radio-map in offline phase, and proposed feasible region narrowing down and particle based algorithm for online tracking. Simulation results show the remarkably performance improvement in terms of positioning accuracy and robustness against RSS noises than the traditional radio-map method.


title: Possible Astrophysical Observables of Quantum Gravity Effects near Black Holes

abstract: Recent implications of results from quantum information theory applied to black holes has led to the confusing conclusions that requires either abandoning the equivalence principle (e.g. the firewall picture), or the no-hair theorem (e.g. the fuzzball picture), or even more unpalatable options. The recent discovery of a pulsar orbiting a black hole opens up new possibilities for tests of theories of gravity. We examine possible observational effects of semiclassical quantum gravity in the vicinity of black holes, as probed by pulsars and event horizon telescope imaging of flares. Pulsar radiation is observable at wavelengths only two orders of magnitude shorter than the Hawking radiation, so precision interferometry of lensed pulsar images may shed light on the quantum gravitational processes and interaction of Hawking radiation with the spacetime near the black hole. This paper discusses the impact on the pulsar radiation interference pattern, which is observable through the modulation index in the foreseeable future, and discusses a possible classical limit of BHC.


title: Bicrossed Products, Matched Pair Deformations and the Factorization Index for Lie Algebras

abstract: For a perfect Lie algebra $\mathfrak{h}$ we classify all Lie algebras containing $\mathfrak{h}$ as a subalgebra of codimension $1$. The automorphism groups of such Lie algebras are fully determined as subgroups of the semidirect product $\mathfrak{h} \ltimes (k^* \times {\rm Aut}_{\rm Lie} (\mathfrak{h}))$. In the non-perfect case the classification of these Lie algebras is a difficult task. Let $\mathfrak{l} (2n+1, k)$ be the Lie algebra with the bracket $[E_i, G] = E_i$, $[G, F_i] = F_i$, for all $i = 1, \dots, n$. We explicitly describe all Lie algebras containing $\mathfrak{l} (2n+1, k)$ as a subalgebra of codimension $1$ by computing all possible bicrossed products $k \bowtie \mathfrak{l} (2n+1, k)$. They are parameterized by a set of matrices ${\rm M}_n (k)^4 \times k^{2n+2}$ which are explicitly determined. Several matched pair deformations of $\mathfrak{l} (2n+1, k)$ are described in order to compute the factorization index of some extensions of the type $k \subset k \bowtie \mathfrak{l} (2n+1, k)$. We provide an example of such extension having an infinite factorization index.


title: Study of the apsidal precession of the Physical Symmetrical Pendulum

abstract: We study the apsidal precession of a Physical Symmetrical Pendulum (Allais' precession) as a generalization of the precession corresponding to the Ideal Spherical Pendulum (Airy's Precession). Based on the Hamilton-Jacobi formalism and using the technics of variation of parameters along with the averaging method, we obtain approximate solutions, in terms of which the motion of both systems admits a simple geometrical description. The method developed in this paper is considerably simpler than the standard one in terms of elliptical functions and the numerical agreement with the exact solutions is excellent. In addition, the present procedure permits to show clearly the origin of the Airy's and Allais' precession, as well as the effect of the spin of the Physical Pendulum on the Allais' precession. Further, the method can be extended to the study of the asymmetrical pendulum in which an exact solution is not possible anymore.


title: Trkalian fields: ray transforms and mini-twistors

abstract: We study X-ray and Divergent beam transforms of Trkalian fields and their relation with Radon transform. We make use of four basic mathematical methods of tomography due to Grangeat, Smith, Tuy and Gelfand-Goncharov for an integral geometric view on them. We also make use of direct approaches which provide a faster but restricted view of the geometry of these transforms. These reduce to well known geometric integral transforms on a sphere of the Radon or the spherical Curl transform in Moses eigenbasis, which are members of an analytic family of integral operators. We also discuss their inversion. The X-ray (also Divergent beam) transform of a Trkalian field is Trkalian. Also the Trkalian subclass of X-ray transforms yields Trkalian fields in the physical space. The Riesz potential of a Trkalian field is proportional to the field. Hence, the spherical mean of the X-ray (also Divergent beam) transform of a Trkalian field over all lines passing through a point yields the field at this point. The pivotal point is the simplification of an intricate quantity: Hilbert transform of the derivative of Radon transform for a Trkalian field in the Moses basis. We also define the X-ray transform of the Riesz potential (of order 2) and Biot-Savart integrals. Then, we discuss a mini-twistor respresentation, presenting a mini-twistor solution for the Trkalian fields equation. This is based on a time-harmonic reduction of wave equation to Helmholtz equation. A Trkalian field is given in terms of a null vector in C3 with an arbitrary function and an exponential factor resulting from this reduction.


title: New approach to finding the maximum number of mutually unbiased bases in $\mathbb{C}^6$

abstract: There has been great interest in finding sets of $m$ mutually unbiased bases which are compatible with a given space $\mathbb{C}^d$, specially in physics due to their interesting applications in quantum information theory. Several general results have been obtained so far, but surprising results may occur for definite $(m,d)$-values. One such case that has remained an open question (the simplest case) is the one regarding the existence of $m=4$ mutually orthogonal bases for $d=6$. In the present work we introduce a new approach to the problem by translating it into an optimization procedure for a given pair $(m,d)$.


title: Extended Armendariz Rings

abstract: In this note we introduce central linear Armendariz rings as a generalization of Armendariz rings and investigate their properties.


title: On A Class of Lifting Modules

abstract: In this paper, we introduce principally $\delta$-lifting modules which are analogous to $\delta$-lifting modules and principally $\delta$-semiperfect modules as a generalization of $\delta$-semiperfect modules and investigate their properties.


title: Rings in which every nilpotent is central

abstract: In this paper, we introduce a class of rings in which every nilpotent element is central. This class of rings generalizes so-called reduced rings. A ring $R$ is called {\it central reduced} if every nilpotent element of $R$ is central. For a ring $R$, we prove that $R$ is central reduced if and only if $R[x_1,x_2,\ldots,x_n]$ is central reduced if and only if $R[[x_1,x_2,\ldots,x_n]]$ is central reduced if and only if $R[x_1,x_1^{-1},x_2,x_2^{-1},\ldots,x_n,x_n^{-1}]$ is central reduced. Moreover, if $R$ is a central reduced ring, then the trivial extension $T(R,R)$ is central Armendariz.


title: Functional renormalization group study of superconductivity in doped Sr$_2$IrO$_4$

abstract: Using functional renormalization group we investigated possible superconductivity in doped Sr$_2$IrO$_4$. In the electron doped case, a $d^*_{x^2-y^2}$-wave superconducting phase is found in a narrow doping region. The pairing is driven by spin fluctuations within the single conduction band. In contrast, for hole doping an $s^*_{\pm}$-wave phase is established, triggered by spin fluctuations within and across the two conduction bands. In all cases there are comparable singlet and triplet components in the pairing function. The Hund's rule coupling reduces (enhances) superconductivity for electron (hole) doping. Our results imply that hole doping is more promising to achieve a higher transition temperature. Experimental perspectives are discussed.


title: Achieving Fully Proportional Representation: Approximability Results

abstract: We study the complexity of (approximate) winner determination under the Monroe and Chamberlin--Courant multiwinner voting rules, which determine the set of representatives by optimizing the total (dis)satisfaction of the voters with their representatives. The total (dis)satisfaction is calculated either as the sum of individual (dis)satisfactions (the utilitarian case) or as the (dis)satisfaction of the worst off voter (the egalitarian case). We provide good approximation algorithms for the satisfaction-based utilitarian versions of the Monroe and Chamberlin--Courant rules, and inapproximability results for the dissatisfaction-based utilitarian versions of them and also for all egalitarian cases. Our algorithms are applicable and particularly appealing when voters submit truncated ballots. We provide experimental evaluation of the algorithms both on real-life preference-aggregation data and on synthetic data. These experiments show that our simple and fast algorithms can in many cases find near-perfect solutions.


title: Application of Computational Physics: Blood Vessel Constrictions and Medical Infuses

abstract: Application of computation in many fields are growing fast in last two decades. Increasing on computation performance helps researchers to understand natural phenomena in many fields of science and technology including in life sciences. Computational fluid dynamic is one of numerical methods which is very popular used to describe those phenomena. In this paper we propose moving particle semi-implicit (MPS) and molecular dynamics (MD) to describe different phenomena in blood vessel. The effect of increasing the blood pressure on vessel wall will be calculate using MD methods, while the two fluid blending dynamics will be discussed using MPS. Result from the first phenomenon shows that around 80% of constriction on blood vessel make blood vessel increase and will start to leak on vessel wall, while from the second phenomenon the result shows the visualization of two fluids mixture (drugs and blood) influenced by ratio of drugs debit to blood debit. Keywords: molecular dynamic, blood vessel, fluid dynamic, moving particle semi implicit.

categories: physics.bio-ph cond-mat.soft physics.med-ph q-bio.TO


title: Isomorphism classes and invariants of low dimensional filiform Leibniz algebras

abstract: In this study, we extend the result on classification of a subclass of filiform Leibniz algebras in low dimensions to dimensions seven and eight based on the technique used by Rakhimov and Bekbaev for classification of subclasses which arise from naturally graded non-Lie filiform Leibniz algebras. The subclass considered here arises from naturally graded filiform Lie algebras. This subclass contains the class of filiform Lie algebras and consequently, by classifying this subclass, we again re-examine the classification of filiform Lie algebras. Our resulting list of filifom Lie algebras is compared with that given by Ancoch\'ea-Bermudez and Goze in 1988 and by G\'{o}mez, Jimenez-Merchan and Khakimdjanov in 1998.


title: Molecular Dynamics Simulation of Soft Grains: Malaria-Infected Red Blood Cells Motion within Obstructed 2-D Capillary Vessel

abstract: Molecular dynamics has been widely used to numerically solve equation of motion of classical many-particle system. It can be used to simulate many systems including biophysics, whose complexity level is determined by the involved elements. Based on this method, a numerical model had been constructed to mimic the behaviour of malaria-infected red blood cells within capillary vessel. The model was governed by three forces namely Coulomb force, normal force, and Stokes force. By utilizing two dimensional four-cells scheme, theoretical observation was carried out to test its capability. Although the parameters were chosen deliberately, all of the quantities were given arbitrary value. Despite this fact, the results were quite satisfactory. Combined with the previous results, it can be said that the proposed model were sufficient enough to mimic the malaria-infected red blood cells motion within obstructed capillary vessel. Keywords: molecular dynamics, two-dimensional model, red-blood cell motion, malaria


title: Polynomial Hamiltonian Systems with Movable Algebraic Singularities

abstract: The singularity structure of solutions of a class of Hamiltonian systems of ordinary differential equations in two dependent variables is studied. It is shown that for any solution, all movable singularities, obtained by analytic continuation along a rectifiable curve, are at most algebraic branch points.


title: Study Of Internal Response Of Epoxy Due To Compressive Load Via Experiment And Simulation Using Abaqus FEA Software

abstract: Epoxy is widely used primarily as a matrix material in the manufacture of polymer matrix composite. Epoxy behavior under compression load has to be understood before the mechanical behavior of polymer matrix composite can be accurately predicted. Simulation model combined with experiment and image analysis have been used to investigate internal response of epoxy polymeric materials subjected to compressive loads. In the experiment, small carbon-based material rods are inserted in the epoxy. The samples are held in one side and subjected to compressive load on the other side. All the samples swell at load sides. Image analysis on the carbon-based rods figures out the internal response, which seems to be isotropy in lateral direction. The results are compared to simulation results. The simulation is conducted using Abaqus FEA software. Similar condition is obtained when a brittle thin material is stuck to the top of simulation model. Keywords: Abaqus FEA, compression, epoxy, polymer matrix composites


title: Static bending and free vibration of cross ply laminated composite plates using NURBS based finite element method and unified formulation

abstract: This paper presents an effective formulation to study the response of laminated composites based on isogeometric approach (IGA) and Carrera unified formulation (CUF). The IGA utilizes the non-uniform rational B-spline (NURBS) functions which allows to construct higher order smooth functions with less computational effort. The static bending and the free vibration of thin and moderately thick laminates plates are studied. The present approach also suffers from shear locking when lower order functions are employed and the shear locking is suppressed by introducing a modification factor. The combination of the IGA with the CUF allows a very accurate prediction of the field variables. The effectiveness of the formulation is demonstrated through numerical examples.


title: On the Homogenization of Geological Fissured Systems With Curved non-periodic Cracks

abstract: We analyze the steady fluid flow in a porous medium containing a network of thin fissures i.e. width $\mathcal{O}(\epsilon)$, where all the cracks are generated by the rigid translation of a continuous piecewise $C^{1}$ functions in a fixed direction. The phenomenon is modeled in mixed variational formulation, using the stationary Darcy's law and setting coefficients of low resistance $\mathcal{O}(\epsilon)$ on the network. The singularities are removed performing asymptotic analysis as $\epsilon \rightarrow 0$ which yields an analogous system hosting only tangential flow in the fissures. Finally the fissures are collapsed into two dimensional manifolds.


title: Analysis of a class of degenerate parabolic equations with saturation mechanisms

abstract: We analyze a family of degenerate parabolic equations with linear growth Lagrangian having the form $u_t=\div (\varphi(u)\psi(\nabla u/u))$. Here $|\psi|\le 1$ and saturates at infinity. We present a simple and natural set of assumptions on the functions $\psi,\varphi$, under which: 1) these equations fall in the framework provided by \cite{ACMEllipticFLDE, ACMMRelat} and hence they are well posed, 2) we can ensure finite propagation speed for these models, 3) a Rankine--Hugoniot analysis on traveling fronts is also performed. On the particular case of $\varphi(u)=u$ we get more detailed information on the spreading rate of compactly supported solutions and some interesting connections with optimal mass transportation theory.


title: Simple predictions from multifield inflationary models

abstract: We explore whether multifield inflationary models make unambiguous predictions for fundamental cosmological observables. Focusing on $N$-quadratic inflation, we numerically evaluate the full perturbation equations for models with 2, 3, and $\mathcal{O}(100)$ fields, using several distinct methods for specifying the initial values of the background fields. All scenarios are highly predictive, with the probability distribution functions of the cosmological observables becoming more sharply peaked as $N$ increases. For $N=100$ fields, 95% of our Monte Carlo samples fall in the ranges $n_s \in (0.9455,0.9534)$; $\alpha \in (-9.741,-7.047)\times 10^{-4}$; $r\in(0.1445,0.1449)$; and $r_\mathrm{iso} \in (0.02137,3.510)\times 10^{-3}$ for the spectral index, running, tensor-to-scalar ratio, and isocurvature-to-adiabatic ratio, respectively. The expected amplitude of isocurvature perturbations grows with $N$, raising the possibility that many-field models may be sensitive to post-inflationary physics and suggesting new avenues for testing these scenarios.