title: Feature selection algorithm based on Catastrophe model to improve the performance of regression analysis

abstract: In this paper we introduce a new feature selection algorithm to remove the irrelevant or redundant features in the data sets. In this algorithm the importance of a feature is based on its fitting to the Catastrophe model. Akaike information crite- rion value is used for ranking the features in the data set. The proposed algorithm is compared with well-known RELIEF feature selection algorithm. Breast Cancer, Parkinson Telemonitoring data and Slice locality data sets are used to evaluate the model.


title: A tale of two modes: Neutrino free-streaming in the early universe

abstract: We present updated constraints on the free-streaming nature of cosmological neutrinos from cosmic microwave background (CMB) power spectra, baryonic acoustic oscillation data, and local measurements of the Hubble constant. Specifically, we consider a Fermi-like four-fermion interaction between massless neutrinos, characterized by an effective coupling constant $ G_{\rm eff}$, and resulting in a neutrino opacity $\dot{\tau}_\nu\propto G_{\rm eff}^2 T_\nu^5$. Using a conservative prior on the parameter $\log_{10}\left(G_{\rm eff} {\rm MeV}^2\right)$, we find a bimodal posterior distribution. The first of these modes is consistent with the standard $\Lambda$CDM cosmology and corresponds to neutrinos decoupling at redshift $z_{\nu,{\rm dec}} > 1.3\times10^5$. The other mode of the posterior, dubbed the "interacting neutrino mode", corresponds to neutrino decoupling occurring within a narrow redshift window centered around $z_{\nu,{\rm dec}}\sim8300$. This mode is characterized by a high value of the effective neutrino coupling constant, together with a lower value of the scalar spectral index and amplitude of fluctuations, and a higher value of the Hubble parameter. Using both a maximum likelihood analysis and the ratio of the two mode's Bayesian evidence, we find the interacting neutrino mode to be statistically disfavored compared to the standard $\Lambda$CDM cosmology. Interestingly, the addition of CMB polarization and direct Hubble constant measurements significantly raises the statistical significance of this secondary mode, indicating that new physics in the neutrino sector could help explain the difference between local measurements of $H_0$, and those inferred from CMB data. A robust consequence of our results is that neutrinos must be free streaming long before the epoch of matter-radiation equality.


title: Black-Hole Growth is Mainly Linked to Host-Galaxy Stellar Mass rather than Star Formation Rate

abstract: We investigate the dependence of black-hole accretion rate (BHAR) on host-galaxy star formation rate (SFR) and stellar mass ($M_*$) in the CANDELS/GOODS-South field in the redshift range of $0.5\leq z < 2.0$. Our sample consists of $\approx 18000$ galaxies, allowing us to probe galaxies with $0.1 \lesssim \mathrm{SFR} \lesssim 100\ M_\odot\ \mathrm{yr}^{-1}$ and/or $10^8 \lesssim M_* \lesssim 10^{11}\ M_{\odot}$. We use sample-mean BHAR to approximate long-term average BHAR. Our sample-mean BHARs are derived from the $Chandra$ Deep Field-South 7 Ms observations, while the SFRs and $M_*$ have been estimated by the CANDELS team through SED fitting. The average BHAR is correlated positively with both SFR and $M_*$, and the BHAR-SFR and BHAR-$M_*$ relations can both be described acceptably by linear models with a slope of unity. However, BHAR appears to be correlated more strongly with $M_*$ than SFR. This result indicates that $M_*$ is the primary host-galaxy property related to black-hole growth, and the apparent BHAR-SFR relation is largely a secondary effect due to the star-forming main sequence. Among our sources, massive galaxies ($M_* \gtrsim 10^{10} M_{\odot}$) have significantly higher BHAR/SFR ratios than less-massive galaxies, indicating the former have higher black-hole fueling efficiency and/or higher SMBH occupation fraction than the latter. Our results can naturally explain the observed proportionality between $M_{\rm BH}$ and $M_*$ for local giant ellipticals, and suggest their $M_{\rm BH}/M_*$ is higher than that of local star-forming galaxies. Among local star-forming galaxies, massive systems might have higher $M_{\rm BH}/M_*$ compared to dwarfs.


title: $R_{D^{(*)}}$ anomaly: A possible hint for natural supersymmetry with $R$-parity violation

abstract: Recently, several $B$-physics experiments have reported an appreciable deviation from the Standard Model (SM) in the tree-level observables $R_{D^{(*)}}$; the combined weighted average now stands at $\approx 4 \sigma$. We first show the anomaly necessarily implies model-independent collider signals of the form $pp \to b \tau \nu$ that should be expediously searched for at ATLAS/CMS as a complementary test of the anomaly. Next we suggest a possible interconnection of the anomaly with the radiative stability of the Standard Model Higgs boson and point to a minimal effective supersymmetric scenario with $R$-parity violation as the underlying cause. We also comment on the possibility of simultaneously explaining the recently reported $R_{K^{(*)}}$ anomaly in this setup.


title: Weyl Metals

abstract: Weyl metal is the first example of a conducting material with a nontrivial electronic structure topology, making it distinct from an ordinary metal. Unlike in insulators, the nontrivial topology is not related to invariants, associated with completely filled bands, but with ones, associated with the Fermi surface. The Fermi surface of a topological metal consists of disconnected sheets, each enclosing a Weyl node, which is a point of contact between two nondegenerate bands. Such a point contact acts as a source of Berry curvature, or a magnetic monopole in momentum space. Its charge, or the flux of the Berry curvature through the enclosing Fermi surface sheet, is a topological invariant. We review the current state of this rapidly growing field, with a focus on bulk transport phenomena in topological metals.


title: Can cluster merger shocks reproduce the luminosity and shape distribution of radio relics?

abstract: Radio relics in galaxy clusters are believed to trace merger shock fronts. If cosmological structure formation determines the luminosity, size and shape distributions of radio relics then merger shocks need to be lighted up in a homogeneous way. We investigate if a mock relic sample, obtained from zoomed galaxy cluster simulations, is able to match the properties of relics measured in the NRAO VLA Sky Survey (NVSS). We compile a list of all radio relics known to date and homogeneously measure their parameters in all NVSS images and apply the same procedure to relics in our simulations. Number counts in the mock relic sample increase more steeply towards lower relic flux densities, suggesting an incompleteness of NVSS in this regime. Overall, we find that NVSS and mock samples show similar properties. However, large simulated relics tend to be somewhat smaller and closer to the cluster centre than observed ones. Besides this, the mock sample reproduces very well-known correlations for radio relics, in particular those relating the radio luminosity with the largest linear size and the X-ray luminosity. We show that these correlations are largely governed by the sensitivity of the NVSS observations. Mock relics show a similar orientation with respect to the direction to the cluster centre as the NVSS sample. Moreover, we find that their maximum radio luminosity roughly correlates with cluster mass, although displaying a large scatter. The overall good agreement between NVSS and the mock sample suggests that properties of radio relics are indeed governed by merger shock fronts, emitting in a homogeneous fashion. Our study demonstrates that the combination of mock observations and data from upcoming radio surveys will allow to shed light on both the origin of radio relics and the nature of the intracluster medium.


title: Fault-Tolerant Quantum Computing in the Pauli or Clifford Frame with Slow Error Diagnostics

abstract: We consider the problem of fault-tolerant quantum computation in the presence of slow error diagnostics, either caused by measurement latencies or slow decoding algorithms. Our scheme offers a few improvements over previously existing solutions, for instance it does not require active error correction and results in a reduced error-correction overhead when error diagnostics is much slower than the gate time. In addition, we adapt our protocol to cases where the underlying error correction strategy chooses the optimal correction amongst all Clifford gates instead of the usual Pauli gates. The resulting Clifford frame protocol is of independent interest as it can increase error thresholds and could find applications in other areas of quantum computation.


title: A scale at 10 MeV, gravitational topological vacuum, and large extra dimensions

abstract: We discuss a possible scale of gravitational origin at around $10$ MeV, or $10^{-12}$ cm, which arises in the MacDowell-Mansouri formalism of gravity due to the topological Gauss-Bonnet term in the action, as pointed out by Bjorken several years ago. A length scale of the same size emerges also in the Kodama solution in gravity, which is known to be closely related to the MacDowell-Mansouri formulation. We particularly draw attention to the intriguing incident that existence of six compact extra dimensions originated from TeV-scale quantum gravity as well points to a length scale of $10^{-12}$ cm, as the compactification scale. The presence of six such extra dimensions is also in remarkable consistency with the MacDowell-Mansouri formalism; it provides a possible explanation for the factor of $\sim10^{120}$ multiplying the Gauss-Bonnet term in the action. We also comment on the relevant implications of such a scale regarding the thermal history of the universe motivated by the fact that it is considerably close to $1-2$ MeV below which the weak interactions freeze out, leading to Big Bang Nucleosynthesis.


title: An opening criterion for dust gaps in protoplanetary discs

abstract: We aim to understand under which conditions a low mass planet can open a gap in viscous dusty protoplanetary discs. For this purpose, we extend the theory of dust radial drift to include the contribution from the tides of an embedded planet and from the gas viscous forces. From this formalism, we derive i) a grain size-dependent criterion for dust gap opening in discs, ii) an estimate of the location of the outer edge of the dust gap and iii) an estimate of the minimum Stokes number above which low-mass planets are able to carve gaps which appear only in the dust disc. These analytical estimates are particularly helpful to appraise the minimum mass of an hypothetical planet carving gaps in discs observed at long wavelengths and high resolution. We validate the theory against 3D SPH simulations of planet-disc interaction in a broad range of dusty protoplanetary discs. We find a remarkable agreement between the theoretical model and the numerical experiments.


title: The Complete Calibration of the Color-Redshift Relation (C3R2) Survey: Survey Overview and Data Release 1

