title: Scheduling Monotone Moldable Jobs in Linear Time

abstract: A moldable job is a job that can be executed on an arbitrary number of processors, and whose processing time depends on the number of processors allotted to it. A moldable job is monotone if its work doesn't decrease for an increasing number of allotted processors. We consider the problem of scheduling monotone moldable jobs to minimize the makespan. We argue that for certain compact input encodings a polynomial algorithm has a running time polynomial in n and log(m), where n is the number of jobs and m is the number of machines. We describe how monotony of jobs can be used to counteract the increased problem complexity that arises from compact encodings, and give tight bounds on the approximability of the problem with compact encoding: it is NP-hard to solve optimally, but admits a PTAS. The main focus of this work are efficient approximation algorithms. We describe different techniques to exploit the monotony of the jobs for better running times, and present a (3/2+{\epsilon})-approximate algorithm whose running time is polynomial in log(m) and 1/{\epsilon}, and only linear in the number n of jobs.


title: A Multiple Data Source Framework for the Identification of Activities of Daily Living Based on Mobile Device Data

abstract: Most mobile devices include motion, magnetic, acoustic, and location sensors. They allow the implementation of a framework for the recognition of Activities of Daily Living (ADL) and its environments, composed by the acquisition, processing, fusion, and classification of data. This study compares different implementations of artificial neural networks, concluding that the obtained results were 85.89% and 100% for the recognition of standard ADL. Additionally, for the identification of standing activities with Deep Neural Networks (DNN) respectively, and 86.50% for the identification of the environments with Feedforward Neural Networks. Numerical results illustrate that the proposed framework can achieve robust performance from the data fusion of off-the-shelf mobile devices.


title: Ejecta Mass Diagnostics of Type Ia Supernovae

abstract: We present one-dimensional non-local thermodynamic equilibrium time-dependent radiative transfer simulations (using CMFGEN) of two sub-Chandrasekhar (sub-$M_{\rm Ch}$), one $M_{\rm Ch}$ and one super-$M_{\rm Ch}$ Type Ia SN ejecta models. Three originate from $M_{\rm Ch}$ delayed detonation models, and the fourth is a sub-$M_{\rm Ch}$ detonation model. Ejecta masses are 1.02, 1.04, 1.40, and 1.70 M$_\odot$, and all models have 0.62 M$_\odot$ of $^{56}{\rm Ni}$. Sub-$M_{\rm Ch}$ model light curves evolve faster, reaching bolometric maximum 2--3 days earlier and having 3--4 days shorter bolometric half light widths. The models vary by $\sim$12 per cent at maximum bolometric luminosity and by 0.17 mag in $B_{\rm max}$. While $\Delta M_{15}(B)$ increases with ejecta mass it only varies by $\sim$5 per cent around 1 mag. Sub-$M_{\rm Ch}$ models are 0.25 mag bluer in $B-R$ at $B_{\rm max}$. Optical spectra share many similarities, but lower mass models exhibit less UV line blanketing during the photospheric phase. At nebular times, significant NIR spectroscopic differences are seen. In particular, emission lines of the Ca II NIR triplet; [S III] $\lambda\lambda$9068,9530; [Ca II] $\lambda\lambda$7291,7324; [Ar III] $\lambda\lambda$7135,7751; and [Ni II] 1.939 $\mu$m are stronger in higher mass models. The [Ni II] 1.939 $\mu$m line is absent in the sub-$M_{\rm Ch}$ detonation model, and provides a valuable potential tool to distinguish sub-$M_{\rm Ch}$ explosions from $M_{\rm Ch}$ explosions. In general, the nebular phase models are too highly ionized. We attribute this to the neglect of clumping and/or the distribution of intermediate mass and iron group elements. The two sub-$M_{\rm Ch}$ models, while exploded by different mechanisms, can be distinguished in the $J$ and $H$ bands at late times (e.g., $+200$ days).


title: DCN+: Mixed Objective and Deep Residual Coattention for Question Answering

abstract: Traditional models for question answering optimize using cross entropy loss, which encourages exact answers at the cost of penalizing nearby or overlapping answers that are sometimes equally accurate. We propose a mixed objective that combines cross entropy loss with self-critical policy learning. The objective uses rewards derived from word overlap to solve the misalignment between evaluation metric and optimization objective. In addition to the mixed objective, we improve dynamic coattention networks (DCN) with a deep residual coattention encoder that is inspired by recent work in deep self-attention and residual networks. Our proposals improve model performance across question types and input lengths, especially for long questions that requires the ability to capture long-term dependencies. On the Stanford Question Answering Dataset, our model achieves state-of-the-art results with 75.1% exact match accuracy and 83.1% F1, while the ensemble obtains 78.9% exact match accuracy and 86.0% F1.


title: Separation of Water and Fat Magnetic Resonance Imaging Signals Using Deep Learning with Convolutional Neural Networks

abstract: Purpose: A new method for magnetic resonance (MR) imaging water-fat separation using a convolutional neural network (ConvNet) and deep learning (DL) is presented. Feasibility of the method with complex and magnitude images is demonstrated with a series of patient studies and accuracy of predicted quantitative values is analyzed. Methods: Water-fat separation of 1200 gradient-echo acquisitions from 90 imaging sessions (normal, acute and chronic myocardial infarction) was performed using a conventional model based method with modeling of R2* and off-resonance and a multi-peak fat spectrum. A U-Net convolutional neural network for calculation of water-only, fat-only, R2* and off-resonance images was trained with 900 gradient-echo Multiple and single-echo complex and magnitude input data algorithms were studied and compared to conventional extended echo modeling. Results: The U-Net ConvNet was easily trained and provided water-fat separation results visually comparable to conventional methods. Myocardial fat deposition in chronic myocardial infarction and intramyocardial hemorrhage in acute myocardial infarction were well visualized in the DL results. Predicted values for R2*, off-resonance, water and fat signal intensities were well correlated with conventional model based water fat separation (R2>=0.97, p<0.001). DL images had a 14% higher signal-to-noise ratio (p<0.001) when compared to the conventional method. Conclusion: Deep learning utilizing ConvNets is a feasible method for MR water-fat separationimaging with complex, magnitude and single echo image data. A trained U-Net can be efficiently used for MR water-fat separation, providing results comparable to conventional model based methods.


title: Beyond Shared Hierarchies: Deep Multitask Learning through Soft Layer Ordering

abstract: Existing deep multitask learning (MTL) approaches align layers shared between tasks in a parallel ordering. Such an organization significantly constricts the types of shared structure that can be learned. The necessity of parallel ordering for deep MTL is first tested by comparing it with permuted ordering of shared layers. The results indicate that a flexible ordering can enable more effective sharing, thus motivating the development of a soft ordering approach, which learns how shared layers are applied in different ways for different tasks. Deep MTL with soft ordering outperforms parallel ordering methods across a series of domains. These results suggest that the power of deep MTL comes from learning highly general building blocks that can be assembled to meet the demands of each task.


title: Spin susceptibility of charge ordered YBa2Cu3Oy across the upper critical field

abstract: The value of the upper critical field Hc2, a fundamental characteristic of the superconducting state, has been subject to strong controversy in high-Tc copper-oxides. Since the issue has been tackled almost exclusively by macroscopic techniques so far, there is a clear need for local-probe measurements. Here, we use 17O NMR to measure the spin susceptibility $\chi_{spin}$ of the CuO2 planes at low temperature in charge ordered YBa2Cu3Oy. We find that $\chi_{spin}$ increases (most likely linearly) with magnetic field H and saturates above field values ranging from 20 to 40 T. This result is consistent with Hc2 values claimed by G. Grissonnanche et al. [Nat. Commun. 5, 3280 (2014)] and with the interpretation that the charge-density-wave (CDW) reduces Hc2 in underdoped YBa2Cu3Oy. Furthermore, the absence of marked deviation in $\chi_{spin}(H)$ at the onset of long-range CDW order indicates that this Hc2 reduction and the Fermi-surface reconstruction are primarily rooted in the short-range CDW order already present in zero field, not in the field-induced long-range CDWorder. Above Hc2, the relatively low values of $\chi_{spin}$ at T=2 K show that the pseudogap is a ground-state property, independent of the superconducting gap.


title: Improved Algorithm for Reconstructing Singular Connection in Multi-Block CFD Applications

abstract: An improved algorithm is proposed for the reconstruction of singular connectivity from the available pairwise connections during preprocessing phase. To evaluate the performance of the algorithm, an in-house computational fluid dynamics (CFD) code is used in which high-order finite-difference method for spatial discretization running on the Tianhe-1A supercomputer is employed. Test cases with a varied amount of mesh points are chosen, and the test results indicate that the improved singular connection reconstruction algorithm can achieve a speedup factor of 1000X or more when compared with the naive search method adopted in the former version of our code. Moreover, the parallel efficiency can benefit from the strategy of local communication based on the algorithm.


title: Multi-Task Learning by Deep Collaboration and Application in Facial Landmark Detection

abstract: Convolutional neural networks (CNNs) have become the most successful approach in many vision-related domains. However, they are limited to domains where data is abundant. Recent works have looked at multi-task learning (MTL) to mitigate data scarcity by leveraging domain-specific information from related tasks. In this paper, we present a novel soft-parameter sharing mechanism for CNNs in a MTL setting, which we refer to as Deep Collaboration. We propose taking into account the notion that task relevance depends on depth by using lateral transformation blocs with skip connections. This allows extracting task-specific features at various depth without sacrificing features relevant to all tasks. We show that CNNs connected with our Deep Collaboration obtain better accuracy on facial landmark detection with related tasks. We finally verify that our approach effectively allows knowledge sharing by showing depth-specific influence of tasks that we know are related.


title: PupilNet v2.0: Convolutional Neural Networks for CPU based real time Robust Pupil Detection

