title: Soft phonon modes in the vicinity of the structural quantum critical point

abstract: The quasi-skutterudite superconductors $A_3T_4$Sn$_{13}$ ($A$=Sr, Ca; $T$=Ir, Rh, Co) are highly tunable featuring a structural quantum critical point. We construct a temperature-lattice constant phase diagram for these isovalent compounds, establishing Ca$_{3}$Rh$_4$Sn$_{13}$ and Ca$_{3}$Co$_4$Sn$_{13}$ as members close to and far away from the structural quantum critical point, respectively. Deconvolution of the lattice specific heat and the electrical resistivity provide an approximate phonon density of states $F(\omega)$ and the electron-phonon transport coupling function $\alpha_{tr}^2F(\omega)$ for Ca$_{3}$Rh$_4$Sn$_{13}$ and Ca$_{3}$Co$_4$Sn$_{13}$, enabling us to investigate the influence of the structural quantum critical point. Our results support the scenario of phonon softening close to the structural quantum critical point, and explain the enhancement of the coupling strength on approaching structural instability.


title: Transfer of magnetic order and anisotropy through epitaxial integration of 3$d$ and 4$f$ spin systems

abstract: Resonant x-ray scattering at the Dy $M_5$ and Ni $L_3$ absorption edges was used to probe the temperature and magnetic field dependence of magnetic order in epitaxial LaNiO$_3$-DyScO$_3$ superlattices. For superlattices with 2 unit cell thick LaNiO$_3$ layers, a commensurate spiral state develops in the Ni spin system below 100 K. Upon cooling below $T_{ind} = 18$ K, Dy-Ni exchange interactions across the LaNiO$_3$-DyScO$_3$ interfaces induce collinear magnetic order of interfacial Dy moments as well as a reorientation of the Ni spins to a direction dictated by the strong magneto-crystalline anisotropy of Dy. This transition is reversible by an external magnetic field of 3 T. Tailored exchange interactions between rare-earth and transition-metal ions thus open up new perspectives for the manipulation of spin structures in metal-oxide heterostructures and devices.


title: Heavy quark form factors in the large $\beta_0$ limit

abstract: Heavy quark form factors are calculated at $\beta_0 \alpha_s \sim 1$ to all orders in $\alpha_s$ at the first order in $1/\beta_0$. The $n_f^2 \alpha_s^3$ terms in the recent results [arXiv:1611.07535] for the vector form factors are confirmed, and $n_f^{L-1} \alpha_s^L$ terms for higher $L$ are predicted.


title: Anisotropic twicing for single particle reconstruction using autocorrelation analysis

abstract: The missing phase problem in X-ray crystallography is commonly solved using the technique of molecular replacement, which borrows phases from a previously solved homologous structure, and appends them to the measured Fourier magnitudes of the diffraction patterns of the unknown structure. More recently, molecular replacement has been proposed for solving the missing orthogonal matrices problem arising in Kam's autocorrelation analysis for single particle reconstruction using X-ray free electron lasers and cryo-EM. In classical molecular replacement, it is common to estimate the magnitudes of the unknown structure as twice the measured magnitudes minus the magnitudes of the homologous structure, a procedure known as `twicing'. Mathematically, this is equivalent to finding an unbiased estimator for a complex-valued scalar. We generalize this scheme for the case of estimating real or complex valued matrices arising in single particle autocorrelation analysis. We name this approach "Anisotropic Twicing" because unlike the scalar case, the unbiased estimator is not obtained by a simple magnitude isotropic correction. We compare the performance of the least squares, twicing and anisotropic twicing estimators on synthetic and experimental datasets. We demonstrate 3D homology modeling in cryo-EM directly from experimental data without iterative refinement or class averaging, for the first time.


title: Absorption and Emission in a Free Electron Lasers

abstract: Multiphoton processes in undulators with plane polarized magnetic field are considered. It is shown that the use of strong magnetic fields in the undulator, for beams with relatively low energy makes it possible to increase substantially the frequencies of the amplified electromagnetic waves without noticeably decreasing the gain. The absorption, emission probabilities and the gain are calculated.


title: A Flexible Framework for Hypothesis Testing in High-dimensions

abstract: Hypothesis testing in the linear regression model is a fundamental statistical problem. We consider linear regression in the high-dimensional regime where the number of parameters exceeds the number of samples ($p> n$). In order to make informative inference, we assume that the model is approximately sparse, that is the effect of covariates on the response can be well approximated by conditioning on a relatively small number of covariates whose identities are unknown. We develop a framework for testing very general hypotheses regarding the model parameters. Our framework encompasses testing whether the parameter lies in a convex cone, testing the signal strength, and testing arbitrary functionals of the parameter. We show that the proposed procedure controls the type I error, and also analyze the power of the procedure. Our numerical experiments confirm our theoretical findings and demonstrate that we control false positive rate (type I error) near the nominal level, and have high power. By duality between hypotheses testing and confidence intervals, the proposed framework can be used to obtain valid confidence intervals for various functionals of the model parameters. For linear functionals, the length of confidence intervals is shown to be minimax rate optimal.

categories: math.ST cs.LG stat.AP stat.ME stat.ML stat.TH


title: Molecular Cloud Formation Via Thermal Instability of Finite Resistive Viscous Radiating Plasma with Finite Larmor Radius Corrections

abstract: The effect of radiative heat-loss function and finite ion Larmor radius (FLR) corrections on the thermal instability of infinite homogeneous viscous plasma has been investigated incorporating the effects of thermal conductivity and finite electrical resistivity for the formation of a molecular cloud. The general dispersion relation is derived using the normal mode analysis method with the help of relevant linearized perturbation equations of the problem. Furthermore the wave propagation along and perpendicular to the direction of external magnetic field has been discussed. Stability of the medium is discussed by applying Routh Hurwitzs criterion and it is found that thermal instability criterion determines the stability of the medium. We find that the presence of radiative heat-loss function and thermal conductivity modify the fundamental criterion of thermal instability into radiatively driven thermal instability criterion. In longitudinal direction FLR corrections, viscosity, magnetic field and finite resistivity have no effect on thermal instability criterion. The presence of radiative heat-loss function and thermal conductivity modify the fundamental thermal instability criterion into radiatively driven thermal instability criterion. Also the FLR corrections modify the growth rate of the Alfven mode. For transverse wave propagation FLR corrections, radiative heat-loss function, magnetic field and thermal conductivity modify the thermal instability criterion. From the curves it is clear that heat-loss function, FLR corrections and viscosity have stabilizing effect, while finite resistivity has destabilizing effect on the thermal modes. Our results show that the FLR corrections and radiative heat-loss functions affect the evolution of interstellar molecular clouds and star formation.


title: Finite dimensional simple modules of deformed current Lie algebras

abstract: The deformed current Lie algebra was introduced by the author to study the representation theory of cyclotomic q-Schur algebras at q=1. In this paper, we classify finite dimensional simple modules of deformed current Lie algebras.


title: Coefficient Estimates for Certain Subclass of Analytic Functions Defined by Subordination

abstract: In this article we determine the coefficient bounds for functions in certain subclasses of analytic functions defined by subordination which are related to the well-known classes of starlike and convex functions. The main results deal with some open problems proposed by Q.H. Xu et al. [20,21]. An application of Jack lemma for certain subclass of starlike functions has been discussed.


title: Suppression of charge noise using barrier control of a singlet-triplet qubit

abstract: It has been recently demonstrated that a singlet-triplet spin qubit in semiconductor double quantum dots can be controlled by changing the height of the potential barrier between the two dots ("barrier control"), which has led to a considerable reduction of charge noises as compared to the traditional tilt control method. In this paper we show, through a molecular-orbital-theoretic calculation of double quantum dots influenced by a charged impurity, that the relative charge noise for a system under the barrier control not only is smaller than that for the tilt control, but actually decreases as a function of an increasing exchange interaction. This is understood as a combined consequence of the greatly suppressed detuning noise when the two dots are symmetrically operated, as well as an enhancement of the inter-dot hopping energy of an electron when the barrier is lowered which in turn reduces the relative charge noise at large exchange interaction values. We have also studied the response of the qubit to charged impurities at different locations, and found that the improvement of barrier control is least for impurities equidistant from the two dots due to the small detuning noise they cause, but is otherwise significant along other directions.


title: Unitary equivalent classes of one-dimensional quantum walks II

abstract: This study investigated the unitary equivalent classes of one-dimensional quantum walks. We determined the unitary equivalent classes of one-dimensional quantum walks, two-phase quantum walks with one defect, complete two-phase quantum walks, one-dimensional quantum walks with one defect and translation-invariant one-dimensional quantum walks. The unitary equivalent classes of one-dimensional quantum walks with initial states were also considered.


title: Smoothed nonparametric two-sample tests

abstract: We propose new smoothed median and the Wilcoxon's rank sum test. As is pointed out by Maesono et al.(2016), some nonparametric discrete tests have a problem with their significance probability. Because of this problem, the selection of the median and the Wilcoxon's test can be biased too, however, we show new smoothed tests are free from the problem. Significance probabilities and local asymptotic powers of the new tests are studied, and we show that they inherit good properties of the discrete tests.


title: On Improving Deep Reinforcement Learning for POMDPs

abstract: Deep Reinforcement Learning (RL) recently emerged as one of the most competitive approaches for learning in sequential decision making problems with fully observable environments, e.g., computer Go. However, very little work has been done in deep RL to handle partially observable environments. We propose a new architecture called Action-specific Deep Recurrent Q-Network (ADRQN) to enhance learning performance in partially observable domains. Actions are encoded by a fully connected layer and coupled with a convolutional observation to form an action-observation pair. The time series of action-observation pairs are then integrated by an LSTM layer that learns latent states based on which a fully connected layer computes Q-values as in conventional Deep Q-Networks (DQNs). We demonstrate the effectiveness of our new architecture in several partially observable domains, including flickering Atari games.


title: Biases in prime factorizations and Liouville functions for arithmetic progressions

