title: One-loop effective LFV $Zl_kl_m$ vertex from heavy neutrinos within the Mass Insertion Approximation

abstract: In this paper we study the effective lepton flavor violating vertex of an electroweak $Z$ gauge boson and two charged leptons with different flavor, $l_k$ and $l_m$, that is generated to one-loop in low scale seesaw models with right handed neutrinos whose masses are heavier than the electroweak scale. We first compute the form factor describing this vertex by using the mass insertion approximation, where the flavor non-diagonal entries of the neutrino Yukawa coupling matrix are the unique origin, to one-loop level, of lepton flavor changing processes with charged leptons in the external legs. Then, by considering the proper large right handed neutrino mass expansion of the form factor, we derive a formula for the $Z l_k l_m$ effective vertex which is very simple and useful for fast phenomenological estimates. In the last part of this work we focus on the phenomenological applications of this vertex for simple and accurate estimates of the $Z \to l_k {\bar l}_m$ decay rates. Concretely, this vertex will allow us to conclude easily on the maximum allowed decay rates by present data in the inverse seesaw model. The found rates are promising, at the reach of future lepton colliders.


title: Effective Field Theories of Post-Newtonian Gravity: A comprehensive review

abstract: [Abridged] This review article presents the progress made over the last decade, since the introduction of effective field theories (EFTs) into post-Newtonian (PN) gravity. These have been put forward in the context of gravitational waves (GWs) from the compact binary inspiral. The mature development of this interdisciplinary field has resulted in significant advances of wide interest to physics at several levels serving various purposes. The field has firmly demonstrated, that seemingly disparate physical domains, such as quantum field theory and classical gravity, are related, and that the EFT framework is a universal one, where it has been proven to supply a robust methodology to boost progress in the development of PN theory. The review is aimed at a broad audience, from general readers new to the field, to specialists and experts in related subjects. The review begins with an overview of the introduction of EFTs into classical gravity and their development. Then, the basic ideas, which form the conceptual foundation of EFTs, are provided, and the strategy of a multi-stage EFT framework, which is utilized for the PN binary inspiral problem, is outlined. The main body of the review is then dedicated to presenting in detail the study of each of the effective theories at each of the intermediate scales in the problem, up to the actual GW observables. The review is concluded with the multiple prospects of building on the progress in the field, and using further modern field theory insights and tools, to specifically address the study of GWs, as well as to broadly expand our fundamental understanding of gauge and gravity theories across the classical and quantum regimes.


title: Superconformal Models for Graphene and Boundary Central Charges

abstract: In the context of boundary conformal field theory, we investigate whether the boundary trace anomaly can depend on marginal directions in the presence of supersymmetry. Recently, it was found that a graphene-like non-supersymmetric conformal field theory with a four-dimensional bulk photon and a three-dimensional boundary electron has two boundary central charges that depend on an exactly marginal direction, namely the gauge coupling. In this work, we supersymmetrize this theory, paying special attention to the boundary terms required by supersymmetry. We study models with 4, 8, and 16 Poincar\'e supercharges in the bulk, half of which are broken by the boundary. In all cases, we find that at all orders in perturbation theory, the gauge coupling is not renormalized, providing strong evidence that these theories are boundary conformal field theories. Moreover, the boundary central charges depend on the coupling. One possible exception to this dependence on marginal directions is that the difference between the two charges is coupling independent at one-loop in the maximally supersymmetric case. In our analysis, a possible boundary Chern-Simons term is incorporated by a bulk $\theta$-term.


title: High-accuracy Geant4 simulation and semi-analytical modeling of nuclear resonance fluorescence

abstract: Nuclear resonance fluorescence (NRF) is a photonuclear interaction that enables highly isotope-specific measurements in both pure and applied physics scenarios. High-accuracy design and analysis of NRF measurements in complex geometries is aided by Monte Carlo simulations of photon physics and transport, motivating Jordan and Warren (2007) to develop the G4NRF codebase for NRF simulation in Geant4. In this work, we enhance the physics accuracy of the G4NRF code and perform improved benchmarking simulations. The NRF cross section calculation in G4NRF, previously a Gaussian approximation, has been replaced with a full numerical integration for improved accuracy in thick-target scenarios. A high-accuracy semi-analytical model of expected NRF count rates in a typical NRF measurement is then constructed and compared against G4NRF simulations for both simple homogeneous and more complex heterogeneous geometries. Agreement between rates predicted by the semi-analytical model and G4NRF simulation is found at a level of ${\sim}1\%$ in simple test cases and ${\sim}3\%$ in more realistic scenarios, improving upon the ${\sim}20\%$ level of the initial benchmarking study and establishing a highly-accurate NRF framework for Geant4.


title: Restructuring Batch Normalization to Accelerate CNN Training

abstract: Batch Normalization (BN) has become a core design block of modern Convolutional Neural Networks (CNNs). A typical modern CNN has a large number of BN layers in its lean and deep architecture. BN requires mean and variance calculations over each mini-batch during training. Therefore, the existing memory access reduction techniques, such as fusing multiple CONV layers, are not effective for accelerating BN due to their inability to optimize mini-batch related calculations during training. To address this increasingly important problem, we propose to restructure BN layers by first splitting a BN layer into two sub-layers (fission) and then combining the first sub-layer with its preceding CONV layer and the second sub-layer with the following activation and CONV layers (fusion). The proposed solution can significantly reduce main-memory accesses while training the latest CNN models, and the experiments on a chip multiprocessor show that the proposed BN restructuring can improve the performance of DenseNet-121 by 25.7%.


title: An efficient quantum circuits optimizing scheme compared with QISKit

abstract: Recently, the development of quantum chips has made great progress-- the number of qubits is increasing and the fidelity is getting higher. However, qubits of these chips are not always fully connected, which sets additional barriers for implementing quantum algorithms and programming quantum programs. In this paper, we introduce a general circuit optimizing scheme, which can efficiently adjust and optimize quantum circuits according to arbitrary given qubits' layout by adding additional quantum gates, exchanging qubits and merging single-qubit gates. Compared with the optimizing algorithm of IBM's QISKit, the quantum gates consumed by our scheme is 74.7%, and the execution time is only 12.9% on average.


title: A Convolutional Neural Network for Aspect Sentiment Classification

abstract: With the development of the Internet, natural language processing (NLP), in which sentiment analysis is an important task, became vital in information processing.Sentiment analysis includes aspect sentiment classification. Aspect sentiment can provide complete and in-depth results with increased attention on aspect-level. Different context words in a sentence influence the sentiment polarity of a sentence variably, and polarity varies based on the different aspects in a sentence. Take the sentence, 'I bought a new camera. The picture quality is amazing but the battery life is too short.'as an example. If the aspect is picture quality, then the expected sentiment polarity is 'positive', if the battery life aspect is considered, then the sentiment polarity should be 'negative'; therefore, aspect is important to consider when we explore aspect sentiment in the sentence. Recurrent neural network (RNN) is regarded as a good model to deal with natural language processing, and RNNs has get good performance on aspect sentiment classification including Target-Dependent LSTM (TD-LSTM) ,Target-Connection LSTM (TC-LSTM) (Tang, 2015a, b), AE-LSTM, AT-LSTM, AEAT-LSTM (Wang et al., 2016).There are also extensive literatures on sentiment classification utilizing convolutional neural network, but there is little literature on aspect sentiment classification using convolutional neural network. In our paper, we develop attention-based input layers in which aspect information is considered by input layer. We then incorporate attention-based input layers into convolutional neural network (CNN) to introduce context words information. In our experiment, incorporating aspect information into CNN improves the latter's aspect sentiment classification performance without using syntactic parser or external sentiment lexicons in a benchmark dataset from Twitter but get better performance compared with other models.


title: Transfer Learning for Clinical Time Series Analysis using Recurrent Neural Networks

abstract: Deep neural networks have shown promising results for various clinical prediction tasks such as diagnosis, mortality prediction, predicting duration of stay in hospital, etc. However, training deep networks -- such as those based on Recurrent Neural Networks (RNNs) -- requires large labeled data, high computational resources, and significant hyperparameter tuning effort. In this work, we investigate as to what extent can transfer learning address these issues when using deep RNNs to model multivariate clinical time series. We consider transferring the knowledge captured in an RNN trained on several source tasks simultaneously using a large labeled dataset to build the model for a target task with limited labeled data. An RNN pre-trained on several tasks provides generic features, which are then used to build simpler linear models for new target tasks without training task-specific RNNs. For evaluation, we train a deep RNN to identify several patient phenotypes on time series from MIMIC-III database, and then use the features extracted using that RNN to build classifiers for identifying previously unseen phenotypes, and also for a seemingly unrelated task of in-hospital mortality. We demonstrate that (i) models trained on features extracted using pre-trained RNN outperform or, in the worst case, perform as well as task-specific RNNs; (ii) the models using features from pre-trained models are more robust to the size of labeled data than task-specific RNNs; and (iii) features extracted using pre-trained RNN are generic enough and perform better than typical statistical hand-crafted features.


title: Mining Periodic Patterns with a MDL Criterion

abstract: The quantity of event logs available is increasing rapidly, be they produced by industrial processes, computing systems, or life tracking, for instance. It is thus important to design effective ways to uncover the information they contain. Because event logs often record repetitive phenomena, mining periodic patterns is especially relevant when considering such data. Indeed, capturing such regularities is instrumental in providing condensed representations of the event sequences. We present an approach for mining periodic patterns from event logs while relying on a Minimum Description Length (MDL) criterion to evaluate candidate patterns. Our goal is to extract a set of patterns that suitably characterises the periodic structure present in the data. We evaluate the interest of our approach on several real-world event log datasets.


title: Primordial Black Holes from the QCD axion