abstract: A key goal of the Stage IV dark energy experiments Euclid, LSST and WFIRST is to measure the growth of structure with cosmic time from weak lensing analysis over large regions of the sky. Weak lensing cosmology will be challenging: in addition to highly accurate galaxy shape measurements, statistically robust and accurate photometric redshift (photo-z) estimates for billions of faint galaxies will be needed in order to reconstruct the three-dimensional matter distribution. Here we present an overview of and initial results from the Complete Calibration of the Color-Redshift Relation (C3R2) survey, designed specifically to calibrate the empirical galaxy color-redshift relation to the Euclid depth. These redshifts will also be important for the calibrations of LSST and WFIRST. The C3R2 survey is obtaining multiplexed observations with Keck (DEIMOS, LRIS, and MOSFIRE), the Gran Telescopio Canarias (GTC; OSIRIS), and the Very Large Telescope (VLT; FORS2 and KMOS) of a targeted sample of galaxies most important for the redshift calibration. We focus spectroscopic efforts on under-sampled regions of galaxy color space identified in previous work in order to minimize the number of spectroscopic redshifts needed to map the color-redshift relation to the required accuracy. Here we present the C3R2 survey strategy and initial results, including the 1283 high confidence redshifts obtained in the 2016A semester and released as Data Release 1.


title: Goodness of fit test under progressive Type-I interval censoring

abstract: In this paper, we propose several statistics for testing uniformity under progressive Type-I interval censoring. We obtain the critical points of these statistics and study the power of the proposed tests against a representative set of alternatives via simulation. Finally, we generalize our methods for continuous and completely specified distributions.


title: Perfect divisibility and 2-divisibility

abstract: A graph $G$ is said to be $2$-divisible if for all (nonempty) induced subgraphs $H$ of $G$, $V(H)$ can be partitioned into two sets $A,B$ such that $\omega(A) < \omega(H)$ and $\omega(B) < \omega(H)$. A graph $G$ is said to be perfectly divisible if for all induced subgraphs $H$ of $G$, $V(H)$ can be partitioned into two sets $A,B$ such that $H[A]$ is perfect and $\omega(B) < \omega(H)$. We prove that if a graph is $(P_5,C_5)$-free, then it is $2$-divisible. We also prove that if a graph is bull-free and either odd-hole-free or $P_5$-free, then it is perfectly divisible.


title: Gate-controlled reversible rectifying behaviour in tunnel contacted atomically-thin MoS$_{2}$ transistor

abstract: Atomically-thin 2D semiconducting materials integrated into van der Waals heterostructures have enabled architectures that hold great promise for next generation nanoelectronics. However, challenges still remain to enable their full acceptance as compliant materials for integration in logic devices. Two key-components to master are the barriers at metal/semiconductor interfaces and the mobility of the semiconducting channel, which endow the building-blocks of ${pn}$ diode and field effect transistor. Here, we have devised a reverted stacking technique to intercalate a wrinkle-free h-BN tunnel layer between MoS$_{2}$ channel and contacting electrodes. Vertical tunnelling of electrons therefore makes it possible to suppress the Schottky barriers and Fermi level pinning, leading to homogeneous gate-control of the channel chemical potential across the bandgap edges. The observed unprecedented features of ambipolar ${pn}$ to ${np}$ diode, which can be reversibly gate tuned, paves the way for future logic applications and high performance switches based on atomically thin semiconducting channel.


title: Solving the Navier-Lame Equation in Cylindrical Coordinates Using the Buchwald Representation: Some Parametric Solutions with Applications

abstract: Using a separable Buchwald representation in cylindrical coordinates, we show how under certain conditions the coupled equations of motion governing the Buchwald potentials can be decoupled and then solved using well-known techniques from the theory of PDEs. Under these conditions, we then construct three parametrized families of particular solutions to the Navier-Lame equation in cylindrical coordinates. In this paper, we specifically construct solutions having 2pi-periodic angular parts. These particular solutions can be directly applied to a fundamental set of linear elastic boundary value problems in cylindrical coordinates and are especially suited to problems involving one or more physical parameters. As an illustrative example, we consider the problem of determining the response of a solid elastic cylinder subjected to a time-harmonic surface pressure that varies sinusoidally along its axis, and we demonstrate how the obtained parametric solutions can be used to efficiently construct an exact solution to this problem. We also briefly consider applications to some related forced-relaxation type problems.


title: Status of the SuperNEMO 0$\nu\beta\beta$ experiment

abstract: SuperNEMO is an ultra-low-background tracker-calorimeter experiment designed to look for the neutrinoless double-beta decay of various isotopes. We present the current state of the experiment's Demonstrator Module, which is currently being installed and and commissioned in the LSM in France.


title: Extreme CO Isotopic Abundances in the ULIRG IRAS 13120-5453: An Extremely Young Starburst or Top-Heavy Initial Mass Function

abstract: We present ALMA $^{12}$CO (J=1-0, 3-2 and 6-5), $^{13}$CO (J=1-0) and C$^{18}$O (J=1-0) observations of the local Ultra Luminous Infrared Galaxy, IRAS 13120-5453 (dubbed "The Yo-yo"). The morphologies of the three isotopic species differ, where $^{13}$CO shows a hole in emission towards the center. We measure integrated brightness temperature line ratios of $^{12}$CO/$^{13}$CO $\geq$ 60 (exceeding 200) and $^{13}$CO/C$^{18}$O $\leq$ 1 in the central region. Assuming optical thin emission, C$^{18}$O is more abundant than $^{13}$CO in several regions. The abundances within the central 500 pc are consistent with enrichment of the ISM via a young starburst ($<$7Myr), a top-heavy initial mass function or a combination of both.


title: Periodic modules and acyclic complexes

abstract: We study the behaviour of modules $M$ that fit into a short exact sequence $0\to M\to C\to M\to 0$, where $C$ belongs to a class of modules $\mathcal C$, the so-called $\mathcal C$-periodic modules. We find a rather general framework to improve and generalize some well-known results of Benson and Goodearl and Simson. In the second part we will combine techniques of hereditary cotorsion pairs and presentation of direct limits, to conclude, among other applications, that if $M$ is any module and $C$ is cotorsion, then $M$ will be also cotorsion. This will lead to some meaningful consequences in the category $\textrm{Ch}(R)$ of unbounded chain complexes and in Gorenstein homological algebra. For example we show that every acyclic complex of cotorsion modules has cotorsion cycles, and more generally, every map $F\to C$ where $C$ is a complex of cotorsion modules and $F$ is an acyclic complex of flat cycles, is null-homotopic. In other words, every complex of cotorsion modules is dg-cotorsion.


title: Revisiting wireless network jamming by SIR-based considerations and Multiband Robust Optimization

abstract: We revisit the mathematical models for wireless network jamming introduced by Commander et al.: we first point out the strong connections with classical wireless network design and then we propose a new model based on the explicit use of signal-to-interference quantities. Moreover, to address the intrinsic uncertain nature of the jamming problem and tackle the peculiar right-hand-side (RHS) uncertainty of the problem, we propose an original robust cutting-plane algorithm drawing inspiration from Multiband Robust Optimization. Finally, we assess the performance of the proposed cutting plane algorithm by experiments on realistic network instances.


title: GUB Covers and Power-Indexed formulations for Wireless Network Design

abstract: We propose a pure 0-1 formulation for the wireless network design problem, i.e. the problem of configuring a set of transmitters to provide service coverage to a set of receivers. In contrast with classical mixed integer formulations, where power emissions are represented by continuous variables, we consider only a finite set of powers values. This has two major advantages: it better fits the usual practice and eliminates the sources of numerical problems which heavily affect continuous models. A crucial ingredient of our approach is an effective basic formulation for the single knapsack problem representing the coverage condition of a receiver. This formulation is based on the GUB cover inequalities introduced by Wolsey (1990) and its core is an extension of the exact formulation of the GUB knapsack polytope with two GUB constraints. This special case corresponds to the very common practical situation where only one major interferer is present. We assess the effectiveness of our formulation by comprehensive computational results over realistic instances of two typical technologies, namely WiMAX and DVB-T.


title: Pseudo-Majoron as Light Mediator of Singlet Scalar Dark Matter

abstract: In the singlet-triplet majoron model of neutrino mass, lepton number is spontaneously broken. If it is also softly broken, then a naturally light pseudoscalar particle $\eta_I$ exists. It may then act as a light mediator for a real singlet scalar $\chi$ with odd dark parity. It is itself unstable, but decays dominantly to two neutrinos through its triplet scalar component, thereby not disturbing the cosmic microwave background (CMB). It also mixes with the standard-model Higgs boson only in one loop, thereby not contributing significantly to the elastic scattering of $\chi$ off nuclei in dark-matter direct-search experiments.


title: Modular Multi-Objective Deep Reinforcement Learning with Decision Values

abstract: In this work we present a method for using Deep Q-Networks (DQNs) in multi-objective environments. Deep Q-Networks provide remarkable performance in single objective problems learning from high-level visual state representations. However, in many scenarios (e.g in robotics, games), the agent needs to pursue multiple objectives simultaneously. We propose an architecture in which separate DQNs are used to control the agent's behaviour with respect to particular objectives. In this architecture we introduce decision values to improve the scalarization of multiple DQNs into a single action. Our architecture enables the decomposition of the agent's behaviour into controllable and replaceable sub-behaviours learned by distinct modules. Moreover, it allows to change the priorities of particular objectives post-learning, while preserving the overall performance of the agent. To evaluate our solution we used a game-like simulator in which an agent - provided with high-level visual input - pursues multiple objectives in a 2D world.


title: Select and Permute: An Improved Online Framework for Scheduling to Minimize Weighted Completion Time

abstract: In this paper, we introduce a new online scheduling framework for minimizing total weighted completion time in a general setting. The framework is inspired by the work of Hall et al. [Mathematics of Operations Research, Vol 22(3):513-544, 1997] and Garg et al. [Proc. of Foundations of Software Technology and Theoretical Computer Science, pp. 96-107, 2007], who show how to convert an offline approximation to an online scheme. Our framework uses two offline approximation algorithms (one for the simpler problem of scheduling without release times, and another for the minimum unscheduled weight problem) to create an online algorithm with provably good competitive ratios. We illustrate multiple applications of this method that yield improved competitive ratios. Our framework gives algorithms with the best or only known competitive ratios for the concurrent open shop, coflow, and concurrent cluster models. We also introduce a randomized variant of our framework based on the ideas of Chakrabarti et al. [Proc. of International Colloquium on Automata, Languages, and Programming, pp. 646-657, 1996] and use it to achieve improved competitive ratios for these same problems.