abstract: We introduce a refinement of the classical Liouville function to primes in arithmetic progressions. Using this, we discover new biases in the appearances of primes in a given arithmetic progression in the prime factorizations of integers. For example, we observe that the primes of the form $4k+1$ tend to appear an even number of times in the prime factorization of a given integer, more so than for primes of the form $4k+3$. We are led to consider variants of P\'olya's conjecture, supported by extensive numerical evidence, and its relation to other conjectures.


title: Memory effects in the velocity relaxation process of the dust particle in dusty plasma

abstract: In this paper, by comparing the time scales associated with the velocity relaxation and correlation time of the random force due to dust charge fluctuations, memory effects in the velocity relaxation of an isolated dust particle exposed to the random force due to dust charge fluctuations are considered, and the velocity relaxation process of the dust particle is considered as a non-Markovian stochastic process. Considering memory effects in the velocity relaxation process of the dust particle yields a retarded friction force, which is introduced by a memory kernel in the fractional Langevin equation. The fluctuation-dissipation theorem for the dust grain is derived from this equation. The mean-square displacement and the velocity autocorrelation function of the dust particle are obtained, and their asymptotic behavior, the dust particle temperature due to charge fluctuations, and the diffusion coefficient are studied in the long-time limit. As an interesting feature, it is found that by considering memory effects in the velocity relaxation process of the dust particle, fluctuating force on the dust particle can cause an anomalous diffusion in a dusty plasma. In this case, the mean-square displacement of the dust grain increases slower than linearly with time, and the velocity autocorrelation function decays as a power-law instead of the exponential decay. Finally, in the Markov limit, these results are in good agreement with those obtained from previous works for Markov (memoryless) process of the velocity relaxation.


title: On novel elastic structures inducing plasmonic resonances with finite frequencies and cloaking due to anomalous localized resonances

abstract: This paper is concerned with the theoretical study of plasmonic resonances for linear elasticity governed by the Lam\'e system in $\mathbb{R}^3$, and their application for cloaking due to anomalous localized resonances. We derive a very general and novel class of elastic structures that can induce plasmonic resonances. It is shown that if either one of the two convexity conditions on the Lam\'e parameters is broken, then we can construct certain plasmon structures that induce resonances. This significantly extends the relevant existing studies in the literature where the violation of both convexity conditions is required. Indeed, the existing plasmonic structures are a particular case of the general structures constructed in our study. Furthermore, we consider the plasmonic resonances within the finite frequency regime, and rigorously verify the quasi-static approximation for diametrically small plasmonic inclusions. Finally, as an application of the newly found structures, we construct a plasmonic device of the core-shell-matrix form that can induce cloaking due to anomalous localized resonance in the quasi-static regime, which also includes the existing study as a special case.


title: Exact Algorithms via Multivariate Subroutines

abstract: We consider the family of $\Phi$-Subset problems, where the input consists of an instance $I$ of size $N$ over a universe $U_I$ of size $n$ and the task is to check whether the universe contains a subset with property $\Phi$ (e.g., $\Phi$ could be the property of being a feedback vertex set for the input graph of size at most $k$). Our main tool is a simple randomized algorithm which solves $\Phi$-Subset in time $(1+b-\frac{1}{c})^n N^{O(1)}$, provided that there is an algorithm for the $\Phi$-Extension problem with running time $b^{n-|X|} c^k N^{O(1)}$. Here, the input for $\Phi$-Extension is an instance $I$ of size $N$ over a universe $U_I$ of size $n$, a subset $X\subseteq U_I$, and an integer $k$, and the task is to check whether there is a set $Y$ with $X\subseteq Y \subseteq U_I$ and $|Y\setminus X|\le k$ with property $\Phi$. We derandomize this algorithm at the cost of increasing the running time by a subexponential factor in $n$, and we adapt it to the enumeration setting where we need to enumerate all subsets of the universe with property $\Phi$. This generalizes the results of Fomin et al. [STOC 2016] who proved the case where $b=1$. As case studies, we use these results to design faster deterministic algorithms for: - checking whether a graph has a feedback vertex set of size at most $k$ - enumerating all minimal feedback vertex sets - enumerating all minimal vertex covers of size at most $k$, and - enumerating all minimal 3-hitting sets. We obtain these results by deriving new $b^{n-|X|} c^k N^{O(1)}$-time algorithms for the corresponding $\Phi$-Extension problems (or enumeration variant). In some cases, this is done by adapting the analysis of an existing algorithm, or in other cases by designing a new algorithm. Our analyses are based on Measure and Conquer, but the value to minimize, $1+b-\frac{1}{c}$, is unconventional and requires non-convex optimization.


title: Performance of the ATLAS Track Reconstruction Algorithms in Dense Environments in LHC Run 2

abstract: With the increase in energy of the Large Hadron Collider to a centre-of-mass energy of 13 TeV for Run 2, events with dense environments, such as in the cores of high-energy jets, became a focus for new physics searches as well as measurements of the Standard Model. These environments are characterized by charged-particle separations of the order of the tracking detectors sensor granularity. Basic track quantities are compared between 3.2 fb$^{-1}$ of data collected by the ATLAS experiment and simulation of proton-proton collisions producing high-transverse-momentum jets at a centre-of-mass energy of 13 TeV. The impact of charged-particle separations and multiplicities on the track reconstruction performance is discussed. The efficiency in the cores of jets with transverse momenta between 200 GeV and 1600 GeV is quantified using a novel, data-driven, method. The method uses the energy loss, dE/dx, to identify pixel clusters originating from two charged particles. Of the charged particles creating these clusters, the measured fraction that fail to be reconstructed is $0.061 \pm 0.006 \textrm{(stat.)} \pm 0.014 \textrm{(syst.)}$ and $0.093 \pm 0.017 \textrm{(stat.)}\pm 0.021 \textrm{(syst.)}$ for jet transverse momenta of 200-400 GeV and 1400-1600 GeV, respectively.


title: Magnetic Proximity Effects in Transition-Metal Dichalcogenides: Converting Excitons

abstract: The two-dimensional character and reduced screening in monolayer transition-metal dichalcogenides (TMDs) lead to the ubiquitous formation of robust excitons with binding energies orders of magnitude larger than in bulk semiconductors. Focusing on neutral excitons, bound electron-hole pairs, that dominate the optical response in TMDs, it is shown that they can provide fingerprints for magnetic proximity effects in magnetic heterostructures. These proximity effects cannot be described by the widely used single-particle description, but instead reveal the possibility of a conversion between optically inactive and active excitons by rotating the magnetization of the magnetic substrate. With recent breakthroughs in fabricating Mo- and W-based magnetic TMD-heterostructures, this emergent optical response can be directly tested experimentally.


title: CSM analyses of $e^+e^- \to t\bar t H, t\bar t Z, t\bar b W$

abstract: We study the modifications of the amplitudes and cross sections of several processes, especially $e^+e^- \to t\bar t H, t\bar t Z, t\bar b W$, generated by Higgs boson and top quark compositeness, in particular within the CSM concept. We illustrate the observable differences that may appear between various, CSM conserving or CSM violating, compositeness possibilities.


title: Other Topics You May Also Agree or Disagree: Modeling Inter-Topic Preferences using Tweets and Matrix Factorization

abstract: We present in this paper our approach for modeling inter-topic preferences of Twitter users: for example, those who agree with the Trans-Pacific Partnership (TPP) also agree with free trade. This kind of knowledge is useful not only for stance detection across multiple topics but also for various real-world applications including public opinion surveys, electoral predictions, electoral campaigns, and online debates. In order to extract users' preferences on Twitter, we design linguistic patterns in which people agree and disagree about specific topics (e.g., "A is completely wrong"). By applying these linguistic patterns to a collection of tweets, we extract statements agreeing and disagreeing with various topics. Inspired by previous work on item recommendation, we formalize the task of modeling inter-topic preferences as matrix factorization: representing users' preferences as a user-topic matrix and mapping both users and topics onto a latent feature space that abstracts the preferences. Our experimental results demonstrate both that our proposed approach is useful in predicting missing preferences of users and that the latent vector representations of topics successfully encode inter-topic preferences.


title: Training L1-Regularized Models with Orthant-Wise Passive Descent Algorithms

abstract: The $L_1$-regularized models are widely used for sparse regression or classification tasks. In this paper, we propose the orthant-wise passive descent algorithm (OPDA) for optimizing $L_1$-regularized models, as an improved substitute of proximal algorithms, which are the standard tools for optimizing the models nowadays. OPDA uses a stochastic variance-reduced gradient (SVRG) to initialize the descent direction, then apply a novel alignment operator to encourage each element keeping the same sign after one iteration of update, so the parameter remains in the same orthant as before. It also explicitly suppresses the magnitude of each element to impose sparsity. The quasi-Newton update can be utilized to incorporate curvature information and accelerate the speed. We prove a linear convergence rate for OPDA on general smooth and strongly-convex loss functions. By conducting experiments on $L_1$-regularized logistic regression and convolutional neural networks, we show that OPDA outperforms state-of-the-art stochastic proximal algorithms, implying a wide range of applications in training sparse models.


title: Joint Hybrid Precoder and Combiner Design for mmWave Spatial Multiplexing Transmission

abstract: Millimeter-wave (mmWave) communications have been considered as a key technology for future 5G wireless networks because of the orders-of-magnitude wider bandwidth than current cellular bands. In this paper, we consider the problem of codebook-based joint analog-digital hybrid precoder and combiner design for spatial multiplexing transmission in a mmWave multiple-input multiple-output (MIMO) system. We propose to jointly select analog precoder and combiner pair for each data stream successively aiming at maximizing the channel gain while suppressing the interference between different data streams. After all analog precoder/combiner pairs have been determined, we can obtain the effective baseband channel. Then, the digital precoder and combiner are computed based on the obtained effective baseband channel to further mitigate the interference and maximize the sum-rate. Simulation results demonstrate that our proposed algorithm exhibits prominent advantages in combating interference between different data streams and offer satisfactory performance improvement compared to the existing codebook-based hybrid beamforming schemes.


title: High-Dimensional Variable Selection and Prediction under Competing Risks with Application to SEER-Medicare Linked Data

