title: Se content $x$ dependence of electron correlation strength in Fe$_{1+y}$Te$_{1-x}$Se$_{x}$

abstract: The iron chalcogenide Fe$_{1+y}$Te$_{1-x}$Se$_{x}$ on the Te-rich side is known to exhibit the strongest electron correlations among the Fe-based superconductors, and is non-superconducting for $x$ < 0.1. In order to understand the origin of such behaviors, we have performed ARPES studies of Fe$_{1+y}$Te$_{1-x}$Se$_{x}$ ($x$ = 0, 0.1, 0.2, and 0.4). The obtained mass renormalization factors for different energy bands are qualitatively consistent with DFT + DMFT calculations. Our results provide evidence for strong orbital dependence of mass renormalization, and systematic data which help us to resolve inconsistencies with other experimental data. The unusually strong orbital dependence of mass renormalization in Te-rich Fe$_{1+y}$Te$_{1-x}$Se$_{x}$ arises from the dominant contribution to the Fermi surface of the $d_{xy}$ band, which is the most strongly correlated and may contribute to the suppression of superconductivity.


title: Divisorial Zariski Decomposition and some properties of full mass currents

abstract: Let $\alpha$ be a big class on a compact K\"ahler manifold. We prove that a decomposition $\alpha=\alpha_1+\alpha_2$ into the sum of a modified nef class $\alpha_1$ and a pseudoeffective class $\alpha_2$ is the divisorial Zariski decomposition of $\alpha$ if and only if $\operatorname{vol}(\alpha)=\operatorname{vol}(\alpha_1)$. We deduce from this result some properties of full mass currents.


abstract: Diffusion of Cu adatoms and dimers on Cu(111) and Ag(111) surfaces is analyzed basing on ab-initio surface potentials. As compared to the simple single adatom jumps that alternate between two states, dimers undergo more complex diffusion process that combines translational and rotational moves. Diffusion coefficients for these processes on both surfaces are calculated and compared with diffusion of single particle. The main parameters of diffusion coefficient: energy barrier and prefactor become temperature dependent when dimer moves jumping between many positions of different values of energy. Small difference in the surface lattice constant for Cu and Ag crystal results in completely different energy landscape for dimer jumps and as an effect different diffusion process. At the same time single adatom diffusion does not change so much.


title: Radial abundance gradients from planetary nebulae at different distances from the galactic plane

abstract: We investigate the variations of the radial O/H abundance gradients from planetary nebulae (PN) located at different distances from the galactic plane. In particular, we determine the abundance gradients at different heights from the plane in order to investigate a possible gradient inversion for the objects at larger distances from the plane. We consider a large sample of PN with known distances, so that the height relative to the galactic plane can be derived, and accurate abundances, so that the gradients can be determined.


title: Real-time evolution of a large-scale relativistic jet

abstract: Context. Astrophysical jets are ubiquitous in the Universe on all scales, but their large-scale dynamics and evolution in time are hard to observe since they usually develop at a very slow pace. Aims. We aim to obtain the first observational proof of the expected large-scale evolution and interaction with the environment in an astrophysical jet. Only jets from microquasars offer a chance to witness the real-time, full-jet evolution within a human lifetime, since they combine a 'short', few parsec length with relativistic velocities. Methods. The methodology of this work is based on a systematic recalibraton of interferometric radio observations of microquasars available in public archives. In particular, radio observations of the microquasar GRS 1758-258 over less than two decades have provided the most striking results. Results. Significant morphological variations in the extended jet structure of GRS 1758-258 are reported here that were previously missed. Its northern radio lobe underwent a major morphological variation that rendered the hotspot undetectable in 2001 and reappeared again in the following years. The reported changes confirm the Galactic nature of the source. We tentatively interpret them in terms of the growth of instabilities in the jet flow. There is also evidence of surrounding cocoon. These results can provide a testbed for models accounting for the evolution of jets and their interaction with the environment.


title: Condensate Nuclei and Magnetic Polarity Reversals in the Sun and Solar-type Stars

abstract: Magnetic field generation in the Sun and solar-type stars is modeled here based on the formation of magnetically polarized condensates (B. M. Mirza, Mod. Phys. Lett. B 28 (2014) 1450148) under the high density and pressure conditions. The correct orders of magnitude for the time period and the energy loss in a solar cycle are deduced, as well as the enhancement in the energy emission observed during solar cycles. It is shown that this feedback magnetic field along with differential rotation is sufficient to generate the toroidal magnetic field in the solar exterior. The model is useful in determining magnetic polarity reversals, energy generation, and reversal times, in solar-type stars.


title: Cumulant expansions for atmospheric flows

abstract: The equations governing atmospheric flows are nonlinear. Consequently, the hierarchy of cumulant equations is not closed. But because atmospheric flows are inhomogeneous and anisotropic, the nonlinearity may manifest itself only weakly through interactions of mean fields with disturbances such as thermals or eddies. In such situations, truncations of the hierarchy of cumulant equations hold promise as a closure strategy. We review how truncations at second order can be used to model and elucidate the dynamics of atmospheric flows. Two examples are considered. First, we study the growth of a dry convective boundary layer, which is heated from below, leading to turbulent upward energy transport and growth of the boundary layer. We demonstrate that a quasilinear truncation of the equations of motion, in which interactions of disturbances among each other are neglected but interactions with mean fields are taken into account, can capture the growth of the convective boundary layer even if it does not capture important turbulent transport terms. Second, we study the evolution of two-dimensional large-scale waves representing waves in Earth's upper atmosphere. We demonstrate that a cumulant expansion truncated at second order (CE2) can capture the evolution of such waves and their nonlinear interaction with the mean flow in some circumstances, for example, when the wave amplitude is small enough or the planetary rotation rate is large enough. However, CE2 fails to capture the flow evolution when nonlinear eddy--eddy interactions in surf zones become important. Higher-order closures can capture these missing interactions. The results point to new ways in which the dynamics of turbulent boundary layers may be represented in climate models, and they illustrate different classes of nonlinear processes that can control wave dissipation and momentum fluxes in the troposphere.

categories: physics.flu-dyn cond-mat.stat-mech physics.ao-ph physics.geo-ph


title: From Floquet to Dicke: quantum spin-Hall insulator interacting with quantum light

abstract: Time-periodic perturbations due to classical electromagnetic fields are useful to engineer the topological properties of matter using the Floquet theory. Here we investigate the effect of quantized electromagnetic fields by focusing on the quantized light-matter interaction on the edge state of a quantum spin-Hall insulator. A Dicke-type superradiant phase transition occurs at arbitrary weak coupling, the electronic spectrum acquires a finite gap and the resulting ground state manifold is topological with Chern number $\pm 1$. When the total number of excitations is conserved, a photocurrent is generated along the edge, being pseudo-quantized as $\omega\ln(1/\omega)$ in the low frequency limit, and decaying as $1/\omega$ for high frequencies with $\omega$ the photon frequency. The photon spectral function exhibits a clean Goldstone mode, a Higgs like collective mode at the optical gap and the polariton continuum.


title: Search for new light gauge bosons in Higgs boson decays to four-lepton final states in $pp$ collisions at $\sqrt{s}=8$ TeV with the ATLAS detector at the LHC

abstract: This paper presents a search for Higgs bosons decaying to four leptons, either electrons or muons, via one or two light exotic gauge bosons $Z_d$, $H\to Z Z_d \to 4\ell$ or $H\to Z_d Z_d \to 4\ell$. The search was performed using $pp$ collision data corresponding to an integrated luminosity of about 20 fb$^{-1}$ at the center-of-mass energy of $\sqrt{s}=8 $TeV recorded with the ATLAS detector at the Large Hadron Collider. The observed data are well described by the Standard Model prediction. Upper bounds on the branching ratio of $H\to Z Z_d \to 4\ell$ and on the kinetic mixing parameter between the $Z_d$ and the Standard Model hypercharge gauge boson are set in the range $(1$--$9)\times10^{-5}$ and $(4$--$17)\times10^{-2}$ respectively, at 95% confidence level assuming the Standard Model branching ratio of $H\to Z Z^* \to 4\ell$, for $Z_d$ masses between 15 and 55 GeV. Upper bounds on the effective mass mixing parameter between the $Z$ and the $Z_d$ are also set using the branching ratio limits in the $H \to ZZ_d\to 4\ell$ search, and are in the range $(1.5$--$8.7)\times 10^{-4}$ for $15<m_{Z_d} < 35$ GeV. Upper bounds on the branching ratio of $H\rightarrow Z_dZ_d \to 4\ell$ and on Higgs portal coupling parameter, controlling the strength of the coupling of the Higgs boson to dark vector bosons, are set in the range $(2$--$3)\times10^{-5}$ and $(1$--$10)\times 10^{-4}$ respectively, at 95% confidence level assuming the Standard Model Higgs boson production cross sections, for $Z_d$ masses between 15 and 60 GeV.


title: Interplay between Detection Strategies and Stochastic Resonance Properties

abstract: We discuss how to exploit stochastic resonance with the methods of statistical theory of decisions. To do so, we evaluate two detection strategies: escape time analysis and strobing. For a standard quartic bistable system with a periodic drive and disturbed by noise, we show that the detection strategies and the physics of the double well are connected, inasmuch as one (the strobing strategy) is based on synchronization, while the other (escape time analysis) is determined by the possibility to accumulate energy in the oscillations. The analysis of the escape times best performs at the frequency of the geometric resonance, while strobing shows a peak of the performances at a special noise level predicted by the stochastic resonance theory. We surmise that the detection properties of the quartic potential are generic for overdamped and underdamped systems, in that the physical nature of resonance decides the competition (in terms of performances) between different detection strategies.


title: Visual Search at Pinterest

abstract: We demonstrate that, with the availability of distributed computation platforms such as Amazon Web Services and open-source tools, it is possible for a small engineering team to build, launch and maintain a cost-effective, large-scale visual search system with widely available tools. We also demonstrate, through a comprehensive set of live experiments at Pinterest, that content recommendation powered by visual search improve user engagement. By sharing our implementation details and the experiences learned from launching a commercial visual search engines from scratch, we hope visual search are more widely incorporated into today's commercial applications.


title: Flexible Queueing Architectures