title: Interactions of Solitary Pulses of E. coli in a One-Dimensional Nutrient Gradient

abstract: We study an anomalous behavior observed in interacting E. coli populations. When two populations of E. coli are placed on opposite ends of a long channel with a supply of nutrient between them, they will travel as pulses toward one another up the nutrient gradient. We present experimental evidence that, counterintuitively, the two pulses will in some cases change direction and be- gin moving away from each other and the nutrient back toward the end of the channel from which they originated. Simulations of the Keller-Segel chemotaxis model reproduce the experimental results. To gain better insight to the phenomenon, we introduce a heuristic approximation to the spatial profile of each population in the Keller-Segel model to derive a system of ordinary differential equations approximating the temporal dynamics of its center of mass and width. This approximate model simplifies analysis of the global dynamics of the bacterial system and allows us to efficiently explore the qualitative behavior changes across variations of parameters, and thereby provides experimentally testable hypotheses about the mechanisms behind the turnaround behavior.


title: A cost-effective isogeometric approach for composite plates based on a stress recovery procedure

abstract: This paper introduces a cost-effective strategy to simulate the behavior of laminated plates by means of isogeometric 3D solid elements. Exploiting the high continuity of spline functions and their properties, a proper out-of-plane stress state is recovered from a coarse displacement solution using a post-processing step based on the enforcement of equilibrium in strong form. Appealing results are obtained and the method is shown to be particularly Peffective on slender composite stacks with a large number of layers.


title: On the techniques for constructing admissible stress fields in model verification: Performances on engineering examples

abstract: Robust global/goal-oriented error estimation is used nowadays to control the approximate finite element solutions obtained from simulation. In the context of Computational Mechanics, the construction of admissible stress fields (\ie stress tensors which verify the equilibrium equations) is required to set up strict and guaranteed error bounds (using residual based error estimators) and plays an important role in the quality of the error estimates. This work focuses on the different procedures used in the calculation of admissible stress fields, which is a crucial and technically complicated point. The three main techniques that currently exist, called the element equilibration technique (EET), the star-patch equilibration technique (SPET), and the element equilibration + star-patch technique (EESPT), are investigated and compared with respect to three different criteria, namely the quality of associated error estimators, computational cost and easiness of practical implementation into commercial finite element codes. The numerical results which are presented focus on industrial problems; they highlight the main advantages and drawbacks of the different methods and show that the behavior of the three estimators, which have the same convergence rate as the exact global error, is consistent. Two- and three-dimensional experiments have been carried out in order to compare the performance and the computational cost of the three different approaches. The analysis of the results reveals that the SPET is more accurate than EET and EESPT methods, but the corresponding computational cost is higher. Overall, the numerical tests prove the interest of the hybrid method EESPT and show that it is a correct compromise between quality of the error estimate, practical implementation and computational cost.


title: Sumo Puff: Tidal Debris or Disturbed Ultra-Diffuse Galaxy?

abstract: We report the discovery of a diffuse stellar cloud with an angular extent $\gtrsim30^{\prime\prime}$, which we term "Sumo Puff", in data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). While we do not have a redshift for this object, it is in close angular proximity to a post-merger galaxy at redshift $z=0.0431$ and is projected within a few virial radii (assuming similar redshifts) of two other ${\sim}L_\star$ galaxies, which we use to bracket a potential redshift range of $0.0055 < z < 0.0431$. The object's light distribution is flat, as characterized by a low Sersic index ($n\sim0.3$). It has a low central $g$-band surface brightness of ${\sim}26.4$ mag arcsec$^{-2}$, large effective radius of ${\sim}13^{\prime\prime}$ (${\sim}11$ kpc at $z=0.0431$ and ${\sim}1.5$ kpc at $z=0.0055$), and an elongated morphology ($b/a\sim0.4$). Its red color ($g-i\sim1$) is consistent with a passively evolving stellar population and similar to the nearby post-merger galaxy, and we may see tidal material connecting Sumo Puff with this galaxy. We offer two possible interpretations for the nature of this object: (1) it is an extreme, galaxy-size tidal feature associated with a recent merger event, or (2) it is a foreground dwarf galaxy with properties consistent with a quenched, disturbed ultra-diffuse galaxy. We present a qualitative comparison with simulations that demonstrates the feasibility of forming a structure similar to this object in a merger event. Follow-up spectroscopy and/or deeper imaging to confirm the presence of the bridge of tidal material will be necessary to reveal the true nature of this object.


title: Magnetar Powered Ordinary Type IIP Supernovae

abstract: We investigate the properties of Type IIP supernovae that are dominantly powered by the rotational kinetic energy of the newly born neutron star. While the spin-down of a magnetar has previously been proposed as a viable energy source in the context of super-luminous supernovae, we show that a similar mechanism could produce both normal and peculiar Type IIP supernova light curves from red supergiant progenitors for a range of initial spin periods and equivalent dipole magnetic field strengths. Although the formation channel for such magnetars in a typical red supergiant progenitor is unknown, it is tantalizing that this proof of concept model is capable of producing ordinary Type IIP lightcurve properties, perhaps implying that rotation rate and magnetic field strength may play important roles in some ordinary looking Type IIP supernova explosions.


title: Shifting the Phase Transition Threshold for Random Graphs and 2-SAT using Degree Constraints

abstract: We show that by restricting the degrees of the vertices of a graph to an arbitrary set \( \Delta \), the threshold point $ \alpha(\Delta) $ of the phase transition for a random graph with $ n $ vertices and $ m = \alpha(\Delta) n $ edges can be either accelerated (e.g., $ \alpha(\Delta) \approx 0.381 $ for $ \Delta = \{0,1,4,5\} $) or postponed (e.g., $ \alpha(\{ 2^0, 2^1, \cdots, 2^k, \cdots \}) \approx 0.795 $) compared to a classical Erd\H{o}s--R\'{e}nyi random graph with $ \alpha(\mathbb Z_{\geq 0}) = \tfrac12 $. In particular, we prove that the probability of graph being nonplanar and the probability of having a complex component, goes from $ 0 $ to $ 1 $ as $ m $ passes $ \alpha(\Delta) n $. We investigate these probabilities and also different graph statistics inside the critical window of transition (diameter, longest path and circumference of a complex component).


title: A hybrid exact-ACO algorithm for the joint scheduling, power and cluster assignment in cooperative wireless networks

abstract: Base station cooperation (BSC) has recently arisen as a promising way to increase the capacity of a wireless network. Implementing BSC adds a new design dimension to the classical wireless network design problem: how to define the subset of base stations (clusters) that coordinate to serve a user. Though the problem of forming clusters has been extensively discussed from a technical point of view, there is still a lack of effective optimization models for its representation and algorithms for its solution. In this work, we make a further step towards filling such gap: 1) we generalize the classical network design problem by adding cooperation as an additional decision dimension; 2) we develop a strong formulation for the resulting problem; 3) we define a new hybrid solution algorithm that combines exact large neighborhood search and ant colony optimization. Finally, we assess the performance of our new model and algorithm on a set of realistic instances of a WiMAX network.


title: A Continuum Poisson-Boltzmann Model for Membrane Channel Proteins

abstract: Membrane proteins constitute a large portion of the human proteome and perform a variety of important functions as membrane receptors, transport proteins, enzymes, signaling proteins, and more. The computational studies of membrane proteins are usually much more complicated than those of globular proteins. Here we propose a new continuum model for Poisson-Boltzmann calculations of membrane channel proteins. Major improvements over the existing continuum slab model are as follows: 1) The location and thickness of the slab model are fine-tuned based on explicit-solvent MD simulations. 2) The highly different accessibility in the membrane and water regions are addressed with a two-step, two-probe grid labeling procedure, and 3) The water pores/channels are automatically identified. The new continuum membrane model is optimized (by adjusting the membrane probe, as well as the slab thickness and center) to best reproduce the distributions of buried water molecules in the membrane region as sampled in explicit water simulations. Our optimization also shows that the widely adopted water probe of 1.4 {\AA} for globular proteins is a very reasonable default value for membrane protein simulations. It gives an overall minimum number of inconsistencies between the continuum and explicit representations of water distributions in membrane channel proteins, at least in the water accessible pore/channel regions that we focus on. Finally, we validate the new membrane model by carrying out binding affinity calculations for a potassium channel, and we observe a good agreement with experiment results.


title: Integrable systems, symmetries and quantization

abstract: These notes correspond to a mini-course given at the Poisson 2016 conference in Geneva. Starting from classical integrable systems in the sense of Liouville, we explore the notion of non-degenerate singularity and expose recent research in connection with semi-toric systems. The quantum and semiclassical counterpart will also be presented, in the viewpoint of the inverse question: from the quantum mechanical spectrum, can you recover the classical system?


title: Scatteract: Automated extraction of data from scatter plots

abstract: Charts are an excellent way to convey patterns and trends in data, but they do not facilitate further modeling of the data or close inspection of individual data points. We present a fully automated system for extracting the numerical values of data points from images of scatter plots. We use deep learning techniques to identify the key components of the chart, and optical character recognition together with robust regression to map from pixels to the coordinate system of the chart. We focus on scatter plots with linear scales, which already have several interesting challenges. Previous work has done fully automatic extraction for other types of charts, but to our knowledge this is the first approach that is fully automatic for scatter plots. Our method performs well, achieving successful data extraction on 89% of the plots in our test set.


title: New bounding techniques for goal-oriented error estimation applied to linear problems

abstract: The paper deals with the accuracy of guaranteed error bounds on outputs of interest computed from approximate methods such as the finite element method. A considerable improvement is introduced for linear problems thanks to new bounding techniques based on Saint-Venant's principle. The main breakthrough of these optimized bounding techniques is the use of properties of homothetic domains which enables to cleverly derive guaranteed and accurate boundings of contributions to the global error estimate over a local region of the domain. Performances of these techniques are illustrated through several numerical experiments.


title: Mass Loss Rates from Coronal Mass Ejections: A Predictive Theoretical Model for Solar-Type Stars

