title: Generalizing the Convolution Operator in Convolutional Neural Networks

abstract: Convolutional neural networks have become a main tool for solving many machine vision and machine learning problems. A major element of these networks is the convolution operator which essentially computes the inner product between a weight vector and the vectorized image patches extracted by sliding a window in the image planes of the previous layer. In this paper, we propose two classes of surrogate functions for the inner product operation inherent in the convolution operator and so attain two generalizations of the convolution operator. The first one is the class of positive definite kernel functions where their application is justified by the kernel trick. The second one is the class of similarity measures defined based on a distance function. We justify this by tracing back to the basic idea behind the neocognitron which is the ancestor of CNNs. Both methods are then further generalized by allowing a monotonically increasing function to be applied subsequently. Like any trainable parameter in a neural network, the template pattern and the parameters of the kernel/distance function are trained with the back-propagation algorithm. As an aside, we use the proposed framework to justify the use of sine activation function in CNNs. Our experiments on the MNIST dataset show that the performance of ordinary CNNs can be achieved by generalized CNNs based on weighted L1/L2 distances, proving the applicability of the proposed generalization of the convolutional neural networks.


title: Remote sensing of forests using discrete return airborne LiDAR

abstract: Airborne discrete return light detection and ranging (LiDAR) point clouds covering forested areas can be processed to segment individual trees and retrieve their morphological attributes. Segmenting individual trees in natural deciduous forests however remained a challenge because of the complex and multi-layered canopy. In this chapter, we present (i) a robust segmentation method that avoids a priori assumptions about the canopy structure, (ii) a vertical canopy stratification procedure that improves segmentation of understory trees, (iii) an occlusion model for estimating the point density of each canopy stratum, and (iv) a distributed computing approach for efficient processing at the forest level. When applied to the University of Kentucky Robinson Forest, the segmentation method detected about 90% of overstory and 47% of understory trees with over-segmentation rates of 14% and 2%. Stratifying the canopy improved the detection rate of understory trees to 68% at the cost of increasing their over-segmentations to 16%. According to our occlusion model, a point density of ~170 pt/m-sqr is needed to segment understory trees as accurately as overstory trees. Lastly, using the distributed approach, we segmented about two million trees in the 7,440-ha forest in 2.5 hours using 192 processors, which is 167 times faster than using a single processor. Keywords: individual tree segmentation, multi-layered stand, vertical canopy stratification, segmentation evaluation, point density, canopy occlusion effect, big data, distributed computing.


title: Guided Co-training for Large-Scale Multi-View Spectral Clustering

abstract: In many real-world applications, we have access to multiple views of the data, each of which characterizes the data from a distinct aspect. Several previous algorithms have demonstrated that one can achieve better clustering accuracy by integrating information from all views appropriately than using only an individual view. Owing to the effectiveness of spectral clustering, many multi-view clustering methods are based on it. Unfortunately, they have limited applicability to large-scale data due to the high computational complexity of spectral clustering. In this work, we propose a novel multi-view spectral clustering method for large-scale data. Our approach is structured under the guided co-training scheme to fuse distinct views, and uses the sampling technique to accelerate spectral clustering. More specifically, we first select $p$ ($\ll n$) landmark points and then approximate the eigen-decomposition accordingly. The augmented view, which is essential to guided co-training process, can then be quickly determined by our method. The proposed algorithm scales linearly with the number of given data. Extensive experiments have been performed and the results support the advantage of our method for handling the large-scale multi-view situation.


title: Unmixing dynamic PET images with variable specific binding kinetics

abstract: To analyze dynamic positron emission tomography (PET) images, various generic multivariate data analysis techniques have been considered in the literature, such as principal component analysis (PCA), independent component analysis (ICA), factor analysis and nonnegative matrix factorization (NMF). Nevertheless, these conventional approaches neglect any possible nonlinear variations in the time activity curves describing the kinetic behavior of tissues with specific binding, which limits their ability to recover a reliable, understandable and interpretable description of the data. This paper proposes an alternative analysis paradigm that accounts for spatial fluctuations in the exchange rate of the tracer between a free compartment and a specifically bound ligand compartment. The method relies on the concept of linear unmixing, usually applied on the hyperspectral domain, which combines NMF with a sum-to-one constraint that ensures an exhaustive description of the mixtures. The spatial variability of the signature corresponding to the specific binding tissue is explicitly modeled through a perturbed component. The performance of the method is assessed on both synthetic and real data and is shown to compete favorably when compared to other conventional analysis methods. The proposed method improved both factor estimation and proportions extraction for specific binding. Modeling the variability of the specific binding factor has a strong potential impact for dynamic PET image analysis.


title: Magnetism from intermetallics and perovskite oxides

abstract: This work has been presented by RJCV to obtain his PhD degree at Fluminense Federal University, in March of 2017. We focused on the synthesis of compounds and then on their magneto-strucutral characterization; mainly due to the interplay of these physical properties. We have prepared intermetallic alloys (including Heusler alloys) and perovskite oxides (manganites and cobaltites); in bulk and nanoparticles. A thorough analysis of the influence of the morphology and crystal structure on the magnetic properties of these compounds is addressed.


title: A comment on the paper Prediction of Kidney Function from Biopsy Images using Convolutional Neural Networks

abstract: This letter presente a comment on the paper Prediction of Kidney Function from Biopsy Images using Convolutional Neural Networks by Ledbetter et al. (2017)


title: Extremely Low Bit Neural Network: Squeeze the Last Bit Out with ADMM

abstract: Although deep learning models are highly effective for various learning tasks, their high computational costs prohibit the deployment to scenarios where either memory or computational resources are limited. In this paper, we focus on compressing and accelerating deep models with network weights represented by very small numbers of bits, referred to as extremely low bit neural network. We model this problem as a discretely constrained optimization problem. Borrowing the idea from Alternating Direction Method of Multipliers (ADMM), we decouple the continuous parameters from the discrete constraints of network, and cast the original hard problem into several subproblems. We propose to solve these subproblems using extragradient and iterative quantization algorithms that lead to considerably faster convergency compared to conventional optimization methods. Extensive experiments on image recognition and object detection verify that the proposed algorithm is more effective than state-of-the-art approaches when coming to extremely low bit neural network.


title: Feature Extraction via Recurrent Random Deep Ensembles and its Application in Gruop-level Happiness Estimation

abstract: This paper presents a novel ensemble framework to extract highly discriminative feature representation of image and its application for group-level happpiness intensity prediction in wild. In order to generate enough diversity of decisions, n convolutional neural networks are trained by bootstrapping the training set and extract n features for each image from them. A recurrent neural network (RNN) is then used to remember which network extracts better feature and generate the final feature representation for one individual image. Several group emotion models (GEM) are used to aggregate face fea- tures in a group and use parameter-optimized support vector regressor (SVR) to get the final results. Through extensive experiments, the great effectiveness of the proposed recurrent random deep ensembles (RRDE) is demonstrated in both structural and decisional ways. The best result yields a 0.55 root-mean-square error (RMSE) on validation set of HAPPEI dataset, significantly better than the baseline of 0.78.


title: Full-Network Embedding in a Multimodal Embedding Pipeline

abstract: The current state-of-the-art for image annotation and image retrieval tasks is obtained through deep neural networks, which combine an image representation and a text representation into a shared embedding space. In this paper we evaluate the impact of using the Full-Network embedding in this setting, replacing the original image representation in a competitive multimodal embedding generation scheme. Unlike the one-layer image embeddings typically used by most approaches, the Full-Network embedding provides a multi-scale representation of images, which results in richer characterizations. To measure the influence of the Full-Network embedding, we evaluate its performance on three different datasets, and compare the results with the original multimodal embedding generation scheme when using a one-layer image embedding, and with the rest of the state-of-the-art. Results for image annotation and image retrieval tasks indicate that the Full-Network embedding is consistently superior to the one-layer embedding. These results motivate the integration of the Full-Network embedding on any multimodal embedding generation scheme, something feasible thanks to the flexibility of the approach.


title: Representation Learning on Large and Small Data

abstract: Deep learning owes its success to three key factors: scale of data, enhanced models to learn representations from data, and scale of computation. This book chapter presented the importance of the data-driven approach to learn good representations from both big data and small data. In terms of big data, it has been widely accepted in the research community that the more data the better for both representation and classification improvement. The question is then how to learn representations from big data, and how to perform representation learning when data is scarce. We addressed the first question by presenting CNN model enhancements in the aspects of representation, optimization, and generalization. To address the small data challenge, we showed transfer representation learning to be effective. Transfer representation learning transfers the learned representation from a source domain where abundant training data is available to a target domain where training data is scarce. Transfer representation learning gave the OM and melanoma diagnosis modules of our XPRIZE Tricorder device (which finished $2^{nd}$ out of $310$ competing teams) a significant boost in diagnosis accuracy.


title: The derivative of the conjugacy for the pair of tent-like maps from an interval into itself

abstract: We consider in this article the properties of the topological conjugacy of the piecewise linear unimodal maps $g:\, [0,\, 1]\rightarrow [0,\, 1]$, all whose kinks belong to the complete pre-image of $0$. We call such maps firm carcass maps. We prove that every firm carcass maps $g_1$ and $g_2$ are topologically conjugated. For the conjugacy $h$ such that $h\circ g_1 = g_2\circ h$ we denote $\{ h_n, n\geq 1\}$ the piecewise linear approximations of $h$, whose graphs connect the points $\{ (x, h(x)),\ g_1^n(x)=0\}$. For any $x\in [0,\, 1]$ we reduce the question about the value of $h'(x)$ to the properties of the sequence $\{h_n'(x),\, n\geq 1\}$. We prove that each conjugacy of firm carcass maps either has the length 2, or is piecewise linear.


title: SAR Target Recognition Using the Multi-aspect-aware Bidirectional LSTM Recurrent Neural Networks