abstract: We study a multi-server model with $n$ flexible servers and $n$ queues, connected through a bipartite graph, where the level of flexibility is captured by the graph's average degree, $d_n$. Applications in content replication in data centers, skill-based routing in call centers, and flexible supply chains are among our main motivations. We focus on the scaling regime where the system size $n$ tends to infinity, while the overall traffic intensity stays fixed. We show that a large capacity region and an asymptotically vanishing queueing delay are simultaneously achievable even under limited flexibility ($d_n \ll n$). Our main results demonstrate that, when $d_n\gg \ln n$, a family of expander-graph-based flexibility architectures has a capacity region that is within a constant factor of the maximum possible, while simultaneously ensuring a diminishing queueing delay for all arrival rate vectors in the capacity region. Our analysis is centered around a new class of virtual-queue-based scheduling policies that rely on dynamically constructed job-to-server assignments on the connectivity graph. For comparison, we also analyze a natural family of modular architectures, which is simpler but has provably weaker performance.


title: A trust-region method for stochastic variational inference with applications to streaming data

abstract: Stochastic variational inference allows for fast posterior inference in complex Bayesian models. However, the algorithm is prone to local optima which can make the quality of the posterior approximation sensitive to the choice of hyperparameters and initialization. We address this problem by replacing the natural gradient step of stochastic varitional inference with a trust-region update. We show that this leads to generally better results and reduced sensitivity to hyperparameters. We also describe a new strategy for variational inference on streaming data and show that here our trust-region method is crucial for getting good performance.


title: Nonsmooth viscosity solutions of elementary symmetric functions of the complex Hessian

abstract: In this paper we prove the existence of nonsmooth viscosity solutions for Dirichlet problems involving elementary symmetric functions of the eigenvalues of the complex Hessian.


title: Graphs with small diameter determined by their $D$-spectra

abstract: Let $G$ be a connected graph with vertex set $V(G)=\{v_{1},v_{2},...,v_{n}\}$. The distance matrix $D(G)=(d_{ij})_{n\times n}$ is the matrix indexed by the vertices of $G,$ where $d_{ij}$ denotes the distance between the vertices $v_{i}$ and $v_{j}$. Suppose that $\lambda_{1}(D)\geq\lambda_{2}(D)\geq\cdots\geq\lambda_{n}(D)$ are the distance spectrum of $G$. The graph $G$ is said to be determined by its $D$-spectrum if with respect to the distance matrix $D(G)$, any graph having the same spectrum as $G$ is isomorphic to $G$. In this paper, we give the distance characteristic polynomial of some graphs with small diameter, and also prove that these graphs are determined by their $D$-spectra.


title: Dispersion for Two Classes of Random Variables: General Theory and Application to Inference of an External Ligand Concentration by a Cell

abstract: We derive expressions for the dispersion for two classes of random variables in Markov processes. Random variables like current and activity pertain to the first class, which is composed by random variables that change whenever a jump in the stochastic trajectory occurs. The second class corresponds to the time the trajectory spends in a state (or cluster of states). While the expression for the first class follows straightforwardly from known results in the literature, we show that a similar formalism can be used to derive an expression for the second class. As an application, we use this formalism to analyze a cellular two-component network estimating an external ligand concentration. The uncertainty related to this external concentration is calculated by monitoring different random variables related to an internal protein. We show that, inter alia, monitoring the time spent in the phosphorylated state of the protein leads to a finite uncertainty only if there is dissipation, whereas the uncertainty obtained from the activity of the transitions of the internal protein can reach the Berg and Purcell limit even in equilibrium.


title: Parity measurement of remote qubits using dispersive coupling and photodetection

abstract: Parity measurement is a key step in many entanglement generation and quantum error correction schemes. We propose a protocol for non-destructive parity measurement of two remote qubits, i.e., macroscopically separated qubits with no direct interaction. The qubits are instead dispersively coupled to separate resonators that radiate to shared photodetectors. The scheme is deterministic in the sense that there is no fundamental bound on the success probability. Compared to previous proposals, our protocol addresses the scenario where number resolving photodetectors are available but the qubit-resonator coupling is time-independent and only dispersive.


title: A Study of Groups through Transversals

abstract: In this note, a necessary and sufficient condition for the normalizer of a core-free subgroup $H$ of a finite group $G$ to be normal in $G$ is obtained. Also, a known result of finite groups is obtained through transversal.


title: QCD Equation of State and Cosmological Parameters in Early Universe

abstract: The time evolution of cosmological parameters in early Universe at the deconfinement transition is studied by an equation of state (EoS) which takes into account the finite baryon density and the background magnetic field. The non perturbative dynamics is described by the Field Correlator Method (FCM) which gives, with a small number of free parameters, a good fit of lattice data. The entire system has two components, i.e. the quark-gluon plasma and the electroweak sector, and the solutions of the Friedmann equation show that the scale factor, $a(t)$, and $H(t)= (1/a)da/dt$ are weakly dependent on the EoS, but the deceleration parameter, $q(t)$, and the jerk, $j(t)$, are strongly modified above the critical temperature $T_c$, corresponding to a critical time $t_c \simeq 20-25 \mu s$. The time evolution of the cosmological parameters suggest that above and around $T_c$ there is a transient state of acceleration typical of a matter dominated Universe; this is entailed by the QCD strong interaction driven by the presence of massive colored objects.


title: {\sigma}-Porosity of the set of strict contractions in a space of non-expansive mappings

abstract: We consider the space of non-expansive mappings on a bounded, closed and convex subset of a Banach space equipped with the metric of uni- form convergence. We show that the set of strict contractions is a {\sigma}-porous subset.


title: Aspects of Quadratic Gravity

abstract: We discuss quadratic gravity where terms quadratic in the curvature tensor are included in the action. After reviewing the corresponding field equations, we analyze in detail the physical propagating modes in some specific backgrounds. First we confirm that the pure $R^2$ theory is indeed ghost free. Then we point out that for flat backgrounds the pure $R^2$ theory propagates only a scalar massless mode and no spin-two tensor mode. However, the latter emerges either by expanding the theory around curved backgrounds like de Sitter or anti-de Sitter, or by changing the long-distance dynamics by introducing the standard Einstein term. In both cases, the theory is modified in the infrared and a propagating graviton is recovered. Hence we recognize a subtle interplay between the UV and IR properties of higher order gravity. We also calculate the corresponding Newton's law for general quadratic curvature theories. Finally, we discuss how quadratic actions may be obtained from a fundamental theory like string- or M-theory. We demonstrate that string theory on non-compact $CY_3$ manifolds, like a line bundle over $\mathbb{CP}^2$, may indeed lead to gravity dynamics determined by a higher curvature action.


title: High-throughput in-situ characterization and modelling of precipitation kinetics in compositionally graded alloys

abstract: The development of new engineering alloy chemistries is a time consuming and iterative process. A necessary step is characterization of the nano/microstructure to provide a link between the processing and properties of each alloy chemistry considered. One approach to accelerate the identification of optimal chemistries is to use samples containing a gradient in composition, ie. combinatorial samples, and to investigate many different chemistries at the same time. However, for engineering alloys, the final properties depend not only on chemistry but also on the path of microstructure development which necessitates characterization of microstructure evolution for each chemistry. In this contribution we demonstrate an approach that allows for the in-situ, nanoscale characterization of the precipitate structures in alloys, as a function of aging time, in combinatorial samples containing a composition gradient. The approach uses small angle x-ray scattering (SAXS) at a synchrotron beamline. The Cu-Co system is used for the proof-of-concept and the combinatorial samples prepared contain a gradient in Co from 0% to 2%. These samples are aged at temperatures between 450{\textdegree}C and 550{\textdegree}C and the precipitate structures (precipitate size, volume fraction and number density) all along the composition gradient are simultaneously monitored as a function of time. This large dataset is used to test the applicability and robustness of a conventional class model for precipitation that considers concurrent nucleation, growth and coarsening and the ability of the model to describe such a large dataset.


title: Synchronising and Non-synchronising Dynamics for a Two-species Aggregation Model

abstract: This paper deals with analysis and numerical simulations of a one-dimensional two-species hyperbolic aggregation model. This model is formed by a system of transport equations with nonlocal velocities, which describes the aggregate dynamics of a two-species population in interaction appearing for instance in bacterial chemotaxis. Blow-up of classical solutions occurs in finite time. This raises the question to define measure-valued solutions for this system. To this aim, we use the duality method developed for transport equations with discontinuous velocity to prove the existence and uniqueness of measure-valued solutions. The proof relies on a stability result. In addition, this approach allows to study the hyperbolic limit of a kinetic chemotaxis model. Moreover, we propose a finite volume numerical scheme whose convergence towards measure-valued solutions is proved. It allows for numerical simulations capturing the behaviour after blow up. Finally, numerical simulations illustrate the complex dynamics of aggregates until the formation of a single aggregate: after blow-up of classical solutions, aggregates of different species are synchronising or nonsynchronising when collide, that is move together or separately, depending on the parameters of the model and masses of species involved.


title: Percolation in suspensions of hard nanoparticles: From spheres to needles

abstract: We investigate geometric percolation and scaling relations in suspensions of nanorods, covering the entire range of aspect ratios from spheres to extremely slender needles. A new version of connectedness percolation theory is introduced and tested against specialized Monte Carlo simulations. The theory accurately predicts percolation thresholds for aspect ratios as low as 10. The percolation threshold for rod-like particles of aspect ratios below 1000 deviates significantly from the inverse aspect ratio scaling prediction, thought to be valid in the limit of infinitely slender rods and often used as a rule of thumb for nano-fibers in composite materials. Hence, most fibers that are currently used as fillers in composite materials cannot be regarded as practically infinitely slender for the purposes of percolation theory. Comparing percolation thresholds of hard rods and new benchmark results for ideal rods, we find that (i) for large aspect ratios, they differ by a factor that is inversely proportional to the connectivity distance between the hard cores, and (ii) they approach the slender rod limit differently.


title: Reactive point processes: A new approach to predicting power failures in underground electrical systems