abstract: We propose a mechanism to generate Primordial Black Holes (PBHs) which is independent of cosmological inflation and occurs slightly below the QCD phase transition. Our setup relies on the collapse of long-lived string-domain wall networks and is naturally realized in QCD axion models with domain wall number $N_{DW}>1$ and Peccei-Quinn symmetry broken after inflation. In our framework, dark matter is mostly composed of axions in the meV mass range along with a small fraction, $\Omega_{\text{PBH}} \gtrsim 10^{-6} \Omega_{\text{CDM}} $ of heavy $M \sim 10^4-10^7 M_\odot$ PBHs. The latter could play a role in alleviating some of the shortcomings of the $\Lambda$CDM model on sub-galactic scales. The scenario has distinct signatures in ongoing axion searches as well as gravitational wave observatories.


title: Berryogenesis: self-induced Berry flux and spontaneous non-equilibrium magnetism

abstract: Spontaneous symmetry breaking is central to the description of interacting phases of matter. Here we reveal a new mechanism through which a driven interacting system subject to a time-reversal symmetric driving field can spontaneously magnetize. We show that the strong internal ac fields of a metal driven close to its plasmon resonance may enable Berryogenesis: the spontaneous generation of a self-induced Bloch band Berry flux. The self-induced Berry flux supports and is sustained by a circulating plasmonic motion, which may arise even for a linearly polarized driving field. This non-equilibrium phase transition occurs above a critical driving amplitude, and may be of either continuous or discontinuous type. Berryogenesis relies on feedback due to interband coherences induced by internal fields, and may readily occur in a wide variety of multiband systems. We anticipate that graphene devices, in particular, provide a natural platform to achieve Berryogenesis and plasmon-mediated spontaneous non-equilibrium magnetization in present-day devices.


title: BSM Hadronic Matrix Elements for $\epsilon'/\epsilon$ and $K\to\pi\pi$ Decays in the Dual QCD Approach

abstract: We calculate for the first time all four-quark hadronic matrix elements of local operators possibly contributing to $K\to\pi\pi$ decays and in particular to the ratio $\epsilon'/\epsilon$ beyond the Standard Model (BSM). To this end we use the Dual QCD (DQCD) approach. In addition to 7 new mirror operators obtained from the SM ones by flipping the chirality, we count 13 BSM four-quark operators of a given chirality linearly independent of each other and of the aforesaid 14 operators for which hadronic matrix elements are already known. We present results in two bases for all these operators, one termed DQCD basis useful for the calculation of the hadronic matrix elements in the DQCD approach and the other called SD basis suited to the short distance renormalization group evolution above the 1~GeV scale. We demonstrate that the pattern of long distance evolution (meson evolution) matches the one of short distance evolution (quark-gluon evolution), a property which to our knowledge cannot be presently achieved in any other analytical framework. The highlights of our paper are chirally enhanced matrix elements of tensor-tensor and scalar-scalar BSM operators. They could thereby explain the emerging $\epsilon'/\epsilon$ anomaly which is strongly indicated within DQCD with some support from lattice QCD. On the other hand we do not expect the BSM operators to be relevant for the $\Delta I=1/2$ rule.


title: Bulk Pumping in 2D Topological Phases

abstract: The notion of topological (Thouless) pumping in topological phases is traditionally associated with Laughlin's pump argument for the quantization of the Hall conductance in two-dimensional (2D) quantum Hall systems. It relies on magnetic flux variations that thread the system of interest without penetrating its bulk, in the spirit of Aharonov-Bohm effects. Here we explore a different paradigm for topological pumping induced, instead, by magnetic flux variations $\delta\chi$ inserted through the bulk of topological phases. We show that $\delta\chi$ generically controls the analog of a topological pump, accompanied by robust physical phenomena. We demonstrate this concept of bulk pumping in two paradigmatic types of 2D topological phases: integer and fractional quantum Hall systems and topological superconductors. We show, in particular, that bulk pumping provides a unifying connection between seemingly distinct physical effects such as density variations described by Streda's formula in quantum Hall phases, and fractional Josephson currents in topological superconductors. We discuss the generalization of bulk pumping to other types of topological phases.

categories: cond-mat.mes-hall cond-mat.quant-gas cond-mat.str-el cond-mat.supr-con quant-ph


title: Einstein's Universe: Cosmological structure formation in numerical relativity

abstract: We perform large-scale cosmological simulations that solve Einstein's equations directly via numerical relativity. Starting with initial conditions sampled from the cosmic microwave background, we track the emergence of a cosmic web without the need for a background cosmology. We measure the backreaction of large-scale structure on the evolution of averaged quantities in a matter-dominated universe. Although our results are preliminary, we find the global backreaction energy density is of order $10^{-8}$ compared to the energy density of matter in our simulations, and is thus unlikely to explain accelerating expansion under our assumptions. Sampling scales above the homogeneity scale of the Universe ($100-180\,h^{-1}$Mpc), in our chosen gauge, we find $2-3\%$ variations in local spatial curvature.


title: Mixing of metals during star cluster formation: statistics and implications for chemical tagging

abstract: Ongoing surveys are in the process of measuring the chemical abundances in large numbers of stars, with the ultimate goal of reconstructing the formation history of the Milky Way using abundances as tracers. However, interpretation of these data requires that we understand the relationship between stellar distributions in chemical and physical space, i.e., how similar in chemical abundance do we expect a pair of stars to be as a function of the distance between their formation sites. We investigate this question by simulating the gravitational collapse of a turbulent molecular cloud extracted from a galaxy-scale simulation, seeded with chemical inhomogeneities with different initial spatial scales. We follow the collapse from galactic scales down to resolutions scales of $\approx 10^{-3}$ pc, and find that, during this process, turbulence mixes the metal patterns, reducing the abundance scatter initially present in the gas by an amount that depends on the initial scale of inhomogeneity of each metal field. However, we find that regardless of the initial spatial structure of the metals at the onset of collapse, the final stellar abundances are highly correlated on distances below a few pc, and nearly uncorrelated on larger distances. Consequently, the star formation process defines a natural size scale of $\sim 1$ pc for chemically-homogenous star clusters, suggesting that any clusters identified as homogenous in chemical space must have formed within $\sim 1$ pc of one another. However, in order to distinguish different star clusters in chemical space, observations across multiple elements will be required, and the elements that are likely to be most efficient at separating distinct clusters in chemical space are those whose correlation length in the ISM is of order tens of pc, comparable to the sizes of individual molecular clouds.


title: Mode changing and giant pulses in the millisecond pulsar PSR B1957+20

abstract: Many radio pulsars have stable pulse profiles, but some exhibit mode changing where the profile switches between two or more quasi-stable modes of emission. So far, these effects had only been seen in relatively slow pulsars, but we show here that the pulse profile of PSR B1957+20, a millisecond pulsar, switches between two modes, with a typical time between mode changes of only $1.7$ s (or $\sim\!1000$ rotations), the shortest observed so far. The two modes differ in both intensity and polarization, with relatively large differences in the interpulse and much more modest ones in the main pulse. We find that the changes in the interpulse precede those in the main pulse by $\sim\!25$ ms, placing an empirical constraint on the timescale over which mode changes occurs. We also find that the properties of the giant pulses emitted by PSR B1957+20 are correlated with the mode of the regular emission: their rate and the rotational phase at which they are emitted both depend on mode. Furthermore, the energy distribution of the giant pulses emitted near the main pulse depends on mode as well. We discuss the ramifications for our understanding of the radio emission mechanisms as well as for pulsar timing experiments.


title: The trouble with Hubble: Local versus global expansion rates in inhomogeneous cosmological simulations with numerical relativity

abstract: In a fully inhomogeneous, anisotropic cosmological simulation performed by solving Einstein's equations with numerical relativity, we find a local measurement of the effective Hubble parameter differs by less than 1\% compared to the global value. This variance is consistent with predictions from Newtonian gravity. We analyse the averaged local expansion rate on scales comparable to Type 1a supernova surveys, and find that local variance cannot resolve the tension between the \citet{riess2018b} and \citet{planck2018a} measurements.


title: Translationally invariant universal classical Hamiltonians

abstract: Spin models are widely studied in the natural sciences, from investigating magnetic materials in condensed matter physics to studying neural networks. Previous work has demonstrated that there exist simple classical spin models that are universal: they can replicate -- in a precise and rigorous sense -- the complete physics of any other classical spin model, to any desired accuracy. However, all previously known universal models break translational invariance. In this paper we show that there exist translationally invariant universal models. Our main result is an explicit construction of a translationally invariant, 2D, nearest-neighbour, universal classical Hamiltonian with a single free parameter. The proof draws on techniques from theoretical computer science, in particular recent complexity theoretic results on tiling problems. Our results imply that there exists a single Hamiltonian which encompasses all classical spin physics, just by tuning a single parameter and varying the size of the lattice. We also prove that our construction is optimal in terms of the number of parameters in the Hamiltonian; there cannot exist a translationally invariant universal Hamiltonian with only the lattice size as a parameter.


title: Formation and destruction of SiS in space

abstract: The presence of SiS in space seems to be restricted to a few selected types of astronomical environments. It is long known to be present in circumstellar envelopes around evolved stars and it has also been detected in a handful of star-forming regions with evidence of outflows, like Sgr\,B2, Orion\,KL and more recently L1157-B1. The kinetics of reactions involving SiS is very poorly known and here we revisit the chemistry of SiS in space by studying some potentially important reactions of formation and destruction of this molecule. We calculated {\it ab initio} potential energy surfaces of the SiOS system and computed rate coefficients in the temperature range 50-2500 K for the reaction of destruction of SiS, in collisions with atomic O, and of its formation, through the reaction between Si and SO. We find that both reactions are rapid, with rate coefficients of a few times 10$^{-10}$ cm$^3$ s$^{-1}$, almost independent of temperature. In the reaction between Si and SO, SiO production is 5-7 times more efficient than SiS formation. The reaction of SiS with O atoms can play an important role in destroying SiS in envelopes around evolved stars. We built a simple chemical model of a postshock gas to study the chemistry of SiS in protostellar outflows and we found that SiS forms with a lower abundance and later than SiO, that SiS is efficiently destroyed through reaction with O, and that the main SiS-forming reactions are Si + SO and Si + SO$_2$.