abstract: Real-time, accurate, and robust pupil detection is an essential prerequisite for pervasive video-based eye-tracking. However, automated pupil detection in realworld scenarios has proven to be an intricate challenge due to fast illumination changes, pupil occlusion, non-centered and off-axis eye recording, as well as physiological eye characteristics. In this paper, we approach this challenge through: I) a convolutional neural network (CNN) running in real time on a single core, II) a novel computational intensive two stage CNN for accuracy improvement, and III) a fast propability distribution based refinement method as a practical alternative to II. We evaluate the proposed approaches against the state-of-the-art pupil detection algorithms, improving the detection rate up to ~9% percent points on average over all data sets (~7% on one CPU core 7ms). This evaluation was performed on over 135,000 images: 94,000 images from the literature, and 41,000 new hand-labeled and challenging images contributed by this work (v1.0).


title: Proving Soundness of Extensional Normal-Form Bisimilarities

abstract: Normal-form bisimilarity is a simple, easy-to-use behavioral equivalence that relates terms in $\lambda$-calculi by decomposing their normal forms into bisimilar subterms. Moreover, it typically allows for powerful up-to techniques, such as bisimulation up to context, which simplify bisimulation proofs even further. However, proving soundness of these relations becomes complicated in the presence of $\eta$-expansion and usually relies on ad hoc proof methods which depend on the language. In this paper we propose a more systematic proof method to show that an extensional normal-form bisimilarity along with its corresponding up to context technique are sound. We illustrate our technique with three calculi: the call-by-value $\lambda$-calculus, the call-by-value $\lambda$-calculus with the delimited-control operators shift and reset, and the call-by-value $\lambda$-calculus with the abortive control operators call/cc and abort. In the first two cases, there was previously no sound up to context technique validating the $\eta$-law, whereas no theory of normal-form bisimulations for a calculus with call/cc and abort has been presented before. Our results have been fully formalized in the Coq proof assistant.


title: Stability of the Yang-Mills heat equation for finite action

abstract: The existence and uniqueness of solutions to the Yang-Mills heat equation over domains in Euclidean three space was proven in a previous paper for initial data lying in the Sobolev space of order one-half, which is the critical Sobolev index for this equation. In the present paper the stability of these solutions will be established. The variational equation, which is only weakly parabolic, and has highly singular coefficients, will be shown to have unique strong solutions up to addition of a vertical solution. Initial data will be taken to be in Sobolev class one-half. The proof relies on an infinitesimal version of the ZDS procedure: one solves first an augmented, strictly parabolic version of the variational equation and then adds to the solution a function which is vertical along the original path. Energy inequalities and Neumann domination techniques will be used to establish apriori initial behavior for solutions.


title: A class of C*-algebraic locally compact quantum groupoids Part I. Motivation and definition

abstract: In this series of papers, we develop the theory of a class of locally compact quantum groupoids, which is motivated by the purely algebraic notion of weak multiplier Hopf algebras. In this Part I, we provide motivation and formulate the definition in the C*-algebra framework. Existence of a certain canonical idempotent element is required and it plays a fundamental role, including the establishment of the coassociativity of the comultiplication. This class contains locally compact quantum groups as a subclass.


title: Inequivalence between gravitational mass and energy due to quantum effects at microscopic and macroscopic levels

abstract: We review recent theoretical results, demonstrating breakdown of the equivalence between active and passive gravitational masses and energy due to quantum effects in General Relativity. In particular, we discuss the simplest composite quantum body - a hydrogen atom - and define its gravitational masses operators. Using Gedanken experiment, we show that the famous Einstein's equation, $E =mc^2$, is broken with small probability for passive gravitational mass of the atom. It is important that the expectation values of both active and passive gravitational masses satisfy the above mentioned equation for stationary quantum states. Nevertheless, we stress that, for quantum superpositions of stationary states in a hydrogen atom, where the expectation values of energy are constant, the expectation values of the masses oscillate in time and, thus, break the Einstein's equation. We briefly discuss experimental possibility to observe the above-mentioned time-dependent oscillations. In this review, we also improve several drawbacks of the original pioneering works.


title: Countering Adversarial Images using Input Transformations

abstract: This paper investigates strategies that defend against adversarial-example attacks on image-classification systems by transforming the inputs before feeding them to the system. Specifically, we study applying image transformations such as bit-depth reduction, JPEG compression, total variance minimization, and image quilting before feeding the image to a convolutional network classifier. Our experiments on ImageNet show that total variance minimization and image quilting are very effective defenses in practice, in particular, when the network is trained on transformed images. The strength of those defenses lies in their non-differentiable nature and their inherent randomness, which makes it difficult for an adversary to circumvent the defenses. Our best defense eliminates 60% of strong gray-box and 90% of strong black-box attacks by a variety of major attack methods


title: Applied Machine Learning to Predict Stress Hotspots I: Face Centered Cubic Materials

abstract: We investigate the formation of stress hotspots in polycrystalline materials under uniaxial tensile deformation by integrating full field crystal plasticity based deformation models and machine learning techniques to gain data driven insights about microstructural properties. Synthetic 3D microstructures are created representing single phase equiaxed microstructures for generic copper alloys. Uniaxial tensile deformation is simulated using a 3-D full-field, image-based Fast Fourier Transform (FFT) technique with rate-sensitive crystal plasticity, to get local micro- mechanical fields (stress and strain rates). Stress hotspots are defined as the grains having stress values above the 90th percentile of the stress distribution. Hotspot neighborhoods are then characterized using metrics that reflect local crystallography, geometry, and connectivity. This data is used to create input feature vectors to train a random forest learning algorithm, which predicts the grains that will become stress hotspots. We are able to achieve an area under the receiving operating characteristic curve (ROC-AUC) of 0.74 for face centered cubic materials modeled on generic copper alloys. The results show the power and the limitations of the machine learning approach applied to the polycrystalline grain networks.


title: Shaping charge excitations in chiral edge states with a time-dependent gate voltage

abstract: We study a coherent conductor supporting a single edge channel in which alternating current pulses are created by local time-dependent gating and sent on a beam-splitter realized by a quantum point contact. The current response to the gate voltage in this setup is intrinsically linear. Based on a fully self-consistent treatment employing a Floquet scattering theory, we analyze the effect of different voltage shapes and frequencies, as well as the role of the gate geometry on the injected signal. In particular, we highlight the impact of frequency-dependent screening on the process of shaping the current signal. The feasibility of creating true single-particle excitations with this method is confirmed by investigating the suppression of excess noise, which is otherwise created by additional electron-hole pair excitations in the current signal.


title: Statistical Modeling of FSO Fronthaul Channel for Drone-based Networks

abstract: We consider a drone-based communication network, where several drones hover above an area and serve as mobile remote radio heads for a large number of mobile users. We assume that the drones employ free space optical (FSO) links for fronthauling of the users' data to a central unit. The main focus of this paper is to quantify the geometric loss of the FSO channel arising from random fluctuation of the position and orientation of the drones. In particular, we derive upper and lower bounds, corresponding approximate expressions, and a closed-form statistical model for the geometric loss. Simulation results validate our derivations and quantify the FSO channel quality as a function of the drone's instability, i.e., the variation of its position and orientation.


title: Dynamical SimRank Search on Time-Varying Networks

abstract: In this article, we study the efficient dynamical computation of all-pairs SimRanks on time-varying graphs. Li {\em et al}.'s approach requires $O(r^4 n^2)$ time and $O(r^2 n^2)$ memory in a graph with $n$ nodes, where $r$ is the target rank of the low-rank SVD. (1) We first consider edge update that does not accompany new node insertions. We show that the SimRank update $\Delta S$ in response to every link update is expressible as a rank-one Sylvester matrix equation. This provides an incremental method requiring $O(Kn^2)$ time and $O(n^2)$ memory in the worst case to update all pairs of similarities for $K$ iterations. (2) To speed up the computation further, we propose a lossless pruning strategy that captures the "affected areas" of $\Delta S$ to eliminate unnecessary retrieval. This reduces the time of the incremental SimRank to $O(K(m+|AFF|))$, where $m$ is the number of edges in the old graph, and $|AFF| (\le n^2)$ is the size of "affected areas" in $\Delta S$, and in practice, $|AFF| \ll n^2$. (3) We also consider edge updates that accompany node insertions, and categorize them into three cases, according to which end of the inserted edge is a new node. For each case, we devise an efficient incremental algorithm that can support new node insertions. (4) We next design an efficient batch incremental method that handles "similar sink edges" simultaneously and eliminates redundant edge updates. (5) To achieve linear memory, we devise a memory-efficient strategy that dynamically updates all pairs of SimRanks column by column in just $O(Kn+m)$ memory, without the need to store all $(n^2)$ pairs of old SimRank scores. Experimental studies on various datasets demonstrate that our solution substantially outperforms the existing incremental SimRank methods, and is faster and more memory-efficient than its competitors on million-scale graphs.


title: Quasi-Euclidean tilings over 2-dimensional Artin groups and their applications

abstract: We describe the structure of quasiflats in two-dimensio\-nal Artin groups. We rely on the notion of metric systolicity developed in our previous work. Using this weak form of non-positive curvature and analyzing in details the combinatorics of tilings of the plane we describe precisely the building blocks for quasiflats in all two-dimensional Artin groups -- atomic sectors. This allows us to provide useful quasi-isometry invariants for such groups -- completions of atomic sectors, stable lines, and the intersection pattern of certain abelian subgroups. These are described combinatorially, in terms of the structure of the graph defining an Artin group. As an important tool, we introduce an analogue of the curve complex in the context of two-dimensional Artin groups -- the intersection graph. We show quasi-isometric invariance of the intersection graph under natural assumptions. As immediate consequences we present a number of results concerning quasi-isometric rigidity for the subclass of CLTTF Artin groups. We give a necessary and sufficient condition for such groups to be strongly rigid (self quasi-isometries are close to automorphisms), we describe quasi-isometry groups, we indicate when quasi-isometries imply isomorphisms for such groups. In particular, there exist many strongly rigid large-type Artin groups. In contrast, none of the right-angled Artin groups are strongly rigid by a previous work of Bestvina, Kleiner and Sageev.


title: Backpropagation through the Void: Optimizing control variates for black-box gradient estimation

abstract: Gradient-based optimization is the foundation of deep learning and reinforcement learning. Even when the mechanism being optimized is unknown or not differentiable, optimization using high-variance or biased gradient estimates is still often the best strategy. We introduce a general framework for learning low-variance, unbiased gradient estimators for black-box functions of random variables. Our method uses gradients of a neural network trained jointly with model parameters or policies, and is applicable in both discrete and continuous settings. We demonstrate this framework for training discrete latent-variable models. We also give an unbiased, action-conditional extension of the advantage actor-critic reinforcement learning algorithm.