abstract: The outstanding pattern recognition performance of deep learning brings new vitality to the synthetic aperture radar (SAR) automatic target recognition (ATR). However, there is a limitation in current deep learning based ATR solution that each learning process only handle one SAR image, namely learning the static scattering information, while missing the space-varying information. It is obvious that multi-aspect joint recognition introduced space-varying scattering information should improve the classification accuracy and robustness. In this paper, a novel multi-aspect-aware method is proposed to achieve this idea through the bidirectional Long Short-Term Memory (LSTM) recurrent neural networks based space-varying scattering information learning. Specifically, we first select different aspect images to generate the multi-aspect space-varying image sequences. Then, the Gabor filter and three-patch local binary pattern (TPLBP) are progressively implemented to extract a comprehensive spatial features, followed by dimensionality reduction with the Multi-layer Perceptron (MLP) network. Finally, we design a bidirectional LSTM recurrent neural network to learn the multi-aspect features with further integrating the softmax classifier to achieve target recognition. Experimental results demonstrate that the proposed method can achieve 99.9% accuracy for 10-class recognition. Besides, its anti-noise and anti-confusion performance are also better than the conventional deep learning based methods.


title: Transformable Biomimetic Liquid Metal Chameleon

abstract: Liquid metal (LM) is of current core interest for a wide variety of newly emerging areas. However, the functional materials thus made so far by LM only could display a single silver-white appearance. Here in this study, the new conceptual colorful LM marbles working like transformable biomimetic chameleons were proposed and fabricated from LM droplets through encasing them with fluorescent nano-particles. We demonstrated that this unique LM marble can be manipulated into various stable magnificent appearances as one desires. And it can also splitt and merge among different colors. Such multifunctional LM chameleon is capable of responding to the outside electric-stimulus and realizing shape transformation and discoloration behaviors as well. Further more, the electric-stimuli has been disclosed to be an easy going way to trigger the release of nano/micro-particles from the LM. The present fluorescent biomimetic liquid metal chameleon is expected to offer important opportunities for diverse unconventional applications, especially in a wide variety of functional smart material and color changeable soft robot areas.


title: YBCO Ceramic Nanofibers obtained by the new technique of Solution Blow Spinning

abstract: This study proposes a novel solution blow spinning technique (SBS) for fabricating YBCO ceramic nanofibers. The precursor solutions were obtained from Y, Ba, and Cu metallic acetates (Ac) and poly(vinyl pyrrolidone) (PVP, Mw = 360,000). Ac:PVP concentrations of 1:1 and 5:1 were tested, resulting in ceramic nanofibers with average diameters of 359 and 375 nm, respectively. X-ray diffraction confirmed the formation of a pure phase of YBa2Cu3O7-x. This is the first study to use SBS for fabricating YBCO nanofibers, and this technique shows promise for obtaining high-quality ceramic materials.


title: Event-plane dependent di-hadron correlations with harmonic $v_n$ subtraction in a hydrodynamic model

abstract: In this work, a hydrodynamic study of the di-hadron azimuthal correlations for the Au+Au collisions at 200 GeV is carried out. The correlations are evaluated using the ZYAM method for the centrality windows as well as the transverse momentum range in accordance with the existing data. Event-plane dependence of the correlation is obtained after the subtraction of contributions from the most dominant harmonic coefficients. In particular, the contribution from the triangular flow, $v_3$, is removed from the proper correlations following the procedure implemented by the STAR collaboration. The resultant structure observed in the correlations was sometimes attributed to the mini-jet dynamics, but the present calculations show that a pure hydrodynamic model gives a reasonable agreement with the main feature of the published data. A brief discussion on the physical content of the present findings is presented.


title: Linguistically Motivated Vocabulary Reduction for Neural Machine Translation from Turkish to English

abstract: The necessity of using a fixed-size word vocabulary in order to control the model complexity in state-of-the-art neural machine translation (NMT) systems is an important bottleneck on performance, especially for morphologically rich languages. Conventional methods that aim to overcome this problem by using sub-word or character-level representations solely rely on statistics and disregard the linguistic properties of words, which leads to interruptions in the word structure and causes semantic and syntactic losses. In this paper, we propose a new vocabulary reduction method for NMT, which can reduce the vocabulary of a given input corpus at any rate while also considering the morphological properties of the language. Our method is based on unsupervised morphology learning and can be, in principle, used for pre-processing any language pair. We also present an alternative word segmentation method based on supervised morphological analysis, which aids us in measuring the accuracy of our model. We evaluate our method in Turkish-to-English NMT task where the input language is morphologically rich and agglutinative. We analyze different representation methods in terms of translation accuracy as well as the semantic and syntactic properties of the generated output. Our method obtains a significant improvement of 2.3 BLEU points over the conventional vocabulary reduction technique, showing that it can provide better accuracy in open vocabulary translation of morphologically rich languages.


title: Sub-Nanometer Channels Embedded in Two-Dimensional Materials

abstract: Two-dimensional (2D) materials are among the most promising candidates for next-generation electronics due to their atomic thinness, allowing for flexible transparent electronics and ultimate length scaling. Thus far, atomically-thin p-n junctions, metal-semiconductor contacts, and metal-insulator barriers have been demonstrated. While 2D materials achieve the thinnest possible devices, precise nanoscale control over the lateral dimensions is also necessary. Here, we report the direct synthesis of sub-nanometer-wide 1D MoS2 channels embedded within WSe2 monolayers, using a dislocation-catalyzed approach. The 1D channels have edges free of misfit dislocations and dangling bonds, forming a coherent interface with the embedding 2D matrix. Periodic dislocation arrays produce 2D superlattices of coherent MoS2 1D channels in WSe2. Using molecular dynamics simulations, we have identified other combinations of 2D materials where 1D channels can also be formed. The electronic band structure of these 1D channels offer the promise of carrier confinement in a direct-gap material and charge separation needed to access the ultimate length scales necessary for future electronic applications.


title: RBF Interpolation with CSRBF of Large Data Sets

abstract: This contribution presents a new analysis of properties of the interpolation using Radial Bases Functions (RBF) related to large data sets interpolation. The RBF application is convenient method for scattered d-dimensional interpolation. The RBF methods lead to a solution of linear system of equations and computational complexity of solution is nearly independent of a dimensionality. However, the RBF methods are usually applied for small data sets with a small span geometric coordinates. This contribution explores properties of the RBF interpolation for large data sets and large span of geometric coordinates of the given data sets with regard to expectable numerical stability of computation.


title: Entropy production for ellipsoidal BGK model of the Boltzmann equation

abstract: The ellipsoidal BGK model (ES-BGK) is a generalized version of the BGK model of the Boltzmann equation designed to yield the correct Prandtl number in the Navier-Stokes limit. In this paper, we make two observations on the entropy production functional of the ES-BGK model. First, we show that the Cercignani type estimate holds for the ES-BGK model in the whole range of parameter $-1/2<\nu<1$. Secondly, we observe that the ellipsoidal relaxation operator satisfies an unexpected sign-definite property. Some implications of these observations are also discussed.


title: Mass transport from the envelope to the disk of V346 Nor: a case study for the luminosity problem in an FUor-type young eruptive star

abstract: A long-standing open issue of the paradigm of low-mass star formation is the luminosity problem: most protostars are less luminous than theoretically predicted. One possible solution is that the accretion process is episodic. FU Ori-type stars (FUors) are thought to be the visible examples for objects in the high accretion state. FUors are often surrounded by massive envelopes, which replenish the disk material and enable the disk to produce accretion outbursts. However, we have insufficient information on the envelope dynamics in FUors, about where and how mass transfer from the envelope to the disk happens. Here we present ALMA observations of the FUor-type star V346 Nor at 1.3 mm continuum and in different CO rotational lines. We mapped the density and velocity structure of its envelope and analyze the results using channel maps, position-velocity diagrams, and spectro-astrometric methods. We found that V346 Nor is surrounded by gaseous material on 10000 au scale in which a prominent outflow cavity is carved. Within the central $\sim$700 au, the circumstellar matter forms a flattened pseudo-disk where material is infalling with conserved angular momentum. Within $\sim$350 au, the velocity profile is more consistent with a disk in Keplerian rotation around a central star of 0.1 $M_{\odot}$. We determined an infall rate from the envelope onto the disk of 6$\times$10$^{-6}\,M_{\odot}$yr$^{-1}$, a factor of few higher than the quiescent accretion rate from the disk onto the star, hinting for a mismatch between the infall and accretion rates as the cause of the eruption.


title: Spatial Kibble-Zurek mechanism through susceptibilities: the inhomogeneous quantum Ising model case

abstract: We study the quantum Ising model in the transverse inhomogeneous magnetic field. Such a system can be approached numerically through exact diagonalization and analytically through the renormalization group techniques. Basic insights into its physics, however, can be obtained by adopting the Kibble-Zurek theory of non-equilibrium phase transitions to description of spatially inhomogeneous systems at equilibrium. We employ all these approaches and focus on derivatives of longitudinal and transverse magnetizations, which have extrema near the critical point. We discuss how these extrema can be used for locating the critical point and for verification of the Kibble-Zurek scaling predictions in the spatial quench.


title: Dynamics of polymers: classic results and recent developments

abstract: In this chapter we review concepts and theories of polymer dynamics. We think of it as an introduction to the topic for scientists specializing in other subfields of statistical mechanics and condensed matter theory, so, for the readers reference, we start with a short review of the equilibrium static properties of polymer systems. Most attention is paid to the dynamics of unentangled polymer systems, where apart from classical Rouse and Zimm models we review some recent scaling and analytical generalizations. The dynamics of systems with entanglements is also briefly reviewed. Special attention is paid to the discussion of comparatively weakly understood topological states of polymer systems and possible approaches to the description of their dynamics.


title: Renormalization scheme dependence of high-Order perturbative QCD predictions