abstract: Reactive point processes (RPPs) are a new statistical model designed for predicting discrete events in time based on past history. RPPs were developed to handle an important problem within the domain of electrical grid reliability: short-term prediction of electrical grid failures ("manhole events"), including outages, fires, explosions and smoking manholes, which can cause threats to public safety and reliability of electrical service in cities. RPPs incorporate self-exciting, self-regulating and saturating components. The self-excitement occurs as a result of a past event, which causes a temporary rise in vulner ability to future events. The self-regulation occurs as a result of an external inspection which temporarily lowers vulnerability to future events. RPPs can saturate when too many events or inspections occur close together, which ensures that the probability of an event stays within a realistic range. Two of the operational challenges for power companies are (i) making continuous-time failure predictions, and (ii) cost/benefit analysis for decision making and proactive maintenance. RPPs are naturally suited for handling both of these challenges. We use the model to predict power-grid failures in Manhattan over a short-term horizon, and to provide a cost/benefit analysis of different proactive maintenance programs.


title: Determining Spin Polarization of Seebeck Coefficients via Anomalous Nernst Effect

abstract: Recently, Seebeck coefficients of ferromagnetic conductors are found to be spin-dependent. However straightforward method of accurately determining its spin polarization is still to be developed. Here, we have derived a linear dependence of anomalous Nernst coefficient on anomalous Hall angle with scaling factor related to spin polarization of Seebeck coefficient, which has been experimentally verified in [Co/Pt]n superlattices. Based on the dependence, we have also evaluated spin polarization of Seebeck coefficient of some ferromagnetic conductors. Besides, we have also found a new mechanism to generate pure spin current from temperature gradient in ferromagnetic/nonmagnetic hybrid system, which could improve efficiency from thermal energy to spin current.


title: Externally triggered coherent two-photon emission from hydrogen molecules

abstract: We report coherent enhancement of two-photon emission from the excited vibrational state of molecular hydrogen triggered by irradiating mid-infrared pulses externally. We previously observed the two-photon emission triggered by the internally generated fourth Stokes photons. By injecting independent mid-infrared pulses externally, it is possible to control experimental parameters and investigate the mechanism in more detail. In this article, we describe the two-photon emission using the external trigger pulses. Its spectrum and dependence on the energy and timing of the trigger pulse are presented along with numerical simulations based on the Maxwell-Bloch equations. The measured number of emitted photons is 6 10^11 photons/pulse and the resulting enhancement factor from the spontaneous emission is more than 10^18. This value is three orders of magnitude higher than that of the previous experiment. External control of emission process is expected to be essential for observation of weaker process of radiative emission of neutrino pair.


title: Empirical Spacings of Unfolded Eigenvalues

abstract: We study random points on the real line generated by the eigenvalues in unitary invariant random matrix ensembles or by more general repulsive particle systems. As the number of points tends to infinity, we prove convergence of the empirical distribution of nearest neighbor spacings. We extend existing results for the spacing distribution in two ways. On the one hand, we believe the empirical distribution to be of more practical relevance than the so far considered expected distribution. On the other hand, we use the unfolding, a non-linear rescaling, which transforms the ensemble such that the density of particles is asymptotically constant. This allows to consider all empirical spacings, where previous results were restricted to a tiny fraction of the particles. Moreover, we prove bounds on the rates of convergence. The main ingredient for the proof, a strong bulk universality result for correlation functions in the unfolded setting including optimal rates, should be of independent interest.


title: Brick polytopes, lattice quotients, and Hopf algebras

abstract: This paper is motivated by the interplay between the Tamari lattice, J.-L. Loday's realization of the associahedron, and J.-L. Loday and M. Ronco's Hopf algebra on binary trees. We show that these constructions extend in the world of acyclic $k$-triangulations, which were already considered as the vertices of V. Pilaud and F. Santos' brick polytopes. We describe combinatorially a natural surjection from the permutations to the acyclic $k$-triangulations. We show that the fibers of this surjection are the classes of the congruence $\equiv^k$ on $\mathfrak{S}_n$ defined as the transitive closure of the rewriting rule $U ac V_1 b_1 \cdots V_k b_k W \equiv^k U ca V_1 b_1 \cdots V_k b_k W$ for letters $a < b_1, \dots, b_k < c$ and words $U, V_1, \dots, V_k, W$ on $[n]$. We then show that the increasing flip order on $k$-triangulations is the lattice quotient of the weak order by this congruence. Moreover, we use this surjection to define a Hopf subalgebra of C. Malvenuto and C. Reutenauer's Hopf algebra on permutations, indexed by acyclic $k$-triangulations, and to describe the product and coproduct in this algebra and its dual in term of combinatorial operations on acyclic $k$-triangulations. Finally, we extend our results in three directions, describing a Cambrian, a tuple, and a Schr\"oder version of these constructions.


title: Mixed timestepping schemes for nonsmooth mechanics with high frequency damping

abstract: This work deals with the integration of nonsmooth flexible multibody systems with impacts and dry friction. We develop a framework which improves a non-impulsive trajectory of state variables by impulsive correction after each time-step if necessary. This correction is automatic and is evaluated on the same kinematic level as the piecewise non-impulsive trajectory. The resulting overall mixed timestepping scheme is consistent with respect to impacts and friction as well as benefits from advantages of the base integration schemes used to calculate the approximation inside the time-step. Therefore, we compare the generalized-$\alpha$ method, the Bathe method and the ED-$\alpha$ method.


title: Intrinsic spin Hall effect in systems with striped spin-orbit coupling

abstract: The Rashba spin-orbit coupling arising from structure inversion asymmetry couples spin and momentum degrees of freedom providing a suitable (and very intensively investigated) environment for spintronic effects and devices. Here we show that in the presence of strong disorder, non-homogeneity in the spin-orbit coupling gives rise to a finite spin Hall conductivity in contrast with the corresponding case of a homogeneous linear spin-orbit coupling. In particular, we examine the inhomogeneity arising from a striped structure for a two-dimensional electron gas, affecting both density and Rashba spin-orbit coupling. We suggest that this situation can be realized at oxide interfaces with periodic top gating.


title: Higgs-$\mu$-$\tau$ Coupling at High and Low Energy Colliders

abstract: There is no tree-level flavor changing neutral current (FCNC) in the standard model (SM) which contains only one Higgs doublet. If more Higgs doublets are introduced for various reasons, the tree level FCNC would be inevitable except extra symmetry was imposed. Therefore FCNC processes are the excellent probes for the physics beyond the SM (BSM). In this paper, we studied the lepton flavor violated (LFV) decay processes $h\rightarrow\mu\tau$ and $\tau\rightarrow\mu\gamma$ induced by Higgs-$\mu$-$\tau$ vertex. For $\tau\rightarrow\mu\gamma$, its branching ratio is also related to the $ht\bar{t}$, $h\tau^+\tau^-$ and $hW^+W^-$ vertices. We categorized the BSM into two scenarios for the Higgs coupling strengths near or away from SM. For the latter scenario, we took the spontaneously broken two Higgs doublet model (Lee model) as an example. We considered the constraints by recent data from LHC and B factories, and found that the measurements gave weak constraints. At LHC Run II, $h\rightarrow\mu\tau$ will be confirmed or set stricter limit on its branching ratio. Accordingly, $\textrm{Br}(\tau\rightarrow\mu\gamma)\lesssim\mathcal{O}(10^{-10}-10^{-8})$ for general chosen parameters. For the positive case, $\tau\rightarrow\mu\gamma$ can be discovered with $\mathcal{O}(10^{10})$ $\tau$ pair samples at SuperB factory, Super $\tau$-charm factory and new Z-factory. The future measurements for $\textrm{Br}(h\rightarrow\mu\tau)$ and $\textrm{Br}(\tau\rightarrow\mu\gamma)$ will be used to distinguish these two scenarios or set strict constraints on the correlations among different Higgs couplings, please see Table II in the text for details.


title: Imaging the cold molecular gas in SDSS J1148 + 5251 at z = 6.4

abstract: We present Karl G. Jansky Very Large Array (VLA) observations of the CO ($J = 2 \rightarrow 1$) line emission towards the $z = 6.419$ quasar SDSS J$114816.64+525150.3$ (J$1148+5251$). The molecular gas is found to be marginally resolved with a major axis of $0.9"$ (consistent with previous size measurements of the CO ($J = 7 \rightarrow 6$) emission). We observe tentative evidence for extended line emission towards the south west on a scale of ~$1.4"$, but this is only detected at $3.3\sigma$ significance and should be confirmed. The position of the molecular emission region is in excellent agreement with previous detections of low frequency radio continuum emission as well as [C ii] line and thermal dust continuum emission. These CO ($J = 2 \rightarrow 1$) observations provide an anchor for the low excitation part of the molecular line SED. We find no evidence for extended low excitation component, neither in the spectral line energy distribution nor the image. We fit a single kinetic gas temperature model of 50 K. We revisit the gas and dynamical masses in light of this new detection of a low order transition of CO, and confirm previous findings that there is no extended reservoir of cold molecular gas in J$1148+5251$, and that the source departs substantially from the low $z$ relationship between black hole mass and bulge mass. Hence, the characteristics of J$1148+5251$ at $z = 6.419$ are very similar to $z$~$2$ quasars, in the lack of a diffuse cold gas reservoir and kpc-size compactness of the star forming region.


title: Domain Wall Formation from Level Crossing in the Axiverse

abstract: We point out that domain wall formation is a more common phenomenon in the Axiverse than previously thought. Level crossing could take place if there is a mixing between axions, and if some of the axions acquire a non-zero mass through non-perturbative effects as the corresponding gauge interactions become strong. The axion potential changes significantly during the level crossing, which affects the axion dynamics in various ways. We find that, if there is a mild hierarchy in the decay constants, the axion starts to run along the valley of the potential, passing through many crests and troughs, until it gets trapped in one of the minima; the {\it axion roulette}. The axion dynamics exhibits a chaotic behavior during the oscillations, and which minimum the axion is finally stabilized is highly sensitive to the initial misalignment angle. Therefore, the axion roulette is considered to be accompanied by domain wall formation. The cosmological domain wall problem can be avoided by introducing a small bias between the vacua. We discuss cosmological implications of the domain wall annihilation for baryogenesis and future gravitational wave experiments.


title: The Age-Redshift Relationship of Old Passive Galaxies