title: Optical spectroscopy of local type-1 AGN LINERs

abstract: The Balmer emission originated in the broad line region (BLR) of active galactic nuclei (AGNs) could be either weak and difficult to detect, or even absent, for low luminosity AGNs, as LINERs. Our goals in this paper are threefold. First, we want to explore the AGN-nature of nearby type-1 LINERs. Second, we aim at deriving a reliable interpretation for the different components of emission lines by studying their kinematics and ionization mechanism. Third, we intend to probe the neutral gas in the nuclei of these LINERs. We study the 22 local (z<0.025) type-1 LINERs from the Palomar Survey, on the basis of optical ground- and space-based long-slit spectroscopic observations taken with TWIN/CAHA and ALFOSC/NOT. Kinematics and fluxes of a set of emission lines, from H{\beta} to [S II], and the NaD doublet in absorption have been modelled and measured, after the subtraction of the underlying starlight. We also use ancillary spectroscopic data from HST/STIS. We found that the broad H{\alpha} component is sometimes elusive in our ground-based spectroscopy whereas it is ubiquitous for space-based data. By combining optical diagnostic diagrams, theoretical models (for AGNs, pAGB-stars and shocks) and the weak/strong-[O I] classification, we exclude the pAGBs-stars scenario in favor of the AGN as the dominant mechanism of ionisation in these LINERs, being shocks however relevant. The kinematical properties of the emission lines may indicate the presence of ionized outflows, preferentially seen in [O I]. However, the neutral gas outflows, diagnosed by NaD, would appear to be less frequent.


title: Fluctuations in non-conformal holographic plasma at criticality

abstract: We continue the study initiated in [arXiv:1708.02252] of the fluctuations of a strongly-coupled non-conformal plasma described holographically by Einstein gravity coupled to a dilaton with an exponential potential. The plasma approaches a critical point of a continuous phase transition in a specific limit, where the metric becomes a linear-dilaton background. This results to an analytic description of the quasi-normal mode spectrum, that can be extended perturbatively in the deviation away from the critical point. In the previous paper we showed that at criticality the quasinormal frequencies coalesce into a branch cut on the real axis. In this paper we give a more extended and complete discussion of these results. We compare in detail the numerical and analytical approximations in order to confirm their validity; we study (numerically and in a WKB approximation) the momentum dependence of the modes, in order to determine the cross-over scale that limits the validity of the hydrodynamic approximation, and which becomes arbitrarily low at the critical point; and we discuss in detail the procedure we use to complete the theory in the UV by gluing a slice of AdS geometry, and the extent to which it should provide a good approximation to a smooth UV-complete situation.


title: Bayesian Source Discrimination in Radio Interferometry

abstract: Methods currently in use for locating and characterising sources in radio interferometry maps are designed for processing images, and require interferometric maps to be preprocessed so as to resemble conventional images. We demonstrate a Bayesian code - BaSC - that takes into account the interferometric visibility data despite working with more computationally manageable image domain data products. This method is better able to discriminate nearby sources than the commonly used SExtractor, and has potential even in more complicated cases. We also demonstrate the correctness of the Bayesian resolving formula for simulated data, and its implications for source discrimination at distances below the full width half maximum of the restoring beam.


title: Cosmological Time Crystals From Einstein-Cubic Gravities

abstract: By including appropriate Riemman cubic invariants, we find that the dynamics of classical time crystals can be straightforwardly realized in Einstein gravity on the FLRW metric. The time reflection symmetry is spontaneously broken in the two vacua with the same scale factor $a$, but opposite $\dot a$. The tunneling from one vacuum to the other provides a robust mechanism for bounce universes; it always occurs for systems with positive energy density. For suitable matter energy-momentum tensor we also construct cyclic universes. Cosmological solutions that resemble the classical time crystals can be constructed in massive gravity.


title: Noise induced Non-Markovianity

abstract: We have studied the non-Markovianity of dichotomously driven spin-boson model in the strong coupling regime in both memory kernel and time convolutionless master equation formulations. A strong correlation between the decay time of the environmental correlation function and the non-zero non-Markovianity is found in the absence of the external noise. Stochastic driving is shown to create strong non-Markovianity when the dynamics of the system without driving is Markovian. Also, exact analytical expressions for the trace distance distinguishability and the non-Markovianity were obtained for the certain range of the parameters that describe the system and its environment.


title: Development in Astronomy and Space Science in Africa

abstract: The development of astronomy and space science in Africa has grown significantly over the past few years. These advancements make the United Nations Sustainable Development Goals more achievable, and open up the possibility of new beneficial collaborations.


title: Disentangling the AGN and Star-Formation connection using XMM-Newton

abstract: There is growing evidence supporting the coeval growth of galaxies and their resident SMBH. Most studies also claim a correlation between the activity of the SMBH and the star-formation of the host galaxy. It is unclear, though, whether this correlation extends to all redz and LX. In this work, we use data from the X-ATLAS and XMM-XXL North fields and compile the largest X-ray sample up to date, to investigate how X-ray AGN affect the star-formation of their host galaxies in a wide redz and luminosity baseline, i.e., 0.03<z<3 and log LX (2-10 keV)= (41-45.5) ergs-1. Our sample consists of 3336 AGN. 1872 X-ray objects have spectroscopic redz. For the remaining sources we calculate photometric redz using TPZ, a machine-learning algorithm. We estimate stellar masses (M*) and Star Formation Rates (SFRs) by applying SED fitting through the CIGALE code, using optical, near-IR and mid-IR photometry. 608 of our sources also have far-IR photometry (Herschel). We use these sources to calibrate the SFR calculations of the rest of our X-ray sample. Our results show a correlation between the LX and the SFR, at all redz and LX spanned by our sample. We also find a dependence of the specific SFR (sSFR) on redz, while there are indications that LX enhances the sSFR even at low redz. We then disentangle the effects of M* and redz on the SFR and study again its dependence on the LX. Towards this end, we use the Schreiber et al. formula to estimate the SFR of main sequence galaxies that have the same stellar mass and redshift as our X-ray AGN. Our analysis reveals that the AGN enhances the star-formation of its host galaxy, when the galaxy lies below the main sequence and quenches the star-formation of the galaxy it lives in, when the host lies above the main sequence. Therefore, the effect of AGN on the SFR of the host galaxy, depends on the location of the galaxy relative to the main sequence.


title: Shortcuts to adiabaticity in Fermi gases

abstract: Shortcuts to adiabaticity (STA) provide an alternative to adiabatic protocols to guide the dynamics of the system of interest without the requirement of slow driving. We report the controlled speedup via STA of the nonadiabatic dynamics of a Fermi gas, both in the non-interacting and strongly coupled, unitary regimes. Friction-free superadiabatic expansion strokes, with no residual excitations in the final state, are demonstrated in the unitary regime by engineering the modulation of the frequencies and aspect ratio of the harmonic trap. STA are also analyzed and implemented in the high-temperature regime, where the shear viscosity plays a pivotal role and the Fermi gas is described by viscous hydrodynamics.


title: A systematic approach to generalisations of General Relativity and their cosmological implications

abstract: A century ago, Einstein formulated his elegant and elaborate theory of General Relativity, which has so far withstood a multitude of empirical tests with remarkable success. Notwithstanding the triumphs of Einstein's theory, the tenacious challenges of modern cosmology and of particle physics have motivated the exploration of further generalised theories of spacetime. Even though Einstein's interpretation of gravity in terms of the curvature of spacetime is commonly adopted, the assignment of geometrical concepts to gravity is ambiguous because General Relativity allows three entirely different, but equivalent approaches of which Einstein's interpretation is only one. From a field-theoretical perspective, however, the construction of a consistent theory for a Lorentz-invariant massless spin-2 particle uniquely leads to General Relativity. Keeping Lorentz invariance then implies that any modification of General Relativity will inevitably introduce additional propagating degrees of freedom into the gravity sector. Adopting this perspective, we will review the recent progress in constructing consistent field theories of gravity based on additional scalar, vector and tensor fields. Within this conceptual framework, we will discuss theories with Galileons, with Lagrange densities as constructed by Horndeski and beyond, extended to DHOST interactions, or containing generalized Proca fields and extensions thereof, or several Proca fields, as well as bigravity theories and scalar-vector-tensor theories. We will review the motivation of their inception, different formulations, and essential results obtained within these classes of theories together with their empirical viability.


title: LaneNet: Real-Time Lane Detection Networks for Autonomous Driving

abstract: Lane detection is to detect lanes on the road and provide the accurate location and shape of each lane. It severs as one of the key techniques to enable modern assisted and autonomous driving systems. However, several unique properties of lanes challenge the detection methods. The lack of distinctive features makes lane detection algorithms tend to be confused by other objects with similar local appearance. Moreover, the inconsistent number of lanes on a road as well as diverse lane line patterns, e.g. solid, broken, single, double, merging, and splitting lines further hamper the performance. In this paper, we propose a deep neural network based method, named LaneNet, to break down the lane detection into two stages: lane edge proposal and lane line localization. Stage one uses a lane edge proposal network for pixel-wise lane edge classification, and the lane line localization network in stage two then detects lane lines based on lane edge proposals. Please note that the goal of our LaneNet is built to detect lane line only, which introduces more difficulties on suppressing the false detections on the similar lane marks on the road like arrows and characters. Despite all the difficulties, our lane detection is shown to be robust to both highway and urban road scenarios method without relying on any assumptions on the lane number or the lane line patterns. The high running speed and low computational cost endow our LaneNet the capability of being deployed on vehicle-based systems. Experiments validate that our LaneNet consistently delivers outstanding performances on real world traffic scenarios.


title: Casimir forces and quantum friction of finite-size atoms in relativistic trajectories