abstract: Conventionally, one adopts typical momentum flow of a physical observable as the renormalization scale for its perturbative QCD (pQCD) approximant. This simple treatment leads to renormalization scheme-and-scale ambiguities due to the renormalization scheme and scale dependence of the strong coupling and the perturbative coefficients do not exactly cancel at any fixed order. It is believed that those ambiguities will be softened by including more higher-order terms. In the paper, to show how the renormalization scheme dependence changes when more loop terms have been included, we discuss the sensitivity of pQCD prediction on the scheme parameters by using the scheme-dependent $\{\beta_{m \geq 2}\}$-terms. We adopt two four-loop examples, $e^+ e^- \to {\rm hadrons}$ and $\tau$ decays into hadrons, for detailed analysis. Our results show that under the conventional scale setting, by including more-and-more loop terms, the scheme dependence of the pQCD prediction cannot be reduced as efficiently as that of the scale dependence. Thus a proper scale-setting approach should be important to reduce the scheme dependence. We observe that the principle of minimum sensitivity could be such a scale-setting approach, which provides a practical way to achieve optimal scheme and scale by requiring the pQCD approximate be independent to the "unphysical" theoretical conventions.


title: Learning Audio - Sheet Music Correspondences for Score Identification and Offline Alignment

abstract: This work addresses the problem of matching short excerpts of audio with their respective counterparts in sheet music images. We show how to employ neural network-based cross-modality embedding spaces for solving the following two sheet music-related tasks: retrieving the correct piece of sheet music from a database when given a music audio as a search query; and aligning an audio recording of a piece with the corresponding images of sheet music. We demonstrate the feasibility of this in experiments on classical piano music by five different composers (Bach, Haydn, Mozart, Beethoven and Chopin), and additionally provide a discussion on why we expect multi-modal neural networks to be a fruitful paradigm for dealing with sheet music and audio at the same time.


title: Few-layer Tellurium: one-dimensional-like layered elementary semiconductor with striking physical properties

abstract: Few-layer Tellurium, an elementary semiconductor, succeeds most of striking physical properties that black phosphorus (BP) offers and could be feasibly synthesized by simple solution-based methods. It is comprised of non-covalently bound parallel Te chains, among which covalent-like feature appears. This feature is, we believe, another demonstration of the previously found covalent-like quasi-bonding (CLQB) where wavefunction hybridization does occur. The strength of this inter-chain CLQB is comparable with that of intra-chain covalent bonding, leading to closed stability of several Te allotropes. It also introduces a tunable bandgap varying from nearly direct 0.31 eV (bulk) to indirect 1.17 eV (2L) and four (two) complex, highly anisotropic and layer-dependent hole (electron) pockets in the first Brillouin zone. It also exhibits an extraordinarily high hole mobility (~10$^5$ cm$^2$/Vs) and strong optical absorption along the non-covalently bound direction, nearly isotropic and layer-dependent optical properties, large ideal strength over 20%, better environmental stability than BP and unusual crossover of force constants for interlayer shear and breathing modes. All these results manifest that the few-layer Te is an extraordinary-high-mobility, high optical absorption, intrinsic-anisotropy, low-cost-fabrication, tunable bandgap, better environmental stability and nearly direct bandgap semiconductor. This "one-dimension-like" few-layer Te, together with other geometrically similar layered materials, may promote the emergence of a new family of layered materials.


title: The Electronics, Trigger and Data Acquisition System for the Liquid Argon Time Projection Chamber of the DarkSide-50 Search for Dark Matter

abstract: The DarkSide-50 experiment at the Laboratori Nazionali del Gran Sasso is a search for dark matter using a dual phase time projection chamber with 50 kg of low radioactivity argon as target. Light signals from interactions in the argon are detected by a system of 38 photo-multiplier tubes (PMTs), 19 above and 19 below the TPC volume inside the argon cryostat. We describe the electronics which processes the signals from the photo-multipliers, the trigger system which identifies events of interest, and the data-acquisition system which records the data for further analysis. The electronics include resistive voltage dividers on the PMTs, custom pre-amplifiers mounted directly on the PMT voltage dividers in the liquid argon, and custom amplifier/discriminators (at room temperature). After amplification, the PMT signals are digitized in CAEN waveform digitizers, and CAEN logic modules are used to construct the trigger, the data acquisition system for the TPC is based on the Fermilab "artdaq" software. The system has been in operation since early 2014.


title: Bearing fault diagnosis under varying working condition based on domain adaptation

abstract: Traditional intelligent fault diagnosis of rolling bearings work well only under a common assumption that the labeled training data (source domain) and unlabeled testing data (target domain) are drawn from the same distribution. When the distribution changes, most fault diagnosis models need to be rebuilt from scratch using newly recollected labeled training data. However, it is expensive or impossible to annotate huge amount of training data to rebuild such new model. Meanwhile, large amounts of labeled training data have not been fully utilized yet, which is apparently a waste of resources. As one of the important research directions of transfer learning, domain adaptation (DA) typically aims at minimizing the differences between distributions of different domains in order to minimize the cross-domain prediction error by taking full advantage of information coming from both source and target domains. In this paper, we present one of the first studies on unsupervised DA in the field of fault diagnosis of rolling bearings under varying working conditions and a novel diagnosis strategy based on unsupervised DA using subspace alignment (SA) is proposed. After processed by unsupervised DA with SA, the distributions of training data and testing data become close and the classifier trained on training data can be used to classify the testing data. Experimental results on the 60 domain adaptation diagnosis problems under varying working condition in Case Western Reserve benchmark data and 12 domain adaptation diagnosis problems under varying working conditions in our new data are given to demonstrate the effectiveness of the proposed method. The proposed methods can effectively distinguish not only bearing faults categories but also fault severities.


title: Penalty Alternating Direction Methods for Mixed-Integer Optimization: A New View on Feasibility Pumps

abstract: Feasibility pumps are highly effective primal heuristics for mixed-integer linear and nonlinear optimization. However, despite their success in practice there are only few works considering their theoretical properties. We show that feasibility pumps can be seen as alternating direction methods applied to special reformulations of the original problem, inheriting the convergence theory of these methods. Moreover, we propose a novel penalty framework that encompasses this alternating direction method, which allows us to refrain from random perturbations that are applied in standard versions of feasibility pumps in case of failure. We present a convergence theory for the new penalty based alternating direction method and compare the new variant of the feasibility pump with existing versions in an extensive numerical study for mixed-integer linear and nonlinear problems.


title: Entropy based fingerprint for local crystalline order

abstract: We introduce a new fingerprint that allows distinguishing between liquid-like and solid-like atomic environments. This fingerprint is based on an approximate expression for the entropy projected on individual atoms. When combined with a local enthalpy, this fingerprint acquires an even finer resolution and it is capable of discriminating between different crystal structures.


title: Quantum Privacy-Preserving Perceptron

abstract: With the extensive applications of machine learning, the issue of private or sensitive data in the training examples becomes more and more serious: during the training process, personal information or habits may be disclosed to unexpected persons or organisations, which can cause serious privacy problems or even financial loss. In this paper, we present a quantum privacy-preserving algorithm for machine learning with perceptron. There are mainly two steps to protect original training examples. Firstly when checking the current classifier, quantum tests are employed to detect data user's possible dishonesty. Secondly when updating the current classifier, private random noise is used to protect the original data. The advantages of our algorithm are: (1) it protects training examples better than the known classical methods; (2) it requires no quantum database and thus is easy to implement.


title: Effects of developmental variability on the dynamics and self-organization of cell populations

abstract: We report experimental and theoretical results on spatiotemporal pattern formation in cell populations, where the parameters vary in space and time due to mechanisms intrinsic to the system, namely Dictyostelium discoideum (D.d.) in the starvation phase. We find that different patterns are formed when the populations are initialized at different developmental stages, or, when populations at different initial developmental stages are mixed. The experimentally observed patterns can be understood with a modified Kessler-Levine model that takes into account the initial spatial heterogeneity of the cell populations and a developmental path introduced by us, i.e., the time dependence of the various biochemical parameters. The dynamics of the parameter agree with known biochemical studies. Most importantly the modified model reproduces not only our results, but also the observations of an independent experiment published earlier. This shows that pattern formation can be used to understand and quantify the temporal evolution of the system parameters.


title: Initial Results from Fitting Resolved Modes using HMI Intensity Observations

abstract: The HMI project recently started processing the continuum intensity images following global helioseismology procedures similar to those used to process the velocity images. The spatial decomposition of these images has produced time series of spherical harmonic coefficients for degrees up to $\ell=300$, using a different apodization than the one used for velocity observations. The first 360 days of observations were processed and made available. In this paper I present initial results from fitting these time series using my state of the art fitting methodology and compare the derived mode characteristics to those estimated using co-eval velocity observations.


title: On the separability of the partial skew groupoid ring

abstract: Given a partial (resp. a global) action $\alpha$ of a connected finite groupoid $G$ on a ring $A$, we determine necessary and sufficient conditions for the partial (resp. global) skew groupoid ring $A\star_{\alpha} G$ to be a separable extension of $A$.


title: Higgs amplitudes from $\mathcal{N}=4$ super Yang-Mills theory

abstract: Higgs plus multi-gluon amplitudes in QCD can be computed in an effective Lagrangian description. In the infinite top-mass limit, an amplitude with a Higgs and $n$ gluons is computed by the form factor of the operator ${\rm Tr}\, F^2$. Up to two loops and for three gluons, its maximally transcendental part is captured entirely by the form factor of the protected stress tensor multiplet operator $\mathcal{T}_2$ in $\mathcal{N}=4$ supersymmetric Yang-Mills theory. The next order correction involves the calculation of the form factor of the higher-dimensional, trilinear operator ${\rm Tr}\, F^3$. We present explicit results at two loops for three gluons, including the subleading transcendental terms derived from a particular descendant of the Konishi operator that contains ${\rm Tr}\, F^3$. These are expressed in terms of a few universal building blocks already identified in earlier calculations. We show that the maximally transcendental part of this quantity, computed in non-supersymmetric Yang-Mills theory, is identical to the form factor of another protected operator, $\mathcal{T}_3$, in the maximally supersymmetric theory. Our results suggest that the maximally transcendental part of Higgs amplitudes in QCD can be entirely computed through $\mathcal{N}=4$ super Yang-Mills.