abstract: We use 32 age measurements of passively evolving galaxies as a function of redshift to test and compare the standard model ($\Lambda$CDM) with the $R_{\rm h}=ct$ Universe. We show that the latter fits the data with a reduced $\chi^2_{\rm dof}=0.435$ for a Hubble constant $H_{0}= 67.2_{-4.0}^{+4.5}$ km $\rm s^{-1}$ $\rm Mpc^{-1}$. By comparison, the optimal flat $\Lambda$CDM model, with two free parameters (including $\Omega_{\rm m}=0.12_{-0.11}^{+0.54}$ and $H_{0}=94.3_{-35.8}^{+32.7}$ km $\rm s^{-1}$ $\rm Mpc^{-1}$), fits the age-\emph{z} data with a reduced $\chi^2_{\rm dof}=0.428$. Based solely on their $\chi^2_{\rm dof}$ values, both models appear to account for the data very well, though the optimized $\Lambda$CDM parameters are only marginally consistent with those of the concordance model ($\Omega_{\rm m}=0.27$ and $H_{0}= 70$ km $\rm s^{-1}$ $\rm Mpc^{-1}$). Fitting the age-$z$ data with the latter results in a reduced $\chi^2_{\rm dof}=0.523$. However, because of the different number of free parameters in these models, selection tools, such as the Akaike, Kullback and Bayes Information Criteria, favour $R_{\rm h}=ct$ over $\Lambda$CDM with a likelihood of $\sim 66.5\%-80.5\%$ versus $\sim 19.5\%-33.5\%$. These results are suggestive, though not yet compelling, given the current limited galaxy age-$z$ sample. We carry out Monte Carlo simulations based on these current age measurements to estimate how large the sample would have to be in order to rule out either model at a $\sim 99.7\%$ confidence level. We find that if the real cosmology is $\Lambda$CDM, a sample of $\sim 45$ galaxy ages would be sufficient to rule out $R_{\rm h}=ct$ at this level of accuracy, while $\sim 350$ galaxy ages would be required to rule out $\Lambda$CDM if the real Universe were instead $R_{\rm h}=ct$.


title: A Generative Model of Natural Texture Surrogates

abstract: Natural images can be viewed as patchworks of different textures, where the local image statistics is roughly stationary within a small neighborhood but otherwise varies from region to region. In order to model this variability, we first applied the parametric texture algorithm of Portilla and Simoncelli to image patches of 64X64 pixels in a large database of natural images such that each image patch is then described by 655 texture parameters which specify certain statistics, such as variances and covariances of wavelet coefficients or coefficient magnitudes within that patch. To model the statistics of these texture parameters, we then developed suitable nonlinear transformations of the parameters that allowed us to fit their joint statistics with a multivariate Gaussian distribution. We find that the first 200 principal components contain more than 99% of the variance and are sufficient to generate textures that are perceptually extremely close to those generated with all 655 components. We demonstrate the usefulness of the model in several ways: (1) We sample ensembles of texture patches that can be directly compared to samples of patches from the natural image database and can to a high degree reproduce their perceptual appearance. (2) We further developed an image compression algorithm which generates surprisingly accurate images at bit rates as low as 0.14 bits/pixel. Finally, (3) We demonstrate how our approach can be used for an efficient and objective evaluation of samples generated with probabilistic models of natural images.


title: An impulsive dynamical systems framework for reset control systems

abstract: Impulsive dynamical systems is a well-established area of dynamical systems theory, and it is used in this work to analyze several basic properties of reset control systems: existence and uniqueness of solutions, and continuous dependence on the initial condition (well-posedness). The work scope is about reset control systems with a linear and time-invariant base system, and a zero-crossing resetting law. A necessary and sufficient condition for existence and uniqueness of solutions, based on the well-posedness of reset instants, is developed. As a result, it is shown that reset control systems (with strictly proper plants) do no have Zeno solutions. It is also shown that full reset and partial reset (with a special structure) always produce well-posed reset instants. Moreover, a definition of continuous dependence on the initial condition is developed, and also a sufficient condition for reset control systems to satisfy that property. Finally, this property is used to analyze sensitivity of reset control systems to sensor noise. This work also includes a number of illustrative examples motivating the key concepts and main results.


title: Lensed: a code for the forward reconstruction of lenses and sources from strong lensing observations

abstract: Robust modelling of strong lensing systems is fundamental to exploit the information they contain about the distribution of matter in galaxies and clusters. In this work, we present Lensed, a new code which performs forward parametric modelling of strong lenses. Lensed takes advantage of a massively parallel ray-tracing kernel to perform the necessary calculations on a modern graphics processing unit (GPU). This makes the precise rendering of the background lensed sources much faster, and allows the simultaneous optimisation of tens of parameters for the selected model. With a single run, the code is able to obtain the full posterior probability distribution for the lens light, the mass distribution and the background source at the same time. Lensed is first tested on mock images which reproduce realistic space-based observations of lensing systems. In this way, we show that it is able to recover unbiased estimates of the lens parameters, even when the sources do not follow exactly the assumed model. Then, we apply it to a subsample of the SLACS lenses, in order to demonstrate its use on real data. The results generally agree with the literature, and highlight the flexibility and robustness of the algorithm.


title: Improved Deep Convolutional Neural Network For Online Handwritten Chinese Character Recognition using Domain-Specific Knowledge

abstract: Deep convolutional neural networks (DCNNs) have achieved great success in various computer vision and pattern recognition applications, including those for handwritten Chinese character recognition (HCCR). However, most current DCNN-based HCCR approaches treat the handwritten sample simply as an image bitmap, ignoring some vital domain-specific information that may be useful but that cannot be learnt by traditional networks. In this paper, we propose an enhancement of the DCNN approach to online HCCR by incorporating a variety of domain-specific knowledge, including deformation, non-linear normalization, imaginary strokes, path signature, and 8-directional features. Our contribution is twofold. First, these domain-specific technologies are investigated and integrated with a DCNN to form a composite network to achieve improved performance. Second, the resulting DCNNs with diversity in their domain knowledge are combined using a hybrid serial-parallel (HSP) strategy. Consequently, we achieve a promising accuracy of 97.20% and 96.87% on CASIA-OLHWDB1.0 and CASIA-OLHWDB1.1, respectively, outperforming the best results previously reported in the literature.


title: MADMM: a generic algorithm for non-smooth optimization on manifolds

abstract: Numerous problems in machine learning are formulated as optimization with manifold constraints. In this paper, we propose the Manifold alternating directions method of multipliers (MADMM), an extension of the classical ADMM scheme for manifold-constrained non-smooth optimization problems and show its application to several challenging problems in dimensionality reduction, data analysis, and manifold learning.


title: Large ground state magnetic moment and magnetocaloric effect in Ni2Mn1.4In0.6

abstract: A large conventional magnetocaloric effect at the second order magnetic transition in cubic Ni2Mn1.4In0.6 Heusler alloy is reported. The isothermal magnetization at 2K shows a huge ground state magnetic moment of about 6.17 {\mu}B/f.u. The theoretical calculations show that the origin of the large magnetic moment in cubic Ni2Mn1.4In0.6 results from the strong ferromagnetic interaction between Mn- Ni and Mn-Mn sublattices. The experimental magnetic moment is in excellent agreement with the moment calculated from the theory. The large magnetic moment gives rise to considerably high adiabatic temperature and entropy changes at the magnetic transition. The present study opens up the possibility to explore cubic Heusler alloys for magnetocaloric applications.


title: Microscopic evaluation of the hypernuclear chart with $\Lambda$-hyperons

abstract: A large number of hypernuclei, where a considerable fraction of nucleons is replaced by strange baryons, and even pure hyperonic species are expected to be bound. Though, the hypernuclear landscape remains largely unknown because of scarce constraints on the $NY$ and $YY$ interactions. We want to estimate the number of potentially bound hypernuclei. In order to evaluate realistic error bars within the theoretical uncertainties associated to the spherical mean-field approach, and the present information from already synthetized hypernuclei on the $N-Y$ and $Y-Y$ channels, we limit ourselves to purely $\Lambda$-hypernuclei, to magic numbers of $\Lambda$'s (for Z $\leq$ 120 and $\Lambda \leq$70), and to even-even-even systems. We consider a density functional approach adjusted to microscopic Bruckner-Hartree-Fock calculations, where the $\Lambda\Lambda$ term is corrected in a phenomenological way, to reproduce present experimental constraints. The number of bound even-even-even $\Lambda$-hypernuclei is estimated to 491680 $\pm$ 34400. This relatively low uncertainty is due to the fact that the well constrained low density and highly unconstrained high density behavior of the energy functional turn out to be largely decoupled. Results in $\Lambda$-hypernuclei appear to be almost independent of the choice for the high-density part of the $\Lambda\Lambda$ interaction. The location of the $\Lambda$-hyperdriplines is also evaluated. Significant deviations from Iron-Nickel elements can be found for $\Lambda$-hypernuclei with the largest binding energy per baryon. Proton, neutron and $\Lambda$-hyperon magicity evolution and triple magic $\Lambda$-hypernuclei are studied. Possible bubbles and haloes effect in $\Lambda$-hypernuclei are also discussed.


title: On effective K\"ahler potential in N=2, d=3 SQED

abstract: We compute the two-loop effective K\"ahler potential in three-dimensional N=2 supersymmetric electrodynamics with Chern-Simons kinetic term for the gauge superfield. The effective action is constructed on the base of background field method with one parametric family of gauges. In such an approach, the quadratic part of quantum action mixes the gauge and matter quantum superfields yielding the complications in the computations of the loop supergraphs. To avoid this obstacle and preserve dependence on the gauge parameter we make a nonlocal change of quantum matter superfields after which the propagator is diagonalized, however the new vertices have appeared. We fix the suitable background and develop the efficient procedure of calculating the two-loop supergraphs with the new vertices. We compute the divergent and finite parts of the superfield effective action, find the two-loop effective K\"ahler potential and show that it does not depend on the gauge parameter.


title: Metric-Independent Volume-Forms in Gravity and Cosmology