abstract: Coronal mass ejections (CMEs) are eruptive events that cause a solar-type star to shed mass and magnetic flux. CMEs tend to occur together with flares, radio storms, and bursts of energetic particles. On the Sun, CME-related mass loss is roughly an order of magnitude less intense than that of the background solar wind. However, on other types of stars, CMEs have been proposed to carry away much more mass and energy than the time-steady wind. Earlier papers have used observed correlations between solar CMEs and flare energies, in combination with stellar flare observations, to estimate stellar CME rates. This paper sidesteps flares and attempts to calibrate a more fundamental correlation between surface-averaged magnetic fluxes and CME properties. For the Sun, there exists a power-law relationship between the magnetic filling factor and the CME kinetic energy flux, and it is generalized for use on other stars. An example prediction of the time evolution of wind/CME mass-loss rates for a solar-mass star is given. A key result is that for ages younger than about 1 Gyr (i.e., activity levels only slightly higher than the present-day Sun), the CME mass loss exceeds that of the time-steady wind. At younger ages, CMEs carry 10 to 100 times more mass than the wind, and such high rates may be powerful enough to dispel circumstellar disks and affect the habitability of nearby planets. The cumulative CME mass lost by the young Sun may have been as much as 1% of a solar mass.


title: Electronic Metamaterials with Tunable Second-order Optical Nonlinearities

abstract: The ability to engineer metamaterials with tunable nonlinear optical properties is crucial for nonlinear optics. Traditionally, metals have been employed to enhance nonlinear optical interactions through field localization. Here, inspired by the electronic properties of materials, we introduce and demonstrate experimentally an asymmetric metal-semiconductor-metal (MSM) metamaterial that exhibits a large and electronically tunable effective second-order optical susceptibility (\c{hi}(2)). The induced \c{hi}(2) originates from the interaction between the third-order optical susceptibility of the semiconductor (\c{hi}(3)) with the engineered internal electric field resulting from the two metals with dissimilar work function at its interfaces. We demonstrate a five times larger second-harmonic intensity from the MSM metamaterial, compared to contributions from its constituents with electrically tunable nonlinear coefficient ranging from 2.8 to 15.6 pm/V.


title: Optical quality assurance of GEM foils

abstract: An analysis software was developed for the high aspect ratio optical scanning system in the Detec- tor Laboratory of the University of Helsinki and the Helsinki Institute of Physics. The system is used e.g. in the quality assurance of the GEM-TPC detectors being developed for the beam diagnostics system of the SuperFRS at future FAIR facility. The software was tested by analyzing five CERN standard GEM foils scanned with the optical scanning system. The measurement uncertainty of the diameter of the GEM holes and the pitch of the hole pattern was found to be 0.5 {\mu}m and 0.3 {\mu}m, respectively. The software design and the performance are discussed. The correlation between the GEM hole size distribution and the corresponding gain variation was studied by comparing them against a detailed gain mapping of a foil and a set of six lower precision control measurements. It can be seen that a qualitative estimation of the behavior of the local variation in gain across the GEM foil can be made based on the measured sizes of the outer and inner holes.


title: Improving Semantic Composition with Offset Inference

abstract: Count-based distributional semantic models suffer from sparsity due to unobserved but plausible co-occurrences in any text collection. This problem is amplified for models like Anchored Packed Trees (APTs), that take the grammatical type of a co-occurrence into account. We therefore introduce a novel form of distributional inference that exploits the rich type structure in APTs and infers missing data by the same mechanism that is used for semantic composition.


title: SREFI: Synthesis of Realistic Example Face Images

abstract: In this paper, we propose a novel face synthesis approach that can generate an arbitrarily large number of synthetic images of both real and synthetic identities. Thus a face image dataset can be expanded in terms of the number of identities represented and the number of images per identity using this approach, without the identity-labeling and privacy complications that come from downloading images from the web. To measure the visual fidelity and uniqueness of the synthetic face images and identities, we conducted face matching experiments with both human participants and a CNN pre-trained on a dataset of 2.6M real face images. To evaluate the stability of these synthetic faces, we trained a CNN model with an augmented dataset containing close to 200,000 synthetic faces. We used a snapshot of this trained CNN to recognize extremely challenging frontal (real) face images. Experiments showed training with the augmented faces boosted the face recognition performance of the CNN.


title: Cataclysm no more: New views on the timing and delivery of lunar impactors

abstract: If properly interpreted, the impact record of the Moon, Earth's nearest neighbour, can be used to gain insights into how the Earth has been influenced by impacting events since its formation ~4.5 billion years (Ga) ago. However, the nature and timing of the lunar impactors - and indeed the lunar impact record itself - are not well understood. Of particular interest are the ages of lunar impact basins and what they tell us about the proposed "lunar cataclysm" and/or the late heavy bombardment (LHB), and how this impact episode may have affected early life on Earth or other planets. Investigations of the lunar impactor population over time have been undertaken and include analyses of orbital data and images; lunar, terrestrial, and other planetary sample data; and dynamical modelling. Here, the existing information regarding the nature of the lunar impact record is reviewed and new interpretations are presented. Importantly, it is demonstrated that most evidence supports a prolonged lunar (and thus, terrestrial) bombardment from ~4.2 to 3.4 Ga and not a cataclysmic spike at ~3.9 Ga. Implications for the conditions required for the origin of life are addressed.


title: Quantum State Tomography with a Single Observable

abstract: Quantum information has been drawing a wealth of research in recent years, shedding light on questions at the heart of quantum mechanics, as well as advancing fields such as complexity theory, cryptography, key distribution, and chemistry. These fundamental and applied aspects of quantum information rely on a crucial issue: the ability to characterize a quantum state from measurements, through a process called Quantum State Tomography (QST). However, QST requires a large number of measurements, each derived from a different physical observable corresponding to a different experimental setup. Unfortunately, changing the setup results in unwanted changes to the data, prolongs the measurement and impairs the assumptions that are always made about the stationarity of the noise. Here, we propose to overcome these drawbacks by performing QST with a single observable. A single observable can often be realized by a single setup, thus considerably reducing the experimental effort. In general, measurements of a single observable do not hold enough information to recover the quantum state. We overcome this lack of information by relying on concepts inspired by Compressed Sensing (CS), exploiting the fact that the sought state - in many applications of quantum information - is close to a pure state (and thus has low rank). Additionally, we increase the system dimension by adding an ancilla that couples to information evolving in the system, thereby providing more measurements, enabling the recovery of the original quantum state from a single-observable measurements. We demonstrate our approach on multi-photon states by recovering structured quantum states from a single observable, in a single experimental setup. We further show how this approach can be used to recover quantum states without number-resolving detectors.


title: Classical capacity per unit cost for quantum channels

abstract: In most communication scenarios, sending a symbol encoded in a quantum state requires spending resources such as energy, which can be quantified by a cost of communication. A standard approach in this context is to quantify the performance of communication protocol by classical capacity, quantifying the maximal amount of information that can be transmitted through a quantum channel per single use of the channel. However, different figures of merit are also possible, and a particularly well-suited one is the classical capacity per unit cost, which quantifies the maximal amount of information that can be transmitted per unit cost. I generalize this concept to account for the quantum nature of the information carriers and communication channels and show that if there exists a state with cost equal to zero, e.g. a vacuum state, the capacity per unit cost can be expressed by a simple formula containing maximization of the relative entropy between two quantum states. This enables me to analyze the behavior of photon information efficiency for general communication tasks and show simple bounds on the capacity per unit cost in terms of quantities familiar from quantum estimation theory. I calculate also the capacity per unit cost for general Gaussian quantum channels.


title: Pairs Trading under Drift Uncertainty and Risk Penalization

abstract: In this work, we study a dynamic portfolio optimization problem related to pairs trading, which is an investment strategy that matches a long position in one security with a short position in another security with similar characteristics. The relationship between pairs, called a spread, is modeled by a Gaussian mean-reverting process whose drift rate is modulated by an unobservable continuous-time, finite-state Markov chain. Using the classical stochastic filtering theory, we reduce this problem with partial information to the one with full information and solve it for the logarithmic utility function, where the terminal wealth is penalized by the riskiness of the portfolio according to the realized volatility of the wealth process. We characterize optimal dollar-neutral strategies as well as optimal value functions under full and partial information and show that the certainty equivalence principle holds for the optimal portfolio strategy. Finally, we provide a numerical analysis for a toy example with a two-state Markov chain.


title: An ab initio Local Exchange Approach to the Calculations of Tc of Dilute Magnetic Semiconductors

abstract: The DFT calculations were performed of densities of states of semiconductor GaAs and magnetic semiconductor Ga0.9375Mn0.0625As. It is obtained that a part of Mn3d- states is hybridized with valence band at Fermi level. The exchange integrals of Anderson impurity model were calculated making use of atomic Hartree-Fock package and angular momentum coupling technique. Theoretical Tc of Ga0.9375Mn0.0625As obtained in the multiscale ab initio method is in reasonable agreement with experiment. The application of Hubbard parameters in DFT calculations using ultrasoft pseudopotentials is discussed.


title: Spin-lattice relaxation of magnetic centers in molecular crystals at low temperature

abstract: We study the spin-phonon relaxation rate of both Kramers and non-Kramers molecular magnets in strongly diluted samples at low temperature. Using the "rotational" contribution to the spin-phonon Hamiltonian, universal formulae for the relaxation rate are obtained. Intriguingly, these formulae are all entirely expressed via measurable or \emph{ab initio} computable physical quantities. Moreover, they are also independent of the energy gaps to excited states involved in the relaxation process. These obtained expressions for direct and Raman processes offer an easy way to determine the lowest limit of the spin-phonon relaxation of any spin system based on magnetic properties of the ground doublet only. In addition, some intriguing properties of Raman process are also found. Particularly, Raman process in Kramers system is found dependent on the magnetic field's orientation but independent of its magnitude, meanwhile the same process in non-Kramers system is significantly reduced out of resonance, i.e. for an applied external field. Interestingly, Raman process is demonstrated to vary as $T^{9}$ for both systems. Application of the theory to a recently investigated cobalt(II) complex shows that it can provide a reasonably good description for the relaxation. Based on these findings, a strategy in developing efficient single-molecule magnets by enhancing the mechanical rigidity of the molecular unit is proposed.


title: On Improving Research Methodology Course at Blekinge Institute of Technology