title: User Environment Detection with Acoustic Sensors Embedded on Mobile Devices for the Recognition of Activities of Daily Living

abstract: The detection of the environment where user is located, is of extreme use for the identification of Activities of Daily Living (ADL). ADL can be identified by use of the sensors available in many off-the-shelf mobile devices, including magnetic and motion, and the environment can be also identified using acoustic sensors. The study presented in this paper is divided in two parts: firstly, we discuss the recognition of the environment using acoustic sensors (i.e., microphone), and secondly, we fuse this information with motion and magnetic sensors (i.e., motion and magnetic sensors) for the recognition of standing activities of daily living. The recognition of the environments and the ADL are performed using pattern recognition techniques, in order to develop a system that includes data acquisition, data processing, data fusion, and artificial intelligence methods. The artificial intelligence methods explored in this study are composed by different types of Artificial Neural Networks (ANN), comparing the different types of ANN and selecting the best methods to implement in the different stages of the system developed. Conclusions point to the use of Deep Neural Networks (DNN) with normalized data for the identification of ADL with 85.89% of accuracy, the use of Feedforward neural networks with non-normalized data for the identification of the environments with 86.50% of accuracy, and the use of DNN with normalized data for the identification of standing activities with 100% of accuracy.


title: The Belyi degree is computable

abstract: We exhibit an algorithm that, given input a curve $X$ over a number field, computes as output the minimal degree of a Belyi map $X \to \mathbb{P}^1$.


title: Accelerated Sparse Subspace Clustering

abstract: State-of-the-art algorithms for sparse subspace clustering perform spectral clustering on a similarity matrix typically obtained by representing each data point as a sparse combination of other points using either basis pursuit (BP) or orthogonal matching pursuit (OMP). BP-based methods are often prohibitive in practice while the performance of OMP-based schemes are unsatisfactory, especially in settings where data points are highly similar. In this paper, we propose a novel algorithm that exploits an accelerated variant of orthogonal least-squares to efficiently find the underlying subspaces. We show that under certain conditions the proposed algorithm returns a subspace-preserving solution. Simulation results illustrate that the proposed method compares favorably with BP-based method in terms of running time while being significantly more accurate than OMP-based schemes.


title: Precise asymptotics of some meeting times arising from the voter model on large random regular graphs

abstract: We consider two independent stationary random walks on large random regular graphs of degree $k\geq 3$ with $N$ vertices. On these graphs, the exponential approximations of the meeting times are known to follow from existing methods and form a basis for the voter model's diffusion approximations. The main result of this note improves the exponential approximations to an explicit form such that the first moments are asymptotically equivalent to $N(k-1)/[2(k-2)]$.


title: Two procedures to flag radio frequency interference in the UV plane

abstract: We present two algorithms to identify and flag radio frequency interference (RFI) in radio interferometric imaging data. The first algorithm utilizes the redundancy of visibilities inside a UV cell in the visibility plane to identify corrupted data, while varying the detection threshold in accordance with the observed reduction in noise with radial UV distance. In the second algorithm, we propose a scheme to detect faint RFI in the visibility time-channel plane of baselines. The efficacy of identifying RFI in the residual visibilities is reduced by the presence of ripples due to inaccurate subtraction of the strongest sources. This can be due to several reasons including primary beam asymmetries and other direction dependent calibration errors. We eliminated these ripples by clipping the corresponding peaks in the associated Fourier plane. RFI was detected in the ripple-free time-channel plane but was flagged in the original visibilities. Application of these two algorithms to 5 different 150 MHz datasets from the GMRT resulted in a reduction in image noise of 20-50% throughout the field along with a reduction in systematics and a corresponding increase in the number of detected sources. However, on comparing the mean flux densities before and after flagging RFI we find a differential change with the fainter sources ($25\sigma <$ S $< 100$ mJy) showing a change of -6% to +1% relative to the stronger sources (S $>$ 100 mJy). We are unable to explain this effect but it could be related to the CLEAN bias known for interferometers.


title: Automata-Guided Hierarchical Reinforcement Learning for Skill Composition

abstract: Skills learned through (deep) reinforcement learning often generalizes poorly across domains and re-training is necessary when presented with a new task. We present a framework that combines techniques in \textit{formal methods} with \textit{reinforcement learning} (RL). The methods we provide allows for convenient specification of tasks with logical expressions, learns hierarchical policies (meta-controller and low-level controllers) with well-defined intrinsic rewards, and construct new skills from existing ones with little to no additional exploration. We evaluate the proposed methods in a simple grid world simulation as well as a more complicated kitchen environment in AI2Thor


title: Weakly coupled conformal gauge theories on the lattice

abstract: Results are reported for the beta-function of weakly coupled conformal gauge theories on the lattice, SU(3) with Nf=14 fundamental and Nf=3 sextet fermions. The models are chosen to be close to the upper end of the conformal window where perturbation theory is reliable hence a fixed point is expected. The study serves as a test of how well lattice methods perform in the weakly coupled conformal cases. We also comment on the 5-loop beta-function of two models close to the lower end of the conformal window, SU(3) with Nf=12 fundamental and Nf=2 sextet fermions.


title: Diffusive Molecular Communications with Reactive Signaling

abstract: This paper focuses on molecular communication (MC) systems where the signaling molecules may participate in a reversible bimolecular reaction in the channel. The motivation for studying these MC systems is that they can realize the concept of constructive and destructive signal superposition, which leads to favorable properties such as inter-symbol interference (ISI) reduction and avoiding environmental contamination due to continuous release of molecules into the channel. This work first derives the maximum likelihood (ML) detector for a binary MC system with reactive signaling molecules under the assumption that the detector has perfect knowledge of the ISI. The performance of this genie-aided ML detector yields an upper bound on the performance of any practical detector. In addition, two suboptimal detectors of different complexity are proposed. The proposed ML detector as well as one of the suboptimal detectors require the channel response (CR) of the considered MC system. Moreover, the CR is needed for the performance evaluation of all proposed detectors. However, analyzing MC with reactive signaling is challenging since the underlying partial differential equations that describe the reaction-diffusion mechanism are coupled and non-linear. Therefore, an algorithm is developed in this paper for efficient computation of the CR to any arbitrary transmit symbol sequence. The accuracy of this algorithm is validated via particle-based simulation. Simulation results using the developed CR algorithm show that the performance of the proposed suboptimal detectors can approach that of the genie- aided ML detector. Moreover, these results show that MC systems with reactive signaling have superior performance relative to those with non-reactive signaling due to the reduction of ISI enabled by the chemical reactions.


title: Explanation of the intriguing quantum phenomenon of the off-resonant cavity mode emission

abstract: We theoretically investigate the unexpected occurrence of an extra emission peak that has been experimentally observed in off-resonant studies of cavity QED systems. Our results within the Markovian master equation approach successfully explain why the central peak arises, and how it reveals that the system is suffering a dynamical phase transition induced by the phonon-mediated coupling. Our findings are in perfect agreement with previous reported experimental results and for the first time the fundamental physics behind this quantum phenomenon is understood.


title: On the internal structure of the current sheet in the pulsar wind

abstract: We investigate the internal structure of the current sheet in the pulsar wind within force-free and two-fluid MHD approximations. Within the force-free approximation we obtain general asymptotic solution of the Grad-Shafranov equation for quasi-spherical pulsar wind up to the second order in small parameter $\varepsilon = (\Omega r/c)^{-1}$. The solution allows an arbitrary latitudinal structure of the radial magnetic field, including that obtained in the numerical simulations of oblique rotators. It is also shown that the shape of the current sheet does not depend on the latitudinal structure. For the internal region of the current sheet outside the fast magnetosonic surface where the force-free approximation is not valid we use two-fluid MHD approximation. Carrying out calculations in the comoving reference frame we succeed in determining intrinsic electric and magnetic fields of a sheet. It allows us to analyze time-dependent effects which were not investigated up to now. In particular, we estimate the efficiency of the particle acceleration inside the sheet. Finally, after investigating the motion of individual particles in the time-dependent current sheet, we find self-consistently the width of the sheet and its time evolution.


title: Recurrence due to periodic multi-soliton fission in the defocusing nonlinear Schrodinger equation

abstract: We address the degree of universality of the Fermi-Pasta-Ulam recurrence induced by multisoliton fission from a harmonic excitation by analysing the case of the semiclassical defocusing nonlinear Schrodinger equation, which models nonlinear wave propagation in a variety of physical settings. Using a suitable Wentzel-Kramers-Brillouin approach to the solution of the associated scattering problem we accurately predict, in full analytical way, the number and the features (amplitude and velocity) of soliton-like excitations emerging post-breaking, as a function of the dispersion smallness parameter. This also permits to predict and analyse the near-recurrences, thereby inferring the universal character of the mechanism originally discovered for the Korteweg-deVries equation. We show, however, that important differences exist between the two models, arising from the different scaling rules obeyed by the soliton velocities.


title: Nonlocal Quantum Effective Actions in Weyl-Flat Spacetimes

abstract: Virtual massless particles in quantum loops lead to nonlocal effects which can have interesting consequences, for example, for primordial magnetogenesis in cosmology or for computing finite $N$ corrections in holography. We describe how the quantum effective actions summarizing these effects can be computed efficiently for Weyl-flat metrics by integrating the Weyl anomaly or, equivalently, the local renormalization group equation. This method relies only on the local Schwinger-DeWitt expansion of the heat kernel and allows for a re-summation of leading large logarithms in situations where the Weyl factor changes by several e-foldings. As an illustration, we obtain the quantum effective action for the Yang-Mills field coupled to massless matter, and the self-interacting massless scalar field. Our action reduces to the nonlocal action obtained using the Barvinsky-Vilkovisky covariant perturbation theory in the regime $R^{2} \ll \nabla^{2} R $ for a typical curvature scale $R$, but has a greater range of validity effectively re-summing the covariant perturbation theory to all orders in curvatures. In particular, it is applicable also in the opposite regime $R^{2} \gg \nabla^{2} R$, which is often of interest in cosmology.