abstract: Competing risk analysis considers event times due to multiple causes, or of more than one event types. Commonly used regression models for such data include 1) cause-specific hazards model, which focuses on modeling one type of event while acknowledging other event types simultaneously; and 2) subdistribution hazards model, which links the covariate effects directly to the cumulative incidence function. Their use and in particular statistical properties in the presence of high-dimensional predictors are largely unexplored. Motivated by an analysis using the linked SEER-Medicare database for the purposes of predicting cancer versus non-cancer mortality for patients with prostate cancer, we study the accuracy of prediction and variable selection of existing statistical learning methods under both models using extensive simulation experiments, including different approaches to choosing penalty parameters in each method. We then apply the optimal approaches to the analysis of the SEER-Medicare data.


title: An artificial neural network to discover Hypervelocity stars: Candidates in Gaia DR1/TGAS

abstract: The paucity of hypervelocity stars (HVSs) known to date has severely hampered their potential to investigate the stellar population of the Galactic Centre and the Galactic Potential. The first Gaia data release gives an opportunity to increase the current sample. The challenge is the disparity between the expected number of hypervelocity stars and that of bound background stars. We have applied a novel data mining algorithm based on machine learning techniques, an artificial neural network, to the Tycho-Gaia astrometric solution (TGAS) catalogue. With no pre-selection of data, we could exclude immediately $\sim 99 \%$ of the stars in the catalogue and find 80 candidates with more than $90\%$ predicted probability to be HVSs, based only on their position, proper motions, and parallax. We have cross-checked our findings with other spectroscopic surveys, determining radial velocities for 30 and spectroscopic distances for 5 candidates. In addition, follow-up observations have been carried out at the Isaac Newton Telescope for 22 stars, for which we obtained radial velocities and distance estimates. We discover 14 stars with a total velocity in the Galactic rest frame > 400 km/s, and 5 of these have a probability $>50\%$ of being unbound from the Milky Way. Tracing back their orbits in different Galactic potential models we find one possible unbound HVS with velocity $\sim$ 520 km/s, 5 bound HVSs, and, notably, 5 runaway stars with median velocity between 400 and 780 km/s. At the moment, uncertainties in the distance estimates and ages are too large to confirm the nature of our candidates by narrowing down their ejection location, and we wait for future Gaia releases to validate the quality of our sample. This test successfully demonstrates the feasibility of our new data mining routine.


title: Isovector properties of the nuclear energy density functional from the quark-meson coupling model

abstract: Background: The Skyrme energy density functional is widely used in mean-field calculations of nuclear structure and dynamics. However, its reliance on phenomenology may compromise its isovector properties and its performance for exotic nuclear systems. Purpose: This work investigates the possibility of removing some phenomenology from the density functional by drawing on the high-energy degrees-of-freedom of the quark-meson coupling (QMC) model. The QMC model has microscopically derived isovector properties and far fewer adjustable parameters. Method: The parameters of the Skyrme interaction are fixed using the energy density functional of the QMC model, to give the Skyrme-QMC (SQMC) parameterisation. Results: Hartree-Fock-Bogoliubov calculations of the Sn, Pb and $N=126$ chains are reported, in which SQMC performs with an accuracy comparable to modern phenomenological functionals. Conclusions: The isovector terms of the energy density functional are significant for the binding energies of neutron-rich nuclei. The isovector dependence of the nuclear spin-orbit interaction must be taken into account for calculations of r-process nucleosynthesis abundances.


title: Solvability of the heat equation with a nonlinear boundary condition

abstract: We obtain necessary conditions and sufficient conditions for the solvability of the heat equation in a half-space of ${\bf R}^N$ with a nonlinear boundary condition. Furthermore, we study the relationship between the life span of the solution and the behavior of the initial function.


title: Hybrid Precoder and Combiner Design with One-Bit Quantized Phase Shifters in mmWave MIMO Systems

abstract: Analog/digital hybrid precoder and combiner have been widely used in millimeter wave (mmWave) multiple-input multiple-output (MIMO) systems due to its energy-efficient and economic superiorities. Infinite resolution of phase shifters (PSs) for the analog beamformer can achieve very close performance compared to the full-digital scheme but will result in high complexity and intensive power consumption. Thus, more cost effective and energy efficient low resolution PSs are typically used in practical mmWave MIMO systems. In this paper, we consider the joint hybrid precoder and combiner design with one-bit quantized PSs in mmWave MIMO systems. We propose to firstly design the analog precoder and combiner pair for each data stream successively, aiming at conditionally maximizing the spectral efficiency. We present a novel binary analog precoder and combiner optimization algorithm under a Rank-1 approximation of the interference-included equivalent channel with lower than quadratic complexity. Then the digital precoder and combiner are computed based on the obtained baseband effective channel to further enhance the spectral efficiency. Simulation results demonstrate that the proposed algorithm outperforms the existing one-bit PSs based hybrid beamforming scheme.


title: Superscattering pattern shaping for radially anisotropic nanowires

abstract: We achieve efficient shaping of superscattering by radially anisotropic nanowires relying on resonant multipolar interferences. It is shown that the radial anisotropy of refractive index can be employed to resonantly overlap electric and magnetic multipoles of various orders, and as a result effective superscattering with different engineered angular patterns can be obtained. We further demonstrate that such superscattering shaping relying on unusual radial anisotropy parameters can be directly realised with isotropic multi-layered nanowires, which may shed new light to many fundamental researches and various applications related to scattering particles.


title: Local discontinuous Galerkin methods for the time tempered fractional diffusion equation

abstract: In this article, we consider discrete schemes for a fractional diffusion equation involving a tempered fractional derivative in time. We present a semi-discrete scheme by using the local discontinuous Galerkin (LDG) discretization in the spatial variables. We prove that the semi-discrete scheme is unconditionally stable in $L^2$ norm and convergence with optimal convergence rate $\mathcal{O}(h^{k+1})$. We develop a class of fully discrete LDG schemes by combining the weighted and shifted Lubich difference operators with respect to the time variable, and establish the error estimates. Finally, numerical experiments are presented to verify the theoretical results.


title: Secure Precise Wireless Transmission with Random-Subcarrier-Selection-based Directional Modulation Transmit Antenna Array

abstract: In this paper, a practical wireless transmission scheme is proposed to transmit confidential messages to the desired user securely and precisely by the joint use of multiple techniques including artificial noise (AN) projection, phase alignment (PA)/beamforming, and random subcarrier selection (RSCS) based on OFDM, and directional modulation (DM), namely RSCS-OFDM-DM. This RSCS-OFDM-DM scheme provides an extremely low-complexity structures for the transmitter and desired receiver and makes the secure and precise wireless transmission realizable in practice. For illegal eavesdroppers, the receive power of confidential messages is so weak that their receivers cannot intercept these confidential messages successfully once it is corrupted by AN. In such a scheme, the design of phase alignment/beamforming vector and AN projection matrix depend intimately on the desired direction angle and distance. It is particularly noted that the use of RSCS leads to a significant outcome that the receive power of confidential messages mainly concentrates on the small neighboring region around the desired receiver and only small fraction of its power leaks out to the remaining large broad regions. This concept is called secure precise transmission. The probability density function of real-time receive signal-to-interference-and-noise ratio (SINR) is derived. Also, the average SINR and its tight upper bound are attained. The approximate closed-form expression for average secrecy rate is derived by analyzing the first-null positions of SINR and clarifying the wiretap region.


title: Carleson measures, BMO spaces and balayages associated to Schrodinger operators

abstract: Let $\L$ be a Schr\"odinger operator of the form $\L=-\Delta+V$ acting on $L^2(\mathbb R^n)$, $n\geq3$, where the nonnegative potential $V$ belongs to the reverse H\"older class $B_q$ for some $q\geq n.$ Let ${\rm BMO}_{{\mathcal{L}}}(\RR)$ denote the BMO space associated to the Schr\"odinger operator $\L$ on $\RR$. In this article we show that for every $f\in {\rm BMO}_{\mathcal{L}}(\RR)$ with compact support, then there exist $g\in L^{\infty}(\RR)$ and a finite Carleson measure $\mu$ such that $$ f(x)=g(x) + S_{\mu, {\mathcal P}}(x) $$ with $\|g\|_{\infty} +\||\mu\||_{c}\leq C \|f\|_{{\rm BMO}_{\mathcal{L}}(\RR)},$ where $$ S_{\mu, {\mathcal P}}=\int_{{\mathbb R}^{n+1}_+} {\mathcal P}_t(x,y) d\mu(y, t), $$ and ${\mathcal P}_t(x,y)$ is the kernel of the Poisson semigroup $\{e^{-t\sqrt{\L}}\}_{t> 0} $ on $L^2(\mathbb R^n)$. Conversely, if $\mu$ is a Carleson measure, then $S_{\mu, {\mathcal P}}$ belongs to the space ${\rm BMO}_{{\mathcal{L}}}(\RR)$. This extends the result for the classical John--Nirenberg BMO space by Carleson \cite{C} (see also \cite{U,GJ,W}) to the BMO setting associated to Schr\"odinger operators.


title: A Strategy for Dynamic Programs: Start over and Muddle through

abstract: In the setting of DynFO, dynamic programs update the stored result of a query whenever the underlying data changes. This update is expressed in terms of first-order logic. We introduce a strategy for constructing dynamic programs that utilises periodic computation of auxiliary data from scratch and the ability to maintain a query for a limited number of change steps. We show that if some program can maintain a query for log n change steps after an AC$^1$-computable initialisation, it can be maintained by a first-order dynamic program as well, i.e., in DynFO. As an application, it is shown that decision and optimisation problems defined by monadic second-order (MSO) formulas are in DynFO, if only change sequences that produce graphs of bounded treewidth are allowed. To establish this result, a Feferman-Vaught-type composition theorem for MSO is established that might be useful in its own right.


title: Iterative Hybrid Precoder and Combiner Design for mmWave MIMO-OFDM Systems

abstract: This paper investigates the problem of hybrid precoder and combiner design for multiple-input multiple-output (MIMO) orthogonal frequency division multiplexing (OFDM) systems operating in millimeter-wave (mmWave) bands. We propose a novel iterative scheme to design the codebook-based analog precoder and combiner in forward and reverse channels. During each iteration, we apply compressive sensing (CS) technology to efficiently estimate the equivalent MIMO-OFDM mmWave channel. Then, the analog precoder or combiner is obtained based on the orthogonal matching pursuit (OMP) algorithm to alleviate the interference between different data streams as well as maximize the spectral efficiency. The digital precoder and combiner are finally obtained based on the effective baseband channel to further enhance the spectral efficiency. Simulation results demonstrate the proposed iterative hybrid precoder and combiner algorithm has significant performance advantages.