abstract: The Research Methodology in Software Engineering and Computer Science (RM) is a compulsory course that must be studied by graduate students at Blekinge Institute of Technology (BTH) prior to undertaking their theses work. The course is focused on teaching research methods and techniques for data collection and analysis in the fields of Computer Science and Software Engineering. It is intended that the course should help students in practically applying appropriate research methods in different courses (in addition to the RM course) including their Master's theses. However, it is believed that there exist deficiencies in the course due to which the course implementation (learning and assessment activities) as well as the performance of different participants (students, teachers, and evaluators) are affected negatively. In this article our aim is to investigate potential deficiencies in the RM course at BTH in order to provide a concrete evidence on the deficiencies faced by students, evaluators, and teachers in the course. Additionally, we suggest recommendations for resolving the identified deficiencies. Our findings gathered through semi-structured interviews with students, teachers, and evaluators in the course are presented in this article. By identifying a total of twenty one deficiencies from different perspectives, we found that there exist critical deficiencies at different levels within the course. Furthermore, in order to overcome the identified deficiencies, we suggest seven recommendations that may be implemented at different levels within the course and the study program. Our suggested recommendations, if implemented, will help in resolving deficiencies in the course, which may lead to achieving an improved teaching and learning in the RM course at BTH.


title: The evolution of CNO isotopes: a new window on cosmic star-formation history and the stellar IMF in the age of ALMA

abstract: We use state-of-the-art chemical models to track the cosmic evolution of the CNO isotopes in the interstellar medium (ISM) of galaxies, yielding powerful constraints on their stellar initial mass function (IMF). We re-assess the relative roles of massive stars, asymptotic giant branch (AGB) stars and novae in the production of rare isotopes such as 13C, 15N, 17O and 18O, along with 12C, 14N and 16O. The CNO isotope yields of super-AGB stars, novae and fast-rotating massive stars are included. Having reproduced the available isotope enrichment data in the solar neighbourhood, and across the Galaxy, and having assessed the sensitivity of our models to the remaining uncertainties, e.g. nova yields and star-formation history, we show that we can meaningfully constrain the stellar IMF in galaxies using C, O and N isotope abundance ratios. In starburst galaxies, where data for multiple isotopologue lines are available, we find compelling new evidence for a top-heavy stellar IMF, with profound implications for their star-formation rates and efficiencies, perhaps also their stellar masses. Neither chemical fractionation nor selective photodissociation can significantly perturb globally-averaged isotopologue abundance ratios away from the corresponding isotope ones, as both these processes will typically affect only small mass fractions of molecular clouds in galaxies. Thus the Atacama Large Millimetre Array now stands ready to probe the stellar IMF, and even the ages of specific starburst events in star-forming galaxies across cosmic time unaffected by the dust obscuration effects that plague optical/near-infrared studies.


title: DFT analysis and FDTD simulation of CH3NH3PbI3-xClx mixed halide perovskite solar cells: role of halide mixing and light trapping technique

abstract: Since perovskite solar cells have attracted a lot of attentions over the past years, the enhancement of their optical absorption and current density are among the basic coming challenges. For this reason, first, we have studied structural and optical properties of organic-inorganic hybrid halide perovskite CH3NH3PbI3 and the compounds doped by chlorine halogen CH3NH3PbI3-xClx in the cubic phase by using density functional theory (DFT). Then, we investigate the light absorption efficiency and optical current density of the single-junction perovskite solar cell CH3NH3PbI3-xClx for the values (x=0,1,2,3) by utilizing optical constants (n,k) resulted from these calculations. The results suggest that increasing the amount of chlorine in the CH3NH3PbI3-xClx compound leads an increase in bandgap energy, as well as a decrease in lattice constants and optical properties like refractive index and extinction coefficient of the structure. Also, the results obtained by simulation express that by taking advantage of light trapping techniques of SiO2, a remarkable increase of light absorption will be achieved to the magnitude of 83.13%, which is noticeable.


title: Liquid-liquid transition revealed by quasi-static cooling of an ultra-viscous metallic liquid

abstract: Temperature-driven polyamorphism has been reported in various supercooled liquids and glasses. The dynamical and structural routes followed by the system during such crossovers are however not universal and appear to be related to intrinsic kinetic properties of the liquid. By combining x-ray photon correlation spectroscopy and high energy x-ray diffraction, we have followed the collective atomic motion and structural changes during quasi-static cooling of the Au49Cu26.9Si16.3Ag5.5Pd2.3 metallic glass-former. Due to this ultra-slow thermal protocol, the glass transition temperature is lowered far enough to reveal a liquid-liquid crossover in the ultra-viscous supercooled phase. This transition is usually hidden by vitrification and leads to a strong liquid characterized by increasing correlation at the medium range.


title: Energy consumption for Shortcuts To Adiabaticity

abstract: Shortcuts to adiabaticity let a system reach the results of a slow adiabatic process in a shorter time. We propose to quantify the "energy cost" of the shortcut by the energy consumption of the system enlarged by including the control device. A mechanical model where the dynamics of the system and control device can be explicitly described illustrates that a broad range of possible values for the consumption are possible, including zero (above the adiabatic energy increment) when friction is negligible and the energy given away as negative power is stored and recovered by perfect regenerative braking.


title: Boson star with particle size effects

abstract: A simple model to study boson stars is to consider these stellar objects as quantum systems of $N$ identical self-gravitating particles within a non-relativistic framework. Some results obtained with point-like particles are recalled as well as the validity limits of this model. Approximate analytical calculations are performed using envelope theory for a truncated Coulomb-like potential simulating a particle size. If the boson mass is sufficiently small, the description of small mass boson stars is possible within non-relativistic formalism. The mass and radius of these stellar objects are strongly dependent on the value of the truncation parameter.


title: Mechanical properties of carbon nanotube reinforced polymer nanocomposites: A coarse-grained model

abstract: In this work, a coarse-grained (CG) model of carbon nanotube (CNT) reinforced polymer matrix composites is developed. A distinguishing feature of the CG model is the ability to capture interactions between polymer chains and nanotubes. The CG potentials for nanotubes and polymer chains are calibrated using the strain energy conservation between CG models and full atomistic systems. The applicability and efficiency of the CG model in predicting the elastic properties of CNT/polymer composites are evaluated through verification processes with molecular simulations. The simulation results reveal that the CG model is able to estimate the mechanical properties of the nanocomposites with high accuracy and low computational cost. The effect of the volume fraction of CNT reinforcements on the Young's modulus of the nanocomposites is investigated. The application of the method in the modeling of large unit cells with randomly distributed CNT reinforcements is examined. The established CG model will enable the simulation of reinforced polymer matrix composites across a wide range of length scales from nano to mesoscale.


title: A population genetic interpretation of GWAS findings for human quantitative traits

abstract: GWAS in humans are revealing the genetic architecture of biomedical and anthropomorphic traits, i.e., the frequencies and effect sizes of variants that contribute to heritable variation in a trait. To interpret these findings, we need to understand how genetic architecture is shaped by basic population genetics processes-notably, by mutation, natural selection and genetic drift. Because many quantitative traits are subject to stabilizing selection and genetic variation that affects one trait often affects many others, we model the genetic architecture of a focal trait that arises under stabilizing selection in a multi-dimensional trait space. We solve the model for the phenotypic distribution and allelic dynamics at steady state and derive robust, closed form solutions for summary statistics of the genetic architecture. Our results provide a simple interpretation for missing heritability and why it varies among traits. They also predict that the distribution of variances contributed by loci identified in GWAS is well approximated by a simple functional form that depends on a single parameter: the expected contribution to genetic variance of a strongly selected site affecting the trait. We test this prediction against the results of GWAS for height and BMI and find that it fits the data well, allowing us to make inferences about the degree of pleiotropy and mutational target size for these traits. Our findings help to explain why the GWAS for height explains more of the heritable variance than similarly-sized GWAS for BMI, and to predict the increase in explained heritability with study sample size. Considering the demographic history of European populations, in which these GWAS were performed, we further find that most of the associations they identified likely involve mutations that arose during the out of Africa bottleneck at sites with selection coefficients around $s=10^{-3}$.


title: Triangular matrix coalgebras and applications

abstract: We study generalized comatrix coalgebras and upper triangular comatrix coalgebras, which are not only a dualization but also an extension of classical generalized matrix algebras. We use these to answer several questions on Noetherian and Artinian type notions in the theory of coalgebras, and to give complete connections between these. We also solve completely the so called finite splitting problem for coalgebras: we show that a coalgebra $C$ has the property that the rational part of every finitely generated left $C^*$-module splits off if and only if $C$ has the form $C=\left(\begin{array}{cc} D & M \\ 0 & E \end{array}\right)$, an upper triangular matrix coalgebra, for a serial coalgebra $D$ whose Ext-quiver is a finite union of cycles, a finite dimensional coalgebra $E$ and a finite dimensional $D$-$E$-bicomodule $M$.


title: Tensile fracture behavior of short carbon nanotube reinforced polymer composites: A coarse-grained model

abstract: Short-fiber-reinforced polymer composites are increasingly used in engineering applications and industrial products owing to their unique combination of superior mechanical properties, and relatively easy and low cost manufacturing process. The mechanical behavior of short carbon nanotube (CNT) polymer composites, however, remains poorly understood due to size and time limitations of experiments and atomistic simulations. To address this issue, the tensile fracture behavior of short CNT reinforced poly (methyl methacrylate) (PMMA) matrix composites is investigated using a coarse-grained (CG) model. The reliability of the CG model is demonstrated by reproducing experimental results on the stress-stain behavior of the polymer material. The effect of the nanotube weight fraction on the mechanical properties, i.e. the Young's modulus, yield strength,tensile strength and critical strain, of the CNT/polymer composites is studied in detail. The dependence of the mechanical properties of the composites on the orientation and length-to-diameter aspect ratio of nanotube reinforcements is also examined.

categories: physics.app-ph cond-mat.mtrl-sci cond-mat.soft physics.comp-ph


title: Continuous monitoring of $\ell_p$ norms in data streams