title: Bayesian model comparison with the Hyv\"arinen score: computation and consistency

abstract: The Bayes factor is a widely used criterion in model comparison and its logarithm is a difference of out-of-sample predictive scores under the logarithmic scoring rule. However, when some of the candidate models involve vague priors on their parameters, the log-Bayes factor features an arbitrary additive constant that hinders its interpretation. As an alternative, we consider model comparison using the Hyv\"arinen score. We propose a method to consistently estimate this score for parametric models, using sequential Monte Carlo methods. We show that this score can be estimated for models with tractable likelihoods as well as nonlinear non-Gaussian state-space models with intractable likelihoods. We prove the asymptotic consistency of this new model selection criterion under strong regularity assumptions in the case of non-nested models, and we provide qualitative insights for the nested case. We also use existing characterizations of proper scoring rules on discrete spaces to extend the Hyv\"arinen score to discrete observations. Our numerical illustrations include L\'evy-driven stochastic volatility models and diffusion models for population dynamics.


title: Pomegranate: fast and flexible probabilistic modeling in python

abstract: We present pomegranate, an open source machine learning package for probabilistic modeling in Python. Probabilistic modeling encompasses a wide range of methods that explicitly describe uncertainty using probability distributions. Three widely used probabilistic models implemented in pomegranate are general mixture models, hidden Markov models, and Bayesian networks. A primary focus of pomegranate is to abstract away the complexities of training models from their definition. This allows users to focus on specifying the correct model for their application instead of being limited by their understanding of the underlying algorithms. An aspect of this focus involves the collection of additive sufficient statistics from data sets as a strategy for training models. This approach trivially enables many useful learning strategies, such as out-of-core learning, minibatch learning, and semi-supervised learning, without requiring the user to consider how to partition data or modify the algorithms to handle these tasks themselves. pomegranate is written in Cython to speed up calculations and releases the global interpreter lock to allow for built-in multithreaded parallelism, making it competitive with---or outperform---other implementations of similar algorithms. This paper presents an overview of the design choices in pomegranate, and how they have enabled complex features to be supported by simple code.


title: Visualizing and Understanding Atari Agents

abstract: While deep reinforcement learning (deep RL) agents are effective at maximizing rewards, it is often unclear what strategies they use to do so. In this paper, we take a step toward explaining deep RL agents through a case study using Atari 2600 environments. In particular, we focus on using saliency maps to understand how an agent learns and executes a policy. We introduce a method for generating useful saliency maps and use it to show 1) what strong agents attend to, 2) whether agents are making decisions for the right or wrong reasons, and 3) how agents evolve during learning. We also test our method on non-expert human subjects and find that it improves their ability to reason about these agents. Overall, our results show that saliency information can provide significant insight into an RL agent's decisions and learning behavior.


title: Segmentation-by-Detection: A Cascade Network for Volumetric Medical Image Segmentation

abstract: We propose an attention mechanism for 3D medical image segmentation. The method, named segmentation-by-detection, is a cascade of a detection module followed by a segmentation module. The detection module enables a region of interest to come to attention and produces a set of object region candidates which are further used as an attention model. Rather than dealing with the entire volume, the segmentation module distills the information from the potential region. This scheme is an efficient solution for volumetric data as it reduces the influence of the surrounding noise which is especially important for medical data with low signal-to-noise ratio. Experimental results on 3D ultrasound data of the femoral head shows superiority of the proposed method when compared with a standard fully convolutional network like the U-Net.


title: Solving the Vlasov equation in two spatial dimensions with the Schr\"odinger method

abstract: We demonstrate that the Vlasov equation describing collisionless self-gravitating matter may be solved with the so-called Schr\"odinger method (ScM). With the ScM, one solves the Schr\"odinger-Poisson system of equations for a complex wave function in d dimensions, rather than the Vlasov equation for a 2d-dimensional phase space density. The ScM also allows calculating the d-dimensional cumulants directly through quasi-local manipulations of the wave function, avoiding the complexity of 2d-dimensional phase space. We perform for the first time a quantitive comparison of the ScM and a conventional Vlasov solver in d=2 dimensions. Our numerical tests were carried out using two types of cold cosmological initial conditions: the classic collapse of a sine wave and those of a gaussian random field as commonly used in cosmological cold dark matter N-body simulations. We compare the first three cumulants, that is, the density, velocity and velocity dispersion, to those obtained by solving the Vlasov equation using the publicly available code ColDICE. We find excellent qualitative and quantitative agreement between these codes, demonstrating the feasibility and advantages of the ScM as an alternative to N-body simulations. We discuss, the emergence of effective vorticity in the ScM through the winding number around the points where the wave function vanishes. As an application we evaluate the background pressure induced by the non-linearity of large scale structure formation, thereby estimating the magnitude of cosmological backreaction. We find that it is negligibly small and has time dependence and magnitude compatible with expectations from the effective field theory of large scale structure.


title: Training GANs with Optimism

abstract: We address the issue of limit cycling behavior in training Generative Adversarial Networks and propose the use of Optimistic Mirror Decent (OMD) for training Wasserstein GANs. Recent theoretical results have shown that optimistic mirror decent (OMD) can enjoy faster regret rates in the context of zero-sum games. WGANs is exactly a context of solving a zero-sum game with simultaneous no-regret dynamics. Moreover, we show that optimistic mirror decent addresses the limit cycling problem in training WGANs. We formally show that in the case of bi-linear zero-sum games the last iterate of OMD dynamics converges to an equilibrium, in contrast to GD dynamics which are bound to cycle. We also portray the huge qualitative difference between GD and OMD dynamics with toy examples, even when GD is modified with many adaptations proposed in the recent literature, such as gradient penalty or momentum. We apply OMD WGAN training to a bioinformatics problem of generating DNA sequences. We observe that models trained with OMD achieve consistently smaller KL divergence with respect to the true underlying distribution, than models trained with GD variants. Finally, we introduce a new algorithm, Optimistic Adam, which is an optimistic variant of Adam. We apply it to WGAN training on CIFAR10 and observe improved performance in terms of inception score as compared to Adam.


title: Sampling and Reconstruction of Graph Signals via Weak Submodularity and Semidefinite Relaxation

abstract: We study the problem of sampling a bandlimited graph signal in the presence of noise, where the objective is to select a node subset of prescribed cardinality that minimizes the signal reconstruction mean squared error (MSE). To that end, we formulate the task at hand as the minimization of MSE subject to binary constraints, and approximate the resulting NP-hard problem via semidefinite programming (SDP) relaxation. Moreover, we provide an alternative formulation based on maximizing a monotone weak submodular function and propose a randomized-greedy algorithm to find a sub-optimal subset. We then derive a worst-case performance guarantee on the MSE returned by the randomized greedy algorithm for general non-stationary graph signals. The efficacy of the proposed methods is illustrated through numerical simulations on synthetic and real-world graphs. Notably, the randomized greedy algorithm yields an order-of-magnitude speedup over state-of-the-art greedy sampling schemes, while incurring only a marginal MSE performance loss.


title: Global Existence of Solutions for a Fractional Caputo Nonlocal Thermistor Problem

abstract: We begin by proving a local existence result for a fractional Caputo nonlocal thermistor problem. Then, additional existence and continuation theorems are obtained, ensuring global existence of solutions.


title: Andreev reflection spectroscopy on Bi$_{2}$X$_{3}$ (X = Se, Te) topological insulators: Implications for the c-axis superconducting proximity effect

abstract: Using Andreev reflection (AR) as an experimental gauge of the superconducting proximity effect (PE), we assess the topological purity of the superconductivity that is induced by the c-axis PE between an s-wave superconductor and the topological insulators Bi$_{2}$X$_{3}$ (X=Se,Te). Point-contact AR spectroscopy is performed with Nb tips on Bi$_{2}$X$_{3}$ single crystals at 4.2 K. Scanning tunneling spectroscopy is also used, to locate the Fermi level $E_F$ relative to the Dirac point in the crystals. The AR data is analyzed with Blonder-Tinkham-Klapwijk theory, taking into account tip-induced spin-orbit coupling, Fermi-surface mismatch, and the co-presence of bulk band and topological surface states at $E_F$. Our results indicate that the superconductivity that can be proximity-induced into Bi$_{2}$X$_{3}$ is predominantly non-topological.


title: The Unreasonable Rigidity of Ulam Sets

abstract: We give a number of results about families of Ulam sets. Generalizing behavior of Ulam sets U(1,n), we prove using an novel model theoretic approach that there is a rigidity phenomenon for Ulam sets U(a,b) as b increases. Based on this, we suggest a natural conjecture, and investigate its potential applications, including a method of proving certain families of Ulam sequences are regular, for which we also provide partial, unconditional, results. Along this same vein, we give an upper bound bound on the density of Ulam sequences U(1,n). Finally, we give classification results for higher dimensional Ulam sets.


title: A multitask deep learning model for real-time deployment in embedded systems

abstract: We propose an approach to Multitask Learning (MTL) to make deep learning models faster and lighter for applications in which multiple tasks need to be solved simultaneously, which is particularly useful in embedded, real-time systems. We develop a multitask model for both Object Detection and Semantic Segmentation and analyze the challenges that appear during its training. Our multitask network is 1.6x faster, lighter and uses less memory than deploying the single-task models in parallel. We conclude that MTL has the potential to give superior performance in exchange of a more complex training process that introduces challenges not present in single-task models.


title: A tableau formula of double Grothendieck polynomials for $321$-avoiding permutations

abstract: In this article, we prove a tableau formula for the double Grothendieck polynomials associated to $321$-avoiding permutations. The proof is based on the compatibility of the formula with the $K$-theoretic divided difference operators.


title: Computing the Lusztig--Vogan Bijection