title: A Framework for Algorithm Stability

abstract: We say that an algorithm is stable if small changes in the input result in small changes in the output. This kind of algorithm stability is particularly relevant when analyzing and visualizing time-varying data. Stability in general plays an important role in a wide variety of areas, such as numerical analysis, machine learning, and topology, but is poorly understood in the context of (combinatorial) algorithms. In this paper we present a framework for analyzing the stability of algorithms. We focus in particular on the tradeoff between the stability of an algorithm and the quality of the solution it computes. Our framework allows for three types of stability analysis with increasing degrees of complexity: event stability, topological stability, and Lipschitz stability. We demonstrate the use of our stability framework by applying it to kinetic Euclidean minimum spanning trees.


title: Multiple solutions for resonant problems of the Robin $p$-Laplacian plus an indefinite potential

abstract: We study a nonlinear boundary value problem driven by the $p$-Laplacian plus an indefinite potential with Robin boundary condition. The reaction term is a Carath\'eodory function which is asymptotically resonant at $\pm\infty$ with respect to a nonprincipal Ljusternik-Schnirelmann eigenvalue. Using variational methods, together with Morse theory and truncation-perturbation techniques, we show that the problem has at least three nontrivial smooth solutions, two of which have a fixed sign.


title: A second-order stochastic maximum principle for generalized mean-field control problem

abstract: In this paper, we study the generalized mean-field stochastic control problem when the usual stochastic maximum principle (SMP) is not applicable due to the singularity of the Hamiltonian function. In this case, we derive a second order SMP. We introduce the adjoint process by the generalized mean-field backward stochastic differential equation. The keys in the proofs are the expansion of the cost functional in terms of a perturbation parameter, and the use of the range theorem for vector-valued measures.


title: Hamiltonian Analysis of Non-Relativistic Non-BPS Dp-brane

abstract: We perform Hamiltonian analysis of non-relativistic non-BPS Dp-brane. We find the constraint structure of this theory and determine corresponding equations of motion. We further discuss property of this theory at the tachyon vacuum.


title: Quantum phase fluctuations and density of states in superconducting nanowires

abstract: We argue that quantum fluctuations of the phase of the order parameter may strongly affect the electron density of states (DOS) in ultrathin superconducting wires. We demonstrate that the effect of such fluctuations is equivalent to that of a quantum dissipative environment formed by sound-like plasma modes propagating along the wire. We derive a non-perturbative expression for the local electron DOS in superconducting nanowires which fully accounts for quantum phase fluctuations. At any non-zero temperature these fluctuations smear out the square-root singularity in DOS near the superconducting gap and generate quasiparticle states at subgap energies. Furthermore, at sufficiently large values of the wire impedance this singularity is suppressed down to $T=0$ in which case DOS tends to zero at subgap energies and exhibits the power-law behavior above the gap. Our predictions can be directly tested in tunneling experiments with superconducting nanowires.


title: Toward Reliable Modeling of S-Nitrosothiol Chemistry: Structure and Properties of Methyl Thionitrite (CH$_3$SNO), an S-Nitrosocysteine Model

abstract: Methyl thionitrite CH$_3$SNO is an important model of S-nitrosated cysteine aminoacid residue (CysNO), a ubiquitous biological S-nitrosothiol (RSNO) involved in numerous physiological processes. Here, we report accurate structure and properties of CH$_3$SNO using accurate ab initio Feller-Peterson-Dixon (FPD) approach. The FPD scheme included CCSD(T)-F12/CBS extrapolated values, as well as corrections for the quadruple coupled cluster excitations $\Delta$(Q), core-valence$\Delta$CV and scalar-relativistic $\Delta$SR effects. The FPD scheme for the energetic parameters also included harmonic zero-point vibrational energy (ZPE) corrected for anharmonicity. The S-N bond length in cis-CH$_3$SNO is calculated as 1.814 {\AA}, and its dissociation energy $D_0=32.4$ kcal/mol in the gas phase. The trans-CH$_3$SNO conformation is 1.2 kcal/mol less stable ($\Delta E_0$) compared to cis-CH$_3$SNO, with a sizeable cis-trans isomerization barrier $\Delta E_0^\ne = 12.7$ kcal/mol. The paradox of the unusually long and weak S-N bond, and hindered rotation along the S-N bond, was rationalized via the detailed analysis of the underlying electronic structure of the -SNO group using Natural Resonance Theory (NRT). After the benchmarking of the density functional theory (DFT) methods against the FPD reference, we recommend mPW2PLYP and mPW2PLYPD double hybrid functionals for calculation of the geometric properties, vibrational frequencies and isomerization barriers of S-nitrosothiols, and PBE0 (PBE0-GD3) hybrid functional for the S-N BDEs. The abovementioned DFT methods are capable of capturing the change in electronic structure and properties of the -SNO fragment, when the CH$_3$SNO molecule is exposed to the influence of physiologically feasible external electric field $F_Z$.


title: Deep Text Classification Can be Fooled

abstract: In this paper, we present an effective method to craft text adversarial samples, revealing one important yet underestimated fact that DNN-based text classifiers are also prone to adversarial sample attack. Specifically, confronted with different adversarial scenarios, the text items that are important for classification are identified by computing the cost gradients of the input (white-box attack) or generating a series of occluded test samples (black-box attack). Based on these items, we design three perturbation strategies, namely insertion, modification, and removal, to generate adversarial samples. The experiment results show that the adversarial samples generated by our method can successfully fool both state-of-the-art character-level and word-level DNN-based text classifiers. The adversarial samples can be perturbed to any desirable classes without compromising their utilities. At the same time, the introduced perturbation is difficult to be perceived.


title: Transmit Filter and Artificial Noise Design for Secure MIMO-OFDM Systems

abstract: Physical layer security has been considered as an important security approach in wireless communications to protect legitimate transmission from passive eavesdroppers. This paper investigates the physical layer security of a wireless multiple-input multiple-output (MIMO) orthogonal frequency division multiplexing (OFDM) communication system in the presence of a multiple-antenna eavesdropper. We first propose a transmit-filter-assisted secure MIMO-OFDM system which can destroy the orthogonality of eavesdropper's signals. Our proposed transmit filter can disturb the reception of eavesdropper while maintaining the quality of legitimate transmission. Then, we propose another artificial noise (AN)-assisted secure MIMO-OFDM system to further improve the security of the legitimate transmission. The time-domain AN signal is designed to disturb the reception of eavesdropper while the legitimate transmission will not be affected. Simulation results are presented to demonstrate the security performance of the proposed transmit filter design and AN-assisted scheme in the MIMO-OFDM system.


title: Phononic frequency comb via three-mode parametric three-wave mixing

abstract: This paper is motivated by the recent demonstration of three-wave mixing based phononic frequency comb. While the previous experiments have shown the existence of three-wave mixing pathway in a system of two-coupled phonon modes, this work demonstrates a similar pathway in a system of three-coupled phonon modes. The paper also presents a number of interesting experimental facts concomitant to the three-mode three-wave mixing based frequency comb observed in a specific micromechanical device. The experimental validation of three-mode three-wave mixing along with the previous demonstration of two-mode three-wave mixing points to the ultimate possibility of multimode frequency combs.


title: Existence theorems for a nonlinear second-order distributional differential equation

abstract: In this work, we are concerned with existence of solutions for a nonlinear second-order distributional differential equation, which contains measure differential equations and stochastic differential equations as special cases. The proof is based on the Leray--Schauder nonlinear alternative and Kurzweil--Henstock--Stieltjes integrals. Meanwhile, examples are worked out to demonstrate that the main results are sharp.


title: Linear foliations of complex spheres I. Chains

abstract: We provide coordinate-free versions of the classical projection Theorem of Marstrand-Kaufman-Mattila. This allows us to generalize this Theorem to the complex setting; in restriction to complex spheres, we obtain further projection Theorems along so-called complex chains.


title: Tsallis entropy and generalized Shannon additivity

abstract: The Tsallis entropy given for a positive parameter $\alpha$ can be considered as a modification of the classical Shannon entropy. For the latter, corresponding to $\alpha=1$, there exist many axiomatic characterizations. One of them based on the well-known Khinchin-Shannon axioms has been simplified several times and adapted to Tsallis entropy, where the axiom of (generalized) Shannon additivity is playing a central role. The main aim of this paper is to discuss this axiom in the context of Tsallis entropy. We show that it is sufficient for characterizing Tsallis entropy with the exceptions of cases $\alpha=1,2$ discussed separately.


title: Topically Driven Neural Language Model

abstract: Language models are typically applied at the sentence level, without access to the broader document context. We present a neural language model that incorporates document context in the form of a topic model-like architecture, thus providing a succinct representation of the broader document context outside of the current sentence. Experiments over a range of datasets demonstrate that our model outperforms a pure sentence-based model in terms of language model perplexity, and leads to topics that are potentially more coherent than those produced by a standard LDA topic model. Our model also has the ability to generate related sentences for a topic, providing another way to interpret topics.


title: Replica Symmetry Breaking in Compressive Sensing

abstract: For noisy compressive sensing systems, the asymptotic distortion with respect to an arbitrary distortion function is determined when a general class of least-square based reconstruction schemes is employed. The sampling matrix is considered to belong to a large ensemble of random matrices including i.i.d. and projector matrices, and the source vector is assumed to be i.i.d. with a desired distribution. We take a statistical mechanical approach by representing the asymptotic distortion as a macroscopic parameter of a spin glass and employing the replica method for the large-system analysis. In contrast to earlier studies, we evaluate the general replica ansatz which includes the RS ansatz as well as RSB. The generality of the solution enables us to study the impact of symmetry breaking. Our numerical investigations depict that for the reconstruction scheme with the "zero-norm" penalty function, the RS fails to predict the asymptotic distortion for relatively large compression rates; however, the one-step RSB ansatz gives a valid prediction of the performance within a larger regime of compression rates.