abstract: Employing alternative spacetime volume-forms (generally-covariant integration measure densities) independent of the pertinent Riemannian spacetime metric have profound impact in general relativity. Although formally appearing as "pure-gauge" dynamical degrees of freedom they trigger a number of remarkable physically important phenomena such as: (i) new mechanism of dynamical generation of cosmological constant; (ii) new type of "quintessential inflation" scenario in cosmology; (iii) non-singular initial "emergent universe" phase of cosmological evolution preceding the inflationary phase; (iv) new mechanism of dynamical spontaneous breakdown of supersymmetry in supergravity; (v) gravitational electrovacuum "bags". We study in some detail the properties, together with their canonical Hamiltonian formulation, of a class of generalized gravity-matter models built with two independent non-Riemannian volume-forms and discuss their implications in cosmology.


title: Random knots using Chebyshev billiard table diagrams

abstract: We use the Chebyshev knot diagram model of Koseleff and Pecker in order to introduce a random knot diagram model by assigning the crossings to be positive or negative uniformly at random. We give a formula for the probability of choosing a knot at random among all knots with bridge index at most 2. Restricted to this class, we define internal and external reduction moves that decrease the number of crossings of the diagram. We make calculations based on our formula, showing the numerics in graphs and providing evidence for our conjecture that the probability of any knot appearing in this model decays to zero as the number of crossings goes to infinity.


title: Optical blocking of electron tunneling into a single self-assembled quantum dot

abstract: Time-resolved resonance fluorescence (RF) is used to analyse electron tunneling between a single self-assembled quantum dot (QD) and an electron reservoir. In equilibrium, the RF intensity reflects the average electron occupation of the QD and exhibits a gate voltage dependence that is given by the Fermi distribution in the reservoir. In the time-resolved signal, however, we find that the relaxation rate for electron tunneling is independent of the occupation in the charge reservoir. Using a master equation approach which includes both the electron tunneling and the optical excitation/recombination, we are able to explain the experimental data by optical blocking, i. e. a reduced electron tunneling rate when the QD is occupied by an exciton.


title: Privacy in the Internet of Things: Threats and Challenges

abstract: The Internet of Things paradigm envisions the pervasive interconnection and cooperation of smart things over the current and future Internet infrastructure. The Internet of Things is, thus, the evolution of the Internet to cover the real-world, enabling many new services that will improve people's everyday lives, spawn new businesses and make buildings, cities and transport smarter. Smart things allow indeed for ubiquitous data collection or tracking, but these useful features are also examples of privacy threats that are already now limiting the success of the Internet of Things vision when not implemented correctly. These threats involve new challenges such as the pervasive privacy-aware management of personal data or methods to control or avoid ubiquitous tracking and profiling. This paper analyzes the privacy issues in the Internet of Things in detail. To this end, we first discuss the evolving features and trends in the Internet of Things with the goal of scrutinizing their privacy implications. Second, we classify and examine privacy threats in this new setting, pointing out the challenges that need to be overcome to ensure that the Internet of Things becomes a reality.


title: The Extreme Risk of Personal Data Breaches & The Erosion of Privacy

abstract: Personal data breaches from organisations, enabling mass identity fraud, constitute an \emph{extreme risk}. This risk worsens daily as an ever-growing amount of personal data are stored by organisations and on-line, and the attack surface surrounding this data becomes larger and harder to secure. Further, breached information is distributed and accumulates in the hands of cyber criminals, thus driving a cumulative erosion of privacy. Statistical modeling of breach data from 2000 through 2015 provides insights into this risk: A current maximum breach size of about 200 million is detected, and is expected to grow by fifty percent over the next five years. The breach sizes are found to be well modeled by an \emph{extremely heavy tailed} truncated Pareto distribution, with tail exponent parameter decreasing linearly from 0.57 in 2007 to 0.37 in 2015. With this current model, given a breach contains above fifty thousand items, there is a ten percent probability of exceeding ten million. A size effect is unearthed where both the frequency and severity of breaches scale with organisation size like $s^{0.6}$. Projections indicate that the total amount of breached information is expected to double from two to four billion items within the next five years, eclipsing the population of users of the Internet. This massive and uncontrolled dissemination of personal identities raises fundamental concerns about privacy.


title: Incommensurate spin density wave at a ferromagnetic quantum critical point in a three-dimensional parabolic semimetal

abstract: We explore the ferromagnetic quantum critical point in a three-dimensional semimetallic system with upward- and downward-dispersing bands touching at the Fermi level. Evaluating the static spin susceptibility to leading order in the coupling between the fermions and the fluctuating ferromagnetic order parameter, we find that the ferromagnetic quantum critical point is masked by an incommensurate, longitudinal spin density wave phase. We first analyze an idealized model which, despite having strong spin-orbit coupling, still possesses O(3) rotational symmetry generated by the total angular momentum operator. In this case, the direction of the incommensurate spin density wave propagation can point anywhere, while the magnetic moment is aligned along the direction of propagation. Including symmetry-allowed anisotropies in the fermion dispersion and the coupling to the order parameter field, however, the ordering wavevector instead breaks a discrete symmetry and aligns along either the [111] or [100] direction, depending on the signs and magnitudes of these two types of anisotropy.


title: Graphene Schottky diodes: an experimental review of the rectifying graphene/semiconductor heterojunction

abstract: In the past decade graphene has been one of the most studied material for several unique and excellent properties. Due to its two dimensional nature, physical and chemical properties and ease of manipulation, graphene offers the possibility of integration with the exiting semiconductor technology for next-generation electronic and sensing devices. In this context, the understanding of the graphene/semiconductor interface is of great importance since it can constitute a versatile standalone device as well as the building-block of more advanced electronic systems. Since graphene was brought to the attention of the scientific community in 2004, the device research has been focused on the more complex graphene transistors, while the graphene/semiconductor junction, despite its importance, has started to be the subject of systematic investigation only recently. As a result, a thorough understanding of the physics and the potentialities of this device is still missing. The studies of the past few years have demonstrated that graphene can form junctions with 3D or 2D semiconducting materials which have rectifying characteristics and behave as excellent Schottky diodes. The main novelty of these devices is the tunable Schottky barrier height, a feature which makes the graphene/semiconductor junction a great platform for the study of interface transport mechanisms as well as for applications in photo-detection, high-speed communications, solar cells, chemical and biological sensing, etc. In this paper, we review the state-of-the art of the research on graphene/semiconductor junctions, the attempts towards a modeling and the most promising applications.


title: Electron localisation in static and time-dependent one-dimensional model systems

abstract: Electron localization is the tendency of an electron in a many-body system to exclude other electrons from its vicinity. Using a new natural measure of localization based on the exact manyelectron wavefunction, we find that localization can vary considerably between different ground-state systems, and can also be strongly disrupted, as a function of time, when a system is driven by an applied electric field. We use our new measure to assess the well-known electron localization function (ELF), both in its approximate single-particle form (often applied within density-functional theory) and its full many-particle form. The full ELF always gives an excellent description of localization, but the approximate ELF fails in time-dependent situations, even when the exact Kohn-Sham orbitals are employed.


title: Anti-magic labeling of regular graphs

abstract: A graph $G=(V,E)$ is antimagic if there is a one-to-one correspondence $f: E \to \{1,2,\ldots, |E|\}$ such that for any two vertices $u,v$, $\sum_{e \in E(u)}f(e) \ne \sum_{e\in E(v)}f(e)$. It is known that bipartite regular graphs are antimagic and non-bipartite regular graphs of odd degree at least three are antimagic. Whether all non-bipartite regular graphs of even degree are antimagic remained an open problem. In this paper, we solve this problem and prove that all even degree regular graphs are antimagic. The paper was submitted to December 2014 by Journal of Graph Theory.


title: Semi-wave and spreading speed for the diffusive competition model with a free boundary

abstract: We determine the asymptotic spreading speed of an invasive species, which invades the territory of a native competitor, governed by a diffusive competition model with a free boundary in a spherically symmetric setting. This free boundary problem was studied recently in [7] but only rough bounds of the spreading speed was obtained there. We show in this paper that the spreading speed has an asymptotic limit (as time goes to infnity), which is given by a certain traveling wave type system of one space dimension, called a semi-wave. This appears to be the first result that gives the precise asymptotic spreading speed for a two species system with free boundaries.


title: Invertible Orientation Scores of 3D Images

abstract: The enhancement and detection of elongated structures in noisy image data is relevant for many biomedical applications. To handle complex crossing structures in 2D images, 2D orientation scores were introduced, which already showed their use in a variety of applications. Here we extend this work to 3D orientation scores. First, we construct the orientation score from a given dataset, which is achieved by an invertible coherent state type of transform. For this transformation we introduce 3D versions of the 2D cake-wavelets, which are complex wavelets that can simultaneously detect oriented structures and oriented edges. For efficient implementation of the different steps in the wavelet creation we use a spherical harmonic transform. Finally, we show some first results of practical applications of 3D orientation scores.


title: Lorenz-Mie theory for 2D scattering and resonance calculations

abstract: This PhD tutorial is concerned with a description of the two-dimensional generalized Lorenz-Mie theory (2D-GLMT), a well-established numerical method used to compute the interaction of light with arrays of cylindrical scatterers. This theory is based on the method of separation of variables and the application of an addition theorem for cylindrical functions. The purpose of this tutorial is to assemble the practical tools necessary to implement the 2D-GLMT method for the computation of scattering by passive scatterers or of resonances in optically active media. The first part contains a derivation of the vector and scalar Helmholtz equations for 2D geometries, starting from Maxwell's equations. Optically active media are included in 2D-GLMT using a recent stationary formulation of the Maxwell-Bloch equations called steady-state ab initio laser theory (SALT), which introduces new classes of solutions useful for resonance computations. Following these preliminaries, a detailed description of 2D-GLMT is presented. The emphasis is placed on the derivation of beam-shape coefficients for scattering computations, as well as the computation of resonant modes using a combination of 2D-GLMT and SALT. The final section contains several numerical examples illustrating the full potential of 2D-GLMT for scattering and resonance computations. These examples, drawn from the literature, include the design of integrated polarization filters and the computation of optical modes of photonic crystal cavities and random lasers.


title: A compactness theorem on Branson's $Q$-curvature equation