abstract: Let $G$ be a connected complex reductive algebraic group with Lie algebra $\mathfrak{g}$. The Lusztig--Vogan bijection relates two bases for the bounded derived category of $G$-equivariant coherent sheaves on the nilpotent cone $\mathcal{N}$ of $\mathfrak{g}$. One basis is indexed by $\Lambda^+$, the set of dominant weights of $G$, and the other by $\Omega$, the set of pairs $(\mathcal{O}, \mathcal{E})$ consisting of a nilpotent orbit $\mathcal{O} \subset \mathcal{N}$ and an irreducible $G$-equivariant vector bundle $\mathcal{E} \rightarrow \mathcal{O}$. The existence of the Lusztig--Vogan bijection $\gamma \colon \Omega \rightarrow \Lambda^+$ was proven by Bezrukavnikov, and an algorithm computing $\gamma$ in type $A$ was given by Achar. Herein we present a combinatorial description of $\gamma$ in type $A$ that subsumes and dramatically simplifies Achar's algorithm.


title: Statistical Inference of Kumaraswamy distribution under imprecise information

abstract: Traditional statistical approaches for estimating the parameters of the Kumaraswamy distribution have dealt with precise information. However, in real world situations, some information about an underlying experimental process might be imprecise and might be represented in the form of fuzzy information. In this paper, we consider the problem of estimating the parameters of a univariate Kumaraswamy distribution with two parameters when the available observations are described by means of fuzzy information. We derive the maximum likelihood estimate of the parameters by using Newton Raphson as well as EM algorithm method. The estimation procedures are discussed in details and compared via Markov Chain Monte Carlo simulations in terms of their average biases and mean squared errors.


title: Erratum: Link prediction in drug-target interactions network using similarity indices

abstract: Background: In silico drug-target interaction (DTI) prediction plays an integral role in drug repositioning: the discovery of new uses for existing drugs. One popular method of drug repositioning is network-based DTI prediction, which uses complex network theory to predict DTIs from a drug-target network. Currently, most network-based DTI prediction is based on machine learning methods such as Restricted Boltzmann Machines (RBM) or Support Vector Machines (SVM). These methods require additional information about the characteristics of drugs, targets and DTIs, such as chemical structure, genome sequence, binding types, causes of interactions, etc., and do not perform satisfactorily when such information is unavailable. We propose a new, alternative method for DTI prediction that makes use of only network topology information attempting to solve this problem. Results: We compare our method for DTI prediction against the well-known RBM approach. We show that when applied to the MATADOR database, our approach based on node neighborhoods yield higher precision for high-ranking predictions than RBM when no information regarding DTI types is available. Conclusion: This demonstrates that approaches purely based on network topology provide a more suitable approach to DTI prediction in the many real-life situations where little or no prior knowledge is available about the characteristics of drugs, targets, or their interactions.


title: Spectral correlations in Anderson insulating wires

abstract: We calculate the spectral level-level correlation function of Anderson insulating wires for all three Wigner-Dyson classes. A measurement of its Fourier transform, the spectral form factor, is within reach of state-of-the-art cold atom quantum quench experiments, and we find good agreement with recent numerical simulations of the latter. Our derivation builds on a representation of the level-level correlation function in terms of a local generating function which may prove useful in other contexts.


title: Searches for supersymmetry

abstract: New and recent results on Supersymmetry searches are shown for the ATLAS and the CMS experiments. Analyses with about 36 fb$^{-1}$ are considered for searches concerning light squarks and gluinos, direct pair production of 3$^{rd}$ generation squarks, electroweak production of charginos, neutralinos, sleptons, R-parity violating scenarios and long-lived particles.


title: Weakly 2-randoms and 1-generics in Scott sets

abstract: Let $S$ be a Scott set, or even an $\omega$-model of $\mathsf{WWKL}$. Then for each $A\in S$, either there is $X \in S$ that is weakly 2-random relative to $A$, or there is $X\in S$ that is 1-generic relative to $A$. It follows that if $A_1,\dots, A_n \in S$ are non-computable, there is $X \in S$ such that each $A_i$ is Turing incomparable with $X$, answering a question of Ku\v{c}era and Slaman. More generally, any $\forall\exists$ sentence in the language of partial orders that holds in $\mathcal D$ also holds in $\mathcal D_S$, where $\mathcal D_S$ is the partial order of Turing degrees of elements of $S$.


title: An effective analysis of the Denjoy rank

abstract: We analyze the descriptive complexity of several $\Pi^1_1$ ranks from classical analysis which are associated to Denjoy integration. We show that $VBG, VBG_\ast, ACG$ and $ACG_\ast$ are $\Pi^1_1$-complete, answering a question of Walsh in case of $ACG_\ast$. Furthermore, we identify the precise descriptive complexity of the set of functions obtainable with at most $\alpha$ steps of the transfinite process of Denjoy totalization: if $|\cdot|$ is the $\Pi^1_1$-rank naturally associated to $VBG, VBG_\ast$ or $ACG_\ast$, and if $\alpha<\omega_1^{ck}$, then $\{F \in C(I): |F| \leq \alpha\}$ is $\Sigma^0_{2\alpha}$-complete. These finer results are an application of the author's previous work on the limsup rank on well-founded trees. Finally, $\{(f,F) \in M(I)\times C(I) : F\in ACG_\ast \text{ and } F'=f \text{ a.e.}\}$ and $\{f \in M(I) : f \text{ is Denjoy integrable}\}$ are $\Pi^1_1$-complete, answering more questions of Walsh.


title: Neural Wikipedian: Generating Textual Summaries from Knowledge Base Triples

abstract: Most people do not interact with Semantic Web data directly. Unless they have the expertise to understand the underlying technology, they need textual or visual interfaces to help them make sense of it. We explore the problem of generating natural language summaries for Semantic Web data. This is non-trivial, especially in an open-domain context. To address this problem, we explore the use of neural networks. Our system encodes the information from a set of triples into a vector of fixed dimensionality and generates a textual summary by conditioning the output on the encoded vector. We train and evaluate our models on two corpora of loosely aligned Wikipedia snippets and DBpedia and Wikidata triples with promising results.


title: Time-resolved ultrafast x-ray scattering from an incoherent electronic mixture

abstract: Time-resolved ultrafast x-ray scattering from photo-excited matter is an emerging method to image ultrafast dynamics in matter with atomic-scale spatial and temporal resolutions. For a correct and rigorous understanding of current and upcoming imaging experiments, we present the theory of time-resolved x-ray scattering from an incoherent electronic mixture using quantum electrodynamical theory of light-matter interaction. We show that the total scattering signal is an incoherent sum of the individual scattering signals arising from different electronic states and therefore heterodyning of the individual signals is not possible for an ensemble of gas-phase photo-excited molecules. We scrutinize the information encoded in the total signal for the experimentally important situation when pulse duration and coherence time of the x-ray pulse are short in comparison to the timescale of the vibrational motion and long in comparison to the timescale of the electronic motion, respectively. Finally, we show that in the case of an electronically excited crystal the total scattering signal imprints the interference of the individual scattering amplitudes associated with different electronic states and heterodyning is possible.

categories: physics.chem-ph cond-mat.mes-hall cond-mat.mtrl-sci physics.atom-ph


title: Bayesian Variable Selection for Multivariate Zero-Inflated Models: Application to Microbiome Count Data

abstract: Microorganisms play critical roles in human health and disease. It is well known that microbes live in diverse communities in which they interact synergistically or antagonistically. Thus for estimating microbial associations with clinical covariates, multivariate statistical models are preferred. Multivariate models allow one to estimate and exploit complex interdependencies among multiple taxa, yielding more powerful tests of exposure or treatment effects than application of taxon-specific univariate analyses. In addition, the analysis of microbial count data requires special attention because data commonly exhibit zero inflation. To meet these needs, we developed a Bayesian variable selection model for multivariate count data with excess zeros that incorporates information on the covariance structure of the outcomes (counts for multiple taxa), while estimating associations with the mean levels of these outcomes. Although there has been a great deal of effort in zero-inflated models for longitudinal data, little attention has been given to high-dimensional multivariate zero-inflated data modeled via a general correlation structure. Through simulation, we compared performance of the proposed method to that of existing univariate approaches, for both the binary and count parts of the model. When outcomes were correlated the proposed variable selection method maintained type I error while boosting the ability to identify true associations in the binary component of the model. For the count part of the model, in some scenarios the the univariate method had higher power than the multivariate approach. This higher power was at a cost of a highly inflated false discovery rate not observed with the proposed multivariate method. We applied the approach to oral microbiome data from the Pediatric HIV/AIDS Cohort Oral Health Study and identified five species (of 44) associated with HIV infection.


title: On the Ristic-Balakrishnan distribution: bivariate extension and characterizations

abstract: Over the last few decades, a significant development has been made towards the augmentation of some well-known lifetime distributions by various strategies. These newly developed models have enjoyed a considerable amount of success in modeling various real life phenomena. Motivated by this, Ristic & Balakrishnan (2012) developed a special class of univariate distributions (see Ristic- Balakrishnan (2012)). Henceforth we call this family of distribution as RB-G family of distributions. The RB-G family has the same parameters of the G distribution plus an additional positive shape parameter a. Several RB-G distribution can be obtained from a specified G distribution. For a = 1, the baseline G distribution is a basic exemplar of the RB-G family with a continuous crossover towards cases with various shapes. In this article we focus our attention on the characterization of this family and discuss some structural properties of the bivariate RB-G family of distributions which are not discussed in detail in Ristic and Balakrishnan (2012).


title: A New Signature Scheme Based on Punctured Reed--Muller Code With Random Insertion

abstract: In this paper, we propose a new signature scheme based on a punctured Reed--Muller (RM) code with random insertion, which improves the Goppa code-based signature scheme developed by Courtois, Finiasz, and Sendrier (CFS). The CFS signature scheme has certain drawbacks in terms of scaling of the parameters and a lack of existential unforgeability under adaptive chosen message attacks (EUF-CMA) security proof. Further, the proposed modified RM code-based signature scheme can use complete decoding, which can be implemented using a recursive decoding method, and thus syndromes for errors larger than the error correctability can be decoded for signing, which improves the probability of successful signing and reduces the signing time. Using the puncturing and insertion methods, the proposed RM code-based signature scheme can avoid some known attacks for RM code-based cryptosystems. The parameters of the proposed signature scheme such as error weight parameter $w$ and the maximum signing trial $N$, can be adjusted in terms of signing time and security level, and it is also proved that the proposed signature scheme achieves EUF-CMA security.


title: Delbruck scattering in small scattering angle region