title: Self-assembly of InAs nanostructures on the sidewalls of GaAs nanowires directed by a Bi surfactant

abstract: Surface energies play a dominant role in the self-assembly of three dimensional (3D) nanostructures. In this letter, we show that using surfactants to modify surface energies can provide a means to externally control nanostructure self-assembly, enabling the synthesis of novel hierarchical nanostructures. We explore Bi as a surfactant in the growth of InAs on the {1-10} sidewall facets of GaAs nanowires. The presence of surface Bi induces the formation of InAs 3D islands by a process resembling the Stranski-Krastanov mechanism, which does not occur in the absence of Bi on these surfaces. The InAs 3D islands nucleate at the corners of the {1-10} facets above a critical shell thickness and then elongate along <110> directions in the plane of the nanowire sidewalls. Exploiting this growth mechanism, we realize a series of novel hierarchical nanostructures, ranging from InAs quantum dots on single {1-10} nanowire facets to zig-zag shaped nanorings completely encircling nanowire cores. Photoluminescence spectroscopy and cathodoluminescence spectral line scans reveal that small surfactant-induced InAs 3D islands behave as optically active quantum dots. This work illustrates how surfactants can provide an unprecedented level of external control over nanostructure self-assembly.


title: A new method of joint nonparametric estimation of probability density and its support

abstract: In this paper we propose a new method of joint nonparametric estimation of probability density and its support. As is well known, nonparametric kernel density estimator has "boundary bias problem" when the support of the population density is not the whole real line. To avoid the unknown boundary effects, our estimator detects the boundary, and eliminates the boundary-bias of the estimator simultaneously. Moreover, we refer an extension to a simple multivariate case, and propose an improved estimator free from the unknown boundary bias.


title: Optimal location of resources for biased movement of species: the 1D case

abstract: In this paper, we investigate an optimal design problem motivated by some issues arising in population dynamics. In a nutshell, we aim at determining the optimal shape of a region occupied by resources for maximizing the survival ability of a species in a given box and we consider the general case of Robin boundary conditions on its boundary. Mathematically, this issue can be modeled with the help of an extremal indefinite weight linear eigenvalue problem. The optimal spatial arrangement is obtained by minimizing the positive principal eigenvalue with respect to the weight, under a L 1 constraint standing for limitation of the total amount of resources. The specificity of such a problem rests upon the presence of nonlinear functions of the weight both in the numerator and denominator of the Rayleigh quotient. By using adapted symmetrization procedures, a well-chosen change of variable, as well as necessary optimality conditions, we completely solve this optimization problem in the unidimensional case by showing first that every minimizer is unimodal and bang-bang. This leads to investigate a finite dimensional optimization problem. This allows to show in particular that every minimizer is (up to additive constants) the characteristic function of three possible domains: an interval that sticks on the boundary of the box, an interval that is symmetrically located at the middle of the box, or, for a precise value of the Robin coefficient, all intervals of a given fixed length.


title: Bohmian Mechanics

abstract: Bohmian mechanics, also known as pilot-wave theory or de Broglie-Bohm theory, is a formulation of quantum mechanics whose fundamental axioms are not about what observers will see if they perform an experiment but about what happens in reality. It is therefore called a "quantum theory without observers," alongside with collapse theories and many-worlds theories and in contrast to orthodox quantum mechanics. It follows from these axioms that in a universe governed by Bohmian mechanics, observers will see outcomes with exactly the probabilities specified by the usual rules of quantum mechanics for empirical predictions. Specifically, Bohmian mechanics asserts that electrons and other elementary particles have a definite position at every time and move according to an equation of motion that is one of the fundamental laws of the theory and involves a wave function that evolves according to the usual Schr\"odinger equation. Bohmian mechanics is named after David Bohm (1917-1992), who was, although not the first to consider this theory, the first to realize (in 1952) that it actually makes correct predictions.


title: Sunspots rotation and magnetic transients associated with flares in NOAA AR 11429

abstract: We analyze sunspots rotation and magnetic transients in NOAA AR 11429 during two X-class (X5.4 and X1.3) flares using the data from the Helioseismic and Magnetic Imager on board the \emph{Solar Dynamics Observatory}. A large leading sunspot with positive magnetic polarity rotated counterclockwise. As expected, the rotation was significantly affected by the two flares. The magnetic transients induced by the flares were clearly evident in the sunspots with negative polarity. They were moving across the sunspots with speed of order $3-7\ \rm km \ s^{-1}$. Furthermore, the trend of magnetic flux evolution of these sunspots exhibited changes associated with the flares. These results may shed light on the understanding of the evolution of sunspots.


title: Structure and energetics of embedded Si patterns in graphene

abstract: Recent experiments have revealed the possibility of precise electron beam manipulation of silicon impurities in graphene. Motivated by these findings and studies on metal surface quantum corrals, the question arises what kind of embedded Si structures are possible within the hexagonal lattice, and how these are limited by the distortion caused by the preference of Si for $sp^{3}$ hybridization. In this work, we study the geometry and stability of elementary Si patterns in graphene, including lines, hexagons, triangles, circles and squares. Due to the size of the required unit cells, to obtain the relaxed geometries we use an empirical bond-order potential as a starting point for density functional theory. Despite some interesting discrepancies, the classical geometries provide an effective route for the simulation of large structures.


title: Solar wind collisional heating

abstract: To properly describe heating in weakly collisional turbulent plasmas such as the solar wind, inter-particle collisions should be taken into account. Collisions can convert ordered energy into heat by means of irreversible relaxation towards the thermal equilibrium. Recently, Pezzi et al. (Phys. Rev. Lett., vol. 116, 2016, p. 145001) showed that the plasma collisionality is enhanced by the presence of fine structures in velocity space. Here, the analysis is extended by directly comparing the effects of the fully nonlinear Landau operator and a linearized Landau operator. By focusing on the relaxation towards the equilibrium of an out of equilibrium distribution function in a homogeneous force-free plasma, here it is pointed out that it is significant to retain nonlinearities in the collisional operator to quantify the importance of collisional effects. Although the presence of several characteristic times associated with the dissipation of different phase space structures is recovered in both the cases of the nonlinear and the linearized operators, the influence of these times is different in the two cases. In the linearized operator case, the recovered characteristic times are systematically larger than in the fully nonlinear operator case, this suggesting that fine velocity structures are dissipated slower if nonlinearities are neglected in the collisional operator.


title: Measurement Matrix Design for Phase Retrieval Based on Mutual Information

abstract: In phase retrieval problems, a signal of interest (SOI) is reconstructed based on the magnitude of a linear transformation of the SOI observed with additive noise. The linear transform is typically referred to as a measurement matrix. Many works on phase retrieval assume that the measurement matrix is a random Gaussian matrix, which, in the noiseless scenario with sufficiently many measurements, guarantees invertability of the transformation between the SOI and the observations, up to an inherent phase ambiguity. However, in many practical applications, the measurement matrix corresponds to an underlying physical setup, and is therefore deterministic, possibly with structural constraints. In this work we study the design of deterministic measurement matrices, based on maximizing the mutual information between the SOI and the observations. We characterize necessary conditions for the optimality of a measurement matrix, and analytically obtain the optimal matrix in the low signal-to-noise ratio regime. Practical methods for designing general measurement matrices and masked Fourier measurements are proposed. Simulation tests demonstrate the performance gain achieved by the proposed techniques compared to random Gaussian measurements for various phase recovery algorithms.


title: Injectivity almost everywhere and mappings with finite distortion in nonlinear elasticity

abstract: We show that a sufficient condition for the weak limit of a sequence of $W^1_q$-homeomorphisms with finite distortion to be almost everywhere injective for $q \geq n-1$, can be stated by means of composition operators. Applying this result, we study nonlinear elasticity problems with respect to these new classes of mappings. Furthermore, we impose loose growth conditions on the stored-energy function for the class of $W^1_n$-homeomorphisms with finite distortion and integrable inner as well as outer distortion coefficients.


title: Quasistationary solutions of scalar fields around collapsing self-interacting boson stars

abstract: There is increasing numerical evidence that scalar fields can form long-lived quasi-bound states around black holes. Recent perturbative and numerical relativity calculations have provided further confirmation in a variety of physical systems, including both static and accreting black holes, and collapsing fermionic stars. In this work we investigate this issue yet again in the context of gravitationally unstable boson stars leading to black hole formation. We build a large sample of spherically symmetric initial models, both stable and unstable, incorporating a self-interaction potential with a quartic term. The three different outcomes of unstable models, namely migration to the stable branch, total dispersion, and collapse to a black hole, are also present for self-interacting boson stars. Our simulations show that for black-hole-forming models, a scalar-field remnant is found outside the black-hole horizon, oscillating at a different frequency than that of the original boson star. This result is in good agreement with recent spherically symmetric simulations of unstable Proca stars collapsing to black holes.


title: On the modelling of shallow turbidity flows

abstract: In this study we investigate shallow turbidity density currents and underflows from mechanical point of view. We propose a simple hyperbolic model for such flows. On one hand, our model is based on very basic conservation principles. On the other hand, the turbulent nature of the flow is also taken into account through the energy dissipation mechanism. Moreover, the mixing with the pure water along with sediments entrainment and deposition processes are considered, which makes the problem dynamically interesting. One of the main advantages of our model is that it requires the specification of only two modeling parameters - the rate of turbulent dissipation and the rate of the pure water entrainment. Consequently, the resulting model turns out to be very simple and self-consistent. This model is validated against several experimental data and several special classes of solutions (such as travelling, self-similar and steady) are constructed. Unsteady simulations show that some special solutions are realized as asymptotic long time states of dynamic trajectories.


title: Paths to collapse for isolated skyrmions in few-monolayer ferromagnetic films