abstract: In insertion-only streaming, one sees a sequence of indices $a_1, a_2, \ldots, a_m\in [n]$. The stream defines a sequence of $m$ frequency vectors $x^{(1)},\ldots,x^{(m)}\in\mathbb{R}^n$ with $(x^{(t)})_i = |\{j : j\in[t], a_j = i\}|$. That is, $x^{(t)}$ is the frequency vector after seeing the first $t$ items in the stream. Much work in the streaming literature focuses on estimating some function $f(x^{(m)})$. Many applications though require obtaining estimates at time $t$ of $f(x^{(t)})$, for every $t\in[m]$. Naively this guarantee is obtained by devising an algorithm with failure probability $\ll 1/m$, then performing a union bound over all stream updates to guarantee that all $m$ estimates are simultaneously accurate with good probability. When $f(x)$ is some $\ell_p$ norm of $x$, recent works have shown that this union bound is wasteful and better space complexity is possible for the continuous monitoring problem, with the strongest known results being for $p=2$ [HTY14, BCIW16, BCINWW17]. In this work, we improve the state of the art for all $0<p<2$, which we obtain via a novel analysis of Indyk's $p$-stable sketch [Indyk06].


title: Band-to-band transitions, selection rules, effective mass and exciton binding energy parameters in monoclinic \beta-Ga2O3

abstract: We employ an eigen polarization model including the description of direction dependent excitonic effects for rendering critical point structures within the dielectric function tensor of monoclinic \beta-Ga2O3 yielding a comprehensive analysis of generalized ellipsometry data obtained from 0.75 eV--9 eV. The eigen polarization model permits complete description of the dielectric response, and we obtain single-electron and excitonic band-to-band transition anisotropic critical point structure model parameters including their polarization eigenvectors within the monoclinic lattice. We compare our experimental analysis with results from density functional theory calculations performed using a recently proposed Gaussian-attenuation-Perdue-Burke-Ernzerhof hybrid density functional, and we present and discuss the order of the fundamental direct band-to-band transitions and their polarization selection rules, the electron and hole effective mass parameters for the three lowest band-to-band transitions, and their exciton binding energy parameters, in excellent agreement with our experimental results. We find that the effective masses for holes are highly anisotropic and correlate with the selection rules for the fundamental band-to-band transitions, where the observed transitions are polarized closely in the direction of the lowest hole effective mass for the valence band participating in the transition.


title: Charge asymmetry dependence of anisotropic flow in pPb and PbPb collisions with CMS experiment

abstract: In nucleus-nucleus collisions, the linear dependence found for the elliptic flow harmonic of both positive or negative charged particles as a function of event charge asymmetry is predicted by the phenomenon known as the Chiral Magnetic Wave (CMW) due to its induced electric quadrupole moment. Here, the event charge asymmetry $A_{\rm ch}$ is defined as $\frac{N_{+}-N_{-}}{N_{+}+N_{-}}$, where $N_{+}$ and $N_{-}$ are the number of positive and negative charged particles, respectively. However, other scenarios are also possible and may provide alternative explanations for the experimental results. New measurements of elliptic ($v_{\rm 2}$) and triangular ($v_{\rm 3}$) flow for positive and negative charged particles as a function of $A_{\rm ch}$ in pPb and PbPb collisions at $\sqrt{s_{_{NN}}} = 5.02~\mathrm{TeV}$ are presented, using data collected by the CMS experiment during the LHC runs 1 and 2. The slopes and intercepts of the charged-dependent $v_{n}$ harmonics vs. $A_{\rm ch}$ are directly compared for pPb and PbPb collisions with similar charged-particle multiplicities, where a strong CMW effect is not expected in very high multiplicity pPb events. Moreover, a comparison is made of the slope parameters between $v_{2}$ and $v_{3}$ harmonics normalized by the inclusive charge particle $v_{n}$ in PbPb collisions as a function of centrality. These results provide a means to discriminate between the CMW and other scenarios such as local charge conservation as possible explanations for the observed charge dependent behavior.


title: A dynamic resource allocation decision model for IT security

abstract: Today, with the continued growth in using information and communication technologies (ICT) for business purposes, business organizations become increasingly dependent on their information systems. Thus, they need to protect them from the different attacks exploiting their vulnerabilities. To do so, the organization has to use security technologies, which may be proactive or reactive ones. Each security technology has a relative cost and addresses specific vulnerabilities. Therefore, the organization has to put in place the appropriate security technologies set that minimizes the information system s vulnerabilities with a minimal cost. This bi objective problem will be considered as a resources allocation problem (RAP) where security technologies represent the resources to be allocated. However, the set of vulnerabilities may change, periodically, with the continual appearance of new ones. Therefore, the security technologies set should be flexible to face these changes, in real time, and the problem becomes a dynamic one. In this paper, we propose a harmony search based algorithm to solve the bi objective dynamic resource allocation decision model. This approach was compared to a genetic algorithm and provided good results.


title: Magneto-Optical Detection of the Spin Hall Effect in Pt and W Thin Films

abstract: Here we report the measurement of the interfacial spin accumulation induced by the spin Hall effect in Pt and W thin films using magneto-optical Kerr microscopy. We show that the Kerr rotation has opposite sign in Pt and W and scales linearly with current density. By comparing the experimental results with ab-initio calculations of the spin Hall and magneto-optical Kerr effects, we quantitatively determine the current-induced spin accumulation at the Pt interface as $5*10^{-12} \mu_B$A$^{-1}$cm$^2$ per atom. From thickness-dependent measurements, we determine the spin diffusion length in a single Pt film to be $11 \pm 3$ nm, which is significantly larger compared to that of Pt adjacent to a magnetic layer.


title: Weighted quasisymmetric enumerator for generalized permutohedra

abstract: We introduce a weighted quasisymmetric enumerator function associated to generalized permutohedra. It refines the Billera, Jia and Reiner quasisymmetric function which also includes the Stanley chromatic symmetric function. Beside that it carries information of face numbers of generalized permutohedra. We consider more systematically the cases of nestohedra and matroid base polytopes.


title: Service Chain (SC) Mapping with Multiple SC Instances in a Wide Area Network

abstract: Network Function Virtualization (NFV) aims to simplify deployment of network services by running Virtual Network Functions (VNFs) on commercial off-the-shelf servers. Service deployment involves placement of VNFs and in-sequence routing of traffic flows through VNFs comprising a Service Chain (SC). The joint VNF placement and traffic routing is usually referred as SC mapping. In a Wide Area Network (WAN), a situation may arise where several traffic flows, generated by many distributed node pairs, require the same SC, one single instance (or occurrence) of that SC might not be enough. SC mapping with multiple SC instances for the same SC turns out to be a very complex problem, since the sequential traversal of VNFs has to be maintained while accounting for traffic flows in various directions. Our study is the first to deal with SC mapping with multiple SC instances to minimize network resource consumption. Exact mathematical modeling of this problem results in a quadratic formulation. We propose a two-phase column-generation-based model and solution in order to get results over large network topologies within reasonable computational times. Using such an approach, we observe that an appropriate choice of only a small set of SC instances can lead to solution very close to the minimum bandwidth consumption.


title: Rattleback dynamics and its reversal time of rotation

abstract: A rattleback is a rigid, semi-elliptic toy which exhibits unintuitive behavior; when it is spun in one direction, it soon begins pitching and stops spinning, then it starts to spin in the opposite direction, but in the other direction, it seems to spin just steadily. This puzzling behavior results from the slight misalignment between the principal axes for the inertia and those for the curvature; the misalignment couples the spinning with the pitching and the rolling oscillations. It has been shown that under the no-slip condition and without dissipation the spin can reverse in both directions, and Garcia and Hubbard obtained the formula for the time required for the spin reversal $t_r$ [Proc. R. Soc. Lond. A 418, 165-197 (1988) ]. In this work, we reformulate the rattleback dynamics in a physically transparent way and reduce it to a three-variable dynamics for spinning, pitching, and rolling. We obtain an expression of the Garcia-Hubbard formula for $t_r$ by a simple product of four factors: (1) the misalignment angle, (2) the difference in the inverses of inertia moment for the two oscillations, (3) that in the radii for the two principal curvatures, and (4) the squared frequency of the oscillation. We perform extensive numerical simulations to examine validity and limitation of the formula, and find that (1) the Garcia-Hubbard formula is good for both spinning directions in the small spin and small oscillation regime, but (2) in the fast spin regime especially for the steady direction, the rattleback may not reverse and shows a rich variety of dynamics including steady spinning, spin wobbling, and chaotic behavior reminiscent of chaos in a dissipative system.


title: Hierarchical Bayesian Data Fusion for Robotic Platform Navigation

abstract: Data fusion has become an active research topic in recent years. Growing computational performance has allowed the use of redundant sensors to measure a single phenomenon. While Bayesian fusion approaches are common in general applications, the computer vision field has largely relegated this approach. Most object following algorithms have gone towards pure machine learning fusion techniques that tend to lack flexibility. Consequently, a more general data fusion scheme is needed. Within this work, a hierarchical Bayesian fusion approach is proposed, which outperforms individual trackers by using redundant measurements. The adaptive framework is achieved by relying on each measurement's local statistics and a global softened majority voting. The proposed approach was validated in a simulated application and two robotic platforms.


title: The neutron skin thickness of $^{48}$Ca from a nonlocal-dispersive-optical-model analysis

abstract: A nonlocal dispersive-optical-model analysis has been carried out for neutrons and protons in $^{48}$Ca. Elastic-scattering angular distributions, total and reaction cross sections, single-particle energies, the neutron and proton numbers, and the charge distribution have been fitted to extract the neutron and proton self-energies both above and below the Fermi energy. From the single-particle propagator resulting from these self-energies we have determined the charge and neutron matter distributions in $^{48}$Ca. A neutron skin of 0.249$\pm$0.023~fm is deduced. The energy dependence of the total neutron cross sections is shown to have strong sensitivity to the skin thickness.


title: Schwarz lemma and Kobayashi metrics for holomorphic and pluriharmonic functions

abstract: The Schwarz lemma as one of the most influential results in complex analysis and it has a great impact to the development of several research fields, such as geometric function theory, hyperbolic geometry, complex dynamical systems, and theory of quasi-conformal mappings. In this note we mainly consider various version of Schwarz lemma and its relatives related to holomorphic functions including several variables.


title: On the complexity of non-orientable Seifert fibre spaces