abstract: Let $(M, g)$ be a closed Riemannian manifold of dimension $5$. Assume that $(M, g)$ is not conformally equivalent to the round sphere. If the scalar curvature $R_g\geq 0$ and the $Q$-curvature $Q_g\geq 0$ on $M$ with $Q_g(p)>0$ for some point $p\in M$, we prove that the set of metrics in the conformal class of $g$ with prescribed constant positive $Q$-curvature is compact in $C^{4, \alpha}$ for any $0 <\alpha < 1$. We also give some estimates for dimension $6$ and $7$.


title: Stable evaluation of Green's functions in cylindrically stratified regions with uniaxial anisotropic layers

abstract: We present a robust algorithm for the computation of electromagnetic fields radiated by point sources (Hertzian dipoles) in cylindrically stratified media where each layer may exhibit material properties (permittivity, permeability, and conductivity) with uniaxial anisotropy. Analytical expressions are obtained based on the spectral representation of the tensor Green's function based on cylindrical Bessel and Hankel eigenfunctions, and extended for layered uniaxial media. Due to the poor scaling of these eigenfunctions for extreme arguments and/or orders, direct numerical evaluation of such expressions can produce numerical instability, i.e., underflow, overflow, and/or round-off errors under finite precision arithmetics. To circumvent these problems, we develop a numerically stable formulation through suitable rescaling of various expressions involved in the computational chain, to yield a robust algorithm for all parameter ranges. Numerical results are presented to illustrate the robustness of the formulation including cases of practical interest to geophysical exploration.


title: The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method

abstract: We consider the zero heat conductivity limit to a contact discontinuity for the mono-dimensional full compressible Navier-Stokes-Fourier system. The method is based on the relative entropy method, and do not assume any smallness conditions on the discontinuity, nor on the $BV$ norm of the initial data. It is proved that for any viscosity $\nu\geq0$, the solution of the compressible Navier-Stokes-Fourier system (with well prepared initial value) converges, when the heat conductivity $\kappa$ tends to zero, to the contact discontinuity solution to the corresponding Euler system. We obtain the decay rate $\kappa^{\frac12}$. It implies that the heat conductivity dominates the dissipation in the regime of the limit to a contact discontinuity. This is the first result, based on the relative entropy, of an asymptotic limit to a discontinuous solutions for a system.


title: Probing gauge-phobic heavy Higgs bosons at high energy hadron colliders

abstract: We study the probe of the gauge-phobic (or nearly gauge-phobic) heavy Higgs bosons (GPHB) at high energy hadron colliders including the 14 TeV LHC and the 50 TeV Super Proton-Proton Collider (SppC). We take the process $pp\to t\bar t t\bar t$, and study it at the hadron level including simulating the jet formation and top quark tagging (with jet substructure). We show that, for a GPHB with $M^{}_H<800$ GeV, $M^{}_H$ can be determined by adjusting the value of $M^{}_H$ in the theoretical $p^{}_T(b_1)$ distribution to fit the observed $p^{}_T(b_1)$ distribution, and the resonance peak can be seen at the SppC for $M^{}_H$=800 GeV and 1 TeV.


title: H2O vapor excitation in dusty AGB envelopes. A PACS view of OH 127.8+0.0

abstract: Context: AGB stars lose a large percentage of their mass in a dust-driven wind. This creates a circumstellar envelope, which can be studied through thermal dust emission and molecular emission lines. In the case of high mass-loss rates, this study is complicated by the high optical depths and the intricate coupling between gas and dust radiative transfer characteristics. An important aspect of the physics of gas-dust interactions is the strong influence of dust on the excitation of several molecules, including H2O. Aims: The dust and gas content of the envelope surrounding the high mass-loss rate OH/IR star OH 127.8+0.0, as traced by Herschel observations, is studied, with a focus on the H2O content and the dust-to-gas ratio. We report detecting a large number of H2O vapor emission lines up to J = 9 in the Herschel data, for which we present the measured line strengths. Methods: The treatments of both gas and dust species are combined using two numerical radiative transfer codes. The method is illustrated for both low and high mass-loss-rate sources. Specifically, we discuss different ways of assessing the dust-to-gas ratio: 1) from the dust thermal emission spectrum and the CO molecular gas line strengths; 2) from the momentum transfer from dust to gas and the measured gas terminal velocity; and 3) from the determination of the required amount of dust to reproduce H2O lines for a given H2O vapor abundance. These three diagnostics probe different zones of the outflow, for the first time allowing an investigation of a possible radial dependence of the dust-to-gas ratio. Results: ... Continued in article.


title: Entanglement entropy and differential entropy for massive flavors

abstract: In this paper we compute the holographic entanglement entropy for massive flavors in the D3-D7 system, for arbitrary mass and various entangling region geometries. We show that the universal terms in the entanglement entropy exactly match those computed in the dual theory using conformal perturbation theory. We derive holographically the universal terms in the entanglement entropy for a CFT perturbed by a relevant operator, up to second order in the coupling; our results are valid for any entangling region geometry. We present a new method for computing the entanglement entropy of any top-down brane probe system using Kaluza-Klein holography and illustrate our results with massive flavors at finite density. Finally we discuss the differential entropy for brane probe systems, emphasising that the differential entropy captures only the effective lower-dimensional Einstein metric rather than the ten-dimensional geometry.


title: Structure-dynamics relation in shaken optical lattices

abstract: Shaken optical lattices permit to coherently modify the tunneling of particles in a controllable manner. We introduce a general relation between the geometry of shaken lattices and their admissible effective dynamics. Using three different examples, we illustrate the symmetries of the emerging tunneling rates. The results provide a clear framework to understand the relation between lattice geometry and accessible dynamics, and a tool to straightforwardly derive truncated effective Hamiltonians on arbitrary geometries.


title: On the renormalization of quasi parton distribution

abstract: Recent developments showed that light-cone parton distributions can be studied by investigating the large momentum limit of the hadronic matrix elements of spacelike correlators, which are known as quasi parton distributions. Like a light-cone parton distribution, a quasi parton distribution also contains ultraviolet divergences and therefore needs renormalization. The renormalization of non-local operators in general is not well understood. However, in the case of quasi quark distribution, the bilinear quark operator with a straight-line gauge link appears to be multiplicatively renormalizable by the quark wave function renormalization in the axial gauge. We first show that the renormalization of the self energy correction to the quasi quark distribution is equivalent to that of the heavy-light quark vector current in heavy quark effective theory at one-loop order. Assuming this equivalence at two-loop order, we then show that the multiplicative renormalizability of the quasi quark distribution is true at two-loop order.


title: Green function for gradient perturbation of unimodal L\'{e}vy processes

abstract: We prove that the Green function of a generator of isotropic unimodal L\'{e}vy processes with the weak lower scaling order bigger than one and the Green function of its gradient perturbations are comparable for bounded smooth open sets if the drift function is from an appropriate Kato class.


title: A Master equation for force distributions in polydisperse frictional particles

abstract: An incremental evolution equation, i.e. a Master equation in statistical mechanics, is introduced for force distributions in polydisperse frictional particle packings. As basic ingredients of the Master equation, the conditional probability distributions of particle overlaps are determined by molecular dynamics simulations. Interestingly, tails of the distributions become much narrower in the case of frictional particles than frictionless particles, implying that correlations of overlaps are strongly reduced by microscopic friction. Comparing different size distributions, we find that the tails are wider for the wider distribution.


title: Path Graphs, Clique Trees, and Flowers

abstract: An \emph{asteroidal triple} is a set of three independent vertices in a graph such that any two vertices in the set are connected by a path which avoids the neighbourhood of the third. A classical result by Lekkerkerker and Boland \cite{6} showed that interval graphs are precisely the chordal graphs that do not have asteroidal triples. Interval graphs are chordal, as are the \emph{directed path graphs} and the \emph{path graphs}. Similar to Lekkerkerker and Boland, Cameron, Ho\'{a}ng, and L\'{e}v\^{e}que \cite{4} gave a characterization of directed path graphs by a "special type" of asteroidal triple, and asked whether or not there was such a characterization for path graphs. We give strong evidence that asteroidal triples alone are insufficient to characterize the family of path graphs, and give a new characterization of path graphs via a forbidden induced subgraph family that we call \emph{sun systems}. Key to our new characterization is the study of \emph{asteroidal sets} in sun systems, which are a natural generalization of asteroidal triples. Our characterization of path graphs by forbidding sun systems also generalizes a characterization of directed path graphs by forbidding odd suns that was given by Chaplick et al.~\cite{9}.


title: Ordered H2O Structures on a Weakly Interacting Surface, A Helium Diffraction Study of H2O/Au(111)

abstract: In this manuscript we report helium atom scattering (HAS) measurements of the structure of the first H2O layer on Au(111). The interaction between H2O and Au(111) is believed to be particularly weak and conflicting evidence from several indirect studies has suggested that water either grows as 3D ice crystals or as an amorphous wetting layer. In contrast, our measurements show that between 110K and 130K, H2O grows as highly commensurate well ordered islands which only partially wet the gold surface. The islands produce a clear (sqrt3Xsqrt3)R30 diffraction pattern and are characterized by a well defined height of ~ 5 Angstrom with respect to the surface gold atoms. These findings provide support for a unique double bilayer model which has recently been suggested for this surface.


title: Damping factors for head-tail modes at strong space charge

abstract: This paper suggests how feedback and Landau damping can be taken into account for transverse oscillations of bunched beam at strong space charge.


title: An analytic recursive method for optimal multiple stopping: Canadization and phase-type fitting

abstract: We study an optimal multiple stopping problem for call-type payoff driven by a spectrally negative Levy process. The stopping times are separated by constant refraction times, and the discount rate can be positive or negative. The computation involves a distribution of the Levy process at a constant horizon and hence the solutions in general cannot be attained analytically. Motivated by the maturity randomization (Canadization) technique by Carr (1998), we approximate the refraction times by independent, identically distributed Erlang random variables. In addition, fitting random jumps to phase-type distributions, our method involves repeated integrations with respect to the resolvent measure written in terms of the scale function of the underlying Levy process. We derive a recursive algorithm to compute the value function in closed form, and sequentially determine the optimal exercise thresholds. A series of numerical examples are provided to compare our analytic formula to results from Monte Carlo simulation.