abstract: Magnetic skyrmions are topological spin configurations in materials with chiral Dzyaloshinskii-Moriya interaction (DMI), that are potentially useful for storing or processing information. To date, DMI has been found in few bulk materials, but can also be induced in atomically thin magnetic films in contact with surfaces with large spin-orbit interactions. Recent experiments have reported that isolated magnetic skyrmions can be stabilized even near room temperature in few-atom thick magnetic layers sandwiched between materials that provide asymmetric spin-orbit coupling. Here we present the minimum-energy path analysis of three distinct mechanisms for the skyrmion collapse, based on ab initio input and the performed atomic-spin simulations. We focus on the stability of a skyrmion in three atomic layers of Co, either epitaxial on the Pt(111) surface, or within a hybrid multilayer where DMI nontrivially varies per monolayer due to competition between different symmetry-breaking from two sides of the Co film. In laterally finite systems, their constrained geometry causes poor thermal stability of the skyrmion toward collapse at the boundary, which we show to be resolved by designing the high-DMI structure within an extended film with lower or no DMI.


title: Lasing on nonlinear localized waves in curved geometry

abstract: The use of geometrical constraints opens many new perspectives in photonics and in fundamental studies of nonlinear waves. By implementing surface structures in vertical cavity surface emitting lasers as manifolds for curved space, we experimentally study the impacts of geometrical constraints on nonlinear wave localization. We observe localized waves pinned to the maximal curvature in an elliptical-ring, and confirm the reduction in the localization length of waves by measuring near and far field patterns, as well as the corresponding dispersion relation. Theoretically, analyses based on a dissipative model with a parabola curve give good agreement remarkably to experimental measurement on the transition from delocalized to localized waves. The introduction of curved geometry allows to control and design lasing modes in the nonlinear regime.


title: Ranking in evolving complex networks