title: A fast finite-difference algorithm for topology optimization of permanent magnets

abstract: We present a finite-difference method for the topology optimization of permanent magnets that is based on the FFT accelerated computation of the stray-field. The presented method employs the density approach for topology optimization and uses an adjoint method for the gradient computation. Comparsion to various state-of-the-art finite-element implementations shows a superior performance and accuracy. Moreover, the presented method is very flexible and easy to implement due to various preexisting FFT stray-field implementations that can be used.


title: Fashioning with Networks: Neural Style Transfer to Design Clothes

abstract: Convolutional Neural Networks have been highly successful in performing a host of computer vision tasks such as object recognition, object detection, image segmentation and texture synthesis. In 2015, Gatys et. al [7] show how the style of a painter can be extracted from an image of the painting and applied to another normal photograph, thus recreating the photo in the style of the painter. The method has been successfully applied to a wide range of images and has since spawned multiple applications and mobile apps. In this paper, the neural style transfer algorithm is applied to fashion so as to synthesize new custom clothes. We construct an approach to personalize and generate new custom clothes based on a users preference and by learning the users fashion choices from a limited set of clothes from their closet. The approach is evaluated by analyzing the generated images of clothes and how well they align with the users fashion style.


title: The complex null string, Galilean conformal algebra and scattering equations

abstract: The scattering equation formalism for scattering amplitudes, and its stringy incarnation, the ambitwistor string, remains a mysterious construction. In this paper, we pursue the study a gauged-unfixed version of the ambitwistor string known as the null string. We explore the following three aspects in detail; its complexification, gauge fixing, and amplitudes. We first study the complexification of the string; the associated symmetries and moduli, and connection to the ambitwistor string. We then look in more details at the leftover symmetry algebra of the string, called Galilean conformal algebra; we study its local and global action and gauge-fixing. We finish by presenting an operator formalism, that we use to compute tree-level scattering amplitudes based on the scattering equations and a one-loop partition function. These results hopefully will open the way to understand conceptual questions related to the loop expansion in these twistor-like string models.


title: Study of Minor Actinides Transmutation in PWR MOX fuel

abstract: The management of long-lived radionuclides in spent fuel is a key issue to achieve the closed nuclear fuel cycle and the sustainable development of nuclear energy. Partitioning-Transmutation is supposed to be an efficient method to treat the long-lived radionuclides in spent fuel. Some Minor Actinides (MAs) have very long half-lives among the radionuclides in the spent fuel. Accordingly, the study of MAs transmutation is a significant work for the post-processing of spent fuel. In the present work, the transmutations in Pressurized Water Reactor (PWR) mixed oxide (MOX) fuel are investigated through the Monte Carlo based code RMC. Two kinds of MAs, $^{237}$Np and five MAs ($^{237}$Np, $^{241}$Am, $^{243}$Am, $^{244}$Cm and $^{245}$Cm) are incorporated homogeneously into the MOX fuel assembly. The transmutation of MAs is simulated with different initial MOX concentrations. The results indicate an overall nice efficiency of transmutation in both initial MOX concentrations, especially for the two kinds of MAs primarily generated in the UOX fuel, $^{237}$Np and $^{241}$Am. In addition, the inclusion of $^{237}$Np in MOX has no large influence for other MAs, while the transmutation efficiency of $^{237}$Np is excellent. The transmutation of MAs in MOX fuel depletion is expected to be a new, efficient nuclear spent fuel management method for the future nuclear power generation.


title: A Relational Event Approach to Modeling Behavioral Dynamics

abstract: This chapter provides an introduction to the analysis of relational event data (i.e., actions, interactions, or other events involving multiple actors that occur over time) within the R/statnet platform. We begin by reviewing the basics of relational event modeling, with an emphasis on models with piecewise constant hazards. We then discuss estimation for dyadic and more general relational event models using the relevent package, with an emphasis on hands-on applications of the methods and interpretation of results. Statnet is a collection of packages for the R statistical computing system that supports the representation, manipulation, visualization, modeling, simulation, and analysis of relational data. Statnet packages are contributed by a team of volunteer developers, and are made freely available under the GNU Public License. These packages are written for the R statistical computing environment, and can be used with any computing platform that supports R (including Windows, Linux, and Mac)."


title: Measurement of cosmic-ray reconstruction efficiencies in the MicroBooNE LArTPC using a small external cosmic-ray counter

abstract: The MicroBooNE detector is a liquid argon time projection chamber at Fermilab designed to study short-baseline neutrino oscillations and neutrino-argon interaction cross-section. Due to its location near the surface, a good understanding of cosmic muons as a source of backgrounds is of fundamental importance for the experiment. We present a method of using an external 0.5 m (L) x 0.5 m (W) muon counter stack, installed above the main detector, to determine the cosmic-ray reconstruction efficiency in MicroBooNE. Data are acquired with this external muon counter stack placed in three different positions, corresponding to cosmic rays intersecting different parts of the detector. The data reconstruction efficiency of tracks in the detector is found to be $\epsilon_{\mathrm{data}}=(97.1\pm0.1~(\mathrm{stat}) \pm 1.4~(\mathrm{sys}))\%$, in good agreement with the Monte Carlo reconstruction efficiency $\epsilon_{\mathrm{MC}} = (97.4\pm0.1)\%$. This analysis represents a small-scale demonstration of the method that can be used with future data coming from a recently installed cosmic-ray tagger system, which will be able to tag $\approx80\%$ of the cosmic rays passing through the MicroBooNE detector.


title: Temporal Hierarchical Clustering

abstract: We study hierarchical clusterings of metric spaces that change over time. This is a natural geometric primitive for the analysis of dynamic data sets. Specifically, we introduce and study the problem of finding a temporally coherent sequence of hierarchical clusterings from a sequence of unlabeled point sets. We encode the clustering objective by embedding each point set into an ultrametric space, which naturally induces a hierarchical clustering of the set of points. We enforce temporal coherence among the embeddings by finding correspondences between successive pairs of ultrametric spaces which exhibit small distortion in the Gromov-Hausdorff sense. We present both upper and lower bounds on the approximability of the resulting optimization problems.


title: Deep Discrete Supervised Hashing

abstract: Hashing has been widely used for large-scale search due to its low storage cost and fast query speed. By using supervised information, supervised hashing can significantly outperform unsupervised hashing. Recently, discrete supervised hashing and deep hashing are two representative progresses in supervised hashing. On one hand, hashing is essentially a discrete optimization problem. Hence, utilizing supervised information to directly guide discrete (binary) coding procedure can avoid sub-optimal solution and improve the accuracy. On the other hand, deep hashing, which integrates deep feature learning and hash-code learning into an end-to-end architecture, can enhance the feedback between feature learning and hash-code learning. The key in discrete supervised hashing is to adopt supervised information to directly guide the discrete coding procedure in hashing. The key in deep hashing is to adopt the supervised information to directly guide the deep feature learning procedure. However, there have not existed works which can use the supervised information to directly guide both discrete coding procedure and deep feature learning procedure in the same framework. In this paper, we propose a novel deep hashing method, called deep discrete supervised hashing (DDSH), to address this problem. DDSH is the first deep hashing method which can utilize supervised information to directly guide both discrete coding procedure and deep feature learning procedure, and thus enhance the feedback between these two important procedures. Experiments on three real datasets show that DDSH can outperform other state-of-the-art baselines, including both discrete hashing and deep hashing baselines, for image retrieval.


title: Common fixed point theorems in C*-algebra-valued b-metric spaces endowed with a graph and applications

abstract: In this paper, we equip a C*-algebra-valued b-metric spaces with a graph G = (V,E) and establish some common fixed point theorems. Also, some examples in support of our main results are provided. Finally, as applications, existence and uniqueness of solution for a type of operator equation and an integral equation are discussed.


title: Characterization of the structure and cross-shore transport properties of a coastal upwelling filament using three-dimensional finite-size Lyapunov exponents

abstract: The three dimensional structure, dynamics and dispersion characteristics of a simulated upwelling filament in the Iberian upwelling system are analyzed using Lagrangian tools. We used a realistic regional simulation of the western Iberian shelf which is concomitant with an in-situ oceanographic campaign that surveyed the area. We compute 3d fields of finite--size Lyapunov exponents (FSLE) from 3d velocity fields and extract the field's ridges to study the spatial distribution and temporal evolution of the Lagrangian Coherent Structures (LCSs) evolving around the filament. We find that the most intense curtain-like LCSs delimit the boundaries of the whole filamentary structure whose general properties match well the observations. The filament interior is characterized by small dispersion of fluid elements. Furthermore, we identify a weak LCS separating the filament into a warmer vein and a colder filament associated with the interaction of a mesoscale eddy with the upwelling front. The cold upwelled water parcels move along the filament conserving their density. The filament itself is characterized by small dispersion of fluid elements in its interior. The comparison of LCSs with potential temperature and salinity gradient fields shows that the outer limits of the filament coincide with regions of large hydrographic gradients, similar to those observed, explaining the isolation of the interior of the filament with the surrounding waters. We conclude that the Lagrangian analysis used in this work is useful in explaining the dynamics of cross-shore exchanges of materials between coastal regions and the open ocean due to mesoscale processes.


title: On the Shapley value of unrooted phylogenetic trees

abstract: The Shapley value, a solution concept from cooperative game theory, has recently been considered for both unrooted and rooted phylogenetic trees. Here, we focus on the Shapley value of unrooted trees and first revisit the so-called split counts of a phylogenetic tree and the Shapley transformation matrix that allows for the calculation of the Shapley value from the edge lengths of a tree. We show that non-isomorphic trees may have permutation-equivalent Shapley transformation matrices and permutation-equivalent null spaces. This implies that estimating the split counts associated with a tree or the Shapley values of its leaves does not suffice to reconstruct the correct tree topology. We then turn to the use of the Shapley value as a prioritization criterion in biodiversity conservation and compare it to a greedy solution concept. Here, we show that for certain phylogenetic trees, the Shapley value may fail as a prioritization criterion, meaning that the diversity spanned by the top $k$ species (ranked by their Shapley values) cannot approximate the total diversity of all $n$ species.