abstract: We study quantum friction and Casimir forces with a full-relativistic formalism for atoms modelled as Unruh-DeWitt detectors in the presence of arbitrary macroscopic objects. We consider the general case of atoms with arbitrary relativistic trajectories in arbitrary quantum states (including coherent superpositions) close to objects that impose arbitrary boundary conditions. Particularizing for conducting plates, we show that, for relative velocities close to the speed of light, the quantum friction diverges while the Casimir force is almost independent of the velocity. Since we include the effect of the finite size of the detector and the finite interaction time, we also obtain quantum friction when the detector is isolated but follows a non-inertial trajectory.

categories: quant-ph cond-mat.mes-hall cond-mat.mtrl-sci gr-qc hep-th


title: Nonlinear Chiral Transport in Dirac Semimetals

abstract: We study the current of chiral charge density in a Dirac semimetal with two Dirac points in momentum space, subjected to an externally applied time dependent electric field and in the presence of a magnetic field. Based on the kinetic equation approach, we find contributions to the chiral charge current, that are proportional to the second power of the electric field and to the first and second powers of the magnetic field, describing the interplay of the chiral anomaly and the drift motion of electrons moving under the action of electric and magnetic fields.


title: On the effect of high order empirical parameters on the nuclear equation of state

abstract: A quantitative knowledge of the nuclear Equation of State (EoS) requires an accurate estimation of the uncertainties on the EoS parameters and their mutual correlations. Such correlations are empirically observed in a large set of EoS models by different authors, but they are not always fully understood. We show that some of these correlations can be interpreted from basic physical constraints imposed on a simple Taylor expansion of the binding energy around saturation density. In particular, we investigate the correlations among the following empirical parameters: the symmetry energy $E_{sym}$, the slope and curvature of the symmetry energy $L_{sym}$ and $K_{sym}$, and the curvature and skewness of the binding energy in symmetric matter $K_{sat}$ and $Q_{sat}$. The uncertainties on these correlations is estimated through analytical modelling as well as a meta-modelling analysis of the EoS subject to physical constraints. We show that a huge dispersion of the correlations among low order empirical parameters is induced by the unknown higher order empirical parameters, such as $K_{sym}$ (second order) and $Q_{sat}$ and $Q_{sym}$ (third order). We also propose an explanation of the reason why $Q_{sat}$ is weakly constrained by present experimental data. We conclude that selective observables on high order parameters, such as $K_{sym}$ and $Q_{sat}$, should be better determined before the present uncertainties on the EoS can be further reduced.


title: Dark Matter, Millicharges, Axion and Scalar Particles, Gauge Bosons, and Other New Physics with LDMX

abstract: The proposed LDMX experiment would provide roughly a meter-long region of instrumented tracking and calorimetry that acts as a beam stop for multi-GeV electrons in which each electron is tagged and its evolution measured. This would offer an unprecedented opportunity to access both collider-invisible and ultra-short lifetime decays of new particles produced in electron (or muon)-nuclear fixed-target collisions. In this paper, we show that the missing momentum channel and displaced decay signals in such an experiment could provide world-leading sensitivity to sub-GeV dark matter, millicharged particles, and visibly or invisibly decaying axions, scalars, dark photons, and a range of other new physics scenarios.


title: Noncommutative dispersion relation and mass-radius relation of white dwarfs

abstract: The equation of state of the electron degenerate gas in a white dwarf is usually treated by employing the ideal dispersion relation. However, the effect of quantum gravity is expected to be inevitably present and when this effect is considered through a non-commutative formulation, the dispersion relation undergoes a substantial modification. In this paper, we take such a modified dispersion relation and find the corresponding equation of state for the degenerate electron gas in white dwarfs. Hence we solve the equation of hydrostatic equilibrium and find that this leads to the possibility of the existence of excessively high values of masses exceeding the Chandrasekhar limit although the quantum gravity effect is taken to be very small. It is only when we impose the additional effect of neutronization that we obtain white dwarfs with masses close to the Chandrasekhar limit with nonzero radii at the neutronization threshold. We demonstrate these results by giving the numerical estimates for the masses and radii of $^4He$ and $^{12}C$, and $^{16}O$ white dwarfs.


title: Recommendation Systems and Self Motivated Users

abstract: Modern recommendation systems rely on the wisdom of the crowd to learn the optimal course of action. This induces an inherent mis-alignment of incentives between the system's objective to learn (explore) and the individual users' objective to take the contemporaneous optimal action (exploit). The design of such systems must account for this and also for additional information available to the users. A prominent, yet simple, example is when agents arrive sequentially and each agent observes the action and reward of his predecessor. We provide an incentive compatible and asymptotically optimal mechanism for that setting. The complexity of the mechanism suggests that the design of such systems for general settings is a challenging task.


title: Detecting parity effect in a superconducting device in the presence of parity switches

abstract: We present a superconducting device showing a clear parity effect in the number of electrons, even when there is, on average, a single nonequilibrium quasiparticle present and the parity of the island switches due to quasiparticles tunneling in and out of the device at rates on the order of 100 Hz. We detect the switching by monitoring in real time the charge state of a superconducting island connected to normal leads by tunnel junctions. The quasiparticles are created by Cooper pairs breaking on the island at a rate of a few kHz. We demonstrate that the pair breaking is caused by the backaction of the single-electron transistor used as a charge detector. With sufficiently low probing currents, our superconducting island is free of quasiparticles 97% of the time.


title: Taelman L-Values for Drinfeld Modules over Tate Algebras

abstract: In the present paper, we analyze Taelman L-values corresponding to Drinfeld modules over Tate algebras of arbitrary rank. Using our results, we also introduce an L-series converging in Tate algebras which can be seen as a generalization of Pellarin L-series.


title: Asymmetric evolutionary games with environmental feedback

abstract: Models in evolutionary game theory traditionally assume symmetric interactions in homogeneous environments. Here, we consider populations evolving in a heterogeneous environment, which consists of patches of different qualities that are occupied by one individual each. The fitness of individuals is not only determined by interactions with others but also by environmental quality. This heterogeneity results in asymmetric interactions where the characteristics of the interaction may depend on an individual's location. Interestingly, in non-varying heterogeneous environments, the long-term dynamics are the same as for symmetric interactions in an average, homogeneous environment. However, introducing environmental feedback between an individual's strategy and the quality of its patch results in rich eco-evolutionary dynamics. Thus, individuals act as ecosystem engineers. The nature of the feedback and the rate of ecological changes can relax or aggravate social dilemmas and promote persistent periodic oscillations of strategy abundance and environmental quality.


title: Transfer with Model Features in Reinforcement Learning

abstract: A key question in Reinforcement Learning is which representation an agent can learn to efficiently reuse knowledge between different tasks. Recently the Successor Representation was shown to have empirical benefits for transferring knowledge between tasks with shared transition dynamics. This paper presents Model Features: a feature representation that clusters behaviourally equivalent states and that is equivalent to a Model-Reduction. Further, we present a Successor Feature model which shows that learning Successor Features is equivalent to learning a Model-Reduction. A novel optimization objective is developed and we provide bounds showing that minimizing this objective results in an increasingly improved approximation of a Model-Reduction. Further, we provide transfer experiments on randomly generated MDPs which vary in their transition and reward functions but approximately preserve behavioural equivalence between states. These results demonstrate that Model Features are suitable for transfer between tasks with varying transition and reward functions.


title: AstroSat-CZTI detection of variable prompt emission polarization in GRB 171010A

abstract: We present spectro-polarimetric analysis of \thisgrb\ using data from \asat, \fermi, and \swift, to provide insights into the physical mechanisms of the prompt radiation and the jet geometry. Prompt emission from \thisgrb\ was very bright (fluence $>10^{-4}$~ergs~cm$^{-2}$) and had a complex structure composed of the superimposition of several pulses. The energy spectra deviate from the typical Band function to show a low energy peak $\sim 15$~keV --- which we interpret as a power-law with two breaks, with a synchrotron origin. Alternately, the prompt spectra can also be interpreted as Comptonized emission, or a blackbody combined with a Band function. Time-resolved analysis confirms the presence of the low energy component, while the peak energy is found to be confined in the range of 100--200~keV. Afterglow emission detected by \fermi-LAT is typical of an external shock model, and we constrain the initial Lorentz factor using the peak time of the emission. \swift-XRT measurements of the afterglow show an indication for a jet break, allowing us to constrain the jet opening angle to $>$ 6$\degr$. Detection of a large number of Compton scattered events by \asat-CZTI provides an opportunity to study hard X-ray polarization of the prompt emission. We find that the burst has high, time-variable polarization, with the emission {\bf have higher polarization} at energies above the peak energy. We discuss all observations in the context of GRB models and polarization arising due to {\bf due to physical or geometric effects:} synchrotron emission from multiple shocks with ordered or random magnetic fields, Poynting flux dominated jet undergoing abrupt magnetic dissipation, sub-photospheric dissipation, a jet consisting of fragmented fireballs, and the Comptonization model.


title: Investigating the role of L1 in automatic pronunciation evaluation of L2 speech

abstract: Automatic pronunciation evaluation plays an important role in pronunciation training and second language education. This field draws heavily on concepts from automatic speech recognition (ASR) to quantify how close the pronunciation of non-native speech is to native-like pronunciation. However, it is known that the formation of accent is related to pronunciation patterns of both the target language (L2) and the speaker's first language (L1). In this paper, we propose to use two native speech acoustic models, one trained on L2 speech and the other trained on L1 speech. We develop two sets of measurements that can be extracted from two acoustic models given accented speech. A new utterance-level feature extraction scheme is used to convert these measurements into a fixed-dimension vector which is used as an input to a statistical model to predict the accentedness of a speaker. On a data set consisting of speakers from 4 different L1 backgrounds, we show that the proposed system yields improved correlation with human evaluators compared to systems only using the L2 acoustic model.


title: Proximal algorithms for large-scale statistical modeling and sensor/actuator selection