abstract: The scattering of photons of x-ray energy off a Coulomb field in very forward scattering region may be thought as the refraction effect due to the Coulomb field. The cross section of the scattering can be computed from the photon bending angle obtained in geometrical optics. Its dependence on the fine structure constant is nonanalytic and softer than would the box diagrams of Delbruck scattering suggest. The refractive scattering is a nonlocal effect occurring over the distance scale of the impact parameter, and therefore may not be described by the single-box diagrams. We suggest that the scattering cross section is a consequence of the sum of the asymptotic series of the multi-box diagrams.


title: Hypoellipticity without loss of derivatives for Fedii's type operators

abstract: We prove that second order linear operators on $\mathbb{R}^{n+m}$ of the form $L(x,y,D_x,D_y) = L_1(x,D_x) + g(x) L_2(y,D_y)$, where $L_1$ and $L_2$ satisfy Morimoto's super-logarithmic estimates and $g$ is smooth, nonnegative, and vanishes only at the origin in $\mathbb{R}^n$ (but to any arbitrary order) are hypoelliptic without loss of derivarives. We also show examples in which our hypotheses are necessary for hypoellipticity.


title: Dynamic quantile linear models: a Bayesian approach

abstract: A new class of models, named dynamic quantile linear models, is presented. It combines dynamic linear models with distribution free quantile regression producing a robust statistical method. Bayesian inference for dynamic quantile linear models can be performed using an efficient Markov chain Monte Carlo algorithm. A fast sequential procedure suited for high-dimensional predictive modeling applications with massive data, in which the generating process is itself changing overtime, is also proposed. The proposed model is evaluated using synthetic and well-known time series data. The model is also applied to predict annual incidence of tuberculosis in Rio de Janeiro state for future years and compared with global strategy targets set by the World Health Organization.


title: Cluster Algebras, Invariant Theory, and Kronecker Coefficients II

abstract: We prove that the semi-invariant ring of the standard representation space of the $l$-flagged $m$-arrow Kronecker quiver is an upper cluster algebra for any $l,m\in \mathbb{N}$. The quiver and cluster are explicitly given. We prove that the quiver with its rigid potential is a polyhedral cluster model. As a consequence, to compute each Kronecker coefficient $g_{\mu,\nu}^\lambda$ with $\lambda$ at most $m$ parts, we only need to count lattice points in at most $m!$ fibre (rational) polytopes inside the ${\rm g}$-vector cone, which is explicitly given.


title: Improving Object Localization with Fitness NMS and Bounded IoU Loss

abstract: We demonstrate that many detection methods are designed to identify only a sufficently accurate bounding box, rather than the best available one. To address this issue we propose a simple and fast modification to the existing methods called Fitness NMS. This method is tested with the DeNet model and obtains a significantly improved MAP at greater localization accuracies without a loss in evaluation rate, and can be used in conjunction with Soft NMS for additional improvements. Next we derive a novel bounding box regression loss based on a set of IoU upper bounds that better matches the goal of IoU maximization while still providing good convergence properties. Following these novelties we investigate RoI clustering schemes for improving evaluation rates for the DeNet wide model variants and provide an analysis of localization performance at various input image dimensions. We obtain a MAP of 33.6%@79Hz and 41.8%@5Hz for MSCOCO and a Titan X (Maxwell). Source code available from: https://github.com/lachlants/denet


title: Deep Neural Networks as Gaussian Processes

abstract: It has long been known that a single-layer fully-connected neural network with an i.i.d. prior over its parameters is equivalent to a Gaussian process (GP), in the limit of infinite network width. This correspondence enables exact Bayesian inference for infinite width neural networks on regression tasks by means of evaluating the corresponding GP. Recently, kernel functions which mimic multi-layer random neural networks have been developed, but only outside of a Bayesian framework. As such, previous work has not identified that these kernels can be used as covariance functions for GPs and allow fully Bayesian prediction with a deep neural network. In this work, we derive the exact equivalence between infinitely wide deep networks and GPs. We further develop a computationally efficient pipeline to compute the covariance function for these GPs. We then use the resulting GPs to perform Bayesian inference for wide deep neural networks on MNIST and CIFAR-10. We observe that trained neural network accuracy approaches that of the corresponding GP with increasing layer width, and that the GP uncertainty is strongly correlated with trained network prediction error. We further find that test performance increases as finite-width trained networks are made wider and more similar to a GP, and thus that GP predictions typically outperform those of finite-width networks. Finally we connect the performance of these GPs to the recent theory of signal propagation in random neural networks.


title: The Mechanism for the Energy Buildup Driving Solar Eruptive Events

abstract: The underlying origin of solar eruptive events (SEEs), ranging from giant coronal mass ejections to small coronal-hole jets, is that the lowest-lying magnetic flux in the Sun's corona undergoes the continual buildup of stress and free energy. This magnetic stress has long been observed as the phenomenon of "filament channels:" strongly sheared magnetic field localized around photospheric polarity inversion lines. However, the mechanism for the stress buildup - the formation of filament channels - is still debated. We present magnetohydrodynamic simulations of a coronal volume that is driven by transient, cellular boundary flows designed to model the processes by which the photosphere drives the corona. The key feature of our simulations is that they accurately preserve magnetic helicity, the topological quantity that is conserved even in the presence of ubiquitous magnetic reconnection. Although small-scale random stress is injected everywhere at the photosphere, driving stochastic reconnection throughout the corona, the net result of the magnetic evolution is a coherent shearing of the lowest-lying field lines. This highly counter-intuitive result - magnetic stress builds up locally rather than spreading out to attain a minimum energy state - explains the formation of filament channels and is the fundamental mechanism underlying SEEs. Furthermore, this mechanism may be relevant to other astrophysical or laboratory plasmas.


title: The Cost of Uncertainty in Curing Epidemics

abstract: Motivated by the study of controlling (curing) epidemics, we consider the spread of an SI process on a known graph, where we have a limited budget to use to transition infected nodes back to the susceptible state (i.e., to cure nodes). Recent work has demonstrated that under perfect and instantaneous information (which nodes are/are not infected), the budget required for curing a graph precisely depends on a combinatorial property called the CutWidth. We show that this assumption is in fact necessary: even a minor degradation of perfect information, e.g., a diagnostic test that is 99% accurate, drastically alters the landscape. Infections that could previously be cured in sublinear time now may require exponential time, or orderwise larger budget to cure. The crux of the issue comes down to a tension not present in the full information case: if a node is suspected (but not certain) to be infected, do we risk wasting our budget to try to cure an uninfected node, or increase our certainty by longer observation, at the risk that the infection spreads further? Our results present fundamental, algorithm-independent bounds that tradeoff budget required vs. uncertainty.


title: Outdoor to Indoor Penetration Loss at 28 GHz for Fixed Wireless Access

abstract: This paper present the results from a 28 GHz channel sounding campaign performed to investigate the effects of outdoor to indoor penetration on the wireless propagation channel characteristics for an urban microcell in a fixed wireless access scenario. The measurements are performed with a real-time channel sounder, which can measure path loss up to 169 dB, and equipped with phased array antennas that allows electrical beam steering for directionally resolved measurements in dynamic environments. Thanks to the short measurement time and the excellent phase stability of the system, we obtain both directional and omnidirectional channel power delay profiles without any delay uncertainty. For outdoor and indoor receiver locations, we compare path loss, delay spreads and angular spreads obtained for two different types of buildings.


title: Dynamic Double Directional Propagation Channel Measurements at 28 GHz

abstract: This paper presents results from the (to our knowledge) first dynamic double-directionally resolved measurement campaign at mm-wave frequencies for an outdoor microcellular scenario. The measurements are performed with USC's real-time channel sounder equipped with phased array antennas that can steer beams electrically in microseconds, allowing directional measurements in dynamic environments. Exploiting the phase coherency of the setup, the multi-path components can be tracked over time to investigate the temporal dependencies of the channel characteristics. We present results for time-varying path-loss, delay spread, mean angles and angular spreads observed at the transmitter (TX) and receiver (RX) in the presence of moving vehicles and pedestrians. Additionally, we investigate excess losses observed due to blockage by vehicles and compare the cases when TX and RX are using fixed beams or when they are capable of adjusting beam directions dynamically.


title: 28 GHz Microcell Measurement Campaign for Residential Environment

abstract: This paper presents results from the (to our knowledge) first double-directionally resolved measurement campaign at mm-wave frequencies in a suburban microcell. The measurements are performed with a real-time channel sounder equipped with phased antenna arrays that allows electrical beam steering in microseconds, and which can measure path-loss of up to 169 dB. Exploiting the phase coherency of the measurements in the different beams, we obtain both directional and omnidirectional channel power delay profiles without any delay uncertainty. We present statistics of channel characteristics such as path-loss, shadowing and delay spread results for line-of-sight and non-line-of-sight cases, as well as sample results for power angular spectrum and extracted multi-path components.


title: On some further properties and application of Weibull-R family of distributions

abstract: In this paper, we provide some new results for the Weibull-R family of distributions (Alzaghal, Ghosh and Alzaatreh (2016)). We derive some new structural properties of the Weibull-R family of distributions. We provide various characterizations of the family via conditional moments, some functions of order statistics and via record values.


title: On $AP_3$ - covering sequences

abstract: Recently, motivated by Stanley sequences, Kiss, S\' andor and Yang introduced a new type sequence: a sequence $A$ of nonnegative integers is called an $AP_k$ - covering sequence if there exists an integer $n_0$ such that if $n > n_0$, then there exist $a_1\in A, \dots , a_{k-1}\in A$, $a_1<a_2<\cdots <a_{k-1}<n$ such that $a_1, \dots , a_{k-1}, n$ form a $k$-term arithmetic progression. They prove that there exists an $AP_3$ - covering sequence $A$ such that $\limsup\limits_{n\to\infty}{A(n)}/{\sqrt n}\le 34$. In this note, we prove that there exists an $AP_3$ - covering sequence $A$ such that $\limsup\limits_{n\to\infty}{A(n)}/{\sqrt n}=\sqrt{15}$.


title: 4-dimensional Riemannian manifolds with a harmonic 2-form of constant length