abstract: In this paper we deal with Seifert fibre spaces, which are compact 3-manifolds admitting a foliation by circles. We give a combinatorial description for these manifolds in all the possible cases: orientable, non-orientable, closed, with boundary. Moreover, we compute a potentially sharp upper bound for their complexity in terms of the invariants of the combinatorial description, extending to the non-orientable case results by Fominykh and Wiest for the orientable case with boundary and by Martelli and Petronio for the closed orientable case.


title: On a theorem of Terzio\u{g}lu

abstract: The theory of compact linear operators acting on a Banach space has such a classical core and is familiar to many. Perhaps lesser known is the characterization theorem of Terzio\u{g}lu for compact maps. In this paper we consider Terzio\u{g}lu's theorem and its consequences. We also give a similar characterization theorem in case where there is an approximation scheme on the Banach space.


title: Non-equilibrium transitions in multiscale systems with a bifurcating slow manifold

abstract: Noise-induced transitions between metastable fixed points in systems evolving on multiple time scales are analyzed in situations where the time scale separation gives rise to a slow manifold with bifurcation. This analysis is performed within the realm of large deviation theory. It is shown that these non-equilibrium transitions make use of a reaction channel created by the bifurcation structure of the slow manifold, leading to vastly increased transition rates. Several examples are used to illustrate these findings, including an insect outbreak model, a system modeling phase separation in the presence of evaporation, and a system modeling transitions in active matter self-assembly. The last example involves a spatially extended system modeled by a stochastic partial differential equation.


title: Complexity Analysis of the Parallel Guided Ejection Search for the Pickup and Delivery Problem with Time Windows

abstract: This paper presents the pessimistic time complexity analysis of the parallel algorithm for minimizing the fleet size in the pickup and delivery problem with time windows. We show how to estimate the pessimistic complexity step by step. This approach can be easily adopted to other parallel algorithms for solving complex transportation problems.


title: Inhomogeneous fluid of penetrable-spheres: application of the random phase approximation

abstract: The focus of the present work is the application of the random phase approximation (RPA), derived for inhomogeneous fluids [Frydel and Ma, Phys. Rev. E 93, 062112 (2016)], to penetrable-spheres. As penetrable-spheres transform into hard-spheres with increasing interactions, they provide an interesting case for exploring the RPA, its shortcomings, and limitations, the weak- versus the strong-coupling limit. Two scenarios taken up by the present study are a one-component and a two-component fluid with symmetric interactions. In the latter case, the mean-field contributions cancel out and any contributions from particle interactions are accounted for by correlations. The accuracy of the RPA for this case is the result of a somewhat lucky cancellation of errors.


title: Distant Supervision for Topic Classification of Tweets in Curated Streams

abstract: We tackle the challenge of topic classification of tweets in the context of analyzing a large collection of curated streams by news outlets and other organizations to deliver relevant content to users. Our approach is novel in applying distant supervision based on semi-automatically identifying curated streams that are topically focused (for example, on politics, entertainment, or sports). These streams provide a source of labeled data to train topic classifiers that can then be applied to categorize tweets from more topically-diffuse streams. Experiments on both noisy labels and human ground-truth judgments demonstrate that our approach yields good topic classifiers essentially "for free", and that topic classifiers trained in this manner are able to dynamically adjust for topic drift as news on Twitter evolves.


title: Spin Density wave instability in a ferromagnet

abstract: Ferromagnetic (FM) and incommensurate spin-density wave (ISDW) states are an unusual set of competing magnetic orders that are seldom observed in the same material without application of a polarizing magnetic field. We report, for the first time, the discovery of an ISDW state that is derived from a FM ground state through a Fermi surface (FS) instability in Fe$_3$Ga$_4$. This was achieved by combining neutron scattering experiments with first principles simulations. Neutron diffraction demonstrates that Fe$_3$Ga$_4$ is in an ISDW state at intermediate temperatures and that there is a conspicuous re-emergence of ferromagnetism above 360 K. First principles calculations show that the ISDW ordering wavevector is in excellent agreement with a prominent nesting condition in the spin-majority FS demonstrating the discovery of a novel instability for FM metals; ISDW formation due to Fermi surface nesting in a spin-polarized Fermi surface.


title: Charged Einstein-\ae ther black holes in $n$-dimensional spacetime

abstract: In this work, we investigate the $n$-dimensional charged static black hole solutions in the Einstein-\ae ther theory. By taking the metric parameter $k$ to be $1,0$, and $-1$, we obtain the spherical, planar, and hyperbolic spacetimes respectively. Three choices of the cosmological constant, $\Lambda>0$, $\Lambda=0$ and $\Lambda<0$, are investigated, which correspond to asymptotically de Sitter, flat and anti-de Sitter spacetimes. The obtained results show the existence of the universal horizon in higher dimensional cases which may trap any particle with arbitrarily large velocity. We analyze the horizon and the surface gravity of 4- and 5-dimensional black holes, and the relations between the above quantities and the electrical charge. It is shown that when the aether coefficient $c_{13}$ or the charge $Q$ increases, the outer Killing horizon shrinks and approaches the universal horizon. Furthermore, the surface gravity decreases and approaches zero in the limit $c_{13}\rightarrow\infty$ or $Q\rightarrow Q_e$, where $Q_e$ is the extreme charge. The main features of the horizon and surface gravity are found to be similar to those in $n=3$ case, but subtle differences are also observed.


title: On Face Segmentation, Face Swapping, and Face Perception

abstract: We show that even when face images are unconstrained and arbitrarily paired, face swapping between them is actually quite simple. To this end, we make the following contributions. (a) Instead of tailoring systems for face segmentation, as others previously proposed, we show that a standard fully convolutional network (FCN) can achieve remarkably fast and accurate segmentations, provided that it is trained on a rich enough example set. For this purpose, we describe novel data collection and generation routines which provide challenging segmented face examples. (b) We use our segmentations to enable robust face swapping under unprecedented conditions. (c) Unlike previous work, our swapping is robust enough to allow for extensive quantitative tests. To this end, we use the Labeled Faces in the Wild (LFW) benchmark and measure the effect of intra- and inter-subject face swapping on recognition. We show that our intra-subject swapped faces remain as recognizable as their sources, testifying to the effectiveness of our method. In line with well known perceptual studies, we show that better face swapping produces less recognizable inter-subject results. This is the first time this effect was quantitatively demonstrated for machine vision systems.


title: $A_{\alpha}$-spectrum of a graph obtained by copies of a rooted graph and applications

abstract: Given a connected graph $R$ on $r$ vertices and a rooted graph $H,$ let $R\{H\}$ be the graph obtained from $r$ copies of $H$ and the graph $R$ by identifying the root of the $i-th$ copy of $H$ with the $i-th$ vertex of $R$. Let $0\leq\alpha\leq1,$ and let \[ A_{\alpha}(G)=\alpha D(G)+(1-\alpha)A(G) \] where $D(G)$ and $A(G)$ are the diagonal matrix of the vertex degrees of $G$ and the adjacency matrix of $G$, respectively. A basic result on the $A_{\alpha}-$ spectrum of $R\{H\}$ is obtained. This result is used to prove that if $H=B_{k}$ is a generalized Bethe tree on $k$ levels, then the eigenvalues of $A_{\alpha}(R\{B_{k}\})$ are the eigenvalues of symmetric tridiagonal matrices of order not exceeding $k$; additionally, the multiplicity of each eigenvalue is determined. Finally, applications to a unicyclic graph are given, including an upper bound on the $\alpha-$ spectral radius in terms of the largest vertex degree and the largest height of the trees obtained by removing the edges of the unique cycle in the graph.


title: Batch-Expansion Training: An Efficient Optimization Framework

abstract: We propose Batch-Expansion Training (BET), a framework for running a batch optimizer on a gradually expanding dataset. As opposed to stochastic approaches, batches do not need to be resampled i.i.d. at every iteration, thus making BET more resource efficient in a distributed setting, and when disk-access is constrained. Moreover, BET can be easily paired with most batch optimizers, does not require any parameter-tuning, and compares favorably to existing stochastic and batch methods. We show that when the batch size grows exponentially with the number of outer iterations, BET achieves optimal $O(1/\epsilon)$ data-access convergence rate for strongly convex objectives. Experiments in parallel and distributed settings show that BET performs better than standard batch and stochastic approaches.


title: Evolution of texture and microstructure during accumulative roll bonding of aluminum AA5086 alloy

abstract: In the present investigation, a strongly bonded strip of an aluminium-magnesium based alloy AA5086 is successfully produced through accumulative roll bonding (ARB). A maximum of up to eight passes has been used for the purpose. Microstructural characterization using electron backscatter diffraction (EBSD) technique indicates the formation of submicron sized (~200-300 nm) subgrains inside the layered microstructure. The material is strongly textured where individual layers possess typical FCC rolling texture components. More than three times enhancement in 0.2% proof stress (PS) has been obtained after 8 passes due to grain refinement and strain hardening


title: A fast second-order implicit difference method for time-space fractional advection-diffusion equation

abstract: In this paper, we consider a fast and second-order implicit difference method for approximation of a class of time-space fractional variable coefficients advection-diffusion equation. To begin with, we construct an implicit difference scheme, based on $L2-1_{\sigma}$ formula [A. A. Alikhanov, A new difference scheme for the time fractional diffusion equation, \emph{J. Comput. Phys.}, 280 (2015)] for the temporal discretization and weighted and shifted Gr\"{u}nwald method for the spatial discretization. Then, unconditional stability of the implicit difference scheme is proved, and we theoretically and numerically show that it converges in the $L_2$-norm with the optimal order $\mathcal{O}(\tau^2 + h^2)$ with time step $\tau$ and mesh size $h$. Secondly, three fast Krylov subspace solvers with suitable circulant preconditioners are designed to solve the discretized linear systems with the Toeplitz matrix. In each iterative step, these methods reduce the memory requirement of the resulting linear equations from $\mathcal{O}(N^2)$ to $\mathcal{O}(N)$ and the computational complexity from $\mathcal{O}(N^3)$ to $\mathcal{O}(N \log N)$, where $N$ is the number of grid nodes. Finally, numerical experiments are carried out to demonstrate that these methods are more practical than the traditional direct solvers of the implicit difference methods, in terms of memory requirement and computational cost.