abstract: Complex networks have emerged as a simple yet powerful framework to represent and analyze a wide range of complex systems. The problem of ranking the nodes and the edges in complex networks is critical for a broad range of real-world problems because it affects how we access online information and products, how success and talent are evaluated in human activities, and how scarce resources are allocated by companies and policymakers, among others. This calls for a deep understanding of how existing ranking algorithms perform, and which are their possible biases that may impair their effectiveness. Well-established ranking algorithms (such as the popular Google's PageRank) are static in nature and, as a consequence, they exhibit important shortcomings when applied to real networks that rapidly evolve in time. The recent advances in the understanding and modeling of evolving networks have enabled the development of a wide and diverse range of ranking algorithms that take the temporal dimension into account. The aim of this review is to survey the existing ranking algorithms, both static and time-aware, and their applications to evolving networks. We emphasize both the impact of network evolution on well-established static algorithms and the benefits from including the temporal dimension for tasks such as prediction of real network traffic, prediction of future links, and identification of highly-significant nodes.


title: Towards Estimating the Upper Bound of Visual-Speech Recognition: The Visual Lip-Reading Feasibility Database

abstract: Speech is the most used communication method between humans and it involves the perception of auditory and visual channels. Automatic speech recognition focuses on interpreting the audio signals, although the video can provide information that is complementary to the audio. Exploiting the visual information, however, has proven challenging. On one hand, researchers have reported that the mapping between phonemes and visemes (visual units) is one-to-many because there are phonemes which are visually similar and indistinguishable between them. On the other hand, it is known that some people are very good lip-readers (e.g: deaf people). We study the limit of visual only speech recognition in controlled conditions. With this goal, we designed a new database in which the speakers are aware of being read and aim to facilitate lip-reading. In the literature, there are discrepancies on whether hearing-impaired people are better lip-readers than normal-hearing people. Then, we analyze if there are differences between the lip-reading abilities of 9 hearing-impaired and 15 normal-hearing people. Finally, human abilities are compared with the performance of a visual automatic speech recognition system. In our tests, hearing-impaired participants outperformed the normal-hearing participants but without reaching statistical significance. Human observers were able to decode 44% of the spoken message. In contrast, the visual only automatic system achieved 20% of word recognition rate. However, if we repeat the comparison in terms of phonemes both obtained very similar recognition rates, just above 50%. This suggests that the gap between human lip-reading and automatic speech-reading might be more related to the use of context than to the ability to interpret mouth appearance.


title: Multi-frequency search for Dark Matter: the role of HESS, CTA, and SKA

abstract: Dark Matter (DM) remains a vital, but elusive, component in our current understanding of the universe. Accordingly, many experimental searches are devoted to uncovering its nature. However, both the existing direct detection methods, and the prominent $\gamma$-ray search with the Fermi Large Area Telescope (Fermi-LAT), are most sensitive to DM particles with masses below 1 TeV, and are significantly less sensitive to the hard spectra produced in annihilation via heavy leptons. The High Energy Stereoscopic System (HESS) has had some success in improving on the Fermi-LAT search for higher mass DM particles, particularly annihilating via heavy lepton states. However, the recent discovery of high J-factor dwarf spheroidal galaxies by the Dark Energy Survey (DES) opens up the possibility of investing more HESS observation time in the search for DM $\gamma$-ray signatures in dwarf galaxies. This work explores the potential of HESS to extend its current limits using these new targets, as well as the future constraints derivable with the up-coming Cherenkov Telescope Array (CTA). These limits are further compared with those we derived at low radio frequencies for the Square Kilometre Array (SKA). Finally, we explore the impact of HESS, CTA, and Fermi-LAT on the phenomenology of the "Madala" boson hypothesized based on anomalies in the data from the Large Hadron Collider (LHC) run 1. The power of these limits from differing frequency bands is suggestive of a highly effective multi-frequency DM hunt strategy making use of both existing and up-coming Southern African telescopes.


title: Airway segmentation from 3D chest CT volumes based on volume of interest using gradient vector flow

abstract: Some lung diseases are related to bronchial airway structures and morphology. Although airway segmentation from chest CT volumes is an important task in the computer-aided diagnosis and surgery assistance systems for the chest, complete 3-D airway structure segmentation is a quite challenging task due to its complex tree-like structure. In this paper, we propose a new airway segmentation method from 3D chest CT volumes based on volume of interests (VOI) using gradient vector flow (GVF). This method segments the bronchial regions by applying the cavity enhancement filter (CEF) to trace the bronchial tree structure from the trachea. It uses the CEF in the VOI to segment each branch. And a tube-likeness function based on GVF and the GVF magnitude map in each VOI are utilized to assist predicting the positions and directions of child branches. By calculating the tube-likeness function based on GVF and the GVF magnitude map, the airway-like candidate structures are identified and their centrelines are extracted. Based on the extracted centrelines, we can detect the branch points of the bifurcations and directions of the airway branches in the next level. At the same time, a leakage detection is performed to avoid the leakage by analysing the pixel information and the shape information of airway candidate regions extracted in the VOI. Finally, we unify all of the extracted bronchial regions to form an integrated airway tree. Preliminary experiments using four cases of chest CT volumes demonstrated that the proposed method can extract more bronchial branches in comparison with other methods.


title: What Can Gamma-rays from Space tell us About the Madala Hypothesis?

abstract: The recent Madala hypothesis, a conjecture that seeks to explain anomalies within Large Hadron Collider (LHC) data (particularly in the transverse momentum of the Higgs boson), is interesting for more than just a statistical hint at unknown and unpredicted physics. This is because the model itself contains additional new particles that may serve as Dark Matter (DM) candidates. These particles interact with the Standard Model via a scalar mediator boson $S$. More interesting still, the conjectured mass range for the DM candidate ($65$ - $100$ GeV) lies within the region of models viable to try explain the recent Galactic Centre (GC) gamma-ray excess seen by Fermi Large Area Telescope (Fermi-LAT) and the High Energy Stereoscopic System (HESS). Therefore, assuming $S$ decays promptly, it should be possible to check what constraints are imposed upon the effective DM annihilation cross-section in the Madala scenario by hunting signatures of $S$ decay that follows DM annihilation within dense astrophysical structures. In order to make use of existing data, we use the Reticulum II dwarf galaxy and the galactic centre gamma-ray excess data sets from Fermi-LAT, and compare these to the consequences of various decay paths for $S$ in the aforementioned environments. We find that, based on this existing data, we can limit $\tau$ lepton, quark, direct gamma-ray, and weak boson channels to levels below the canonical relic cross-section. This allows us to set new limits on the branching ratios of $S$ decay, which can rule out a Higgs-like decay branching for $S$, in the case where the Madala DM candidate is assumed to comprise all DM.


title: Reconstructing Detailed Line Profiles of Lamellar Gratings from GISAXS Patterns with a Maxwell Solver

abstract: Laterally periodic nanostructures were investigated with grazing incidence small angle X-ray scattering (GISAXS) by using the diffraction patterns to reconstruct the surface shape. To model visible light scattering, rigorous calculations of the near and far field by numerically solving Maxwell's equations with a finite-element method are well established. The application of this technique to X-rays is still challenging, due to the discrepancy between incident wavelength and finite-element size. This drawback vanishes for GISAXS due to the small angles of incidence, the conical scattering geometry and the periodicity of the surface structures, which allows a rigorous computation of the diffraction efficiencies with sufficient numerical precision. To develop dimensional metrology tools based on GISAXS, lamellar gratings with line widths down to 55 nm were produced by state-of-the-art e-beam lithography and then etched into silicon. The high surface sensitivity of GISAXS in conjunction with a Maxwell solver allows a detailed reconstruction of the grating line shape also for thick, non-homogeneous substrates. The reconstructed geometrical line shape models are statistically validated by applying a Markov chain Monte Carlo (MCMC) sampling technique which reveals that GISAXS is able to reconstruct critical parameters like the widths of the lines with sub-nm uncertainty.


title: Observation of $e^{+}e^{-} \to \eta h_{c}$ at center-of-mass energies from 4.085 to 4.600 GeV

abstract: We observe for the first time the process $e^{+}e^{-} \rightarrow \eta h_c$ with data collected by the BESIII experiment. Significant signals are observed at the center-of-mass energy $\sqrt{s}=4.226$ GeV, and the Born cross section is measured to be $(9.5^{+2.2}_{-2.0} \pm 2.7)$ pb. Evidence for $\eta h_c$ is observed at $\sqrt{s}=4.358$ GeV with a Born cross section of $(10.0^{+3.1}_{-2.7} \pm 2.6)$ pb, and upper limits on the production cross section at other center-of-mass energies between 4.085 and 4.600 GeV are determined.


title: Optimal Decentralized Economical-sharing Scheme in Islanded AC Microgrids with Cascaded Inverters

abstract: To address the economical dispatch problem without communications in islanded AC microgrids consisting of cascaded inverters, this paper proposes an optimal decentralized economical-sharing scheme. In proposed scheme, optimal sharing function of the current is applied to generate the reference voltages. And the frequency is used to drive all distributed generators (DGs) synchronize operation in microgrids. When the microgrid is in steady state, DGs share a single common frequency and current in terms of the proposed scheme. Thus the potential advantages of simplicity and decentralized manner are retained. The AC microgrid model has been developed through simulations and experiments to verify the effectiveness and performance of the proposed scheme.


title: Automatic Viseme Vocabulary Construction to Enhance Continuous Lip-reading

abstract: Speech is the most common communication method between humans and involves the perception of both auditory and visual channels. Automatic speech recognition focuses on interpreting the audio signals, but it has been demonstrated that video can provide information that is complementary to the audio. Thus, the study of automatic lip-reading is important and is still an open problem. One of the key challenges is the definition of the visual elementary units (the visemes) and their vocabulary. Many researchers have analyzed the importance of the phoneme to viseme mapping and have proposed viseme vocabularies with lengths between 11 and 15 visemes. These viseme vocabularies have usually been manually defined by their linguistic properties and in some cases using decision trees or clustering techniques. In this work, we focus on the automatic construction of an optimal viseme vocabulary based on the association of phonemes with similar appearance. To this end, we construct an automatic system that uses local appearance descriptors to extract the main characteristics of the mouth region and HMMs to model the statistic relations of both viseme and phoneme sequences. To compare the performance of the system different descriptors (PCA, DCT and SIFT) are analyzed. We test our system in a Spanish corpus of continuous speech. Our results indicate that we are able to recognize approximately 58% of the visemes, 47% of the phonemes and 23% of the words in a continuous speech scenario and that the optimal viseme vocabulary for Spanish is composed by 20 visemes.


title: Spatially-Flat Robertson-Walker models with combined CDM and stiff matter sources and the corresponding thermodynamics

abstract: In the present paper, we are considering a spatially-flat Friedmann-Robertson-Walker cosmological model, fueled with stiff matter and dust, treated as non-interacting ideal fluid sources. By solving the corresponding Friedmann equation with a non-zero cosmological constant, we are deriving the scale function and the fundamental cosmological parameters. Within a thermodynamic approach, the general form of the Equation of State is obtained, together with the explicit dependence of the energy density and pressure on temperature.


title: Filmor Theorem for integers

abstract: Fillmore Theorem says that if $A$ is a nonscalar matrix of order $n$ over a field $\mathbb{F}$ and $\gamma_1,\ldots,\gamma_n\in \mathbb{F}$ are such that $\gamma_1+\cdots+\gamma_n=\text{tr} \, A$, then there is a matrix $B$ similar to $A$ with diagonal $(\gamma_1,\ldots,\gamma_n)$. Fillmore proof works by induction on the size of $A$ and implicitly provides an algorithm to construct $B$. We develop an explicit and extremely simple algorithm that finish in two steps (two similarities), and with its help we extend Fillmore Theorem to integers (if $A$ is integer then we can require to $B$ to be integer).


title: Cleaning graphene: comparing heat treatments in air and in vacuum

abstract: Surface impurities and contamination often seriously degrade the properties of two-dimensional materials such as graphene. To remove contamination, thermal annealing is commonly used. We present a comparative analysis of annealing treatments in air and in vacuum, both ex situ and "pre-situ", where an ultra-high vacuum treatment chamber is directly connected to an aberration-corrected scanning transmission electron microscope. While ex situ treatments do remove contamination, it is challenging to obtain atomically clean surfaces after ambient transfer. However, pre-situ cleaning with radiative or laser heating appears reliable and well suited to clean graphene without undue damage to its structure.


title: Bose-Einstein condensation and liquid-gas phase transition in alpha-matter

abstract: Systems of Bose particles with both repulsive and attractive interactions are studied using the Skyrme-like mean-field model. The phase diagram of such systems exhibits two special lines in the chemical potential-temperature plane: one line which represents the first-order liquid-gas phase transition with the critical end point, and another line which represents the onset of Bose-Einstein condensation. The calculations are made for strongly-interacting matter composed of alpha particles. The phase diagram of this matter is qualitatively similar to that observed for the atomic He4 liquid. The sensitivity of the results to the model parameters is studied. For weak interaction coupling the critical point is located at the Bose-condensation line.


title: Graviton fluctuations erase the cosmological constant

abstract: Graviton fluctuations induce strong non-perturbative infrared renormalization effects for the cosmological constant. The functional renormalization flow drives a positive cosmological constant towards zero, solving the cosmological constant problem without the need to tune parameters. We propose a simple computation of the graviton contribution to the flow of the effective potential for scalar fields. Within variable gravity we find that the potential increases asymptotically at most quadratically with the scalar field. With effective Planck mass proportional to the scalar field, the solutions of the derived cosmological equations lead to an asymptotically vanishing cosmological "constant" in the infinite future, providing for dynamical dark energy in the present cosmological epoch. Beyond a solution of the cosmological constant problem, our simplified computation also entails a sizeable positive graviton-induced anomalous dimension for the quartic Higgs coupling in the ultraviolet regime, substantiating the successful prediction of the Higgs boson mass within the asymptotic safety scenario for quantum gravity.


title: Implementation of BVD (boundary variation diminishing) algorithm in simulations of compressible multiphase flows

abstract: We present in this work a new reconstruction scheme, so-called MUSCL-THINC-BVD scheme, to solve the five-equation model for interfacial two phase flows. This scheme employs the traditional shock capturing MUSCL (Monotone Upstream-centered Schemes for Conservation Law) scheme as well as the interface sharpening THINC (Tangent of Hyperbola for INterface Capturing) scheme as two building-blocks of spatial reconstruction using the BVD (boundary variation diminishing) principle that minimizes the variations (jumps) of the reconstructed variables at cell boundaries, and thus effectively reduces the numerical dissipations in numerical solutions. The MUSCL-THINC-BVD scheme is implemented to all state variables and volume fraction, which realizes the consistency among volume fraction and other physical variables. Benchmark tests are carried out to verify the capability of the present method in capturing the material interface as a well-defined sharp jump in volume fraction, as well as significant improvement in solution quality. The proposed scheme is a simple and effective method of practical significance for simulating compressible interfacial multiphase flows.


title: Derivations, Automorphisms, and Representations of Complex $\omega$-Lie Algebras

abstract: Let $(\mathfrak{g},\omega)$ be a finite-dimensional non-Lie complex $\omega$-Lie algebra. We study the derivation algebra $Der(\mathfrak{g})$ and the automorphism group $Aut(\mathfrak{g})$ of $(\mathfrak{g},\omega)$. We introduce the notions of $\omega$-derivations and $\omega$-automorphisms of $(\mathfrak{g},\omega)$ which naturally preserve the bilinear form $\omega$. We show that the set $Der_{\omega}(\mathfrak{g})$ of all $\omega$-derivations is a Lie subalgebra of $Der(\mathfrak{g})$ and the set $Aut_{\omega}(\mathfrak{g})$ of all $\omega$-automorphisms is a subgroup of $Aut(\mathfrak{g})$. For any 3-dimensional and 4-dimensional nontrivial $\omega$-Lie algebra $\mathfrak{g}$, we compute $Der(\mathfrak{g})$ and $Aut(\mathfrak{g})$ explicitly, and study some Lie group properties of $Aut(\mathfrak{g})$. We also study representation theory of $\omega$-Lie algebras. We show that all 3-dimensional nontrivial $\omega$-Lie algebras are multiplicative, as well as we provide a 4-dimensional example of $\omega$-Lie algebra that is not multiplicative. Finally, we show that any irreducible representation of the simple $\omega$-Lie algebra $C_{\alpha}(\alpha\neq 0,-1)$ is 1-dimensional.


title: Statistical Methods for the Neutrino Mass Hierarchy

abstract: In the next decade several experiments will attempt to determine the neutrino mass hierarchy, i.e. the sign of $\Delta m_{31}^2$. In the last years it was noticed that the two hierarchies are disjoint hypotheses and, for this reason, Wilks' theorem cannot be applied: this means that $\Delta\chi^2=\chi^2_{IH}-\chi^2_{NH}$ does not follow a one-degree-of-freedom chi-square distribution. It was proven that, under certain assumptions, it follows instead a Gaussian distribution with $\sigma=2\sqrt{\mu}$. I will present several possible definition of sensitivity and review the approaches proposed in the literature, both within the Bayesian and the frequentist framework, examining advantages and disadvantages and discussing how they should be modified if the conditions for Gaussianity are not fulfilled. I will also discuss the possibility of introducing a new pull parameter in order to avoid the issue related to the non-nested hypotheses and the differences between marginalization and minimization, showing under which conditions the two procedures yield the same $\Delta\chi^2$.


title: A remark on the negative activation volume for defects in solids

abstract: The evolution of the aging process of glassy materials quenched from temperatures above their glass transition temperature $T_g$, when plotting the relaxed enthalpy versus the decrease in volume, leads to a slope comparable to the isothermal compressibility close to $T_g$. This empirical result was explained in an earlier publication (V. Katsika-Tsigourakou, G. E. Zardas J. Non-Cryst. Solids 356 (2010) 179-180) by means of a thermodynamical model. Here, we show that the same model enables the explanation of the rare cases of negative defect activation in solids.


title: The loss surface of deep and wide neural networks

abstract: While the optimization problem behind deep neural networks is highly non-convex, it is frequently observed in practice that training deep networks seems possible without getting stuck in suboptimal points. It has been argued that this is the case as all local minima are close to being globally optimal. We show that this is (almost) true, in fact almost all local minima are globally optimal, for a fully connected network with squared loss and analytic activation function given that the number of hidden units of one layer of the network is larger than the number of training points and the network structure from this layer on is pyramidal.


title: On the exponent of exponential convergence of the $p$-version FEM spaces

abstract: We study the exponent of the exponential rate of convergence in terms of the number of degrees of freedom for various non-standard {$p$-version} finite element spaces employing reduced cardinality basis. More specifically, we show that serendipity finite element methods and discontinuous Galerkin finite element methods with total degree $\mathcal{P}_p$ basis have a faster exponential convergence with respect to the number of degrees of freedom than their counterparts employing the tensor product $\mathcal{Q}_p$ basis for quadrilateral/hexahedral elements, for piecewise analytic problems under $p$-refinement. The above results are proven by using a new $p$-optimal error bound for the $L^2$-orthogonal projection onto the total degree $\mathcal{P}_p$ basis, and for the $H^1$-projection onto the serendipity finite element space over tensor product elements with dimension $d\geq2$. These new $p$-optimal error bounds lead to a larger exponent of the exponential rate of convergence with respect to the number of degrees of freedom. Moreover, these results show that part of the basis functions in $\mathcal{Q}_p$ basis {plays} no roles in achieving the $hp$-optimal error bound in the Sobolev space. The sharpness of theoretical results is also verified by a series of numerical examples.


title: Study of the composition dependence of the ionic conductivity of LiH_{x}D_{1-x} alloys

abstract: Recent observations have shown that at a certain composition of solid solutions AgCl$_x$Br$_{1-x}$, AgCl-CdCl$_2$, AgBr-CdBr$_2$, of ionic crystals, the electrical conductivity exhibits a value appreciably larger than that of the end constituents. Here, we investigate the electrical conductivity $\sigma$ of LiH$_x$D$_{1-x}$ solid solutions -which are of prominent importance in fuel hydrogen storage applications- and find that in the whole composition range no maximum is likely to occur in the $\sigma$ versus x dependence.


title: Cluster algebras of finite type via a Coxeter element and Demazure Crystals of type B,C,D

abstract: For a classical group $G$ and a Coxeter element $c$ of the Weyl group, it is known that the coordinate ring $\mathbb{C}[G^{e,c^2}]$ of the double Bruhat cell $G^{e,c^2}:=B\cap B_-c^2B_-$ has a structure of cluster algebra of finite type, where $B$ and $B_-$ are opposite Borel subgroups. In this article, we consider the case $G$ is of type ${\rm B}_r$, ${\rm C}_r$ or ${\rm D}_r$ and describe all the cluster variables in $\mathbb{C}[G^{e,c^2}]$ as monomial realizations of certain Demazure crystals.


title: Epsilon-shapes: characterizing, detecting and thickening thin features in geometric models

abstract: We focus on the analysis of planar shapes and solid objects having thin features and propose a new mathematical model to characterize them. Based on our model, that we call an epsilon-shape, we show how thin parts can be effectively and efficiently detected by an algorithm, and propose a novel approach to thicken these features while leaving all the other parts of the shape unchanged. When compared with state-of-the-art solutions, our proposal proves to be particularly flexible, efficient and stable, and does not require any unintuitive parameter to fine-tune the process. Furthermore, our method is able to detect thin features both in the object and in its complement, thus providing a useful tool to detect thin cavities and narrow channels. We discuss the importance of this kind of analysis in the design of robust structures and in the creation of geometry to be fabricated with modern additive manufacturing technology.


title: On Maximizing Sensor Network Lifetime by Energy Balancing

abstract: Many physical systems, such as water/electricity distribution networks, are monitored by battery-powered Wireless Sensor Networks (WSNs). Since battery replacement of sensor nodes is generally difficult, long-term monitoring can be only achieved if the operation of the WSN nodes contributes to a long WSN lifetime. Two prominent techniques to long WSN lifetime are i) optimal sensor activation and ii) efficient data gathering and forwarding based on compressive sensing. These techniques are feasible only if the activated sensor nodes establish a connected communication network (connectivity constraint), and satisfy a compressive sensing decoding constraint (cardinality constraint). These two constraints make the problem of maximizing network lifetime via sensor node activation and compressive sensing NP-hard. To overcome this difficulty, an alternative approach that iteratively solves energy balancing problems is proposed. However, understanding whether maximizing network lifetime and energy balancing problems are aligned objectives is a fundamental open issue. The analysis reveals that the two optimization problems give different solutions, but the difference between the lifetime achieved by the energy balancing approach and the maximum lifetime is small when the initial energy at sensor nodes is significantly larger than the energy consumed for a single transmission. The lifetime achieved by the energy balancing is asymptotically optimal, and that the achievable network lifetime is at least $50$\% of the optimum. Analysis and numerical simulations quantify the efficiency of the proposed energy balancing approach.