abstract: Several problems in modeling and control of stochastically-driven dynamical systems can be cast as regularized semi-definite programs. We examine two such representative problems and show that they can be formulated in a similar manner. The first, in statistical modeling, seeks to reconcile observed statistics by suitably and minimally perturbing prior dynamics. The second seeks to optimally select a subset of available sensors and actuators for control purposes. To address modeling and control of large-scale systems we develop a unified algorithmic framework using proximal methods. Our customized algorithms exploit problem structure and allow handling statistical modeling, as well as sensor and actuator selection, for substantially larger scales than what is amenable to current general-purpose solvers. We establish linear convergence of the proximal gradient algorithm, draw contrast between the proposed proximal algorithms and alternating direction method of multipliers, and provide examples that illustrate the merits and effectiveness of our framework.


title: The radius of analyticity for solutions to a problem in epitaxial growth on the torus

abstract: A certain model for epitaxial film growth has recently attracted attention, with the existence of small global solutions having being proved in both the case of the n-dimensional torus and free space. We address a regularity question for these solutions, showing that in the case of the torus, the solutions become analytic at any positive time, with the radius of analyticity growing linearly for all time. As other authors have, we take the Laplacian of the initial data to be in the Wiener algebra, and we find an explicit smallness condition on the size of the data. Our particular condition on the torus is that the Laplacian of the initial data should have norm less than 1/4 in the Wiener algebra.


title: Sparse Compression of Expected Solution Operators

abstract: We show that the expected solution operator of prototypical linear elliptic partial differential operators with random coefficients is well approximated by a computable sparse matrix. This result is based on a random localized orthogonal multiresolution decomposition of the solution space that allows both the sparse approximate inversion of the random operator represented in this basis as well as its stochastic averaging. The approximate expected solution operator can be interpreted in terms of classical Haar wavelets. When combined with a suitable sampling approach for the expectation, this construction leads to an efficient method for computing a sparse representation of the expected solution operator.


title: Electric field modulated topological magnetoelectric effect in Bi$_2$Se$_3$

abstract: Topological insulators have been predicted to exhibit a variety of interesting phenomena including a quantized magnetoelectric response and novel spintronics effects due to spin textures on their surfaces. However, experimental observation of these phenomena has proved difficult due to the finite bulk carrier density which may overwhelm the intrinsic topological responses that are expressed at the surface. Here, we demonstrate a novel ionic gel gating technique to tune the chemical potential of Bi$_{2}$Se$_{3}$ thin films while simultaneously performing THz spectroscopy. We can tune the carrier concentration by an order of magnitude and shift the Fermi energy, E$_{F} $ to as low as $\simeq$ 10 meV above the Dirac point. At high bias voltage and magnetic field, we observe a quantized Faraday angle consistent with the topological magnetoelectric effect that can be tuned by ionic gel gating through a number of plateau states


title: Anomalies Of Phase Diagrams And Physical Properties Of Antiferrodistortive Perovskide Oxides

abstract: The influence of rotomagnetic (RM), rotoelectric (RE) and magnetoelectric (ME) coupling on phase diagram and properties of antiferrodistortive (AFD) perovskite oxides was reviewed. The main examples we consider in the review are typical AFD perovkites, such as incipient ferroelectrics EuTiO3, SrTiO3, EuxSr1-xTiO3, multiferroic BiFeO3 and Bi1-xRxFeO3 (x = La, Nd). The strong influence of RM, RE and ME couplings on the physical properties and phase diagrams including antiferromagnetic (AFM), ferroelectric (FE) and structural AFD phases has been revealed in the framework of Landau-Ginzburg-Devonshire (LGD) theory, as well as the prediction of novel (double and triple) multiferroic phases has been demonstrated. In the review we are especially focused on (a) the possibility to induce FM (FE) phase in EuTiO3 (as well as in other paraelectric AFM oxides) by the application of an electric (magnetic) field due to the ME coupling; (b) the analysis of the size effects and novel phases in EuxSr1-xTiO3 nanosystems, where the LGD predicts the presence of the triple AFD-FE-FM(AFM) phase at low temperatures; (c) the appearance of improper spontaneous polarization and pyroelectricity in the vicinity of antiphase domain boundaries, structural twin walls, surfaces and interphases in the AFD phase of non-ferroelectric SrTiO3 induced by the flexoelectricity and rotostriction; (d) the occurrence of low symmetry monoclinic phase with in-plane FE polarization in thin strained EuxSr1-xTiO3 films and its stabilization over wide temperature range by AFD oxygen octahedron tilts due to flexoelectric and rotostriction coupling; (e) discussion of a surprisingly strong size-induced increase of AFM transition temperature caused by the joint action of RM coupling with elastic stress accumulated in the intergrain spaces of BiFeO3 dense ceramics.


title: Formulas for Chebotarev densities of Galois extensions of number fields

abstract: We generalize the Chebotarev density formulas of Dawsey (2017) and Alladi (1977) to the setting of arbitrary finite Galois extensions of number fields $L/K$. In particular, if $C \subset G = \textrm{Gal}(L/K)$ is a conjugacy class, then we establish that the Chebotarev density is the following limit of partial sums of ideals of $K$: \[ -\lim_{X\rightarrow\infty} \sum_{\substack{2\leq N(I)\leq X \\ I \in S(L/K; C)}} \frac{\mu_K(I)}{N(I)} = \frac{|C|}{|G|}, \] where $\mu_K(I)$ denotes the generalized M\"obius function and $S(L/K;C)$ is the set of ideals $I\subset \mathcal{O}_K$ such that $I$ has a unique prime divisor $\mathfrak{p}$ of minimal norm and the Artin symbol $\left[\frac{L/K}{\mathfrak{p}}\right]$ is $C$. To obtain this formula, we generalize several results from classical analytic number theory, as well as Alladi's concept of duality for minimal and maximal prime divisors, to the setting of ideals in number fields.


title: Global Transition-based Non-projective Dependency Parsing

abstract: Shi, Huang, and Lee (2017) obtained state-of-the-art results for English and Chinese dependency parsing by combining dynamic-programming implementations of transition-based dependency parsers with a minimal set of bidirectional LSTM features. However, their results were limited to projective parsing. In this paper, we extend their approach to support non-projectivity by providing the first practical implementation of the MH_4 algorithm, an $O(n^4)$ mildly nonprojective dynamic-programming parser with very high coverage on non-projective treebanks. To make MH_4 compatible with minimal transition-based feature sets, we introduce a transition-based interpretation of it in which parser items are mapped to sequences of transitions. We thus obtain the first implementation of global decoding for non-projective transition-based parsing, and demonstrate empirically that it is more effective than its projective counterpart in parsing a number of highly non-projective languages


title: Power-law entropies for continuous systems and generalized operations

abstract: We present our view in a standing debate about the definition and meaning of power-law entropies for continuous systems. Our suggestion is that such arguments should take into account the generalized operations of addition and multiplication induced by the power-law entropies' composition properties. To be concrete, we highlight our view using the case of the $q$- also known as "Tsallis", entropic functionals.


title: Shannon entropy for intuitionistic fuzzy information

abstract: The paper presents an extension of Shannon fuzzy entropy for intuitionistic fuzzy one. Firstly, we presented a new formula for calculating the distance and similarity of intuitionistic fuzzy information. Then, we constructed measures for information features like score, certainty and uncertainty. Also, a new concept was introduced, namely escort fuzzy information. Then, using the escort fuzzy information, Shannon's formula for intuitionistic fuzzy information was obtained. It should be underlined that Shannon's entropy for intuitionistic fuzzy information verifies the four defining conditions of intuitionistic fuzzy uncertainty. The measures of its two components were also identified: fuzziness (ambiguity) and incompleteness (ignorance).


title: Significant acceleration of development by automating quality assurance of a medical particle accelerator safety system using a formal language driven test stand

abstract: At the Centre for Proton Therapy at the Paul Scherrer Institute cancer patients are treated with a fixed beamline and in two gantries for ocular and non-ocular malignancies, respectively. For the installation of a third gantry a new patient safety system (PaSS) was developed and is sequentially being rolled out to update the existing areas. The aim of PaSS is to interrupt the treatment whenever any sub-system detects a hazardous condition. To ensure correct treatment delivery, this system needs to be thoroughly tested as part of the regular quality assurance (QA) protocols as well as after any upgrade. In the legacy safety systems, unit testing required an extensive use of resources: two weeks of work per area in the laboratory in addition to QA beam time. In order to significantly reduce the time, an automated PaSS test stand for unit testing was developed based on a PXI chassis with virtually unlimited IOs that are synchronously stimulated or sampled at 1 MHz. It can emulate the rest of the facility using adapters to connect each type of interface. With it PaSS can be tested under arbitrary conditions. A VHDL-based formal language was developed to describe stimuli, expected behaviour and specific measurements, interpreted by a LabView runtime environment. This article describes the tools and methodology being applied for unit testing and QA release tests for the new PaSS. It shows how automation and formalization made possible an increase in test coverage while significantly cutting down the laboratory testing time and facility's beam usage.


title: Graph functionality

abstract: Let $G=(V,E)$ be a graph and $A$ its adjacency matrix. We say that a vertex $y \in V$ is a function of vertices $x_1, \ldots, x_k \in V$ if there exists a Boolean function $f$ of $k$ variables such that for any vertex $z \in V - \{y, x_1, \ldots, x_k\}$, $A(y,z)=f(A(x_1,z),\ldots,A(x_k,z))$. The functionality $fun(y)$ of vertex $y$ is the minimum $k$ such that $y$ is a function of $k$ vertices. The functionality $fun(G)$ of the graph $G$ is $\max\limits_H\min\limits_{y\in V(H)}fun(y)$, where the maximum is taken over all induced subgraphs $H$ of $G$. In the present paper, we show that functionality generalizes simultaneously several other graph parameters, such as degeneracy or clique-width, by proving that bounded degeneracy or bounded clique-width imply bounded functionality. Moreover, we show that this generalization is proper by revealing classes of graphs of unbounded degeneracy and clique-width, where functionality is bounded by a constant. This includes permutation graphs, unit interval graphs and line graphs. We also observe that bounded functionality implies bounded VC-dimension, i.e. graphs of bounded VC-dimension extend graphs of bounded functionality, and this extension is also proper.