abstract: It was shown by Seaman that if a compact, oriented 4-dimensional riemannian manifold (M, g) of positive sectional curvature admits a harmonic 2-form of constant length, its intersection form is definite and such a harmonic form is unique up to constant multiples. In this paper, we show that such a manifold is diffeomorphic to $\mathbb{CP}_{2}$ with a slightly weaker curvature hypothesis and there is an infinite dimensional moduli space of such metrics near the Fubini-Study metric on $\mathbb{CP}_{2}$. We discuss some of conditions which can be added in order to get the Fubini-Study metric up to diffeomorphisms and rescaling.


title: On a problem of Nathanson

abstract: A set $A$ of nonnegative integers is an asymptotic basis of order $h$ if every sufficiently large integer can be represented as the sum of $h$ integers (not necessarily distinct) of $A$. An asymptotic basis $A$ of order $h$ is minimal if no proper subset of $A$ is an asymptotic basis of order $h$. In this paper, we resolve a problem of Nathanson on minimal asymptotic bases of order $h$.


title: The number of spanning trees in circulant graphs, its arithmetic properties and asymptotic

abstract: In this paper, we develop a new method to produce explicit formulas for the number $\tau(n)$ of spanning trees in the undirected circulant graphs $C_{n}(s_1,s_2,\ldots,s_k)$ and $C_{2n}(s_1,s_2,\ldots,s_k,n).$ Also, we prove that in both cases the number of spanning trees can be represented in the form $\tau(n)=p \,n \,a(n)^2,$ where $a(n)$ is an integer sequence and $p$ is a prescribed natural number depending on the parity of $n.$ Finally, we find an asymptotic formula for $\tau(n)$ through the Mahler measure of the associated Laurent polynomial $L(z)=2k-\sum\limits_{i=1}^k(z^{s_i}+z^{-s_i}).$


title: On the Lang-Trotter conjecture for two elliptic curves

abstract: Following Lang and Trotter we describe a probabilistic model that predicts the distribution of primes $p$ with given Frobenius traces at $p$ for two fixed elliptic curves over $\mathbb{Q}$. In addition, we propose explicit Euler product representations for the constant in the predicted asymptotic formula and describe in detail the universal component of this constant. A new feature is that in some cases the $\ell$-adic limits determining the $\ell$-factors of the universal constant, unlike the Lang-Trotter conjecture for a single elliptic curve, do not stabilize. We also prove the conjecture on average over a family of elliptic curves following the work of David, Koukoulopoulos, and Smith.


title: Bandwidth selection for nonparametric modal regression

abstract: In the context of estimating local modes of a conditional density based on kernel density estimators, we show that existing bandwidth selection methods developed for kernel density estimation are unsuitable for mode estimation. We propose two methods to select bandwidths tailored for mode estimation in the regression setting. Numerical studies using synthetic data and a real-life data set are carried out to demonstrate the performance of the proposed methods in comparison with several well received bandwidth selection methods for density estimation.


title: SDSS-IV MaNGA: Constraints on the Conditions for Star Formation in Galaxy Discs