title: Comparison of Vacuum Static Quadrupolar Metrics

abstract: We investigate the properties of static and axisymmetric vacuum solutions of Einstein equations which generalize the Schwarzschild spherically symmetric solution to include a quadrupole parameter. We test all the solutions with respect to elementary and asymptotic flatness and curvature regularity. Analyzing their multipole structure, according to the relativistic invariant Geroch definition, we show that all of them are equivalent up to the level of the quadrupole. We conclude that the $q-$metric, a variant of the Zipoy-Voorhees metric, is the simplest generalization of the Schwarzschild metric, containing a quadrupole parameter.


title: Asynchronous Distributed Variational Gaussian Processes for Regression

abstract: Gaussian processes (GPs) are powerful non-parametric function estimators. However, their applications are largely limited by the expensive computational cost of the inference procedures. Existing stochastic or distributed synchronous variational inferences, although have alleviated this issue by scaling up GPs to millions of samples, are still far from satisfactory for real-world large applications, where the data sizes are often orders of magnitudes larger, say, billions. To solve this problem, we propose ADVGP, the first Asynchronous Distributed Variational Gaussian Process inference for regression, on the recent large-scale machine learning platform, PARAMETERSERVER. ADVGP uses a novel, flexible variational framework based on a weight space augmentation, and implements the highly efficient, asynchronous proximal gradient optimization. While maintaining comparable or better predictive performance, ADVGP greatly improves upon the efficiency of the existing variational methods. With ADVGP, we effortlessly scale up GP regression to a real-world application with billions of samples and demonstrate an excellent, superior prediction accuracy to the popular linear models.


title: Parity anomaly in four dimensions

abstract: In an analogy to the odd-dimensional case we define the parity anomaly as the part of the one-loop effective action for fermions associated with spectral asymmetry of the Dirac operator. This quantity is computed directly on four-dimensional manifolds with boundary and related to the Chern-Simons current on the boundary. Despite a quite unusual Chern-Simons level obtained, the action is gauge invariant and passes all consistency checks.


title: Circumcentering the Douglas--Rachford method

abstract: We introduce and study a geometric modification of the Douglas-Rach\-ford method called the Circumcentered-Douglas-Rachford method. This method iterates by taking the intersection of bisectors of reflection steps for solving certain classes of feasibility problems. The convergence analysis is established for best approximation problems involving two (affine) subspaces and both our theoretical and numerical results compare favorably to the original Douglas-Rachford method. Under suitable conditions, it is shown that the linear rate of convergence of the Circumcentered-Douglas-Rachford method is at least the cosine of the Friedrichs angle between the (affine) subspaces, which is known to be the sharp rate for the Douglas-Rachford method. We also present a preliminary discussion on the Circumcentered-Douglas-Rachford method applied to the many set case and to examples featuring non-affine convex sets.


title: Towards Distributed Machine Learning in Shared Clusters: A Dynamically-Partitioned Approach

abstract: Many cluster management systems (CMSs) have been proposed to share a single cluster with multiple distributed computing systems. However, none of the existing approaches can handle distributed machine learning (ML) workloads given the following criteria: high resource utilization, fair resource allocation and low sharing overhead. To solve this problem, we propose a new CMS named Dorm, incorporating a dynamically-partitioned cluster management mechanism and an utilization-fairness optimizer. Specifically, Dorm uses the container-based virtualization technique to partition a cluster, runs one application per partition, and can dynamically resize each partition at application runtime for resource efficiency and fairness. Each application directly launches its tasks on the assigned partition without petitioning for resources frequently, so Dorm imposes flat sharing overhead. Extensive performance evaluations showed that Dorm could simultaneously increase the resource utilization by a factor of up to 2.32, reduce the fairness loss by a factor of up to 1.52, and speed up popular distributed ML applications by a factor of up to 2.72, compared to existing approaches. Dorm's sharing overhead is less than 5% in most cases.


title: Hints on the gradual re-sizing of the torus in AGN by decomposing IRS/Spitzer spectra

abstract: Several authors have claimed that the less luminous active galactic nuclei (AGN) are not capable of sustaining the dusty torus structure. Thus, a gradual re-sizing of the torus is expected when the AGN luminosity decreases. Our aim is to confront mid-infrared observations of local AGN of different luminosities with this scenario. We decomposed about ~100 IRS/Spitzer spectra of LLAGN and powerful Seyferts in order to decontaminate the torus component from other contributors. We have used the affinity propagation (AP) method to cluster the data into five groups within the sample according to torus contribution to the 5-15 um range (Ctorus) and bolometric luminosity. The AP groups show a progressively higher torus contribution and an increase of the bolometric luminosity, from Group 1 (Ctorus~ 0% and logLbol ~ 41) and up to Group 5 (Ctorus ~80% and log(Lbol) ~44). We have fitted the average spectra of each of the AP groups to clumpy models. The torus is no longer present in Group 1, supporting the disappearance at low-luminosities. We were able to fit the average spectra for the torus component in Groups 3 (Ctorus~ 40% and log(Lbol)~ 42.6), 4 (Ctorus~ 60% and log(Lbol)~ 43.7), and 5 to Clumpy torus models. We did not find a good fitting to Clumpy torus models for Group 2 (Ctorus~ 18% and log(Lbol)~ 42). This might suggest a different configuration and/or composition of the clouds for Group 2, which is consistent with a different gas content seen in Groups 1, 2, and 3, according to the detections of H2 molecular lines. Groups 3, 4, and 5 show a trend to decrease of the width of the torus (which yields to a likely decrease of the geometrical covering factor), although we cannot confirm it with the present data. Finally, Groups 3, 4, and 5 show an increase on the outer radius of the torus for higher luminosities, consistent with a re-sizing of the torus according to the AGN luminosity.


title: Multiscale Simulations of Defect Dipole-Enhanced Electromechanical Coupling at Dilute Defect Concentrations

abstract: The role of defects in solids of mixed ionic-covalent bonds such as ferroelectric oxides is complex. Current understanding of defects on ferroelectric properties at the single-defect level remains mostly at the empirical level, and the detailed atomistic mechanisms for many defect-mediated polarization-switching processes have not been convincingly revealed quantum mechanically. We simulate the polarization-electric field (P-E) and strain-electric field ({$\epsilon$}-E) hysteresis loops for BaTiO3 in the presence of generic defect dipoles with large-scale molecular dynamics and provide a detailed atomistic picture of the defect dipole-enhanced electromechanical coupling. We develop a general first-principles-based atomistic model, enabling a quantitative understanding of the relationship between macroscopic ferroelectric properties and dipolar impurities of different orientations, concentrations, and dipole moments. We find that the collective orientation of dipolar defects relative to the external field is the key microscopic structure feature that strongly affects materials hardening/softening and electromechanical coupling. We show that a small concentration (0.1 at.%) of defect dipoles dramatically improves electromechanical response. This offers the opportunity to improve the performance of inexpensive polycrystalline ferroelectric ceramics through defect dipole engineering for a range of applications including piezoelectric sensors, actuators, and transducers.


title: Decoupled field integral equations for electromagnetic scattering from homogeneous penetrable obstacles

abstract: We present a new method for the analysis of electromagnetic scattering from homogeneous penetrable bodies. Our approach is based on a reformulation of the governing Maxwell equations in terms of two uncoupled vector Helmholtz systems: one for the electric feld and one for the magnetic field. This permits the derivation of resonance-free Fredholm equations of the second kind that are stable at all frequencies, insensitive to the genus of the scatterers, and invertible for all passive materials including those with negative permittivities or permeabilities. We refer to these as decoupled field integral equations.


title: Testing Network Structure Using Relations Between Small Subgraph Probabilities

abstract: We study the problem of testing for structure in networks using relations between the observed frequencies of small subgraphs. We consider the statistics \begin{align*} T_3 & =(\text{edge frequency})^3 - \text{triangle frequency}\\ T_2 & =3(\text{edge frequency})^2(1-\text{edge frequency}) - \text{V-shape frequency} \end{align*} and prove a central limit theorem for $(T_2, T_3)$ under an Erd\H{o}s-R\'{e}nyi null model. We then analyze the power of the associated $\chi^2$ test statistic under a general class of alternative models. In particular, when the alternative is a $k$-community stochastic block model, with $k$ unknown, the power of the test approaches one. Moreover, the signal-to-noise ratio required is strictly weaker than that required for community detection. We also study the relation with other statistics over three-node subgraphs, and analyze the error under two natural algorithms for sampling small subgraphs. Together, our results show how global structural characteristics of networks can be inferred from local subgraph frequencies, without requiring the global community structure to be explicitly estimated.


title: Robust, Deep and Inductive Anomaly Detection

abstract: PCA is a classical statistical technique whose simplicity and maturity has seen it find widespread use as an anomaly detection technique. However, it is limited in this regard by being sensitive to gross perturbations of the input, and by seeking a linear subspace that captures normal behaviour. The first issue has been dealt with by robust PCA, a variant of PCA that explicitly allows for some data points to be arbitrarily corrupted, however, this does not resolve the second issue, and indeed introduces the new issue that one can no longer inductively find anomalies on a test set. This paper addresses both issues in a single model, the robust autoencoder. This method learns a nonlinear subspace that captures the majority of data points, while allowing for some data to have arbitrary corruption. The model is simple to train and leverages recent advances in the optimisation of deep neural networks. Experiments on a range of real-world datasets highlight the model's effectiveness.


title: Reducibility of quantum harmonic oscillator on $ R^d$ with differential and quasi-periodic in time potential

abstract: We improve the results by Gr\'ebert and Paturel in \cite{GP} and prove that a linear Schr\"odinger equation on $R^d$ with harmonic potential $|x|^2$ and small $t$-quasiperiodic potential as $$ {\rm i}u_t - \Delta u+|x|^2u+\varepsilon V(\omega t,x)u=0, \ (t,x)\in R\times R^d $$ reduces to an autonomous system for most values of the frequency vector $\omega\in R^n$. The new point is that the potential $V(\theta,\cdot )$ is only in ${\mathcal{C}^{\beta}}(T^n, \mathcal{H}^{s}(R^d))$ with $\beta$ large enough. As a consequence any solution of such a linear PDE is almost periodic in time and remains bounded in some suitable Sobolev norms.