title: Three- and four-electron integrals involving Gaussian geminals: fundamental integrals, upper bounds and recurrence relations

abstract: We report the three main ingredients to calculate three- and four-electron integrals over Gaussian basis functions involving Gaussian geminal operators: fundamental integrals, upper bounds, and recurrence relations. In particular, we consider the three- and four-electron integrals that may arise in explicitly-correlated F12 methods. A straightforward method to obtain the fundamental integrals is given. We derive vertical, transfer and horizontal recurrence relations to build up angular momentum over the centers. Strong, simple and scaling-consistent upper bounds are also reported. This latest ingredient allows to compute only the $\order{N^2}$ significant three- and four-electron integrals, avoiding the computation of the very large number of negligible integrals.


title: The columnar-disorder phase boundary in a mixture of hard squares and dimers

abstract: A mixture of hard squares, dimers and vacancies on a square lattice is known to undergo a transition from a low-density disordered phase to high-density columnar ordered phase. Along the fully packed square-dimer line, the system undergoes an Kosterliz-Thouless type transition to a phase with power law correlations. We estimate the phase boundary separating the ordered and disordered phases by calculating the interfacial tension between two differently ordered phases within two different approximation schemes. The analytically obtained phase boundary is in good agreement with Monte Carlo simulations.


title: Canonical RDEs and general semimartingales as rough paths

abstract: In the spirit of Marcus canonical stochastic differential equations, we study a similar notion of rough differential equations (RDEs), notably dropping the assumption of continuity prevalent in the rough path literature. A new metric is exhibited in which the solution map is a continuous function of the driving rough path and a so-called path function, which directly models the effect of the jump on the system. In a second part, we show that general multidimensional semimartingales admit canonically defined rough path lifts. An extension of L\'epingle's BDG inequality to this setting is given, and in turn leads to a number of novel limit theorems for semimartingale driven differential equations, both in law and in probability, conveniently phrased via Kurtz-Protter's uniformly-controlled-variations (UCV) condition. A number of examples illustrate the scope of our results.


title: Anomalous lattice dynamics in a sigma-Fe60V40 alloy: Mossbauer spectroscopic study

abstract: Lattice dynamics in a sigma-Fe60V40 compound, which shows a re-entrant magnetism and orders ferromagnetic ally at ca. 170K, was investigated with the Mossbauer spectroscopy in the temperature interval of 5-300 K. Two relevant spectral parameters viz. the average center shift, <CS>, and the relative recoil-free fraction, f/fo, were explored. The former yielded the Debye temperature, T_D1, and the mean-square velocity of vibration, <v^2>, while the latter T_D2 and the mean-square amplitude of vibrations, <x^2>. Significant differences in the lattice-dynamical behaviors in the magnetic and paramagnetic phases were revealed. In particular, the values of T_D were notably lower and those of f/fo greatly higher in the former. This anomalous result has likely its origin in a remarkably high inharmonic contribution to the vibrations found for the ground magnetic state (spin-glass). Especially anomalous behavior vs. temperature exhibits <x^2> where four well defined ranges could have been identified and ascribed to the paramagnetic, ferromagnetic and two spin-glass phases. Linear correlations between <v^2>-<x^2> were found within each of the four ranges. They enabled determination of force constants, hence a change of the potential energy, Ep, in each of the ranges. The total change of Ep=30 meV while the corresponding one of the kinetic energy, determined from the knowledge of <v^2>, was Ek=21 meV. The lack of balance between Ep and Ek follows from the anharmonic lattice-dynamical behavior observed in the spin-glass state. The results give a strong evidence that magnetism can significantly affect the lattice dynamics.


title: Optimizing Automata Learning via Monads

abstract: Automata learning has been successfully applied in the verification of hardware and software. The size of the automaton model learned is a bottleneck for scalability, and hence optimizations that enable learning of compact representations are important. This paper exploits monads, both as a mathematical structure and a programming construct, to design, prove correct, and implement a wide class of such optimizations. The former perspective on monads allows us to develop a new algorithm and accompanying correctness proofs, building upon a general framework for automata learning based on category theory. The new algorithm is parametric on a monad, which provides a rich algebraic structure to capture non-determinism and other side-effects. We show that our approach allows us to uniformly capture existing algorithms, develop new ones, and add optimizations. The latter perspective allows us to effortlessly translate the theory into practice: we provide a Haskell library implementing our general framework, and we show experimental results for two specific instances: non-deterministic and weighted automata.


title: Nonequilibrium interactions between ideal polymers and a repulsive surface

abstract: We use Newtonian and overdamped Langevin dynamics to study long flexible polymers dragged by an external force at a constant velocity $v$. The work $W$ by that force depends on the initial state of the polymer and the details of the process. Jarzynski equality can be used to relate the non-equilibrium work distribution $P(W)$ obtained from repeated experiments to equilibrium free energy difference $\Delta F$ between the initial and final states. We use the power law dependence of the geometrical and dynamical characteristics of the polymer on the number of monomers $N$ to suggest the existence of a critical velocity $v_c(N)$, such that for $v<v_c$ the reconstruction of $\Delta F$ is an easy task, while for $v$ significantly exceeding $v_c$ it becomes practically impossible. We demonstrate the existence of such $v_c$ analytically for ideal polymer in free space and numerically for a polymer being dragged away from a repulsive wall. Our results suggest that the distribution of the dissipated work $W_{\rm d}=W-\Delta F$ in properly scaled variables approaches a limiting shape for large $N$.


title: Local $h$-vectors of Quasi-Geometric and Barycentric Subdivisions

abstract: In this paper, we answer two questions on local $h$-vectors, which were asked by Athanasiadis. First, we characterize all possible local $h$-vectors of quasi-geometric subdivisions of a simplex. Second, we prove that the local $\gamma$-vector of the barycentric subdivision of any CW-regular subdivision of a simplex is nonnegative. Along the way, we derive a new recurrence formula for the derangement polynomials.