title: Understanding and Accelerating Particle-Based Variational Inference

abstract: Particle-based variational inference methods (ParVIs) have gained attention in the Bayesian inference literature, for their capacity to yield flexible and accurate approximations. We explore ParVIs from the perspective of Wasserstein gradient flows, and make both theoretical and practical contributions. We unify various finite-particle approximations that existing ParVIs use, and recognize that the approximation is essentially a compulsory smoothing treatment, in either of two equivalent forms. This novel understanding reveals the assumptions and relations of existing ParVIs, and also inspires new ParVIs. We propose an acceleration framework and a principled bandwidth-selection method for general ParVIs; these are based on the developed theory and leverage the geometry of the Wasserstein space. Experimental results show the improved convergence by the acceleration framework and enhanced sample accuracy by the bandwidth-selection method.


title: Massively-Parallel Break Detection for Satellite Data

abstract: The field of remote sensing is nowadays faced with huge amounts of data. While this offers a variety of exciting research opportunities, it also yields significant challenges regarding both computation time and space requirements. In practice, the sheer data volumes render existing approaches too slow for processing and analyzing all the available data. This work aims at accelerating BFAST, one of the state-of-the-art methods for break detection given satellite image time series. In particular, we propose a massively-parallel implementation for BFAST that can effectively make use of modern parallel compute devices such as GPUs. Our experimental evaluation shows that the proposed GPU implementation is up to four orders of magnitude faster than the existing publicly available implementation and up to ten times faster than a corresponding multi-threaded CPU execution. The dramatic decrease in running time renders the analysis of significantly larger datasets possible in seconds or minutes instead of hours or days. We demonstrate the practical benefits of our implementations given both artificial and real datasets.


title: Magnetic field-free circularly polarized thermal emission from chiral metasurface

abstract: Thermal radiation from bulk disorderly placed nonresonant emitters is incoherent, broadband and isotropic. In an external magnetic field the thermal radiation from any source is circularly polarized. Here we propose a thermal radiation source which emits circularly polarized radiation and which is not placed in a magnetic field. The thermal source consists of a slab waveguide with etched chiral metasurface. Due to the absence of a mirror symmetry of the metasurface, the thermally generated electromagnetic waves become circularly polarized. In this letter we discuss the origin of this phenomenon in details. Using the Fourier modal method we analyze the eigenmodes of the structure and the emissivity spectra. We demonstrate that the degree of circular polarization in an optimized structure can be as high as 0.87.


title: Composition of rational functions: state-space realization and applications

abstract: We define two versions of compositions of matrix-valued rational functions of appropriate sizes and whenever analytic at infinity, offer a set of formulas for the corresponding state-space realization, in terms of the realizations of the original functions. Focusing on positive real functions, the first composition is applied to electrical circuits theory along with introducing a connection to networks of feedback loops. The second composition is applied to Stieltjes functions.


title: Moduli space singularities for $3d$ $\mathcal{N} = 4$ circular quiver gauge theories

abstract: The singularity structure of the Coulomb and Higgs branches of good $3d$ $\mathcal{N}=4$ circular quiver gauge theories (CQGTs) with unitary gauge groups is studied. The central method employed is the Kraft--Procesi transition. CQGTs are described as a generalisation of a class of linear quivers. This class degenerates into the familiar class $T_{\rho}^{\sigma}(SU(N))$ in the linear case, however the circular case does not have the degeneracy and so the class of CQGTs contains many more theories and much more structure. We describe a collection of good, unitary, CQGTs from which the entire class can be found using Kraft--Procesi transitions. The singularity structure of a general member of this collection is fully determined, encompassing the singularity structure of a generic CQGT. Higher-level Hasse diagrams are introduced in order to write the results compactly. In higher-level Hasse diagrams, single nodes represent lattices of nilpotent orbit Hasse diagrams and edges represent traversing structure between lattices. The results generalise the case of linear quiver moduli spaces which are known to be nilpotent varieties of $\mathfrak{sl}_n$.


title: Stellar masses, metallicity gradients and suppressed star formation revealed in a new sample of absorption selected galaxies

abstract: Context. Absorbing galaxies are selected via the detection of characteristic absorption lines which their gas-rich media imprint in the spectra of distant light-beacons. The proximity of the typically faint foreground absorbing galaxies to bright background sources makes it challenging to robustly identify these in emission, and hence to characterise their relation to the general galaxy population. Aims. We search for emission to confirm and characterise ten galaxies hosting damped, metal-rich quasar absorbers at redshift z < 1. Methods. We identify the absorbing galaxies by matching spectroscopic absorption -and emission redshifts and from projected separations. Combining emission-line diagnostics with existing absorption spectroscopy and photometry of quasar-fields hosting metal-rich, damped absorbers, we compare our new detections with reference samples and place them on scaling relations. Results. We spectroscopically confirm seven galaxies harbouring damped absorbers (a 70% success-rate). Our results conform to the emerging picture that neutral gas on scales of tens of kpc in galaxies is what causes the characteristic Hi absorption. Our key results are: (I) Absorbing galaxies with $\log _{10} [M_\star ~(M_\odot)] \gtrsim 10$ have star formation rates that are lower than predicted for the main sequence of star formation. (II) The distribution of impact parameter with Hi column density and with absorption-metallicity for absorbing galaxies at $z\sim 2-3$ extends to $z\sim 0.7$ and to lower Hi column densities. (III) A robust mean metallicity gradient of $\langle \Gamma \rangle = 0.022 \pm 0.001~[dex~kpc^{-1}]$. (IV) By correcting absorption metallicities for $\langle \Gamma \rangle$ and imposing a truncation-radius at $12~\mathrm{kpc}$, absorbing galaxies fall on top of predicted mass-metallicity relations, with a statistically significant decrease in scatter.


title: Option Pricing with Heavy-Tailed Distributions of Logarithmic Returns

abstract: A growing body of literature suggests that heavy tailed distributions represent an adequate model for the observations of log returns of stocks. Motivated by these findings, here we develop a discrete time framework for pricing of European options. Probability density functions of log returns for different periods are conveniently taken to be convolutions of the Student's t-distribution with three degrees of freedom. The supports of these distributions are truncated in order to obtain finite values for the options. Within this framework, options with different strikes and maturities for one stock rely on a single parameter -- the standard deviation of the Student's t-distribution for unit period. We provide a study which shows that the distribution support width has weak influence on the option prices for certain range of values of the width. It is furthermore shown that such family of truncated distributions approximately satisfies the no-arbitrage principle and the put-call parity. The relevance of the pricing procedure is empirically verified by obtaining remarkably good match of the numerically computed values by our scheme to real market data.


title: Thermal and transport properties of U3Si2

abstract: We have studied U_3Si_2 by means of the heat capacity, electrical resistivity, Seebeck and Hall effects, and thermal conductivity in the temperature range 2-300 K and in magnetic fields up to 9 T. All the results obtained point to delocalized nature of 5f-electrons in this material. The low temperature heat capacity is enhanced (gamma_el ~ 150 mJ/mol-K2) and shows an upturn in Cp/T (T), characteristic of spin fluctuations. The thermal conductivity of U3Si2 is ~8.5 W/m-K at room temperature and we show that the electronic part dominates heat transport above 300 K as expected for a metallic system, although the lattice contribution cannot be completely neglected.


title: Path integral for quantum Mabuchi K-energy

abstract: We construct a path integral based on the coupling of the Liouville action and the Mabuchi K-energy on a one-dimensional complex manifold. To the best of our knowledge this is the first rigorous construction of such an object and this is done by means of probabilistic tools. Both functionals play an important role respectively in Riemannian geometry (in the case of surfaces) and K\"ahler geometry. As an output, we obtain a path integral whose Weyl anomaly displays the standard Liouville anomaly plus an additional K-energy term. Motivations come from theoretical physics where these type of path integrals arise as a model for fluctuating metrics on surfaces when coupling (small) massive perturbations of conformal field theories to quantum gravity as advocated by A. Bilal, F. Ferrari, S. Klevtsov and S. Zelditch. Interestingly, our computations show that quantum corrections perturb the classical Mabuchi K-energy and produce a quantum Mabuchi K-energy: this type of correction is reminiscent of the quantum Liouville theory. Our construction is probabilistic and relies on a variant of Gaussian multiplicative chaos (GMC), the Derivative GMC (DGMC for short). The technical backbone of our construction consists in two estimates on (derivative and standard) GMC which are of independent interest in probability theory. Firstly, we show that these DGMC random variables possess negative exponential moments and secondly we derive optimal small deviations estimates for the GMC associated with a recentered Gaussian Free Field.


title: Learning Personalized Representation for Inverse Problems in Medical Imaging Using Deep Neural Network

abstract: Recently deep neural networks have been widely and successfully applied in computer vision tasks and attracted growing interests in medical imaging. One barrier for the application of deep neural networks to medical imaging is the need of large amounts of prior training pairs, which is not always feasible in clinical practice. In this work we propose a personalized representation learning framework where no prior training pairs are needed, but only the patient's own prior images. The representation is expressed using a deep neural network with the patient's prior images as network input. We then applied this novel image representation to inverse problems in medical imaging in which the original inverse problem was formulated as a constraint optimization problem and solved using the alternating direction method of multipliers (ADMM) algorithm. Anatomically guided brain positron emission tomography (PET) image reconstruction and image denoising were employed as examples to demonstrate the effectiveness of the proposed framework. Quantification results based on simulation and real datasets show that the proposed personalized representation framework outperform other widely adopted methods.


title: Stable topological edge states in a non-Hermitian four-band model

abstract: We introduce a one dimensional non-Hermitian four band tight binding lattice system. We find stable topological edge states protected by particle-hole and parity-time symmetries. We show that topological phase appears in the system. We discuss that the system can also be used as a topological laser if the gain and loss positions in a unit cell are judiciously arranged.


title: The SPectral Image Typer (SPIT)