abstract: Regions of disc galaxies with widespread star formation tend to be both gravitationally unstable and self-shielded against ionizing radiation, whereas extended outer discs with little or no star formation tend to be stable and unshielded on average. We explore what drives the transition between these two regimes, specifically whether discs first meet the conditions for self-shielding (parameterized by dust optical depth, $\tau$) or gravitational instability (parameterized by a modified version of Toomre's instability parameters, $Q_{\rm thermal}$, which quantifies the stability of a gas disc that is thermally supported at $T=10^4$ K). We first introduce a new metric formed by the product of these quantities, $Q_{\rm thermal}\tau$, which indicates whether the conditions for disk instability or self-shielding are easier to meet in a given region of a galaxy, and we discuss how $Q_{\rm thermal}\tau$ can be constrained even in the absence of direct gas information. We then analyse a sample of 13 galaxies with resolved gas measurements and find that on average galaxies will reach the threshold for disk instabilities ($Q_{\rm thermal}<1$) before reaching the threshold for self-shielding ($\tau>1$). Using integral field spectroscopic observations of a sample of 236 galaxies from the MaNGA survey, we find that the value of $Q_{\rm thermal}\tau$ in star-forming discs is consistent with similar behavior. These results support a scenario where disc fragmentation and collapse occurs before self-shielding, suggesting that gravitational instabilities are the primary condition for widespread star formation in galaxy discs. Our results support similar conclusions based on recent galaxy simulations.


title: Keyword-based Query Comprehending via Multiple Optimized-Demand Augmentation

abstract: In this paper, we consider the problem of machine reading task when the questions are in the form of keywords, rather than natural language. In recent years, researchers have achieved significant success on machine reading comprehension tasks, such as SQuAD and TriviaQA. These datasets provide a natural language question sentence and a pre-selected passage, and the goal is to answer the question according to the passage. However, in the situation of interacting with machines by means of text, people are more likely to raise a query in form of several keywords rather than a complete sentence. The keyword-based query comprehension is a new challenge, because small variations to a question may completely change its semantical information, thus yield different answers. In this paper, we propose a novel neural network system that consists a Demand Optimization Model based on a passage-attention neural machine translation and a Reader Model that can find the answer given the optimized question. The Demand Optimization Model optimizes the original query and output multiple reconstructed questions, then the Reader Model takes the new questions as input and locate the answers from the passage. To make predictions robust, an evaluation mechanism will score the reconstructed questions so the final answer strike a good balance between the quality of both the Demand Optimization Model and the Reader Model. Experimental results on several datasets show that our framework significantly improves multiple strong baselines on this challenging task.


title: Diophantine equations involving Euler's totient function

abstract: In this paper, we consider the equations involving Euler's totient function $\phi$ and Lucas type sequences. In particular, we prove that the equation $\phi (x^m-y^m)=x^n-y^n$ has no solutions in positive integers $x, y, m, n$ except for the trivial solutions $(x, y, m , n)=(a+1, a, 1, 1)$, where $a$ is a positive integer, and the equation $\phi ((x^m-y^m)/(x-y))=(x^n-y^n)/(x-y)$ has no solutions in positive integers $x, y, m, n$ except for the trivial solutions $(x, y, m , n)=(a, b, 1, 1)$, where $a, b$ are integers with $a>b\ge 1$.


title: Finding all Maximal Area Parallelograms in a Convex Polygon

abstract: Polygon inclusion problems have been studied extensively in geometric optimization. In this paper, we consider the variant of computing the maximum area parallelograms (MAPs) and all the locally maximal area parallelograms (LMAPs) in a given convex polygon. By proving and utilizing several structural properties of the LMAPs, we compute all of them (including all the MAPs) in $O(n^2)$ time, where $n$ denotes the number of edges of the given polygon. In addition, we prove that the LMAPs interleave each other and thus the number of LMAPs is $O(n)$. We discuss applications of our result to, among others, the problem of computing the maximum area centrally-symmetric convex body inside a convex polygon, and the simplest case of the Heilbronn triangle problem.


title: Textured heterogeneity in square artificial spin ice

abstract: We report evidence of spontaneous formation of a heterogeneous network of superdomains in two-dimensional square artificial spin ice nanostructures in externally applied magnetic fields. The magnetic heterogeneity is locally disordered but has a zig-zag texture at longer length scales. Resonant coherent soft-x-ray scattering off such textures give rise to unique internal structure in Bragg peaks. Our result shows that the macroscopic magnetic texture is derived from the microscopic structure of the Dirac strings.


title: Supernova forecast with strong lensing

abstract: In the coming LSST era, we will observe $\mathcal{O}(100)$ of lensed supernovae (SNe). In this paper, we investigate possibility for predicting time and sky position of a supernova using strong lensing. We find that it will be possible to predict the time and position of the fourth image of SNe which produce four images by strong lensing, with combined information from the three previous images. It is useful to perform multi-messenger observations of the very early phase of supernova explosions including the shock breakout.


title: On the denominators of harmonic numbers

abstract: Let $H_n$ be the $n$-th harmonic number and let $v_n$ be its denominator. It is well known that $v_n$ is even for every integer $n\ge 2$. In this paper, we study the properties of $v_n$. One of our results is: the set of positive integers $n$ such that $v_n$ is divisible by the least common multiple of $1, 2, \cdots, \lfloor {n^{1/4}}\rfloor $ has density one. In particular, for any positive integer $m$, the set of positive integers $n$ such that $v_n$ is divisible by $m$ has density one.


title: Surface Imaging of Proxima b and Other Exoplanets: Topography, Biosignatures, and Artificial Mega-Structures

abstract: Seeing oceans, continents, quasi-static weather, and other surface features on exoplanets may allow us to detect and characterize life outside the solar system. The Proxima b planet resides within the stellar habitable zone allowing for liquid water on its surface, and it may be Earth-like. However, even the largest planned telescopes will not be able to resolve its surface features directly. Here, we demonstrate an inversion technique to image indirectly exoplanet surfaces using observed unresolved reflected light variations over the course of the exoplanets orbital and axial rotation: ExoPlanet Surface Imaging (EPSI). We show that the reflected light curve contains enough information to detect both longitudinal and latitudinal structures and to map exoplanet surface features. We demonstrate this using examples of Solar system planets and moons as well as simulated planets with Earth-like life and artificial megastructures. We also describe how it is possible to infer the planet and orbit geometry from light curves. In particular, we show how albedo maps of Proxima b can be successfully reconstructed for tidally locked, resonance, and unlocked axial and orbital rotation. Such albedo maps obtained in different wavelength passbands can provide "photographic" views of distant exoplanets. We estimate the signal-to-noise ratio necessary for successful inversions and analyse telescope and detector requirements necessary for the first surface images of Proxima b and other nearby exoplanets.


title: On additive representation functions

abstract: Let $A$ be an infinite set of natural numbers. For $n\in \mathbb{N}$, let $r(A, n)$ denote the number of solutions of the equation $n=a+b$ with $a, b\in A, a\le b$. Let $|A(x)|$ be the number of integers in $A$ which are less than or equal to $x$. In this paper, we prove that, if $r(A, n)\not= 1$ for all sufficiently large integers $n$, then $|A(x)|> \frac 12 (\log x/\log\log x)^2$ for all sufficiently large $x$.


title: Single-mode instability in standing-wave lasers: The quantum cascade laser as a self-pumped parametric oscillator

abstract: We report the observation of a clear single-mode instability threshold in continuous-wave Fabry-Perot quantum cascade lasers (QCLs). The instability is characterized by the appearance of sidebands separated by tens of free spectral ranges (FSR) from the first lasing mode, at a pump current not much higher than the lasing threshold. As the current is increased, higher-order sidebands appear that preserve the initial spacing, and the spectra are suggestive of harmonically phase-locked waveforms. We present a theory of the instability that applies to all homogeneously broadened standing-wave lasers. The low instability threshold and the large sideband spacing can be explained by the combination of an unclamped, incoherent Lorentzian gain due to the population grating, and a coherent parametric gain caused by temporal population pulsations that changes the spectral gain line shape. The parametric term suppresses the gain of sidebands whose separation is much smaller than the reciprocal gain recovery time, while enhancing the gain of more distant sidebands. The large gain recovery frequency of the QCL compared to the FSR is essential to observe this parametric effect, which is responsible for the multiple-FSR sideband separation. We predict that by tuning the strength of the incoherent gain contribution, for example by engineering the modal overlap factors and the carrier diffusion, both amplitude-modulated (AM) or frequency-modulated emission can be achieved from QCLs. We provide initial evidence of an AM waveform emitted by a QCL with highly asymmetric facet reflectivities, thereby opening a promising route to ultrashort pulse generation in the mid-infrared. Together, the experiments and theory clarify a deep connection between parametric oscillation in optically pumped microresonators and the single-mode instability of lasers, tying together literature from the last 60 years.


title: Two-dimensional MoS$_2$ electromechanical actuators

abstract: We investigate electromechanical properties of two-dimensional MoS$_2$ monolayers in the 1H, 1T, and 1T$^\prime$ structures as a function of charge doping by using density functional theory. We find isotropic elastic moduli in the 1H and 1T structures, while the 1T$^\prime$ structure exhibits an anisotropic elastic modulus. Moreover, the 1T structure is shown to have a negative Poisson's ratio, while Poisson's ratios of the 1H and 1T$^\prime$ are positive. By charge doping, the monolayer MoS$_2$ shows a reversibly strain and work density per cycle ranging from -0.68% to 2.67% and from 4.4 to 36.9 MJ/m$^3$, respectively, making them suitable for applications in electromechanical actuators. Stress generated is also examined in this work and we find that 1T and 1T$^\prime$ MoS$_2$ monolayers relatively have better performance than 1H MoS$_2$ monolayer. We argue that such excellent electromechanical performance originate from the electrical conductivity of the metallic 1T and semimetallic 1T$^\prime$ structures high Young's modulus of about $150-200$ GPa.


title: A generalized concatenation construction for q-ary 1-perfect codes

abstract: We consider perfect 1-error correcting codes over a finite field with $q$ elements (briefly $q$-ary 1-perfect codes). In this paper, a generalized concatenation construction for $q$-ary 1-perfect codes is presented that allows us to construct $q$-ary 1-perfect codes of length $(q - 1)nm + n + m$ from the given $q$-ary 1-perfect codes of length $n =(q^{s_1} - 1) / (q - 1)$ and $m = (q^{s_2} - 1) / (q - 1)$, where $s_1, s_2$ are natural numbers not less than two. This construction allows us to also construct $q$-ary codes with parameters $(q^{s_1 + s_2}, q^{q^{s_1 + s_2} - (s_1 + s_2) - 1}, 3)_q$ and can be regarded as a $q$-ary analogue of the well-known Phelps construction.


title: Two-loop binding corrections to the electron gyromagnetic factor

abstract: We compute corrections to the gyromagnetic factor of an electron bound in a hydrogen-like ion at order alpha^2(Z alpha)^5. This result removes a major uncertainty in predictions for silicon and carbon ions, used to determine the atomic mass of the electron.


title: Modeling synthetic spectra for transiting extrasolar giant planets: detectability of H$_2$S and PH$_3$ with JWST

abstract: JWST$'$s large aperture and wide wavelength coverage will enable it to collect the highest quality transit spectra observed so far. For exoplanetary atmospheres we expect to retrieve the abundance of the most abundant molecules, such as H$_2$O, CO, and CH$_4$. Other molecules, such as H$_2$S and PH$_3$, have been observed in Jupiter and Saturn but their chemistry and detectability in strongly irradiated planets is highly unknown. In this paper, we make the first effort to study their spectral features in solar composition atmospheres, and evaluate their detectability with JWST. We model the chemistry of phosphorus and sulfur in solar composition atmospheres. Our model includes the effect of vertical transport. Photochemistry effects are not included in our calculations. Using the abundance profiles, we model the JWST transmission and emission spectra for K=6.8 G-type star and for planets with cloud-free solar composition atmospheres. We find PH$_3$ is detectable at 3 sigma from transmission spectra of the simulated atmosphere with $T_{\rm eq}$ $<$ 500K using the NIRCam LW grism F444W mode with a total observing time of 28.8 hrs. H$_2$S is detectable at 3 sigma in the transmission and emission spectra for the simulated planet with $T_{\rm eq}$ $>$ 1500K using the NIRCam LW grism F322W2 mode with a total observing time of 24.0 hrs. Our results specifically highlight the importance of including H$_2$S for future abundances retrieval with JWST. The presence of clouds and hazes challenges the detections of PH$_3$ and H$_2$S, but H$_2$S features are still expected to be present in the emission spectra.


title: Grand Fujii-Fujii-Nakamoto operator inequality dealing with operator order and operator chaotic order

abstract: In this paper, we shall prove that a grand Fujii-Fujii-Nakamoto operator inequality implies operator order and operator chaotic order under different conditions.


title: Functional Connectomics from Data: Probabilistic Graphical Models for Neuronal Network of C. elegans

abstract: We propose a data-driven approach to represent neuronal network dynamics as a Probabilistic Graphical Model (PGM). Our approach learns the PGM structure by employing dimension reduction to network response dynamics evoked by stimuli applied to each neuron separately. The outcome model captures how stimuli propagate through the network and thus represents functional dependencies between neurons, i.e., functional connectome. The benefit of using a PGM as the functional connectome is that posterior inference can be done efficiently and circumvent the complexities in direct inference of response pathways in dynamic neuronal networks. In particular, posterior inference reveals the relations between known stimuli and downstream neurons or allows to query which stimuli are associated with downstream neurons. For validation and as an example for our approach we apply our methodology to a model of Caenorhabiditis elegans nervous system which structure and dynamics are well-studied. From its dynamical model we collect time series of the network response and use singular value decomposition to obtain a low-dimensional projection of the time series data. We then extract dominant patterns in each data matrix to get pairwise dependency information and create a graphical model for the full somatic nervous system. The PGM enables us to obtain and verify underlying neuronal pathways dominant for known behavioral scenarios and to detect possible pathways for novel scenarios.


title: Domino tilings of the expanded Aztec diamond

abstract: The expanded Aztec diamond is a generalized version of the Aztec diamond, with an arbitrary number of long columns and long rows in the middle. In this paper, we count the number of domino tilings of the expanded Aztec diamond. The exact number of domino tilings is given by recurrence relations of state matrices by virtue of the state matrix recursion algorithm, recently developed by the author to solve various two-dimensional regular lattice model enumeration problems.


title: Stability of Minkowski space and polyhomogeneity of the metric

abstract: We study the nonlinear stability of the $(3+1)$-dimensional Minkowski spacetime as a solution of the Einstein vacuum equation. Similarly to our previous work on the stability of cosmological black holes, we construct the solution of the nonlinear initial value problem using an iteration scheme in which we solve a linearized equation globally at each step; we use a generalized harmonic gauge and implement constraint damping to fix the geometry of null infinity. The linear analysis is largely based on energy and vector field methods originating in work by Klainerman. The weak null condition of Lindblad and Rodnianski arises naturally as a nilpotent coupling of certain metric components in a linear model operator at null infinity. Upon compactifying $\mathbb{R}^4$ to a manifold with corners, with boundary hypersurfaces corresponding to spacelike, null, and timelike infinity, we show, using the framework of Melrose's b-analysis, that polyhomogeneous initial data produce a polyhomogeneous spacetime metric. Finally, we relate the Bondi mass to a logarithmic term in the expansion of the metric at null infinity and prove the Bondi mass loss formula.


title: Hidden Antipolar Order Parameter and Entangled Neel-Type Charged Domain Walls in Hybrid Improper Ferroelectrics

abstract: Hybrid improper ferroelectricity (HIF) denotes a new class of polar instability by the mixture of two octahedral-distortion modes and can feature the coexistence of abundant head-to-head and tail-to-tail polar domains, of which the domain walls tend to be charged due to the respective screening charges with an opposite sign. However, no such coexisting carriers are available in the materials. Using group-theoretical, microscopic, and spectroscopic analyses, we established the existence of hidden antipolar order parameter in model HIF (Ca,Sr)3Ti2O7 by the condensation of a weak, previously unnoticed antipolar lattice instability, turning the order-parameter spaces to be multicomponent with the distinct polar-antipolar intertwining and accompanied formation of Neel-type twin-like antipolar domain walls (few nm) between the head-to-head and tail-to-tail domains. The finite-width Neel walls and correlated domain topology inherently lift the polar divergences between the domains, casting an emergent exemplification of charged domain-wall screening by an antipolar ingredient. Comparisons to topological defects in improper-ferroelectrics hexagonal manganites were discussed.


title: Stochastic Modeling and Forecast of Hydrological Contributions in the Colombian Electric System

abstract: In this paper as show that hydrological contributions in the colombian electrical system during the period between 2004 and 2016 have a periodic dynamic, with fundamental periods that are repeated every three years and that tend to oscillate around a long term average, in this context, such contributions can be characterized by mean reversion stochastic processes with periodic functional tendency. The objective of this paper is modeling and forecast the dynamic behavioral of hydrological contributions in the colombian electric system. A description of climate and hydrology in Colombia is carried out, as well as an analysis of periodicity and distributional properties of the data. A Gaussian estimation is performed, which allows to find all the constant and functional parameters from the historical data with daily frequency. The forecasts of the hydrological contributions are graphically illustrated for a period of three years and two alternatives to construct forecast fringes are proposed.