title: Polynomial mixing under a certain stationary Euler flow

abstract: We study the mixing properties of a scalar $\rho$ advected by a certain incompressible velocity field $u$ on the two dimensional unit ball, which is a stationary radial solution of the Euler equation. The scalar $\rho$ solves the continuity equation with velocity field $u$ and we can measure the degree of mixedness of~$\rho$ with two different scales commonly used in this setting, namely the geometric and the functional mixing scale. We develop a physical space approach well adapted for the quantitative analysis of the decay in time of the geometric mixing scale, which turns out to be polynomial for a large class of initial data. This extends previous results for the functional mixing scale, based on the explicit expression for the solution in Fourier variable, results that are also partially recovered by our approach.


title: Conformal and covariant Z4 formulation of the Einstein equations: strongly hyperbolic first-order reduction and solution with discontinuous Galerkin schemes

abstract: We present a strongly hyperbolic first-order formulation of the Einstein equations based on the conformal and covariant Z4 system (CCZ4) with constraint-violation damping, which we refer to as FO-CCZ4. As CCZ4, this formulation combines the advantages of a conformal and traceless formulation, with the suppression of constraint violations given by the damping terms, but being first order in time and space, it is particularly suited for a discontinuous Galerkin (DG) implementation. The strongly hyperbolic first-order formulation has been obtained by making careful use of first and second-order ordering constraints. A proof of strong hyperbolicity is given for a selected choice of standard gauges via an analytical computation of the entire eigenstructure of the FO-CCZ4 system. The resulting governing partial differential equations system is written in non-conservative form and requires the evolution of 58 unknowns. A key feature of our formulation is that the first-order CCZ4 system decouples into a set of pure ordinary differential equations and a reduced hyperbolic system of partial differential equations that contains only linearly degenerate fields. We implement FO-CCZ4 in a high-order path-conservative arbitrary-high-order-method-using-derivatives (ADER)-DG scheme with adaptive mesh refinement and local time-stepping, supplemented with a third-order ADER-WENO subcell finite-volume limiter in order to deal with singularities arising with black holes. We validate the correctness of the formulation through a series of standard tests in vacuum, performed in one, two and three spatial dimensions, and also present preliminary results on the evolution of binary black-hole systems. To the best of our knowledge, these are the first successful three-dimensional simulations of moving punctures carried out with high-order DG schemes using a first-order formulation of the Einstein equations.


title: Dimension towers of SICs. I. Aligned SICs and embedded tight frames

abstract: Algebraic number theory relates SIC-POVMs in dimension $d>3$ to those in dimension $d(d-2)$. We define a SIC in dimension $d(d-2)$ to be aligned to a SIC in dimension $d$ if and only if the squares of the overlap phases in dimension $d$ appear as a subset of the overlap phases in dimension $d(d-2)$ in a specified way. We give 19 (mostly numerical) examples of aligned SICs. We conjecture that given any SIC in dimension $d$ there exists an aligned SIC in dimension $d(d-2)$. In all our examples the aligned SIC has lower dimensional equiangular tight frames embedded in it. If $d$ is odd so that a natural tensor product structure exists, we prove that the individual vectors in the aligned SIC have a very special entanglement structure, and the existence of the embedded tight frames follows as a theorem. If $d-2$ is an odd prime number we prove that a complete set of mutually unbiased bases can be obtained by reducing an aligned SIC to this dimension.


title: JH$\im$-suboperator pairs with application to invariant approximation by using C-class functions

abstract: In this paper, by using C-class functions some results and common fixed point theorems are established for generalized JH-operator pairs of Sintu- navarat and Kumam (Journal of Inequalities and Applications, 67 (2011), 10 pages, doi:10.1186/1029-242X-2011-67) . Also, a new class of non-commuting self-mappings as JH$\im$-suboperator pairs are introduced. Final, as applications, several invariant approximation results are discussed


title: PyPSA: Python for Power System Analysis

abstract: Python for Power System Analysis (PyPSA) is a free software toolbox for simulating and optimising modern electrical power systems over multiple periods. PyPSA includes models for conventional generators with unit commitment, variable renewable generation, storage units, coupling to other energy sectors, and mixed alternating and direct current networks. It is designed to be easily extensible and to scale well with large networks and long time series. In this paper the basic functionality of PyPSA is described, including the formulation of the full power flow equations and the multi-period optimisation of operation and investment with linear power flow equations. PyPSA is positioned in the existing free software landscape as a bridge between traditional power flow analysis tools for steady-state analysis and full multi-period energy system models. The functionality is demonstrated on two open datasets of the transmission system in Germany (based on SciGRID) and Europe (based on GridKit).


title: On the ambiguity in the notion of transverse traceless modes of gravitational waves