title: Surpassing the Path-Limited Resolution of a Fourier Transform Spectrometer with Frequency Combs

abstract: Fourier transform spectroscopy based on incoherent light sources is a well-established tool in research fields from molecular spectroscopy and atmospheric monitoring to material science and biophysics. It provides broadband molecular spectra and information about the molecular structure and composition of absorptive media. However, the spectral resolution is fundamentally limited by the maximum delay range ({\Delta}$_{max}$) of the interferometer, so acquisition of high-resolution spectra implies long measurement times and large instrument size. We overcome this limit by combining the Fourier transform spectrometer with an optical frequency comb and measuring the intensities of individual comb lines by precisely matching the {\Delta}$_{max}$ to the comb line spacing. This allows measurements of absorption lines narrower than the nominal (optical path-limited) resolution without ringing effects from the instrumental lineshape and reduces the acquisition time and interferometer length by orders of magnitude.


title: Type Ia supernova Hubble diagram with near-infrared and optical observations

abstract: We main goal of this paper is to test whether the NIR peak magnitudes of SNe Ia could be accurately estimated with only a single observation obtained close to maximum light, provided the time of B band maximum and the optical stretch parameter are known. We obtained multi-epoch UBVRI and single-epoch J and H photometric observations of 16 SNe Ia in the redshift range z=0.037-0.183, doubling the leverage of the current SN Ia NIR Hubble diagram and the number of SNe beyond redshift 0.04. This sample was analyzed together with 102 NIR and 458 optical light curves (LCs) of normal SNe Ia from the literature. The analysis of 45 well-sampled NIR LCs shows that a single template accurately describes them if its time axis is stretched with the optical stretch parameter. This allows us to estimate the NIR peak magnitudes even with one observation obtained within 10 days from B-band maximum. We find that the NIR Hubble residuals show weak correlation with DM_15 and E(B-V), and for the first time we report a possible dependence on the J_max-H_max color. The intrinsic NIR luminosity scatter of SNe Ia is estimated to be around 0.10 mag, which is smaller than what can be derived for a similarly heterogeneous sample at optical wavelengths. In conclusion, we find that SNe Ia are at least as good standard candles in the NIR as in the optical. We showed that it is feasible to extended the NIR SN Ia Hubble diagram to z=0.2 with very modest sampling of the NIR LCs, if complemented by well-sampled optical LCs. Our results suggest that the most efficient way to extend the NIR Hubble diagram to high redshift would be to obtain a single observation close to the NIR maximum. (abridged)


title: SKIRT: the design of a suite of input models for Monte Carlo radiative transfer simulations

abstract: The Monte Carlo method is the most popular technique to perform radiative transfer simulations in a general 3D geometry. The algorithms behind and acceleration techniques for Monte Carlo radiative transfer are discussed extensively in the literature, and many different Monte Carlo codes are publicly available. On the contrary, the design of a suite of components that can be used for the distribution of sources and sinks in radiative transfer codes has received very little attention. The availability of such models, with different degrees of complexity, has many benefits. For example, they can serve as toy models to test new physical ingredients, or as parameterised models for inverse radiative transfer fitting. For 3D Monte Carlo codes, this requires algorithms to efficiently generate random positions from 3D density distributions. We describe the design of a flexible suite of components for the Monte Carlo radiative transfer code SKIRT. The design is based on a combination of basic building blocks (which can be either analytical toy models or numerical models defined on grids or a set of particles) and the extensive use of decorators that combine and alter these building blocks to more complex structures. For a number of decorators, e.g. those that add spiral structure or clumpiness, we provide a detailed description of the algorithms that can be used to generate random positions. Advantages of this decorator-based design include code transparency, the avoidance of code duplication, and an increase in code maintainability. Moreover, since decorators can be chained without problems, very complex models can easily be constructed out of simple building blocks. Finally, based on a number of test simulations, we demonstrate that our design using customised random position generators is superior to a simpler design based on a generic black-box random position generator.


title: Results of a campaign to observe outbursts of the dwarf nova CSS 121005:212625+201948

abstract: A monitoring programme of CSS 121005:212625+201948 covering nearly two observing seasons has shown that it is a typical SU UMa dwarf nova, but it has one of the shortest supercycles of its class, at 66.9(6) d. The superoutbursts are interspersed with 3 to 7 short duration (~2 days) normal outbursts each of which are separated by a mean interval of 11 days, but can be as short as 2 days. The most intensively studied superoutburst was that of 2014 November, which lasted 14 days and had an outburst amplitude of >4.8 magnitudes, reaching magnitude 15.7 at its brightest. Time resolved photometry revealed superhumps with a peak-to-peak amplitude of 0.2 magnitudes, later declining to 0.1 magnitude. The superhump period was Psh = 0.08838(18) d. Time resolved photometry was conducted during several other superoutbursts, which gave broadly similar results.


title: Bayesian Trend Filtering

abstract: We develop a fully Bayesian hierarchical model for trend filtering, itself a new development in nonparametric, univariate regression. The framework more broadly applies to the generalized lasso, but focus is on Bayesian trend filtering. We compare two shrinkage priors, double exponential and generalized double Pareto. A simulation study, comparing Bayesian trend filtering to the original formulation and a number of other popular methods shows our method to improve estimation error while maintaining if not improving coverage probability. Two time series data sets demonstrate Bayesian trend filtering's robustness to possible violations of its assumptions.


title: Estimates on the amplitude of the first Dirichlet eigenvector in discrete frameworks

abstract: Consider a finite absorbing Markov generator, irreducible on the non-absorbing states. Perron-Frobenius theory ensures the existence of a corresponding positive eigenvector $\varphi$. The goal of the paper is to give bounds on the amplitude $\max \varphi/\min\varphi$. Two approaches are proposed: one using a path method and the other one, restricted to the reversible situation, based on spectral estimates. The latter approach is extended to denumerable birth and death processes absorbing at 0 for which infinity is an entrance boundary. The interest of estimating the ratio is the reduction of the quantitative study of convergence to quasi-stationarity to the convergence to equilibrium of related ergodic processes, as seen in [7].


title: A Category Theory of Communication Theory

abstract: A theory of how agents can come to understand a language is presented. If understanding a sentence $\alpha$ is to associate an operator with $\alpha$ that transforms the representational state of the agent as intended by the sender, then coming to know a language involves coming to know the operators that correspond to the meaning of any sentence. This involves a higher order operator that operates on the possible transformations that operate on the representational capacity of the agent. We formalize these constructs using concepts and diagrams analogous to category theory.


title: Effect of O3 on the atmospheric temperature structure of early Mars

abstract: Ozone is an important radiative trace gas in the Earth's atmosphere. The presence of ozone can significantly influence the thermal structure of an atmosphere, and by this e.g. cloud formation. Photochemical studies suggest that ozone can form in carbon dioxide-rich atmospheres. We investigate the effect of ozone on the temperature structure of simulated early Martian atmospheres. With a 1D radiative-convective model, we calculate temperature-pressure profiles for a 1 bar carbon dioxide atmosphere. Ozone profiles are fixed, parameterized profiles. We vary the location of the ozone layer maximum and the concentration at this maximum. The maximum is placed at different pressure levels in the upper and middle atmosphere (1-10 mbar). Results suggest that the impact of ozone on surface temperatures is relatively small. However, the planetary albedo significantly decreases at large ozone concentrations. Throughout the middle and upper atmospheres, temperatures increase upon introducing ozone due to strong UV absorption. This heating of the middle atmosphere strongly reduces the zone of carbon dioxide condensation, hence the potential formation of carbon dioxide clouds. For high ozone concentrations, the formation of carbon dioxide clouds is inhibited in the entire atmosphere. In addition, due to the heating of the middle atmosphere, the cold trap is located at increasingly higher pressures when increasing ozone. This leads to wetter stratospheres hence might increase water loss rates on early Mars. However, increased stratospheric H2O would lead to more HOx, which could efficiently destroy ozone. This result emphasizes the need for consistent climate-chemistry calculations to assess the feedback between temperature structure, water content and ozone chemistry. Furthermore, convection is inhibited at high ozone amounts, leading to a stably stratified atmosphere.


title: Simultaneous X-ray and optical spectroscopy of the Oef supergiant lambda Cep

abstract: Probing the structures of stellar winds is of prime importance for the understanding of massive stars. Based on their optical spectral morphology and variability, the stars of the Oef class have been suggested to feature large-scale structures in their wind. High-resolution X-ray spectroscopy and time-series of X-ray observations of presumably-single O-type stars can help us understand the physics of their stellar winds. We have collected XMM-Newton observations and coordinated optical spectroscopy of the O6Ief star lambda Cep to study its X-ray and optical variability and to analyse its high-resolution X-ray spectrum. We investigate the line profile variability of the He II 4686 and H-alpha emission lines in our time series of optical spectra, including a search for periodicities. We further discuss the variability of the broadband X-ray flux and analyse the high-resolution spectrum of lambda Cep using line-by-line fits as well as a code designed to fit the full high-resolution X-ray spectrum consistently. During our observing campaign, the He II 4686 line varies on a timescale of ~18 hours. On the contrary, the H-alpha line profile displays a modulation on a timescale of 4.1 days which is likely the rotation period of the star. The X-ray flux varies on time-scales of days and could in fact be modulated by the same 4.1 days period as H-alpha, although both variations are shifted in phase. The high-resolution X-ray spectrum reveals broad and skewed emission lines as expected for the X-ray emission from a distribution of wind-embedded shocks. Most of the X-ray emission arises within less than 2R* above the photosphere.


title: The O(3,2) Symmetry derivable from the Poincar\'e Sphere

abstract: Henri Poincar\'e formulated the mathematics of the Lorentz transformations, known as the Poincar\'e group. He also formulated the Poincar\'e sphere for polarization optics. It is noted that his sphere contains the symmetry of the Lorentz group applicable to the momentum-energy four-vector of a particle in the Lorentz-covariant world. Since the particle mass is a Lorentz-invariant quantity, the Lorentz group does not allow its variations. However, the Poincar\'e sphere contains the symmetry corresponding to the mass variation, leading to the $O(3,2)$ symmetry. An illustrative calculation is given.


title: Polly's Polyhedral Scheduling in the Presence of Reductions