abstract: We present the Spectral Image Typer (SPIT), a convolutional neural network (CNN) built to classify spectral images. In contrast to traditional, rules-based algorithms which rely on meta data provided with the image (e.g. header cards), SPIT is trained solely on the image data. We have trained SPIT on 2,004 human-classified images taken with the Kast spectrometer at Lick Observatory with types of Bias, Arc, Flat, Science and Standard. We include several pre-processing steps (scaling, trimming) motivated by human practice and also expanded the training set to balance between image type and increase diversity. The algorithm achieved an accuracy of 98.7% on the held-out validation set and an accuracy of 98.7% on the test set of images. We then adopt a slightly modified classification scheme to improve robustness at a modestly reduced cost in accuracy (98.2%). The majority of mis-classifications are Science frames with very faint sources confused with Arc images (e.g. faint emission-line galaxies) or Science frames with very bright sources confused with Standard stars. These are errors that even a well-trained human is prone to make. Future work will increase the training set from Kast, will include additional optical and near-IR instruments, and may expand the CNN architecture complexity. We are now incorporating SPIT in the PYPIT data reduction pipeline (DRP) and are willing to facilitate its inclusion in other DRPs.


title: Random cherry graphs

abstract: Due to the popularity of randomly evolving graph processes, there exists a randomized version of many recursively defined graph models. This is also the case with the cherry tree, which was introduced by Buksz\'ar and Pr\'ekopa to improve Bonferroni type upper bounds on the probability of the union of random events. Here we consider a substantially extended random analogue of that model, embedding it into a general time dependent branching process.


title: Seq2RDF: An end-to-end application for deriving Triples from Natural Language Text

abstract: We present an end-to-end approach that takes unstructured textual input and generates structured output compliant with a given vocabulary. Inspired by recent successes in neural machine translation, we treat the triples within a given knowledge graph as an independent graph language and propose an encoder-decoder framework with an attention mechanism that leverages knowledge graph embeddings. Our model learns the mapping from natural language text to triple representation in the form of subject-predicate-object using the selected knowledge graph vocabulary. Experiments on three different data sets show that we achieve competitive F1-Measures over the baselines using our simple yet effective approach. A demo video is included.


title: Life-times of quantum resonances through the Geometrical Phase Propagator Approach

abstract: We employ the recently introduced Geometric Phase Propagator Approach (GPPA) [Phys. Rev. A85, 062110 (2012)] to develop an improved perturbative scheme for the calculation of life times in driven quantum systems. This incorporates a resummation of the contributions of virtual processes starting and ending at the same state in the considered time interval. The proposed procedure allows for a strict determination of the conditions leading to finite life times in a general driven quantum system by isolating the resummed terms in the perturbative expansion contributing to their generation. To illustrate how the derived conditions apply in practice, we consider the effect of driving in a system with purely discrete energy spectrum, as well as in a system for which the eigenvalue spectrum contains a continuous part. We show that in the first case, when the driving contains a dense set of frequencies acting as a noise to the system, the corresponding bound states acquire a finite life time. When the energy spectrum contains also a continuum set of eigenvalues then the bound states, due to the driving, couple to the continuum and become quasi-bound resonances. The benchmark of this change is the appearance of a Fano-type peak in the associated transmission profile. In both cases the corresponding life-time can be efficiently estimated within the reformulated GPPA approach.


title: Monte Carlo algorithm based on internal bridging moves for the atomistic simulation of thiophene oligomers and polymers

abstract: We introduce a powerful Monte Carlo (MC) algorithm for the atomistic simulation of bulk models of oligo- and poly-thiophenes by redesigning MC moves originally developed for considerably simpler polymer structures and architectures, such as linear and branched polyethylene, to account for the ring structure of the thiophene monomer. Elementary MC moves implemented include bias reptation of an end thiophene ring, flip of an internal thiophene ring, rotation of an end thiophene ring, concerted rotation of three thiophene rings, rigid translation of an entire molecule, rotation of an entire molecule and volume fluctuation. In the implementation of all moves we assume that thiophene ring atoms remain rigid and strictly co-planar; on the other hand, inter-ring torsion and bond bending angles remain fully flexible subject to suitable potential energy functions. Test simulations with the new algorithm of an important thiophene oligomer, {\alpha}-sexithiophene ({\alpha}-6T), at a high enough temperature (above its isotropic-to-nematic phase transition) using a new united atom model specifically developed for the purpose of this work provide predictions for the volumetric, conformational and structural properties that are remarkably close to those obtained from detailed atomistic Molecular Dynamics (MD) simulations using an all-atom model. The new algorithm is particularly promising for exploring the rich (and largely unexplored) phase behavior and nanoscale ordering of very long (also more complex) thiophene-based polymers which cannot be addressed by conventional MD methods due to the extremely long relaxation times characterizing chain dynamics in these systems.


title: Renormalization in Condensed Matter: Fermionic Systems - from Mathematics to Materials

abstract: Renormalization plays an important role in the theoretically and mathematically careful analysis of models in condensed-matter physics. I review selected results about correlated-fermion systems, ranging from mathematical theorems to applications in models relevant for materials science, such as the prediction of equilibrium phases of systems with competing ordering tendencies, and quantum criticality.


title: Fisher information analysis of list-mode SPECT emission data for joint estimation of activity and attenuation distribution

abstract: The potential to perform attenuation and scatter compensation (ASC) in single-photon emission computed tomography (SPECT) imaging using only the SPECT emission data is highly significant. In this context, attenuation in SPECT is primarily due to Compton scattering, where the probability of Compton scatter is proportional to the attenuation coefficient of the tissue and the energy of the scattered photon and the scattering angle are related. Given this premise, we investigate whether the SPECT scattered-photon data acquired in list-mode (LM) format and including the energy information can be used to estimate the attenuation map. For this purpose, we propose a Fisher-information-based method that yields the Cramer-Rao bound (CRB) for the task of jointly estimating the activity/attenuation distribution using only the SPECT emission data. The proposed method is applied to analyze the information content of SPECT LM emission data in a 2D SPECT system using computational studies with digital phantoms for different photon-count levels. The results show that scattered photons contain information to estimate the attenuation coefficients. An increase in the number of detected photons leads to lower CRB for both the attenuation and activity coefficients. Also, the CRB obtained for the attenuation and activity coefficients is typically much lower than the true value of these coefficients. Further, processing the emission data in LM format yields a lower CRB in comparison to binning data. Finally, we observe that systems with better energy resolution yield a lower CRB for the attenuation coefficient. Overall, the results provide strong evidence that LM SPECT emission data, including the scattered photons, contains information to jointly estimate the activity and attenuation coefficients.


title: Mechanisms of Calcium Leak from Cardiac Sarcoplasmic Reticulum Revealed by Statistical Mechanics

abstract: Heart muscle contraction is normally activated by a synchronized Ca release from sarcoplasmic reticulum (SR), a major intracellular Ca store. However, under abnormal conditions Ca leaks from the SR, decreasing heart contraction amplitude and increasing risk of life-threatening arrhythmia. The mechanisms and regimes of SR operation generating the abnormal Ca leak remain unclear. Here we employed both numerical and analytical modeling to get mechanistic insights into the emergent Ca leak phenomenon. Our numerical simulations using a detailed realistic model of Ca release unit (CRU) reveal sharp transitions resulting in Ca leak. The emergence of leak is closely mapped mathematically to the Ising model from statistical mechanics. The system steady-state behavior is determined by two aggregate parameters: the analogues of magnetic field ($h$) and the inverse temperature ($\beta$) in the Ising model, for which we have explicit formulas in terms of SR Ca and release channel opening/closing rates. The classification of leak regimes takes the shape of a phase $\beta$-$h$ diagram, with the regime boundaries occurring at $h$=0 and a critical value of $\beta$ ($\beta*$) which we estimate using a classical Ising model and mean field theory. Our theory predicts that a synchronized Ca leak will occur when $h$>0 and $\beta>\beta*$ and a disordered leak occurs when $\beta<\beta*$ and $h$ is not too negative. The disorder leak is distinguished from synchronized leak (in long-lasting sparks) by larger Peierls contour lengths, an output parameter reflecting degree of disorder. Thus, in addition to our detailed numerical model approach we also offer an instantaneous computational tool using analytical formulas of the Ising model for respective RyR parameters and SR Ca load that describe and classify phase transitions and leak emergence.


title: FluidSim: modular, object-oriented Python package for high-performance CFD simulations

abstract: The Python package fluidsim is introduced in this article as an extensible framework for Computational Fluid Mechanics (CFD) solvers. It is developed as a part of FluidDyn project (Augier et al., 2018), an effort to promote open-source and open-science collaboration within fluid mechanics community and intended for both educational as well as research purposes. Solvers in fluidsim are scalable, High-Performance Computing (HPC) codes which are powered under the hood by the rich, scientific Python ecosystem and the Application Programming Interfaces (API) provided by fluiddyn and fluidfft packages (Mohanan et al., 2018). The present article describes the design aspects of fluidsim, viz. use of Python as the main language; focus on the ease of use, reuse and maintenance of the code without compromising performance. The implementation details including optimization methods, modular organization of features and object-oriented approach of using classes to implement solvers are also briefly explained. Currently, fluidsim includes solvers for a variety of physical problems using different numerical methods (including finite-difference methods). However, this metapaper shall dwell only on the implementation and performance of its pseudo-spectral solvers, in particular the two- and three-dimensional Navier-Stokes solvers. We investigate the performance and scalability of fluidsim in a state of the art HPC cluster. Three similar pseudo-spectral CFD codes based on Python (Dedalus, SpectralDNS) and Fortran (NS3D) are presented and qualitatively and quantitatively compared to fluidsim. The source code is hosted at Bitbucket as a Mercurial repository bitbucket.org/fluiddyn/fluidsim and the documentation generated using Sphinx can be read online at fluidsim.readthedocs.io.


title: The Systems of Computer Mathematics in the Cloud-Based Learning Environment of Educational Institutions