abstract: Somewhat surprisingly, in many of the widely used monographs and review articles the term Transverse-Traceless modes of linearized gravitational waves is used to denote two entirely different notions. These treatments generally begin with a decomposition of the metric perturbation that is local in the momentum space (and hence non-local in physical space), and denote the resulting transverse traceless modes by $h_{ab}^{\TT}$. However, while discussing gravitational waves emitted by an isolated system --typically in a later section-- the relevant modes are extracted using a `projection operator' that is local in physical space. These modes are also called transverse-traceless and again labeled $h_{ab}^{\TT}$, implying that this is just a reformulation of the previous notion. But the two notions are conceptually distinct and the difference persists even in the asymptotic region. We show that this confusion arises already in Maxwell theory that is often discussed as a prelude to the gravitational case. Finally, we discuss why the distinction has nonetheless remained largely unnoticed, and also point out that there are some important physical effects where only one of the notions gives the correct answer.


title: A matrix Bougerol identity and the Hua-Pickrell measures

abstract: We prove a Hermitian matrix version of Bougerol's identity. Moreover, we construct the Hua-Pickrell measures on Hermitian matrices, as stochastic integrals with respect to a drifting Hermitian Brownian motion and with an integrand involving a conjugation by an independent, matrix analogue of the exponential of a complex Brownian motion with drift.


title: Robust Private Information Retrieval on Coded Data

abstract: We consider the problem of designing PIR scheme on coded data when certain nodes are unresponsive. We provide the construction of $\nu$-robust PIR schemes that can tolerate up to $\nu$ unresponsive nodes. These schemes are adaptive and universally optimal in the sense of achieving (asymptotically) optimal download cost for any number of unresponsive nodes up to $\nu$.


title: A breakthrough in Speech emotion recognition using Deep Retinal Convolution Neural Networks

abstract: Speech emotion recognition (SER) is to study the formation and change of speaker's emotional state from the speech signal perspective, so as to make the interaction between human and computer more intelligent. SER is a challenging task that has encountered the problem of less training data and low prediction accuracy. Here we propose a data augmentation algorithm based on the imaging principle of the retina and convex lens, to acquire the different sizes of spectrogram and increase the amount of training data by changing the distance between the spectrogram and the convex lens. Meanwhile, with the help of deep learning to get the high-level features, we propose the Deep Retinal Convolution Neural Networks (DRCNNs) for SER and achieve the average accuracy over 99%. The experimental results indicate that DRCNNs outperforms the previous studies in terms of both the number of emotions and the accuracy of recognition. Predictably, our results will dramatically improve human-computer interaction.


title: Bounce statistics for rational lattice paths

abstract: Given two relatively prime positive integers $\alpha$ and $\beta$, we consider simple lattice paths (with unit East and unit North steps) from $(0,0)$ to $(\alpha k,\beta k)$, and enumerate them by their left and right bounces with respect to the line $y=\frac{\beta}{\alpha} x$. We give the corresponding multivariate generating functions for all such paths as well as for subclasses of paths that start and end with a prescribed step. For illustration purposes, we discuss the case $\beta=1$ and express some of our functions in terms of the Fuss-Catalan generating function $c_\alpha(x)$.


title: Stability of ALE Ricci-flat manifolds under Ricci flow

abstract: We prove that if an ALE Ricci-flat manifold $(M,g)$ is linearly stable and integrable, it is dynamically stable under Ricci flow, i.e. any Ricci flow starting close to g exists for all time and converges modulo diffeomorphism to an ALE Ricci-flat metric close to $g$. By adapting Tian's approach in the closed case, we show that integrability holds for ALE Calabi-Yau manifolds which implies that they are dynamically stable.


title: Regularization techniques for fine-tuning in neural machine translation

abstract: We investigate techniques for supervised domain adaptation for neural machine translation where an existing model trained on a large out-of-domain dataset is adapted to a small in-domain dataset. In this scenario, overfitting is a major challenge. We investigate a number of techniques to reduce overfitting and improve transfer learning, including regularization techniques such as dropout and L2-regularization towards an out-of-domain prior. In addition, we introduce tuneout, a novel regularization technique inspired by dropout. We apply these techniques, alone and in combination, to neural machine translation, obtaining improvements on IWSLT datasets for English->German and English->Russian. We also investigate the amounts of in-domain training data needed for domain adaptation in NMT, and find a logarithmic relationship between the amount of training data and gain in BLEU score.


title: A Family of High-order Gas-kinetic Schemes and Its Comparison with Riemann Solver Based High-order Methods

abstract: Most high order computational fluid dynamics (CFD) methods for compressible flows are based on Riemann solver for the flux evaluation and Runge-Kutta (RK) time stepping technique for temporal accuracy. The main advantage of this kind of approach is the easy implementation and stability enhancement by introducing more middle stages. However, the nth-order time accuracy needs no less than n stages for the RK method, which is very time and memory consuming for a high order method. On the other hand, the multi-stage multi-derivative (MSMD) method can be used to achieve the same order of time accuracy using less middle stages, once the time derivatives of the flux function is used.The gas kinetic scheme (GKS) provides such a time accurate evolution model. By combining the second-order or third-order GKS flux functions with the MSMD technique, a family of high order gas kinetic methods can be constructed. As an extension of the previous 2-stage 4th-order GKS, the 5th-order schemes with 2 and 3 stages will be developed in this paper. Based on the same 5th-order WENO reconstruction, the performance of gas kinetic schemes from the 2nd- to the 5th-order time accurate methods will be evaluated. The results show that the 5th-order scheme can achieve the theoretical order of accuracy for the Euler equations, and present accurate Navier-Stokes solutions as well due to the coupling of inviscid and viscous terms in the GKS formulation. In comparison with Riemann solver based 5th-order RK method, the high order GKS has advantages in terms of efficiency, accuracy, and robustness, for all test cases. It provides a promising direction for the development of high-order CFD methods for the computation of complex flows, such as turbulence and acoustics with shock interactions.


title: Random integral operators related to the point processes

abstract: In the article we study properties of the random integral operator in $L_2(\mathbb{R})$ whose kernel is obtained as a convolution of Gaussian density with a stationary point process.


title: Color superfluidity of neutral ultra-cold fermions in the presence of color-flip and color-orbit fields

abstract: We describe how color superfluidity is modified in the presence of color-flip and color-orbit fields in the context of ultra-cold atoms, and discuss connections between this problem and that of color superconductivity in quantum chromodynamics. We consider s-wave contact interactions between different colors, and we identify superfluid phases, with five being nodal and one being fully gapped. When our system is described in a mixed color basis, the superfluid order parameter tensor is characterized by six independent components with explicit momentum dependence induced by color-orbit coupling. The nodal superfluid phases are topological in nature, and the low temperature phase diagram of color-flip field versus interaction parameter exhibits a pentacritical point, where all five nodal color superfluid phases converge. These results are in sharp contrast to the case of zero color-flip and color-orbit fields, where the system has perfect U(3) symmetry and possess a superfluid phase that is characterized by fully gapped quasiparticle excitations with a single complex order parameter with no momentum dependence and by inert unpaired fermions representing a non-superfluid component. Furthermore, we analyse the order parameter tensor in a total pseudo-spin basis, investigate its momentum dependence in the singlet, triplet and quintet sectors, and compare the results with the simpler case of spin-1/2 fermions in the presence of spin-flip and spin-orbit fields. Finally, we analyse in detail spectroscopic properties of color superfluids in the presence of color-flip and color-orbit fields, such as the quasiparticle excitation spectrum, momentum distribution, and density of states to help characterize all the encountered topological quantum phases, which can be realized in fermionic isotopes of Lithium, Potassium and Ytterbium atoms with three internal states trapped.


title: Surprise in simplicity: an unusual spectral evolution of a single pulse GRB 151006A

abstract: We present a detailed analysis of GRB 151006A, the first GRB detected by Astrosat CZT Imager (CZTI). We study the long term spectral evolution by exploiting the capabilities of \emph{Fermi} and \emph{Swift} satellites at different phases, which is complemented by the polarization measurement with the CZTI. While the light curve of the GRB in different energy bands show a simple pulse profile, the spectrum shows an unusual evolution. The first phase exhibits a hard-to-soft (HTS) evolution until $\sim16-20$\,s, followed by a sudden increase in the spectral peak reaching a few MeV. Such a dramatic change in the spectral evolution in case of a single pulse burst is reported for the first time. This is captured by all models we used namely, Band function, Blackbody+Band and two blackbodies+power law. Interestingly, the \emph{Fermi} Large Area Telescope (LAT) also detects its first photon ($>100$\,MeV) during this time. This new injection of energy may be associated with either the beginning of afterglow phase, or a second hard pulse of the prompt emission itself which, however, is not seen in the otherwise smooth pulse profile. By constructing Bayesian blocks and studying the hardness evolution we find a good evidence for a second hard pulse. The \emph{Swift} data at late epochs ($>T_{90}$ of the GRB) also shows a significant spectral evolution consistent with the early second phase. The CZTI data (100--350\,keV), though having low significance ($1\sigma$), show high values of polarization in the two epochs ($77\%$ to $94\%$), in agreement with our interpretation.


title: Quaternionic Arithmetic Lattices of Rank 2 and a Fake Quadric in Characteristic 2

abstract: We construct a torsion-free arithmetic lattice in $\mathrm{PGL}_2(\mathbb{F}_2(\!(t)\!))\times\mathrm{PGL}_2(\mathbb{F}_2(\!(t)\!))$ arising from a quaternion algebra over $\mathbb{F}_2(z)$. It is the fundamental group of a square complex with universal covering $T_3\times T_3$, a product of trees with constant valency $3$, which has minimal Euler characteristic. Furthermore, our lattice gives rise to a fake quadric over $\mathbb{F}_2(\!(t)\!)$ by means of non-archimedean uniformization.


title: A Framework for Super-Resolution of Scalable Video via Sparse Reconstruction of Residual Frames

abstract: This paper introduces a framework for super-resolution of scalable video based on compressive sensing and sparse representation of residual frames in reconnaissance and surveillance applications. We exploit efficient compressive sampling and sparse reconstruction algorithms to super-resolve the video sequence with respect to different compression rates. We use the sparsity of residual information in residual frames as the key point in devising our framework. Moreover, a controlling factor as the compressibility threshold to control the complexity-performance trade-off is defined. Numerical experiments confirm the efficiency of the proposed framework in terms of the compression rate as well as the quality of reconstructed video sequence in terms of PSNR measure. The framework leads to a more efficient compression rate and higher video quality compared to other state-of-the-art algorithms considering performance-complexity trade-offs.


title: 2004 TT357: A potential contact binary in the Trans-Neptunian belt

abstract: We report photometric observations of the trans-Neptunian object 2004~TT$_{357}$ obtained in 2015 and 2017 using the 4.3~m Lowell's Discovery Channel Telescope. We derive a rotational period of 7.79$\pm$0.01~h and a peak-to-peak lightcurve amplitude of 0.76$\pm$0.03~mag. 2004 TT$_{357}$ displays a large variability that can be explained by a very elongated single object or can be due to a contact/close binary. The most likely scenario is that 2004 TT$_{357}$ is a contact binary. If it is in hydrostatic equilibrium, we find that the lightcurve can be explained by a system with a mass ratio q$_{min}$=0.45$\pm$0.05, and a density of $\rho_{min}$=2g cm$^{-3}$, or less likely a system with q$_{max}$=0.8$\pm$0.05, and $\rho_{max}$=5g cm$^{-3}$. Considering a single triaxial ellipsoid in hydrostatic equilibrium, we derive a lower limit to the density of 0.78g cm$^{-3}$, and an elongation (a/b) of 2.01 assuming an equatorial view. From Hubble Space Telescope data, we report no resolved companion orbiting 2004 TT$_{357}$. Despite an expected high fraction of contact binaries in the trans-Neptunian belt, 2001 QG$_{298}$ is the unique confirmed contact binary in the trans-Neptunian belt, and 2004 TT$_{357}$ is only the second candidate to this class of systems, with 2003 SQ$_{317}$.


title: Quantifying Coherence and Entanglement via Simple Measurements

abstract: Coherence and entanglement are fundamental properties of quantum systems, promising to power the near future quantum computers, sensors and simulators. Yet, their experimental detection is challenging, usually requiring full reconstruction of the system state. We show that one can extract quantitative bounds to the relative entropy of coherence and the coherent information, coherence and entanglement quantifiers respectively, by a limited number of purity measurements. The scheme is readily implementable with current technology to verify quantum computations in large scale registers, without carrying out expensive state tomography.


title: Gravitational lensing by wave dark matter halos

abstract: Wave Dark Matter (WaveDM) has recently gained attention as a viable candidate to account for the dark matter content of the Universe. In this paper we explore the extent to which dark matter halos in this model, and under what conditions, are able to reproduce strong lensing systems. First, we analytically explore the lensing properties of the model -- finding that a pure WaveDM density profile, a soliton profile, produces a weaker lensing effect than other similar cored profiles. Then we analyze models with a soliton embedded in an NFW profile, as has been found in numerical simulations of structure formation. We use a benchmark model with a boson mass of $m_a=10^{-22}{\rm eV}$, for which we see that there is a bi-modality in the contribution of the external NFW part of the profile, and actually some of the free parameters associated with it are not well constrained. We find that for configurations with boson masses $10^{-23}$ -- $10^{-22}{\rm eV}$, a range of masses preferred by dwarf galaxy kinematics, the soliton profile alone can fit the data but its size is incompatible with the luminous extent of the lens galaxies. Likewise, boson masses of the order of $10^{-21}{\rm eV}$, which would be consistent with Lyman-$\alpha$ constraints and consist of more compact soliton configurations, necessarily require the NFW part in order to reproduce the observed Einstein radii. We then conclude that lens systems impose a conservative lower bound $m_a > 10^{-24}$ and that the NFW envelope around the soliton must be present to satisfy the observational requirements.


title: Debugging Transactions and Tracking their Provenance with Reenactment

abstract: Debugging transactions and understanding their execution are of immense importance for developing OLAP applications, to trace causes of errors in production systems, and to audit the operations of a database. However, debugging transactions is hard for several reasons: 1) after the execution of a transaction, its input is no longer available for debugging, 2) internal states of a transaction are typically not accessible, and 3) the execution of a transaction may be affected by concurrently running transactions. We present a debugger for transactions that enables non-invasive, post-mortem debugging of transactions with provenance tracking and supports what-if scenarios (changes to transaction code or data). Using reenactment, a declarative replay technique we have developed, a transaction is replayed over the state of the DB seen by its original execution including all its interactions with concurrently executed transactions from the history. Importantly, our approach uses the temporal database and audit logging capabilities available in many DBMS and does not require any modifications to the underlying database system nor transactional workload.


title: Production of sub-gigabar pressures by a hyper-velocity impact in the collider using laser-induced cavity pressure acceleration

abstract: Production of high dynamic pressure using a strong shock wave is a topic of high relevance for high energy density physics, inertial confinement fusion and materials science. Although the pressures in the multi-megabar range can be produced by the shocks generated with a large variety of methods, the higher pressures, in the sub-gigabar or gigabar range, are achievable only with nuclear explosions or laser-driven shocks. However, the laser-to-shock energy conversion efficiency in the laser-based methods currently applied is low and, as a result, multi-kilojoule multi-beam lasers have to be used to produce such extremely high pressures. In this paper, the generation of high-pressure shocks in the newly proposed collider in which the projectile impacting a solid target is driven by the laser-induced cavity pressure acceleration (LICPA) mechanism is investigated using two-dimensional hydrodynamic simulations. A special attention is paid to the dependence of shock parameters and the laser-to-shock energy conversion efficiency on the impacted target material and the laser driver energy. It has been found that both in case of low-density and high-density solid targets the shock pressures in the sub-gigabar range can be produced in the LICPA-based collider with the laser energy of only a few hundred of joules, and the laser-to-shock energy conversion efficiency can reach values of 10 to 20 percent, by an order of magnitude higher than the conversion efficiencies achieved with other laser-based methods used so far.


title: Effects of the selection function on metallicity trends in spectroscopic surveys of the Milky Way

abstract: We investigate here the effect of the selection function on the metallicity distribution function (MDF) and on the vertical metallicity gradient by studying similar lines of sight using four different spectroscopic surveys (APOGEE, LAMOST, RAVE, and Gaia-ESO), which have different targeting strategies and therefore different selection functions. We use common fields between the spectroscopic surveys of APOGEE, LAMOST, RAVE (ALR) and APOGEE, RAVE, Gaia-ESO (AGR) and use two stellar population synthesis models, GALAXIA and TRILEGAL, to create mock fields for each survey. We apply the selection function in the form of colour and magnitude cuts of the respective survey to the mock fields to replicate the observed source sample. We make a basic comparison between the models to check which best reproduces the observed sample distribution. We carry out a quantitative comparison between the synthetic MDF from the mock catalogues using both models to understand the effect of the selection function on the MDF and on the vertical metallicity gradient. Using both models, we find a negligible effect of the selection function on the MDF for APOGEE, LAMOST, and RAVE. We find a negligible selection function effect on the vertical metallicity gradients as well, though GALAXIA and TRILEGAL have steeper and shallower slopes, respectively, than the observed gradient. After applying correction terms on the metallicities of RAVE and LAMOST with respect to our reference APOGEE sample, our observed vertical metallicity gradients between the four surveys are consistent within 1-sigma. We also find consistent gradient for the combined sample of all surveys in ALR and AGR. We estimated a mean vertical metallicity gradient of -0.241+/-0.028 dex kpc-1. There is a significant scatter in the estimated gradients in the literature, but our estimates are within their ranges.


title: Learning Neural Network Classifiers with Low Model Complexity

abstract: Modern neural network architectures for large-scale learning tasks have substantially higher model complexities, which makes understanding, visualizing and training these architectures difficult. Recent contributions to deep learning techniques have focused on architectural modifications to improve parameter efficiency and performance. In this paper, we derive a continuous and differentiable error functional for a neural network that minimizes its empirical error as well as a measure of the model complexity. The latter measure is obtained by deriving a differentiable upper bound on the Vapnik-Chervonenkis (VC) dimension of the classifier layer of a class of deep networks. Using standard backpropagation, we realize a training rule that tries to minimize the error on training samples, while improving generalization by keeping the model complexity low. We demonstrate the effectiveness of our formulation (the Low Complexity Neural Network - LCNN) across several deep learning algorithms, and a variety of large benchmark datasets. We show that hidden layer neurons in the resultant networks learn features that are crisp, and in the case of image datasets, quantitatively sharper. Our proposed approach yields benefits across a wide range of architectures, in comparison to and in conjunction with methods such as Dropout and Batch Normalization, and our results strongly suggest that deep learning techniques can benefit from model complexity control methods such as the LCNN learning rule.


title: Two-species boson mixture on a ring: A group theoretic approach to the quantum dynamics of low-energy excitations

abstract: We investigate the weak excitations of a system made up of two condensates trapped in a Bose-Hubbard ring and coupled by an interspecies repulsive interaction. Our approach, based on the Bogoliubov approximation scheme, shows that one can reduce the problem Hamiltonian to the sum of sub-Hamiltonians $\hat{H}_k$, each one associated to momentum modes $\pm k$. Each $\hat{H}_k$ is then recognized to be an element of a dynamical algebra. This uncommon and remarkable property allows us to present a straightforward diagonalization scheme, to find constants of motion, to highlight the significant microscopic processes, and to compute their time evolution. The proposed solution scheme is applied to a simple but still very interesting closed circuit, the trimer. The dynamics of low-energy excitations, corresponding to weakly-populated vortices, is investigated considering different choices of the initial conditions, and the angular-momentum transfer between the two condensates is evidenced. Finally, the condition for which the spectral collapse and dynamical instability are observed is derived analytically.


title: SN 2014C: VLBI Images of a Supernova Interacting with a Circumstellar Shell

abstract: We report on VLBI measurements of supernova 2014C at several epochs between $t = 384$ and 1057 days after the explosion. SN 2014C was an unusual supernova that initially had Type Ib optical spectrum, but after $t = 130$ d it developed a Type IIn spectrum with prominent H$\alpha$ lines, suggesting the onset of strong circumstellar interaction. Our first VLBI observation was at $t = 384$ d, and we find that the outer radius of SN 2014C was $(6.40 \pm 0.26) \times 10^{16}$ cm (for a distance of 15.1 Mpc), implying an average expansion velocity of $19300 \pm 790$ \kms\ up to that time. At our last epoch, SN 2014C was moderately resolved and shows an approximately circular outline but with an enhancement of the brightness on the W side. The outer radius of the radio emission at $t = 1057$ d is $(14.9 \pm 0.6) \times 10^{16}$ cm. We find that the expansion between $t = 384$ and 1057 d is well described by a constant velocity expansion with $v = 13600 \pm 650$ \kms. SN 2014C had clearly been substantially decelerated by $t = 384$ d. Our measurements are compatible with a scenario where the expanding shock impacted upon a shell of dense circumstellar material during the first year, as suggested by the observations at other wavelengths, but had progressed through the dense shell by the time of the VLBI observations.


title: Towards an Homological Generalization of the Direct Summand Theorem

abstract: We present a more general (parametric-) homological characterization of the Direct Summand Theorem. Specifically, we state two new conjectures: the Socle-Parameter conjecture (SPC) in its weak and strong forms. We give a proof for the week form by showing that it is equivalent to the Direct Summand Conjecture (DSC), now known to be true after the work of Y. Andr\'{e}, based on Scholze's theory of perfectoids. Furthermore, we prove the SPC in its strong form for the case when the multiplicity of the parameters is smaller or equal than two. Finally, we present a new proof of the DSC in the equicharacteristic case, based on the techniques thus developed.


title: Superinjective Simplicial Maps of the Two-sided Curve Complexes on Nonorientable Surfaces

abstract: Let $N$ be a compact, connected, nonorientable surface of genus $g$ with $n$ boundary components with $g \geq 5$, $n \geq 0$. Let $\mathcal{T}(N)$ be the two-sided curve complex of $N$. If $\lambda :\mathcal{T}(N) \rightarrow \mathcal{T}(N)$ is a superinjective simplicial map, then there exists a homeomorphism $h : N \rightarrow N$ unique up to isotopy such that $H(\alpha) = \lambda(\alpha)$ for every vertex $\alpha$ in $\mathcal{T}(N)$ where $H=[h]$.


title: Deep Convolutional Framelet Denosing for Low-Dose CT via Wavelet Residual Network

abstract: Model based iterative reconstruction (MBIR) algorithms for low-dose X-ray CT are computationally expensive. To address this problem, we recently proposed a deep convolutional neural network (CNN) for low-dose X-ray CT and won the second place in 2016 AAPM Low-Dose CT Grand Challenge. However, some of the texture were not fully recovered. To address this problem, here we propose a novel framelet-based denoising algorithm using wavelet residual network which synergistically combines the expressive power of deep learning and the performance guarantee from the framelet-based denoising algorithms. The new algorithms were inspired by the recent interpretation of the deep convolutional neural network (CNN) as a cascaded convolution framelet signal representation. Extensive experimental results confirm that the proposed networks have significantly improved performance and preserves the detail texture of the original images.


title: An Analysis of the Twitter Discussion on the 2016 Austrian Presidential Elections

abstract: In this paper, we provide a systematic analysis of the Twitter discussion on the 2016 Austrian presidential elections. In particular, we extracted and analyzed a data-set consisting of 343645 Twitter messages related to the 2016 Austrian presidential elections. Our analysis combines methods from network science, sentiment analysis, as well as bot detection. Among other things, we found that: a) the winner of the election (Alexander Van der Bellen) was considerably more popular and influential on Twitter than his opponent, b) the Twitter followers of Van der Bellen substantially participated in the spread of misinformation about him, c) there was a clear polarization in terms of the sentiments spread by Twitter followers of the two presidential candidates, d) the in-degree and out-degree distributions of the underlying communication network are heavy-tailed, and e) compared to other recent events, such as the 2016 Brexit referendum or the 2016 US presidential elections, only a very small number of bots participated in the Twitter discussion on the 2016 Austrian presidential election.


title: A Bilinear Equalizer for Massive MIMO Systems

abstract: We present a novel approach for low-complexity equalizer design well-suited for cellular massive MIMO systems. Our design allows to exploit the channel structure in terms of covariance matrices to improve the performance in the face of pilot-contamination, while basically keeping the complexity of a matched filter. This is achieved by restricting the equalizer to functions which are bilinear in the received data signals and the observations from a training phase. The proposed design generalizes several previous approaches to equalizer design for massive MIMO. We show by asymptotic analysis that with the proposed design the achievable rate grows without bound for growing numbers of antennas even in the presence of pilot-contamination. We demonstrate with numerical results that the proposed design is competitive with more complex approaches in a practical cellular setup.


title: Formal Analysis of Continuous-time Systems using Fourier Transform

abstract: To study the dynamical behaviour of the engineering and physical systems, we often need to capture their continuous behaviour, which is modeled using differential equations, and perform the frequency-domain analysis of these systems. Traditionally, Fourier transform methods are used to perform this frequency domain analysis using paper-and-pencil based analytical techniques or computer simulations. However, both of these methods are error prone and thus are not suitable for analyzing systems used in safety-critical domains, like medicine and transportation. In order to provide an accurate alternative, we propose to use higher-order-logic theorem proving to conduct the frequency domain analysis of these systems. For this purpose, the paper presents a higher-order-logic formalization of Fourier transform using the HOL-Light theorem prover. In particular, we use the higher-order-logic based formalizations of differential, integral, transcendental and topological theories of multivariable calculus to formally define Fourier transform and reason about the correctness of its classical properties, such as existence, linearity, time shifting, frequency shifting, modulation, time scaling, time reversal and differentiation in time domain, and its relationships with Fourier Cosine, Fourier Sine and Laplace transforms. We use our proposed formalization for the formal verification of the frequency response of a generic n-order linear system, an audio equalizer and a MEMs accelerometer, using the HOL-Light theorem prover.


title: Haffner 16 Redux: Revisiting a Young Cluster in the Outer Galaxy

abstract: Images and spectra recorded with the Gemini Multi-Object Spectrograph on Gemini South are used to investigate the stellar content of the open cluster Haffner 16. The (i', g'-i') color-magnitude diagram (CMD) constructed from these data samples the cluster main sequence (MS) and 5 magnitudes of the pre-MS (PMS). The isochrones do not track the PMS on the CMD, in the sense that the PMS has a shallower slope on the CMD than predicted by the models. Still, a dip in star counts that is associated with the relaxation of PMS stars onto the MS is identified near i' = 17. The depth and brightness of this feature - as well as the morphology of the cluster MS on the CMD - are matched by models with a slightly sub-solar metallicity that have an age ~ 20 Myr and a distance modulus of 12.3 +/- 0.2. A light profile of Haffner 16 is constructed in the WISE W1 filter which suggests that the cluster is surrounded by a diffuse stellar halo. Spectra are presented of candidate cluster MS and PMS stars selected according to location on the CMD. The spectra show characteristics that are suggestive of a sub-solar metallicity. Halpha emission is common among objects on the PMS locus on the CMD near i' = 18. It is suggested that the location of the Haffner 16 PMS on the CMD is affected by large-scale cool spot activity, likely induced by rapid stellar rotation.


title: A "converse" stability condition is necessary for a compact higher order scheme on non-uniform meshes

abstract: The stability bounds and error estimates for a compact higher order Numerov-Crank-Nicolson scheme on non-uniform space meshes for the 1D time-dependent Schr\"odinger equation have been recently derived. This analysis has been done in $L^2$ and $H^1$ mesh norms and used the non-standard "converse" condition $h_\omega\leq c_0\tau$, where $h_\omega$ is the mean space step, $\tau$ is the time step and $c_0>0$. Now we prove that such condition is necessary for some families of non-uniform meshes and any space norm. Also numerical results show unacceptably wrong behavior of numerical solutions (their dramatic mass non-conservation) when this condition is violated.


title: Laser induced strong-field ionization gas jet tomography

abstract: We introduce a novel in-situ strong field ionization tomography approach for characterizing the spatial density distribution of gas jets. We show that for typical intensities in high harmonic generation experiments, the strong field ionization mechanism used in our approach provides an improvement in the resolution close to factor of 2 (resolving about 8 times smaller voxel volume), when compared to linear/single-photon imaging modalities. We find, that while the depth of scan in linear tomography is limited by resolution loss due to the divergence of the driving laser beam, in the proposed approach the depth of focus is localized due to the inherent physical nature of strong-field interaction and discuss implications of these findings. We explore key aspects of the proposed method and compare it with commonly used single- and multi-photon imaging mechanisms. The proposed method will be particularly useful for strong field and attosecond science experiments.

categories: physics.optics physics.app-ph physics.atm-clus physics.atom-ph physics.chem-ph


title: Effective Fluid Description of the Dark Universe

abstract: We propose an effective anisotropic fluid description for a generic infrared-modified theory of gravity. In our framework, the additional component of the acceleration, commonly attributed to dark matter, is explained as a radial pressure generated by the reaction of the dark energy fluid to the presence of baryonic matter. Using quite general assumptions, and a microscopic description of the fluid in terms of a Bose-Einstein condensate of gravitons, we find the static, spherically symmetric solution for the metric in terms of the Misner-Sharp mass function and the fluid pressure. At galactic scales, we correctly reproduce the leading MOND-like $\log(r)$ and subleading $(1/r)\,\log(r)$ terms in the weak-field expansion of the potential. Our description also predicts a tiny (of order $10^{-6}$ for a typical spiral galaxy) Machian modification of the Newtonian potential at galactic scales, which is controlled by the cosmological acceleration.


title: Multi-chain models of Conserved Lattice Gas

abstract: Conserved lattice gas (CLG) models in one dimension exhibit absorbing state phase transition (APT) with simple integer exponents $\beta=1=\nu=\eta$ whereas the same on a ladder belong to directed percolation (DP)universality. We conjecture that additional stochasticity in particle transfer is a relevant perturbation and its presence on a ladder force the APT to be in DP class. To substantiate this we introduce a class of restricted conserved lattice gas models on a multi-chain system ($M\times L$ square lattice with periodic boundary condition in both directions), where particles which have exactly one vacant neighbor are active and they move deterministically to the neighboring vacant site. We show that for odd number of chains , in the thermodynamic limit $L \to \infty,$ these models exhibit APT at $\rho_c= \frac{1}{2}(1+\frac1M)$ with $\beta =1.$ On the other hand, for even-chain systems transition occurs at $\rho_c=\frac12$ with $\beta=1,2$ for $M=2,4$ respectively, and $\beta= 3$ for $M\ge6.$ We illustrate this unusual critical behavior analytically using a transfer matrix method.


title: Influence of sample preparation on the transformation of low-density to high-density amorphous ice: An explanation based on the Potential Energy Landscape

abstract: Experiments and computer simulations of the transformations of amorphous ices display different behavior depending on sample preparation methods, and on the rates of change of temperature and pressure to which samples are subjected. In addition to these factors, simulation results also depend strongly on the chosen water model. Using computer simulations of the ST2 water model, we study how the sharpness of the compression-induced transition from low-density amorphous ice (LDA) to high-density amorphous ice (HDA) is influenced by the preparation of LDA. By studying LDA samples prepared using widely different procedures, we find that the sharpness of the LDA-to- HDA transformation is correlated with the depth of the initial LDA sample in the potential energy landscape (PEL), as characterized by the inherent structure energy. Our results show that the complex phenomenology of the amorphous ices reported in experiments and computer simulations can be understood and predicted in a unified way from a knowledge of the PEL of the system.


title: Gaussian Process-Based Model Predictive Control of Blood Glucose for Patients with Type 1 Diabetes Mellitus

abstract: The insulin sensitivity (IS) of the human body changes with a circadian rhythm. This adds to the time-varying feature of the glucose metabolism process and places challenges on the blood glucose (BG) control of patients with Type 1 Diabetes Mellitus. This paper presents a Model Predictive Controller that takes the periodic IS into account, in order to enhance BG control. The future effect of the IS is predicted using a machine learning technique, namely, a customized Gaussian Process (GP), based on historical training data. The training data for the GP is continuously updated during closed-loop control, which enables the control scheme to learn and adapt to intra-individual and inter-individual changes of the circadian IS rhythm. The necessary state information is provided by an Unscented Kalman Filter. The closed-loop performance of the proposed control scheme is evaluated for different scenarios (including fasting, announced meals and skipped meals) through in silico studies on simulation models of G\"ottingen Minipigs.


title: Switching dynamics in cholesteric liquid crystal emulsions

abstract: In this work we numerically study the switching dynamics of a 2D cholesteric emulsion droplet immersed in an isotropic fluid under an electric field, which is either uniform or rotating with constant speed. The overall dynamics depend strongly on the magnitude and on the direction (with respect to the cholesteric axis) of the applied field, on the anchoring of the director at the droplet surface and on the elasticity. If the surface anchoring is homeotropic and a uniform field is parallel to the cholesteric axis, the director undergoes deep elastic deformations and the droplet typically gets stuck into metastable states which are rich in topological defects. When the surface anchoring is tangential, the effects due to the electric field are overall less dramatic, as a small number of topological defects form at equilibrium. The application of the field perpendicular to the cholesteric axis usually has negligible effects on the defect dynamics. The presence of a rotating electric field of varying frequency fosters the rotation of the defects and of the droplet as well, typically at lower speed than that of the field, due to the inertia of the liquid crystal. If the surface anchoring is homeotropic a periodic motion is found. Our results represent a first step to understand the dynamical response of a cholesteric droplet under an electric field and its possible application in designing novel liquid crystal-based devices.


title: Linear Boltzmann dynamics in a strip with large reflective obstacles: stationary state and residence time

abstract: The presence of obstacles modify the way in which particles diffuse. In cells, for instance, it is observed that, due to the presence of macromolecules playing the role of obstacles, the mean square displacement ofbiomolecules scales as a power law with exponent smaller than one. On the other hand, different situations in grain and pedestrian dynamics in which the presence of an obstacle accelerate the dynamics are known. We focus on the time, called residence time, needed by particles to cross a strip assuming that the dynamics inside the strip follows the linear Boltzmann dynamics. We find that the residence time is not monotonic with the sizeand the location of the obstacles, since the obstacle can force those particles that eventually cross the strip to spend a smaller time in the strip itself. We focus on the case of a rectangular strip with two open sides and two reflective sides and we consider reflective obstaclea into the strip.