abstract: The polyhedral model provides a powerful mathematical abstraction to enable effective optimization of loop nests with respect to a given optimization goal, e.g., exploiting parallelism. Unexploited reduction properties are a frequent reason for polyhedral optimizers to assume parallelism prohibiting dependences. To our knowledge, no polyhedral loop optimizer available in any production compiler provides support for reductions. In this paper, we show that leveraging the parallelism of reductions can lead to a significant performance increase. We give a precise, dependence based, definition of reductions and discuss ways to extend polyhedral optimization to exploit the associativity and commutativity of reduction computations. We have implemented a reduction-enabled scheduling approach in the Polly polyhedral optimizer and evaluate it on the standard Polybench 3.2 benchmark suite. We were able to detect and model all 52 arithmetic reductions and achieve speedups up to 2.21$\times$ on a quad core machine by exploiting the multidimensional reduction in the BiCG benchmark.


title: Exploring multimodal data fusion through joint decompositions with flexible couplings

abstract: A Bayesian framework is proposed to define flexible coupling models for joint tensor decompositions of multiple data sets. Under this framework, a natural formulation of the data fusion problem is to cast it in terms of a joint maximum a posteriori (MAP) estimator. Data driven scenarios of joint posterior distributions are provided, including general Gaussian priors and non Gaussian coupling priors. We present and discuss implementation issues of algorithms used to obtain the joint MAP estimator. We also show how this framework can be adapted to tackle the problem of joint decompositions of large datasets. In the case of a conditional Gaussian coupling with a linear transformation, we give theoretical bounds on the data fusion performance using the Bayesian Cramer-Rao bound. Simulations are reported for hybrid coupling models ranging from simple additive Gaussian models, to Gamma-type models with positive variables and to the coupling of data sets which are inherently of different size due to different resolution of the measurement devices.


title: Comparison of the magneto-Peltier and magneto-Seebeck effects in magnetic tunnel junctions

abstract: Understanding heat generation and transport processes in a magnetic tunnel junction (MTJ) is a significant step towards improving its application in current memory devices. Recent work has experimentally demonstrated the magneto-Seebeck effect in MTJs, where the Seebeck coefficient of the junction varies as the magnetic configuration changes from a parallel (P) to an anti-parallel (AP) configuration. Here we report the study on its as-yet-unexplored reciprocal effect, the magneto-Peltier effect, where the heat flow carried by the tunneling electrons is altered by changing the magnetic configuration of the MTJ. The magneto-Peltier signal that reflects the change in the temperature difference across the junction between the P and AP configurations scales linearly with the applied current in the small bias but is greatly enhanced in the large bias regime, due to higher-order Joule heating mechanisms. By carefully extracting the linear response which reflects the magneto-Peltier effect, and comparing it with the magneto-Seebeck measurements performed on the same device, we observe results consistent with Onsager reciprocity. We estimate a magneto-Peltier coefficient of 13.4 mV in the linear regime using a three-dimensional thermoelectric model. Our result opens up the possibility of programmable thermoelectric devices based on the Peltier effect in MTJs.


title: On entropy production in nonequilibrium systems

abstract: In this paper we discuss the meaning of the Schnakenberg formula for entropy production in non-equilibrium systems. To this end we consider a non-equilibrium system as part of a larger isolated system which includes the environment. We prove that the Schnakenberg formula provides only a lower bound to the actual entropy production in the environment. This is also demonstrated in the simplest example of a three-state clock model.


title: Model structure on differential graded commutative algebras over the ring of differential operators

abstract: We construct a cofibrantly generated model structure on the category of differential non-negatively graded quasi-coherent commutative $D_X$-algebras, where $D_X$ is the sheaf of differential operators of a smooth afine algebraic variety X. The paper contains an extensive appendix on D-modules, sheaves versus global sections, some more technical model categorical issues, as well as on relative Sullivan algebras. This article is the first of a series of works -located at the interface of homotopical algebra, algebraic geometry, and mathematical physics - on a derived D-geometric approach to the BV-formalism.


title: Multiple scattering of elastic waves by pinned dislocation segments in a continuum

abstract: The coherent propagation of elastic waves in a solid filled with a random distribution of pinned dislocation segments is studied to all orders in perturbation theory. It is shown that, within the independent scattering approximation, the perturbation series that generates the mass operator is a geometric series that can thus be formally summed. A divergent quantity is shown to be renormalizable to zero at low frequencies. At higher frequencies said quantity can be expressed in terms of a cut-off with dimensions of length, related to the dislocation length, and physical quantities can be computed in terms of two parameters, to be determined by experiment. The approach used in this problem is compared and contrasted with the scattering of de Broglie waves by delta-function potentials as described by the Schr\"odinger equation.


title: Stratified Patient Appointment Scheduling for Community-based Chronic Disease Management Programs

abstract: Disease management programs have emerged as a cost-effective approach to treat chronic diseases. Appointment adherence is critical to the success of such programs; missed appointment are costly, resulting in reduced resource utilization and worsening of patients' health states. The time of an appointment is one of the factors that impacts adherence. We investigate the benefits, in terms of improved adherence, of incorporating patients' time-of-day preferences during appointment schedule creation and, thus, ultimately, on population health outcomes. Through an extensive computational study, we demonstrate, more generally, the usefulness of patient stratification in appointment scheduling in the environment that motivates our research, an asthma management program offered in Chicago. We find that capturing patient characteristics in appointment scheduling, especially their time preferences, leads to substantial improvements in community health outcomes. We also identify settings in which simple, easy-to-use policies can produce schedules that are comparable in quality to those obtained with an optimization-based approach.


title: Vorticity evolution in a rigid pipe of circular cross-section

abstract: In this paper, we show that the spatio-temporal evolution of incompressible flows in a long circular pipe can be described by vorticity dynamics. The principal techniques to obtain solutions are similar to those used for flows in the whole space. As the consideration of the Navier-Stokes equations is given in a cylindrical co-ordinates system, two aspects of complication arise. One is the interaction of the velocity components in the radial and azimuthal directions, due to the fictitious centrifugal force in the equations of motion. The rate of the vorticity production at the pipe wall depends on the initial data at entry, and hence is unknown a priori; it must be determined as part of the solution. The vorticity solution obtained defines an intricate flow-field of multitudinous degrees of freedom. As the Reynolds number increases, the analytical solution predicts vorticity-scale proliferations in succession. For sufficiently large initial data, pipe flows are of a turbulent nature. The solution of the governing equations is globally regular and does not bifurcate in space or in time. It is asserted that laminar-turbulent transition is a dynamic process inbred in the non-linearity. The presence of exogenous disturbances, due to imperfect test environments or purpose-made artificial forcing, distorts the course of the intrinsic transition. The flow structures observed by Reynolds (1883) and others can be synthesised and elucidated in light of the current theory.


title: Inhomogeneous incompressible viscous flows with slowly varying initial data

abstract: The purpose of this paper is to provide a large class of initial data which generates global smooth solution of the 3-D inhomogeneous incompressible Navier-Stokes system in the whole space~$\R^3$. This class of data is based on functions which vary slowly in one direction. The idea is that 2-D inhomogeneous Navier-Stokes system with large data is globally well-posedness and we construct the 3-D approximate solutions by the 2-D solutions with a parameter. One of the key point of this study is the investigation of the time decay properties of the solutions to the 2-D inhomogeneous Navier-Stokes system. We obtained the same optimal decay estimates as the solutions of 2-D homogeneous Navier-Stokes system.


title: Performance-Complexity Analysis for MAC ML-based Decoding with User Selection

abstract: This work explores the rate-reliability-complexity limits of the quasi-static K-user multiple access channel (MAC), with or without feedback. Using high-SNR asymptotics, the work first derives bounds on the computational resources required to achieve near-optimal (ML-based) decoding performance. It then bounds the (reduced) complexity needed to achieve any (including suboptimal) diversity-multiplexing performance tradeoff (DMT) performance, and finally bounds the same complexity, in the presence of feedback-aided user selection. This latter effort reveals the ability of a few bits of feedback not only to improve performance, but also to reduce complexity. In this context, the analysis reveals the interesting finding that proper calibration of user selection can allow for near-optimal ML-based decoding, with complexity that need not scale exponentially in the total number of codeword bits. The derived bounds constitute the best known performance-vs-complexity behavior to date for ML-based MAC decoding, as well as a first exploration of the complexity-feedback-performance interdependencies in multiuser settings.


title: Linear Codes from a Generic Construction

abstract: A generic construction of linear codes over finite fields has recently received a lot of attention, and many one-weight, two-weight and three-weight codes with good error correcting capability have been produced with this generic approach. The first objective of this paper is to establish relationships among some classes of linear codes obtained with this approach, so that the parameters of some classes of linear codes can be derived from those of other classes with known parameters. In this way, linear codes with new parameters will be derived. The second is to present a class of three-weight binary codes and consider their applications in secret sharing.


title: Chiral degeneracies and Fermi-surface Chern numbers in bcc Fe

abstract: The degeneracies in the spinor bandstructure of bcc Fe are studied from first principles. We find numerous isolated band touchings carrying chiral charges of magnitude one (Weyl points) or two (double-Weyl nodes), as well as nonchiral degeneracy loops (nodal rings). Some degeneracies are located on symmetry lines or planes in the Brillouin zone and others at generic low-symmetry points, realizing all possible scenarios consistent with the magnetic point group. We clarify the general theory relating the chiral band touchings to the Chern numbers of the Fermi sheets enclosing them, and use this approach to determine the Chern numbers on the Fermi surface of bcc Fe. Although most Fermi sheets enclose Weyl nodes, in almost all cases the net enclosed charge vanishes for symmetry reasons, resulting in a vanishing Chern number. The exceptions are two inversion-symmetric electron pockets along the symmetry line $\Delta$ parallel to the magnetization. Each of them surrounds a single Weyl point, leading to Chern numbers of $\pm 1$. These small topological pockets contribute a sizable amount to the nonquantized part of the intrinsic anomalous Hall conductivity, proportional to their reciprocal-space separation. Variation of the Fermi level (or other system parameters) may lead to a touching event between Fermi sheets, accompanied by a transfer of Chern number between them.