abstract: The article highlights the promising ways of providing access to the mathematical software in higher educational institutions. It is emphasized that the cloud computing services implementation is the actual trend of modern ICT pedagogical systems development. The analysis and evaluation of existing ex-perience and educational research of different types of mathematical software packages use are proposed. The methodological guidance on the organization of the cloud-based learning environment in higher educational institutions using the systems of computer mathematics (SCM) is outlined. The advantages and disadvantages of different cloud service models of access to SCM are de-scribed. The cloud-based learning component with the use of Maxima system is described and evaluated.


title: Direct Uncertainty Prediction for Medical Second Opinions

abstract: The issue of disagreements amongst human experts is a ubiquitous one in both machine learning and medicine. In medicine, this often corresponds to doctor disagreements on a patient diagnosis. In this work, we show that machine learning models can be trained to give uncertainty scores to data instances that might result in high expert disagreements. In particular, they can identify patient cases that would benefit most from a medical second opinion. Our central methodological finding is that Direct Uncertainty Prediction (DUP), training a model to predict an uncertainty score directly from the raw patient features, works better than Uncertainty Via Classification, the two-step process of training a classifier and postprocessing the output distribution to give an uncertainty score. We show this both with a theoretical result, and on extensive evaluations on a large scale medical imaging application.


title: Mapping the research software sustainability space

abstract: A growing number of largely uncoordinated initiatives focus on research software sustainability. A comprehensive mapping of the research software sustainability space can help identify gaps in their efforts, track results, and avoid duplication of work. To this end, this paper suggests enhancing an existing schematic of activities in research software sustainability, and formalizing it in a directed graph model. Such a model can be further used to define a classification schema which, applied to research results in the field, can drive the identification of past activities and the planning of future efforts.


title: Semi-infinite cohomology and Kazhdan-Lusztig equivalence at positive level

abstract: A positive level Kazhdan-Lusztig functor is defined using Arkhipov-Gaitsgory duality for affine Lie algebras. The functor sends objects in the DG category of G(O)-equivariant positive level affine Lie algebra modules to objects in the DG category of modules over Lusztig's quantum group at a root of unity. We prove that the semi-infinite cohomology functor for positive level modules factors through the Kazhdan-Lusztig functor at positive level and the quantum group cohomology functor with respect to the positive part of Lusztig's quantum group.


title: BOHB: Robust and Efficient Hyperparameter Optimization at Scale

abstract: Modern deep learning methods are very sensitive to many hyperparameters, and, due to the long training times of state-of-the-art models, vanilla Bayesian hyperparameter optimization is typically computationally infeasible. On the other hand, bandit-based configuration evaluation approaches based on random search lack guidance and do not converge to the best configurations as quickly. Here, we propose to combine the benefits of both Bayesian optimization and bandit-based methods, in order to achieve the best of both worlds: strong anytime performance and fast convergence to optimal configurations. We propose a new practical state-of-the-art hyperparameter optimization method, which consistently outperforms both Bayesian optimization and Hyperband on a wide range of problem types, including high-dimensional toy functions, support vector machines, feed-forward neural networks, Bayesian neural networks, deep reinforcement learning, and convolutional neural networks. Our method is robust and versatile, while at the same time being conceptually simple and easy to implement.


title: FluidFFT: common API (C++ and Python) for Fast Fourier Transform HPC libraries

abstract: The Python package fluidfft provides a common Python API for performing Fast Fourier Transforms (FFT) in sequential, in parallel and on GPU with different FFT libraries (FFTW, P3DFFT, PFFT, cuFFT). fluidfft is a comprehensive FFT framework which allows Python users to easily and efficiently perform FFT and the associated tasks, such as as computing linear operators and energy spectra. We describe the architecture of the package composed of C++ and Cython FFT classes, Python "operator" classes and Pythran functions. The package supplies utilities to easily test itself and benchmark the different FFT solutions for a particular case and on a particular machine. We present a performance scaling analysis on three different computing clusters and a microbenchmark showing that fluidfft is an interesting solution to write efficient Python applications using FFT.


title: Disentangling basal and accrued height-for-age for cross-population comparisons

abstract: Objectives: Current standards for comparing stunting across human populations assume a universal model of child growth. Such comparisons ignore population differences that are independent of deprivation and health outcomes. This paper partitions variation in height-for-age that is specifically associated with deprivation and health outcomes to provide a basis for cross-population comparisons. Materials & Methods: Using a multi-level model with a sigmoid relationship of resources and growth, we partition variation in height-for-age z-scores (HAZ) from 1,522,564 children across 70 countries into two components: 1) "accrued HAZ" shaped by environmental inputs (e.g., undernutrition, infectious disease, inadequate sanitation, poverty), and 2) a country-specific "basal HAZ" independent of such inputs. We validate these components against population-level infant mortality rates, and assess how these basal differences may affect cross-population comparisons of stunting. Results: Basal HAZ differs reliably across countries (range of 1.5 SD) and is independent of measures of infant mortality. By contrast, accrued HAZ captures stunting as impaired growth due to deprivation and is more closely associated with infant mortality than observed HAZ. Ranking populations by accrued HAZ suggest that populations in West Africa and the Caribbean suffer much greater levels of stunting than suggested by observed HAZ. Discussion: Current universal standards may dramatically underestimate stunting in populations with taller basal HAZ. Relying on observed HAZ rather than accrued HAZ may also lead to inappropriate cross-population comparisons, such as concluding that Haitian children enjoy better conditions for growth than do Indian or Guatemalan children.


title: The Quest for Gravity Agent: from Newton to Einstein and Feynman

abstract: Detection of gravitational waves by LIGO-Virgo observatories disclosed the agent of gravity (gravitons) and together with information paradox of black hole event horizon points to paradigm shift from non-quantum geometry to quantum gravity physics. Historical, philosophical, mathematical, physical and astronomical arguments are presented for demonstration that in the beginning of the 21st century we are on the verge of the forthcoming change in understanding of gravity nature, which will bring new ideas in performing gravitational experiments and interpretation of astrophysical phenomena.


title: Uncertainty Quantification of Electronic and Photonic ICs with Non-Gaussian Correlated Process Variations

abstract: Since the invention of generalized polynomial chaos in 2002, uncertainty quantification has impacted many engineering fields, including variation-aware design automation of integrated circuits and integrated photonics. Due to the fast convergence rate, the generalized polynomial chaos expansion has achieved orders-of-magnitude speedup than Monte Carlo in many applications. However, almost all existing generalized polynomial chaos methods have a strong assumption: the uncertain parameters are mutually independent or Gaussian correlated. This assumption rarely holds in many realistic applications, and it has been a long-standing challenge for both theorists and practitioners. This paper propose a rigorous and efficient solution to address the challenge of non-Gaussian correlation. We first extend generalized polynomial chaos, and propose a class of smooth basis functions to efficiently handle non-Gaussian correlations. Then, we consider high-dimensional parameters, and develop a scalable tensor method to compute the proposed basis functions. Finally, we develop a sparse solver with adaptive sample selections to solve high-dimensional uncertainty quantification problems. We validate our theory and algorithm by electronic and photonic ICs with 19 to 57 non-Gaussian correlated variation parameters. The results show that our approach outperforms Monte Carlo by $2500\times$ to $3000\times$ in terms of efficiency. Moreover, our method can accurately predict the output density functions with multiple peaks caused by non-Gaussian correlations, which is hard to handle by existing methods. Based on the results in this paper, many novel uncertainty quantification algorithms can be developed and can be further applied to a broad range of engineering domains.


title: Synthetic contrast enhancement in cardiac CT with Deep Learning

abstract: In Europe the 20% of the CT scans cover the thoracic region. The acquired images contain information about the cardiovascular system that often remains latent due to the lack of contrast in the cardiac area. On the other hand, the contrast enhanced computed tomography (CECT) represents an imaging technique that allows to easily assess the cardiac chambers volumes and the contrast dynamics. With this work we aim to face the problem of extraction and presentation of these latent information, using a deep learning approach with convolutional neural networks. Starting from the extraction of relevant features from the image without contrast medium, we try to re-map them on features typical of CECT, to synthesize an image characterized by an attenuation in the cardiac chambers as if a virtually iodine contrast medium was injected. The purposes are to guarantee an estimation of the left cardiac chambers volume and to perform an evaluation of the contrast dynamics. Our approach is based on a deconvolutional network trained on a set of 120 patients who underwent both CT acquisitions in the same contrastographic arterial phase and the same cardiac phase. To ensure a reliable predicted CECT image, in terms of values and morphology, a custom loss function is defined by combining an error function to find a pixel-wise correspondence, which takes into account the similarity in term of Hounsfield units between the input and output images and by a cross-entropy computed on the binarized versions of the synthesized and of the real CECT image. The proposed method is finally tested on 20 subjects.


title: An Efficient Algorithm for Automatic Structure Optimization in X-ray Standing-Wave Experiments

abstract: X-ray standing-wave photoemission experiments involving multilayered samples are emerging as unique probes of the buried interfaces that are ubiquitous in current device and materials research. Such data require for their analysis a structure optimization process comparing experiment to theory that is not straightforward. In this work, we present a new computer program for optimizing the analysis of standing-wave data, called SWOPT, that automates this trial-and-error optimization process. The program includes an algorithm that has been developed for computationally expensive problems: so-called black-box simulation optimizations. It also includes a more efficient version of the Yang X-ray Optics Program (YXRO) [Yang, S.-H., Gray, A.X., Kaiser, A.M., Mun, B.S., Sell, B.C., Kortright, J.B., Fadley, C.S., J. Appl. Phys. 113, 1 (2013)] which is about an order of magnitude faster than the original version. Human interaction is not required during optimization. We tested our optimization algorithm on real and hypothetical problems and show that it finds better solutions significantly faster than a random search approach. The total optimization time ranges, depending on the sample structure, from minutes to a few hours on a modern laptop computer, and can be up to 100x faster than a corresponding manual optimization. These speeds make the SWOPT program a valuable tool for realtime analyses of data during synchrotron experiments.