title: Numerical study of the chiral $\mathbb{Z}_3$ quantum phase transition in one spatial dimension

abstract: Recent experiments on a one-dimensional chain of trapped alkali atoms [arXiv:1707.04344] have observed a quantum transition associated with the onset of period-3 ordering of pumped Rydberg states. This spontaneous $\mathbb{Z}_3$ symmetry breaking is described by a constrained model of hard-core bosons proposed by Fendley $et\, \,al.$ [arXiv:cond-mat/0309438]. By symmetry arguments, the transition is expected to be in the universality class of the $\mathbb{Z}_3$ chiral clock model with parameters preserving both time-reversal and spatial-inversion symmetries. We study the nature of the order-disorder transition in these models, and numerically calculate its critical exponents with exact diagonalization and density-matrix renormalization group techniques. We use finite-size scaling to determine the dynamical critical exponent $z$ and the correlation length exponent $\nu$. Our analysis presents the only known instance of a strongly-coupled transition between gapped states with $z \ne 1$, implying an underlying nonconformal critical field theory.

categories: cond-mat.str-el cond-mat.quant-gas cond-mat.stat-mech physics.atom-ph


title: The Gaia-ESO survey: impact of extra-mixing on C- and N-abundances of giant stars

abstract: The GES survey using FLAMES at the VLT has obtained high-resolution UVES spectra for a large number of giant stars, allowing a determination of the abundances of the key chemical elements C and N at their surface. The surface abundances of these chemical species are well-known to change in stars during their evolution on the red giant branch after the first dredge-up episod, as a result of extra-mixing phenomena. We investigate the effects of thermohaline mixing on C and N abundances using the first comparison between the GES [C/N] determinations with simulations of the observed fields using a model of stellar population synthesis. We explore the effects of thermohaline mixing on the chemical properties of giants through stellar evolutionary models computed with the stellar evolution code STAREVOL. We include these stellar evolution models in the Besan\c{c}on Galaxy model to simulate the [C/N] distributions determined from the UVES spectra of the GES and compare them with the observations. Theoretical predictions including the effect of thermohaline mixing are in good agreement with the observations. However, the field stars in the GES with C and N-abundance measurements have a metallicity close to solar, where the efficiency of thermohaline mixing is not very large. The C and N abundances derived by the GES in open and globular clusters clearly show the impact of thermohaline mixing at low-metallicity, allowing to explain the [C/N] ratio observed in lower-mass and older giant stars. Using independent observations of carbon isotopic ratio in clump field stars and open clusters, we also confirm that thermohaline mixing should be taken into account to explain the behavior of 12C/13C ratio as a function of stellar age. Overall the current model including thermohaline mixing is able to reproduce very well the C- and N-abundances over the whole metallicity range investigated by the GES data.


title: What do gravitons say about (unimodular) gravity?

abstract: We revisit the problem of constraining the weak field limit of the gravitational lagrangian from S-matrix properties. From unitarity and Lorentz invariance of the S-matrix of massless gravitons, we derive on-shell gauge invariance to consist on the transverse part of the linearised diffeomorphisms group. Moreover, by looking to the interaction between sources, we conclude that there exist only two possible lagrangians that lead to a well-defined covariant interaction, corresponding to the weak field limits of General Relativity and Unimodular Gravity. Additionally, this result confirms the equivalence of the S-matrix of both theories around flat space-time.


title: Exact Simulation of the Extrema of Stable Processes

abstract: We exhibit an exact simulation algorithm for the supremum of a stable process over a finite time interval using dominated coupling from the past (DCFTP). We establish a novel perpetuity equation for the supremum (via the representation of the concave majorants of L\'evy processes) and apply it to construct a Markov chain in the DCFTP algorithm. We prove that the number of steps taken backwards in time before the coalescence is detected is finite. We analyse numerically the performance of the algorithm (the code, written in Julia 1.0, is available on GitHub).


title: Statistical properties of paired fixed fields

abstract: The initial conditions of cosmological simulations are commonly drawn from a Gaussian ensemble. The limited number of modes inside a simulation volume gives rise to statistical fluctuations known as \textit{sample variance}, limiting the accuracy of simulation predictions. Fixed fields offer an alternative initialization strategy; they have the same power spectrum as standard Gaussian fields but without intrinsic amplitude scatter at linear order. Paired fixed fields consists of two fixed fields with opposite phases that cancel phase correlations which otherwise induce second-order scatter in the non-linear power spectrum. We study the statistical properties of those fields for 19 different quantities at different redshifts through a large set of 600 N-body and 506 state-of-the-art magneto-hydrodynamic simulations covering a wide range of scales, mass and spatial resolutions. We find that paired fixed simulations do not introduce a bias on any of the examined quantities. We quantify the statistical improvement brought by these simulations, over standard ones, on different power spectra such as matter, halos, CDM, gas, stars, black-holes and magnetic fields, finding that they can reduce their variance by factors as large as $10^6$. We quantify the improvement achieved by fixing and by pairing, showing that sample variance in some quantities can be highly suppressed by pairing after fixing. Paired fixed simulations do not change the scatter in quantities such as the probability distribution function of matter density, or the halo, void or stellar mass functions. We argue that procedures aiming at reducing the sample variance of those quantities are unlikely to work. Our results show that paired fixed simulations do not affect either mean relations or scatter of galaxy properties, and suggest that the information embedded in 1-pt statistics is highly complementary to that in clustering.


title: Observational Signature of the Logarithmic Terms in the Soft Graviton Theorem

abstract: We show that the recently discovered logarithmic terms in the soft graviton theorem induce a late time component in the gravitational wave-form that falls off as inverse power of time, producing a tail term to the linear memory effect.


title: Focal Visual-Text Attention for Visual Question Answering

abstract: Recent insights on language and vision with neural networks have been successfully applied to simple single-image visual question answering. However, to tackle real-life question answering problems on multimedia collections such as personal photos, we have to look at whole collections with sequences of photos or videos. When answering questions from a large collection, a natural problem is to identify snippets to support the answer. In this paper, we describe a novel neural network called Focal Visual-Text Attention network (FVTA) for collective reasoning in visual question answering, where both visual and text sequence information such as images and text metadata are presented. FVTA introduces an end-to-end approach that makes use of a hierarchical process to dynamically determine what media and what time to focus on in the sequential data to answer the question. FVTA can not only answer the questions well but also provides the justifications which the system results are based upon to get the answers. FVTA achieves state-of-the-art performance on the MemexQA dataset and competitive results on the MovieQA dataset.


title: Mass Hierarchies and Dynamical Field Range

abstract: Several swampland conjectures suggest that there is a critical field range beyond which the effective field theory (EFT) description breaks down in quantum gravity. In applications of these conjectures, however, the field range of interest is the field space distance traced by the physical trajectory that solves the equations of motion. We refer to this field space distance as the dynamical field range. We show that in the absence of a mass hierarchy between the light and heavy fields, the trajectory of the light field does not, in general, follow a geodesic in field space. Then, stabilizing the heavy fields at the minimum of their potential does not accurately describe the dynamics of the light field in general. A mass hierarchy can delay the breakdown of the EFT, and extend the effective field range. We illustrate these subtleties of multi-field dynamics with axions in Type II string compactifications.


title: EEG-GAN: Generative adversarial networks for electroencephalograhic (EEG) brain signals

abstract: Generative adversarial networks (GANs) are recently highly successful in generative applications involving images and start being applied to time series data. Here we describe EEG-GAN as a framework to generate electroencephalographic (EEG) brain signals. We introduce a modification to the improved training of Wasserstein GANs to stabilize training and investigate a range of architectural choices critical for time series generation (most notably up- and down-sampling). For evaluation we consider and compare different metrics such as Inception score, Frechet inception distance and sliced Wasserstein distance, together showing that our EEG-GAN framework generated naturalistic EEG examples. It thus opens up a range of new generative application scenarios in the neuroscientific and neurological context, such as data augmentation in brain-computer interfacing tasks, EEG super-sampling, or restoration of corrupted data segments. The possibility to generate signals of a certain class and/or with specific properties may also open a new avenue for research into the underlying structure of brain signals.


title: When Freeze-out Precedes Freeze-in: Sub-TeV Fermion Triplet Dark Matter with Radiative Neutrino Mass

abstract: We propose a minimal predictive scenario for dark matter and radiative neutrino mass where the relic abundance of dark matter is generated from a hybrid setup comprising of both thermal freeze-out as well as non-thermal freeze-in mechanisms. Considering three copies of fermion triplets and one scalar doublet, odd under an unbroken $\mathbb{Z}_2$ symmetry, to be responsible for radiative origin of neutrino mass, we consider the lightest fermion triplet as a dark matter candidate which remains under-abundant in the sub-TeV regime from usual thermal freeze-out. Late decay of the $\mathbb{Z}_2$-odd scalar doublet into dark matter serves as the non-thermal (freeze-in) contribution which not only fills the thermal dark matter deficit, but also constrains the mother particle's parameter space so that the correct relic abundance of dark matter is generated. Apart from showing interesting differences from the purely freeze-out and purely freeze-in dark matter scenarios, the model remains testable through disappearing charge track signatures at colliders, observable direct and indirect detection rates for dark matter and prediction of almost vanishing lightest neutrino mass.


title: Chains in CR geometry as geodesics of a Kropina metric

abstract: With the help of a generalization of the Fermat principle in general relativity, we show that chains in CR geometry are geodesics of a certain Kropina metric constructed from the CR structure. We study the projective equivalence of Kropina metrics and show that if the kernel distributions of the corresponding 1-forms are non-integrable then two projectively equivalent metrics are trivially projectively equivalent. As an application, we show that sufficiently many chains determine the CR structure up to conjugacy, generalizing and reproving the main result of [J.-H. Cheng, 1988]. The correspondence between geodesics of the Kropina metric and chains allows us to use the methods of metric geometry and the calculus of variations to study chains. We use these methods to re-prove the result of [H. Jacobowitz, 1985] that locally any two points of a strictly pseudoconvex CR manifolds can be joined by a chain. Finally, we generalize this result to the global setting by showing that any two points of a connected compact strictly pseudoconvex CR manifold which admits a pseudo-Einstein contact form with positive Tanaka-Webster scalar curvature can be joined by a chain.


title: Growth strategy determines network performance

abstract: The interplay between structure and function is crucial in determining some emerging properties of many natural systems. Here we use an adaptive neural network model inspired in observations of synaptic pruning that couples activity and topological dynamics and reproduces experimental temporal profiles of synaptic density, including an initial transient period of relatively high synaptic connectivity. Using a simplified framework, we prove that the existence of this transient is critical in providing ordered stationary states that have the property of being able to store stable memories. In fact, there is a discontinuous phase transition between the ordered memory phase and a disordered one as a function of the initial transient synaptic density. We also show that intermediate synaptic density values are optimal in order to obtain these stable memory states with a minimum energy consumption, and that ultimately it is the transient heterogeneity in the network what determines the stationary state. Our results here could explain why the pruning curves observed in actual brain areas present their characteristic temporal profiles and, eventually, anomalies such as autism and schizophrenia associated, respectively, with a deficit or an excess of pruning.


title: An explicit analysis of the entropic penalty in linear programming

abstract: Solving linear programs by using entropic penalization has recently attracted new interest in the optimization community, since this strategy forms the basis for the fastest-known algorithms for the optimal transport problem, with many applications in modern large-scale machine learning. Crucial to these applications has been an analysis of how quickly solutions to the penalized program approach true optima to the original linear program. More than 20 years ago, Cominetti and San Mart\'in showed that this convergence is exponentially fast; however, their proof is asymptotic and does not give any indication of how accurately the entropic program approximates the original program for any particular choice of the penalization parameter. We close this long-standing gap in the literature regarding entropic penalization by giving a new proof of the exponential convergence, valid for any linear program. Our proof is non-asymptotic, yields explicit constants, and has the virtue of being extremely simple. We provide matching lower bounds and show that the entropic approach does not lead to a near-linear time approximation scheme for the linear assignment problem.


title: A tunable quantum dissipator for active resonator reset in circuit QED

abstract: We propose a method for fast, deterministic resonator reset based on tunable dissipative modes. The dissipator is based on a Josephson junction with relatively low quality factor. When the dissipator is tuned into resonance with a high quality microwave resonator, resonator photons are absorbed by the dissipator at a rate orders of magnitude faster than the resonator relaxation rate. We determine the optimal parameters for realization of the tunable dissipator, and examine application of the dissipator to removing spurious photon population in the qubit readout resonator in circuit quantum electrodynamics. We show that even in the nonlinear large photon occupation regime, this enhanced resonator decay rate can be attained by appropriate modulation of the dissipator frequency.


title: Particle Swarm Optimization based search for gravitational waves from compact binary coalescences: performance improvements

abstract: While a fully-coherent all-sky search is known to be optimal for detecting signals from compact binary coalescences (CBCs), its high computational cost has limited current searches to less sensitive coincidence-based schemes. For a network of first generation GW detectors, it has been demonstrated that Particle Swarm Optimization (PSO) can reduce the computational cost of this search, in terms of the number of likelihood evaluations, by a factor of $\approx 10$ compared to a grid-based optimizer. Here, we extend the PSO-based search to a network of second generation detectors and present further substantial improvements in its performance by adopting the local-best variant of PSO and an effective strategy for tuning its configuration parameters. It is shown that a PSO-based search is viable over the entire binary mass range relevant to second generation detectors at realistic signal strengths.


title: Supersymmetry-guided method for mode selection and optimization in coupled systems

abstract: Single-mode operation of coupled systems such as optical-fiber bundles, lattices of photonic waveguides, or laser arrays requires an efficient method to suppress unwanted super-modes. Here, we propose a systematic supersymmetry-based approach to selectively eliminate modes of such systems by decreasing their lifetime relative to the lifetime of the mode of interest. The proposed method allows to explore the opto-geometric parameters of the coupled system and to maximize the relative lifetime of a selected mode. We report a ten-fold increase in the relative lifetime of the fundamental modes of large one-dimensional coupled arrays in comparison to simple 'head-to-tail' coupling geometries. The ability to select multiple supported modes in one- and two-dimensional arrays is also demonstrated.


title: Evolutionary Innovations and Where to Find Them: Routes to Open-Ended Evolution in Natural and Artificial Systems

abstract: This paper presents a high-level conceptual framework to help orient the discussion and implementation of open-endedness in evolutionary systems. Drawing upon earlier work by Banzhaf et al., three different kinds of open-endedness are identified: exploratory, expansive, and transformational. These are characterised in terms of their relationship to the search space of phenotypic behaviours. A formalism is introduced to describe three key processes required for an evolutionary process: the generation of a phenotype from a genetic description, the evaluation of that phenotype, and the reproduction with variation of individuals according to their evaluation. The distinction is made between intrinsic and extrinsic implementations of these processes. A discussion then investigates how various interactions between these processes, and their modes of implementation, can lead to open-endedness. However, an important contribution of the paper is the demonstration that these considerations relate to exploratory open-endedness only. Conditions for the implementation of the more interesting kinds of open-endedness - expansive and transformational - are also discussed, emphasizing factors such as multiple domains of behaviour, transdomain bridges, and non-additive compositional systems. These factors relate not to the generic evolutionary properties of individuals and populations, but rather to the nature of the building blocks out of which individual organisms are constructed, and the laws and properties of the environment in which they exist. The paper ends with suggestions of how the framework can be used to categorise and compare the open-ended evolutionary potential of different systems, how it might guide the design of systems with greater capacity for open-ended evolution, and how it might be further improved.


title: Infinite lattice models by expansion with a non-Gaussian initial approximation

abstract: Recently, a convergent series employing a non-Gaussian initial approximation was constructed and shown to be an effective computational tool for the finite size lattice models with a polynomial interaction. Here we numerically examine the applicability of the convergent series method to models defined on infinite lattices. The comparison of the convergent series computations and the infinite lattice extrapolations of the Monte Carlo simulations reveals an agreement between two approaches.


title: Enabling Cooperative IoT Security via Software Defined Networks (SDN)

abstract: Internet of Things (IoT) is becoming an increasingly attractive target for cybercriminals. We observe that many attacks to IoTs are launched in a collusive way, such as brute-force hacking usernames and passwords, to target at a particular victim. However, most of the time our defending mechanisms to such kind of attacks are carried out individually and independently, which leads to ineffective and weak defense. To this end, we propose to leverage Software Defined Networks (SDN) to enable cooperative security for legacy IP-based IoT devices. SDN decouples control plane and data plane, and can help bridge the knowledge divided between the application and network layers. In this paper, we discuss the IoT security problems and challenges, and present an SDN-based architecture to enable IoT security in a cooperative manner. Furthermore, we implemented a platform that can quickly share the attacking information with peer controllers and block the attacks. We carried out our experiments in both virtual and physical SDN environments with OpenFlow switches. Our evaluation results show that both environments can scale well to handle attacks, but hardware implementation is much more efficient than a virtual one.


title: The Fullerene-like Structure of Superheavy Element Z=120 (Greinerium)-a tribute to Walter Greiner

abstract: We review the basic ideas and some theoretical models behind the concept of fullerene-like structures made of $\alpha$-particles. The possibility of such a peculiar nuclear shape developing in a double magic superheavy nucleus with $Z$=120 was mentioned for the first time in the literature by Walter Greiner and then constantly defended by him. We provide estimates of the energy of such metastable states within the liquid-drop model. In the second part of this paper we discuss a simple model rooted in the nontopological soliton model consisting of a complex scalar field describing a finite system of $\alpha$ bosons with non-linear self-interactions coupled to the electromagnetic field. We demonstrate that this model predicts density depletion in the central region of a soliton-like structure for a range of model parameters.


title: An arithmetic topos for integer matrices

abstract: We study the topos of sets equipped with an action of the monoid of regular $2 \times 2$ matrices over the integers. In particular, we show that the topos-theoretic points are given by the double quotient $\left. GL_2(\hat{\mathbb{Z}}) ~\middle\backslash~ M_2(\mathbb{A}_f)~\middle/~GL_2(\mathbb{Q})\right.$, so they classify the groups $\mathbb{Z}^2 \subseteq A \subseteq \mathbb{Q}^2$ up to isomorphism. We determine the topos automorphisms and then point out the relation with Conway's big picture and the work of Connes and Consani on the Arithmetic Site. As an application to number theory, we show that classifying extensions of $\mathbb{Q}$ by $\mathbb{Z}$ up to isomorphism relates to Goormaghtigh conjecture.


title: High-Dimensional Econometrics and Regularized GMM

abstract: This chapter presents key concepts and theoretical results for analyzing estimation and inference in high-dimensional models. High-dimensional models are characterized by having a number of unknown parameters that is not vanishingly small relative to the sample size. We first present results in a framework where estimators of parameters of interest may be represented directly as approximate means. Within this context, we review fundamental results including high-dimensional central limit theorems, bootstrap approximation of high-dimensional limit distributions, and moderate deviation theory. We also review key concepts underlying inference when many parameters are of interest such as multiple testing with family-wise error rate or false discovery rate control. We then turn to a general high-dimensional minimum distance framework with a special focus on generalized method of moments problems where we present results for estimation and inference about model parameters. The presented results cover a wide array of econometric applications, and we discuss several leading special cases including high-dimensional linear regression and linear instrumental variables models to illustrate the general results.


title: Series representations in spaces of vector-valued functions via Schauder decompositions

abstract: It is a classical result that every $\mathbb{C}$-valued holomorphic function has a local power series representation. This even remains true for holomorphic functions with values in a locally complete locally convex Hausdorff space $E$ over $\mathbb{C}$. Motivated by this example we try to answer the following question. Let $E$ be a locally convex Hausdorff space over a field $\mathbb{K}$, $\mathcal{F}(\Omega)$ be a locally convex Hausdorff space of $\mathbb{K}$-valued functions on a set $\Omega$ and $\mathcal{F}(\Omega,E)$ be an $E$-valued counterpart of $\mathcal{F}(\Omega)$ (where the term $E$-valued counterpart needs clarification itself). For which spaces is it possible to lift series representations of elements of $\mathcal{F}(\Omega)$ to elements of $\mathcal{F}(\Omega,E)$? We derive sufficient conditions for the answer to be affirmative using Schauder decompositions which are applicable for many classical spaces of functions $\mathcal{F}(\Omega)$ having an equicontinuous Schauder basis.


title: Resolutions of nilpotent orbit closures via Coulomb branches of 3-dimensional N=4 theories

abstract: The Coulomb branches of certain 3-dimensional N=4 quiver gauge theories are closures of nilpotent orbits of classical or exceptional algebras. The monopole formula, as Hilbert series of the associated Coulomb branch chiral ring, has been successful in describing the singular hyper-K\"ahler structure. By means of the monopole formula with background charges for flavour symmetries, which realises real mass deformations, we study the resolution properties of all (characteristic) height two nilpotent orbits. As a result, the monopole formula correctly reproduces (i) the existence of a symplectic resolution, (ii) the form of the symplectic resolution, and (iii) the Mukai flops in the case of multiple resolutions. Moreover, the (characteristic) height two nilpotent orbit closures are resolved by cotangent bundles of Hermitian symmetric spaces and the unitary Coulomb branch quiver realisations exhaust all the possibilities.


title: Global analysis of charge exchange meson production at high energies

abstract: Many experiments that are conducted to study the hadron spectrum rely on peripheral resonance production. Hereby, the rapidity gap allows the process to be viewed as an independent fragmen- tation of the beam and the target, with the beam fragmentation dominated by production and decays of meson resonances. We test this separation by determining the kinematic regimes that are dominated by factorizable contributions, indicating the most favorable regions to perform this kind of experiments. In doing so, we use a Regge model to analyze the available world data of charge exchange meson production with beam momentum above 5 GeV in the laboratory frame, that are not dominated by either pion or Pomeron exchanges. We determine the Regge residues and point out the kinematic regimes which are dominated by factorizable contributions.


title: Construction of a Triple-GEM Detector Using Commercially Manufactured Large GEM Foils

abstract: Many experiments are currently using or proposing to use large area GEM foils in their detectors, which is creating a need for commercially available GEM foils. Currently CERN is the only main distributor of large GEM foils, however with the growing interest in GEM technology keeping up with the increasing demand for GEMs will be difficult. We present here an update on the assembly and testing of triple-GEM tracking detectors utilizing single-masked $40 \times 40$ cm$^2$ commercial GEM foils produced by Tech-Etch. The triple-GEM detectors will allow us to characterize the overall quality of these Tech-Etch foils through gain, efficiency, and energy resolution measurements. This will be done by constructing four single-mask triple-GEM detectors, using foils manufactured by Tech-Etch, which follow the design used by the STAR Forward GEM Tracker (FGT). The stack is formed by gluing the foils to the frames and then gluing the frames together. The stack also includes a Tech-Etch produced high voltage foil and a 2D $r-\phi$ readout foil. While one of the four triple-GEM detectors will be built identically to the STAR FGT, the other three will investigate ways in which to further decrease the material budget and increase the efficiency of the detector by incorporating perforated Kapton spacer rings rather than G10 spacing grids to reduce the dead area of the detector.


title: Expanding CO Shells in the Orion A Molecular Cloud

abstract: We present the discovery of expanding spherical shells around low to intermediate-mass young stars in the Orion A giant molecular cloud using observations of $^{12}$CO (1-0) and $^{13}$CO (1-0) from the Nobeyama Radio Observatory 45-meter telescope. The shells have radii from 0.05 to 0.85 pc and expand outward at 0.8 to 5 km/s. The total energy in the expanding shells is comparable to protostellar outflows in the region. Together, shells and outflows inject enough energy and momentum to maintain the cloud turbulence. The mass-loss rates required to power the observed shells are two to three orders of magnitude higher than predicted for line-driven stellar winds from intermediate-mass stars. This discrepancy may be resolved by invoking accretion-driven wind variability. We describe in detail several shells in this paper and present the full sample in the online journal.


title: Large and Persistent Photoconductivity due to Hole-Hole Correlation in CdS

abstract: Large and persistent photoconductivity (LPPC) in semiconductors is due to the trapping of photo-generated minority carriers at crystal defects. Theory has suggested that anion vacancies in II-VI semiconductors are responsible for LPPC due to negative-U behavior, whereby two minority carriers become kinetically trapped by lattice relaxation following photo-excitation. By performing a detailed analysis of photoconductivity in CdS, we provide experimental support for this negative-U model of LPPC. We also show that LPPC is correlated with sulfur deficiency. We use this understanding to vary the photoconductivity of CdS films over nine orders of magnitude, and vary the LPPC characteristic decay time from seconds to 10^4 seconds, by controlling the activities of Cd^{2+} and S^{2-} ions during chemical bath deposition. We suggest a screening method to identify other materials with long-lived, non-equilibrium, photo-excited states based on the results of ground-state calculations of atomic rearrangements following defect redox reactions, with a conceptual connection to polaron formation.


title: Secrecy Outage Analysis of Mixed RF-FSO Downlink SWIPT Systems

abstract: We analyze a secure dual-hop mixed radio frequency-free space optical (RF-FSO) downlink simultaneous wireless information and power transfer (SWIPT) systems. The FSO link and all RF links experience Gamma-Gamma, independent and identical Nakagami-$m$ fading, respectively. We analyze the effects of atmospheric turbulence, pointing error, detection technology, path loss, and energy harvesting on secrecy performance. Signal-to-noise ratios at both legitimate and illegitimate receivers are not independent since they are both simultaneously influenced by the FSO link. We derive the closed-form expression of the secrecy outage probability (SOP) as well as the asymptotic result for SOP when signal-to-noise ratios at relay and legitimate destinations tend to infinity. Monte-Carlo simulations are performed to verify the accuracy of our analysis. The results show that the secrecy diversity order (SDO) depends on the fading parameter of the relay-destination link and the number of the destination's antennas. Additionally, the SDO also depends on the fading parameters, the pointing error parameter, and the detection type of the FSO link.


title: Performance Evaluation of Deep Learning Networks for Semantic Segmentation of Traffic Stereo-Pair Images

abstract: Semantic image segmentation is one the most demanding task, especially for analysis of traffic conditions for self-driving cars. Here the results of application of several deep learning architectures (PSPNet and ICNet) for semantic image segmentation of traffic stereo-pair images are presented. The images from Cityscapes dataset and custom urban images were analyzed as to the segmentation accuracy and image inference time. For the models pre-trained on Cityscapes dataset, the inference time was equal in the limits of standard deviation, but the segmentation accuracy was different for various cities and stereo channels even. The distributions of accuracy (mean intersection over union - mIoU) values for each city and channel are asymmetric, long-tailed, and have many extreme outliers, especially for PSPNet network in comparison to ICNet network. Some statistical properties of these distributions (skewness, kurtosis) allow us to distinguish these two networks and open the question about relations between architecture of deep learning networks and statistical distribution of the predicted results (mIoU here). The results obtained demonstrated the different sensitivity of these networks to: (1) the local street view peculiarities in different cities that should be taken into account during the targeted fine tuning the models before their practical applications, (2) the right and left data channels in stereo-pairs. For both networks, the difference in the predicted results (mIoU here) for the right and left data channels in stereo-pairs is out of the limits of statistical error in relation to mIoU values. It means that the traffic stereo pairs can be effectively used not only for depth calculations (as it is usually used), but also as an additional data channel that can provide much more information about scene objects than simple duplication of the same street view images.


title: Mean dimension and an embedding theorem for real flows

abstract: We develop mean dimension theory for $\mathbb{R}$-flows. We obtain fundamental properties and examples and prove an embedding theorem: Any real flow $(X,\mathbb{R})$ of mean dimension strictly less than $r$ admits an extension $(Y,\mathbb{R})$ whose mean dimension is equal to that of $(X,\mathbb{R})$ and such that $(Y,\mathbb{R})$ can be embedded in the $\mathbb{R}$-shift on the compact function space $\{f\in C(\mathbb{R},[-1,1])|\;\mathrm{supp}(\hat{f})\subset [-r,r]\}$, where $\hat{f}$ is the Fourier transform of $f$ considered as a tempered distribution. These canonical embedding spaces appeared previously as a tool in embedding results for $\mathbb{Z}$-actions.


title: Classical Nonrelativistic Effective Field Theories for a Real Scalar Field

abstract: A classical nonrelativistic effective field theory for a real Lorentz-scalar field $\phi$ is most conveniently formulated in terms of a complex scalar field $\psi$. There have been two derivations of effective Lagrangians for the complex field $\psi$ in which terms in the effective potential were determined to order $(\psi^* \psi)^4$. We point out an error in each of the effective Lagrangians. After correcting the errors, we demonstrate the equivalence of the two effective Lagrangians by verifying that they both reproduce $T$-matrix elements of the relativistic real scalar field theory and by also constructing a redefinition of the complex field $\psi$ that transforms terms in one effective Lagrangian into the corresponding terms of the other.


title: MRPC: An R package for accurate inference of causal graphs

abstract: We present MRPC, an R package that learns causal graphs with improved accuracy over existing packages, such as pcalg and bnlearn. Our algorithm builds on the powerful PC algorithm, the canonical algorithm in computer science for learning directed acyclic graphs. The improvement in accuracy results from online control of the false discovery rate (FDR) that reduces false positive edges, a more accurate approach to identifying v-structures (i.e., $T_1 \rightarrow T_2 \leftarrow T_3$), and robust estimation of the correlation matrix among nodes. For genomic data that contain genotypes and gene expression for each sample, MRPC incorporates the principle of Mendelian randomization to orient the edges. Our package can be applied to continuous and discrete data.


title: Shear viscosity and wall slip behavior of dense suspensions of polydisperse particles

abstract: Shear viscosity and wall slip of dense suspensions of a silicone polymer incorporated with polydisperse particles were investigated. Three types of particles with low aspect ratios were used to achieve a relatively high maximum packing fraction of 0.86. Such a high maximum packing fraction allowed the preparation of suspensions with a wide range of solid volume fractions in the range of 0.62 to 0.82. The wall slip velocities of the suspensions in steady torsional and capillary flows were characterized and determined to be fully consistent with the mechanism of apparent slip layer formation at the wall. Upon wall slip corrections it was found that at shear stresses which are significantly above the yield stress the relative shear viscosity of the suspensions obeys well the Krieger-Dougherty relationship that links the relative shear viscosity behavior of dense suspensions solely to the ratio of the volume fraction over the maximum packing fraction of solids. However, at lower shear stresses that are in the vicinity of the yield stresses the relative shear viscosity becomes functions of both the ratio of the volume fraction over the maximum packing fraction and the shear stress. It is clearly demonstrated that without wall slip analysis the accurate characterization of the relative shear viscosity of dense suspensions is not possible.


title: Soft Charges and Electric-Magnetic Duality

abstract: The main focus of this work is to study magnetic soft charges of the four dimensional Maxwell theory. Imposing appropriate asymptotic falloff conditions, we compute the electric and magnetic soft charges and their algebra both at spatial and at null infinity. While the commutator of two electric or two magnetic soft charges vanish, the electric and magnetic soft charges satisfy a complex $U(1)$ current algebra. This current algebra through Sugawara construction yields two $U(1)$ Kac-Moody algebras. We repeat the charge analysis in the electric-magnetic duality-symmetric Maxwell theory and construct the duality-symmetric phase space where the electric and magnetic soft charges generate the respective boundary gauge transformations. We show that the generator of the electric-magnetic duality and the electric and magnetic soft charges form infinite copies of $iso(2)$ algebra. Moreover, we study the algebra of charges associated with the global Poincar\'e symmetry of the background Minkowski spacetime and the soft charges. We discuss physical meaning and implication of our charges and their algebra.


title: Swimming of peritrichous bacteria is enabled by an elastohydrodynamic instability

abstract: Peritrichously-flagellated bacteria, such as Escherichia coli, self-propel in fluids by using specialised motors to rotate multiple helical filaments. The rotation of each motor is transmitted to a short flexible segment called the hook which in turn transmits it to a flagellar filament, enabling swimming of the whole cell. Since multiple motors are spatially distributed on the body of the organism, one would expect the propulsive forces from the filaments to push against each other leading to negligible swimming. We use a combination of computations and theory to show that the swimming of multi-flagellated bacteria is enabled by an elastohydrodynamic bending instability occurring for hooks more flexible than a critical threshold. Using past measurements of hook bending stiffness, we demonstrate how the design of real bacteria allows them to be safely on the side of this instability that promotes systematic swimming.


title: Regular Rotating MOG Dark Compact Object

abstract: A regular rotating MOG (modified gravity) compact object is derived that reduces to the Kerr black hole when the parameter $\alpha=0$. Physical consequences of the dark compact object, which is regular everywhere in spacetime for $\alpha > \alpha_{\rm crit}=0.674$ and is a rotating Kerr-MOG black hole for $\alpha < \alpha_{\rm crit}$ are investigated.


title: Heavy Flavor Azimuthal Correlations in Cold Nuclear Matter

abstract: Background: It has been proposed that the azimuthal distributions of heavy flavor quark-antiquark pairs may be modified in the medium of a heavy-ion collision. Purpose: This work tests this proposition through next-to-leading order (NLO) calculations of the azimuthal distribution, $d\sigma/d\phi$, including transverse momentum broadening, employing $<k_T^2>$ and fragmentation in exclusive $Q \bar Q$ pair production. While these studies were done for $p+p$, $p + \bar p$ and $p+$Pb collisions, understanding azimuthal angle correlations between heavy quarks in these smaller, colder systems is important for their interpretation in heavy-ion collisions. Methods: First, single inclusive $p_T$ distributions calculated with the exclusive HVQMNR code are compared to those calculated in the fixed-order next-to-leading logarithm approach. Next the azimuthal distributions are calculated and sensitivities to $<k_T^2>$, $p_T$ cut, and rapidity are studied at $\sqrt{s} = 7$ TeV. Finally, calculations are compared to $Q \bar Q$ data in elementary $p+p$ and $p + \bar p$ collisions at $\sqrt{s} = 7$ TeV and 1.96 TeV as well as to the nuclear modification factor $R_{p {\rm Pb}}(p_T)$ in $p+$Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV measured by ALICE. Results: The low $p_T$ ($p_T < 10$ GeV) azimuthal distributions are very sensitive to the $k_T$ broadening and rather insensitive to the fragmentation function. The NLO contributions can result in an enhancement at $\phi \sim 0$ absent any other effects. Agreement with the data was found to be good. Conclusions: The NLO calculations, assuming collinear factorization and introducing $k_T$ broadening, result in significant modifications of the azimuthal distribution at low $p_T$ which must be taken into account in calculations of these distributions in heavy-ion collisions.


title: A New Theoretical Evaluation Framework for Satisfaction Equilibria in Wireless Networks

abstract: In this paper, a theoretical evaluation framework regarding the \textit{Satisfaction Equilibrium (SE)} in wireless communication networks is introduced and examined. To study these equilibria operation points, we coin some new concepts, namely the \textit{Valued Satisfaction Equilibrium}, the \textit{Price of Efficiency} and the \textit{Max Price of Satisfaction}, which can be used for measuring the efficiency of the obtained equilibria solutions. The aforementioned framework is analyzed and evaluated in a wireless communication environment under the presence of the Gaussian Interference channel (GIC). Within this setting, a non-cooperative game among the users is studied, where users aim in a selfish manner to meet their Quality of Service (QoS) prerequisite. However instead of maximizing the QoS which is generally energy costly, we evangelize that better energy-efficiency is achieved by targeting satisfactory QoS levels only. The sufficient and necessary conditions that lead to the \textit{Satisfaction Equilibrium} are provided for the two-user case and the \textit{Efficient Satisfaction Equilibrium (ESE)} is determined, where the users satisfy their QoS constraints with the lowest possible cost. Moreover, specific measures for evaluating the efficiency of various satisfaction equilibria, in a formal and quantitative manner, expressing the tradeoff with respect to the achieved utility or a given objective function and corresponding cost, are defined and analyzed.


title: Achieving Data Dissemination with Security using FIWARE and Intel Software Guard Extensions (SGX)

abstract: The Internet of Things (IoT) field has gained much attention from industry and academia, being the main subject for numerous research and development projects. Frequently, the dense amount of generated data from IoT applications is sent to a cloud service, that is responsible for processing and storage. Many of these applications demand security and privacy for their data because of their sensitive nature. This is specially true when such data must be processed in entities hosted in public clouds, where the environment in which applications run may not be trusted. Some concerns are then raised since it is not trivial to provide the needed protection for these sensitive data. We present a solution that considers the security components of FIWARE and the Intel SGX capabilities. FIWARE is a platform created to support the development of Smart Applications, including IoT systems, and SGX is the Intel solution for Trusted Execution Environment (TEE). We propose a new component for key management that, together with other FIWARE components, can be used to provide privacy, confidentiality, and integrity guarantees for IoT data. A case study illustrates how this proposed solution can be employed in a realistic scenario, which allows the dissemination of sensitive data through public clouds without risking privacy issues. The results of the experiments provide evidence that our approach does not harm scalability or availability of the system. In addition, it presents acceptable memory costs when considering the benefit of the privacy guarantees achieved.


title: Y-Net: A deep Convolutional Neural Network for Polyp Detection

abstract: Colorectal polyps are important precursors to colon cancer, the third most common cause of cancer mortality for both men and women. It is a disease where early detection is of crucial importance. Colonoscopy is commonly used for early detection of cancer and precancerous pathology. It is a demanding procedure requiring significant amount of time from specialized physicians and nurses, in addition to a significant miss-rates of polyps by specialists. Automated polyp detection in colonoscopy videos has been demonstrated to be a promising way to handle this problem. {However, polyps detection is a challenging problem due to the availability of limited amount of training data and large appearance variations of polyps. To handle this problem, we propose a novel deep learning method Y-Net that consists of two encoder networks with a decoder network. Our proposed Y-Net method} relies on efficient use of pre-trained and un-trained models with novel sum-skip-concatenation operations. Each of the encoders are trained with encoder specific learning rate along the decoder. Compared with the previous methods employing hand-crafted features or 2-D/3-D convolutional neural network, our approach outperforms state-of-the-art methods for polyp detection with 7.3% F1-score and 13% recall improvement.


title: Cohomological localization for transverse Lie algebra actions on Riemannian foliations

abstract: We prove localization and integration formulas for the equivariant basic cohomology of Riemannian foliations. As a corollary we obtain a Duistermaat-Heckman theorem for transversely symplectic foliations.


title: Tilted shear-free axially symmetric fluids

abstract: We carry on a systematic study of the physical properties of axially symmetric fluid distributions, which appear to be geodesic, shear--free, irrotational, non--dissipative and purely electric, for the comoving congruence of observers, from the point of view of the tilted congruence. The vanishing of the magnetic part of the Weyl tensor for the comoving congruence of observers, suggests that no gravitational radiation is produced during the evolution of the system. Instead, the magnetic part of the Weyl tensor as measured by tilted observers is non vanishing (as well as the shear, the four--acceleration, the vorticity and the dissipation), giving rise to a flux of gravitational radiation that can be characterized through the super--Poynting vector. This result strengthens further the relevance of the role of observers in the description of a physical system. An explanation of this dual interpretation in terms of the information theory, is provided.


title: Probabilistic Deep Learning using Random Sum-Product Networks

abstract: The need for consistent treatment of uncertainty has recently triggered increased interest in probabilistic deep learning methods. However, most current approaches have severe limitations when it comes to inference, since many of these models do not even permit to evaluate exact data likelihoods. Sum-product networks (SPNs), on the other hand, are an excellent architecture in that regard, as they allow to efficiently evaluate likelihoods, as well as arbitrary marginalization and conditioning tasks. Nevertheless, SPNs have not been fully explored as serious deep learning models, likely due to their special structural requirements, which complicate learning. In this paper, we make a drastic simplification and use random SPN structures which are trained in a "classical deep learning manner", i.e. employing automatic differentiation, SGD, and GPU support. The resulting models, called RAT-SPNs, yield prediction results comparable to deep neural networks, while still being interpretable as generative model and maintaining well-calibrated uncertainties. This property makes them highly robust under missing input features and enables them to naturally detect outliers and peculiar samples.


title: Adversarial Scene Editing: Automatic Object Removal from Weak Supervision

abstract: While great progress has been made recently in automatic image manipulation, it has been limited to object centric images like faces or structured scene datasets. In this work, we take a step towards general scene-level image editing by developing an automatic interaction-free object removal model. Our model learns to find and remove objects from general scene images using image-level labels and unpaired data in a generative adversarial network (GAN) framework. We achieve this with two key contributions: a two-stage editor architecture consisting of a mask generator and image in-painter that co-operate to remove objects, and a novel GAN based prior for the mask generator that allows us to flexibly incorporate knowledge about object shapes. We experimentally show on two datasets that our method effectively removes a wide variety of objects using weak supervision only


title: A note on domination in intersecting linear systems

abstract: A linear system is a pair $(P,\mathcal{L})$ where $\mathcal{L}$ is a family of subsets on a ground finite set $P$ such that $|l\cap l^\prime|\leq 1$, for every $l,l^\prime \in \mathcal{L}$. The elements of $P$ and $\mathcal{L}$ are called points and lines, respectively, and the linear system is called intersecting if any pair of lines intersect in exactly one point. A subset $D$ of points of a linear system $(P,\mathcal{L})$ is a dominating set of $(P,\mathcal{L})$ if for every $u\in P\setminus D$ there exists $v\in D$ such that $u,v\in l$, for some $l\in\mathcal{L}$. The cardinality of a minimum dominating set of a linear system $(P,\mathcal{L})$ is called domination number of $(P,\mathcal{L})$, denoted by $\gamma(P,\mathcal{L})$. On the other hand, a subset $R$ of lines of a linear system $(P,\mathcal{L})$ is a 2-packing if any three elements of $R$ have not a common point (are triplewise disjoint). The cardinality of a maximum 2-packing of a linear system $(P,\mathcal{L})$ is called 2-packing number of $(P,\mathcal{L})$, denoted by $\nu_2(P,\mathcal{L})$. It is know for intersecting linear systems $(P,\mathcal{L})$ of rank $r$ it satisfies $\gamma(P,\mathcal{L})\leq r-1$. In this note we prove, if $q$ is an even prime power and $(P,\mathcal{L})$ is an intersecting linear system of rank $q+2$ satisfying $\gamma(P,\mathcal{L})=q+1$, then this linear system can be constructed from a spanning $(q+1)$-uniform intersecting linear subsystem $(P^\prime,\mathcal{L}^\prime)$ of the projective plane of order $q$ satisfying $\tau(P^\prime,\mathcal{L}^\prime)=\nu_2(P^\prime,\mathcal{L}^\prime)-1=q+1$.


title: The WIYN One Degree Imager in 2018: An Extended 30-Detector Focal Plane

abstract: We report on the upgraded One Degree Imager (ODI) at the WIYN 3.5 meter telescope at the Kitt Peak Observatory after the focal plane was expanded by an additional seventeen detectors in spring 2015. The now thirty Orthogonal Transfer Array CCD detectors provide a total field of view of 40' x 48' on the sky. The newly added detectors underwent a design revision to mitigate reduced charge transfer efficiency under low light conditions. We discuss the performance of the focal plane and challenges in the photometric calibration of the wide field of view, helped by the addition of telescope baffles. In a parallel project, we upgraded the instrument's three filter arm mechanisms, where a degrading worm-gear mechanism was replaced by a chain drive that is operating faster and with high reliability. Three more filters, a u' band and two narrow band filters were added to the instrument's complement, with two additional narrow band filters currently in procurement (including an H-alpha filter). We review the lessons learned during nearly three years of operating the instrument in the observatory environment and discuss infrastructure upgrades that were driven by ODI's needs.


title: Error scaling of large-eddy simulation in the outer region of wall-bounded turbulence

abstract: We study the error scaling properties of large-eddy simulation (LES) in the outer region of wall-bounded turbulence at moderately high Reynolds numbers. In order to avoid the additional complexity of wall-modeling, we perform LES of turbulent channel flows in which the no-slip condition at the wall is replaced by a Neumann condition supplying the exact mean wall-stress. The statistics investigated are the mean velocity profile, turbulence intensities, and kinetic energy spectra. The errors follow $(\Delta/L)^{\alpha}Re_\tau^{\gamma}$, where $\Delta$ is the characteristic grid resolution, $Re_\tau$ is the friction Reynolds number, and $L$ is the meaningful length-scale to normalize $\Delta$ in order to collapse the errors across the wall-normal distance. We show that $\Delta$ can be expressed as the $L_2$-norm of the grid vector and that $L$ is well represented by the ratio of the friction velocity and mean shear. The exponent $\alpha$ is estimated from theoretical arguments for each statistical quantity of interest and shown to roughly match the values computed by numerical simulations. For the mean profile and kinetic energy spectra, $\alpha\approx1$, whereas the turbulence intensities converge at a slower rate $\alpha<1$. The exponent $\gamma$ is approximately $0$, i.e. the LES solution is independent of the Reynolds number. The expected behavior of the turbulence intensities at high Reynolds numbers is also derived and shown to agree with the classic log-layer profiles for grid resolutions lying within the inertial range. Further examination of the LES turbulence intensities and spectra reveals that both quantities resemble their filtered counterparts from direct numerical simulation (DNS) data, but that the mechanism responsible for this similarity is related to the balance between the input power and dissipation rather than to filtering.


title: Parameter estimation and treatment optimization in a stochastic model for immunotherapy of cancer

abstract: Adoptive Cell Transfer therapy of cancer is currently in full development and mathematical modeling is playing a critical role in this area. We study a stochastic model developed by Baar et al. in 2015 for modeling immunotherapy against melanoma skin cancer. First, we estimate the parameters of the deterministic limit of the model based on biological data of tumor growth in mice. A Nonlinear Mixed Effects Model is estimated by the Stochastic Approximation Expectation Maximization algorithm. With the estimated parameters, we head back to the stochastic model and calculate the probability that the T cells all get exhausted during the treatment. We show that for some relevant parameter values, an early relapse is due to stochastic fluctuations (complete T cells exhaustion) with a non negligible probability. Then, focusing on the relapse related to the T cell exhaustion, we propose to optimize the treatment plan (treatment doses and restimulation times) by minimizing the T cell exhaustion probability in the parameter estimation ranges.


title: Selecting Reduced Models in the Cross-Entropy Method

abstract: This paper deals with the estimation of rare event probabilities using importance sampling (IS), where an optimal proposal distribution is computed with the cross-entropy (CE) method. Although, IS optimized with the CE method leads to an efficient reduction of the estimator variance, this approach remains unaffordable for problems where the repeated evaluation of the score function represents a too intensive computational effort. This is often the case for score functions related to the solution of a partial differential equation (PDE) with random inputs. This work proposes to alleviate computation by the parsimonious use of a hierarchy of score function approximations in the CE optimization process. The score function approximation is obtained by selecting the surrogate of lowest dimensionality, whose accuracy guarantees to pass the current CE optimization stage. The selection relies on certified upper bounds on the error norm. An asymptotic analysis provides some theoretical guarantees on the efficiency and convergence of the proposed algorithm. Numerical results demonstrate the gain brought by the method in the context of pollution alerts and a system modeled by a PDE.


title: Transverse localization of transmission eigenchannels

abstract: Transmission eigenchannels are building blocks of coherent wave transport in diffusive media, and selective excitation of individual eigenchannels can lead to diverse transport behavior. An essential yet poorly understood property is the transverse spatial profile of each eigenchannel, which is critical for coupling into and out of it. Here, we discover that the transmission eigenchannels of a disordered slab possess localized incident and outgoing profiles, even in the diffusive regime far from Anderson localization. Such transverse localization arises from a combination of reciprocity, local coupling of spatial modes, and nonlocal correlations of scattered waves. Experimentally, we observe signatures of such localization despite finite illumination area. Our results reveal the intrinsic characteristics of transmission eigenchannels in the open slab geometry, commonly used for applications in imaging and energy transfer through turbid media.


title: A Framework for the construction of upper bounds on the number of affine linear regions of ReLU feed-forward neural networks

abstract: We present a framework to derive upper bounds on the number of regions that feed-forward neural networks with ReLU activation functions are affine linear on. It is based on an inductive analysis that keeps track of the number of such regions per dimensionality of their images within the layers. More precisely, the information about the number regions per dimensionality is pushed through the layers starting with one region of the input dimension of the neural network and using a recursion based on an analysis of how many regions per output dimensionality a subsequent layer with a certain width can induce on an input region with a given dimensionality. The final bound on the number of regions depends on the number and widths of the layers of the neural network and on some additional parameters that were used for the recursion. It is stated in terms of the $L1$-norm of the last column of a product of matrices and provides a unifying treatment of several previously known bounds: Depending on the choice of the recursion parameters that determine these matrices, it is possible to obtain the bounds from Mont\'{u}far (2014), (2017) and Serra et. al. (2017) as special cases. For the latter, which is the strongest of these bounds, the formulation in terms of matrices provides new insight. In particular, by using explicit formulas for a Jordan-like decomposition of the involved matrices, we achieve new tighter results for the asymptotic setting, where the number of layers of the same fixed width tends to infinity.


title: Intuitive Analyses via Drift Theory

abstract: Drift theory is an intuitive tool for reasoning about random processes: It allows turning expected stepwise changes into expected first-hitting times. While drift theory is used extensively by the community studying randomized search heuristics, it has seen hardly any applications outside of this field, in spite of many research questions that can be formulated as first-hitting times. We state the most useful drift theorems and demonstrate their use for various randomized processes, including the coupon collector process, winning streaks, approximating vertex cover, and a random sorting algorithm. We also consider processes without expected stepwise change and give theorems based on drift theory applicable in such scenarios. We use these theorems for the analysis of the gambler's ruin process, for a coloring algorithm, for an algorithm for 2-SAT, and for a version of the Moran process without bias. A final tool we present is a tight theorem for processes on finite state spaces, which we apply to the Moran process. We aim to enable the reader to apply drift theory in their own research to derive accessible proofs and to teach it as a simple tool for the analysis of random processes.


title: The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: a tomographic measurement of structure growth and expansion rate from anisotropic galaxy clustering in Fourier space

abstract: We perform a tomographic structure growth and expansion rate analysis using the monopole, quadrupole and hexadecapole of the redshift-space galaxy power spectrum derived from the Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 combined sample, which covers the redshift range of $0.20<z<0.75$. By allowing for overlap between neighbouring redshift slices in order to extract information on the light-cone, we successfully obtain joint BAO and RSD constraints with a precision of $2-3\%$ for $D_A$, $3-10\%$ for $H$ and $9-12\%$ for $f\sigma_8$ with a redshift resolution of $\Delta z\sim0.04$. Our measurement is consistent with that presented in arXiv:1709.05173, where the analysis is performed in configuration space. We apply our measurement to constrain the $f(R)$ gravity model, and find that the 95\% CL upper limit of ${\rm log_{10}}B_0$ can be reduced by 11\% by our tomographic BAO and RSD measurement.


title: Least gradient problem with respect to a non-strictly convex norm

abstract: We study the planar least gradient problem with respect to an anisotropic norm $\phi$ for continuous boundary data. We prove existence of minimizers for strictly convex domains $\Omega$. Furthermore, we inspect the issue of uniqueness and regularity of minimizers only in terms of the modes of convexity of $\phi$ and $\Omega$. The results are independent from the regularity of $\phi$.


title: Some analysis on a fractional differential equation involving a noncontinuous right-hand side

abstract: By developing new techniques we establish local existence and uniqueness theorems for an initial value problem involving a nonlinear equation in the sense of Riemann-Liouville fractional derivative in the case that the nonlinear function on the right hand side of the equation is not continuous on $[0,T]\times\mathbb{R}.$


title: Cuscuton kinks and branes

abstract: In this paper, we study a peculiar model for the scalar field. We add the cuscuton term in a standard model and investigate how this inclusion modifies the usual behavior of kinks. We find the first order equations and calculate the energy density and the total energy of the system. Also, we investigate the linear stability of the model, which is governed by a Sturm-Liouville eigenvalue equation that can be transformed in an equation of the Shcr\"odinger type. The model is also investigated in the braneworld scenario, where a first order formalism is also obtained and the linear stability is investigated.


title: Dark Markets with Multiple Assets: Segmentation, Asymptotic Stability, and Equilibrium Prices

abstract: We study a generalization of the model of a dark market due to Duffie-G\^arleanu- Pedersen [6]. Our market is segmented and involves multiple assets. We show that this market has a unique asymptotically stable equilibrium. In order to establish this result, we use a novel approach inspired by a theory due to McKenzie and Hawkins-Simon. Moreover, we obtain a closed form solution for the price of each asset at which investors trade at equilibrium. We conduct a comparative statics analysis which shows, among other sensitivities, how equilibrium prices respond to the level of interactions between investors.


title: Non-Topological Majorana Zero Modes in Inhomogeneous Spin Ladders

abstract: We show that the coupling of homogeneous Heisenberg spin-1/2 ladders in different phases leads to the formation of interfacial zero energy Majorana bound states. Unlike Majorana bound states at the interfaces of topological quantum wires, these states are void of topological protection and generally susceptible to local perturbations of the host spin system. However, a key message of our work is that in practice they show a high degree of resilience over wide parameter ranges which may make them interesting candidates for applications.


title: Deployment of a Blockchain-Based Self-Sovereign Identity

abstract: Digital identity is unsolved: after many years of research there is still no trusted communication over the Internet. To provide identity within the context of mutual distrust, this paper presents a blockchain-based digital identity solution. Without depending upon a single trusted third party, the proposed solution achieves passport-level legally valid identity. This solution for making identities Self-Sovereign, builds on a generic provable claim model for which attestations of truth from third parties need to be collected. The claim model is then shown to be both blockchain structure and proof method agnostic. Four different implementations in support of these two claim model properties are shown to offer sub-second performance for claim creation and claim verification. Through the properties of Self-Sovereign Identity, legally valid status and acceptable performance, our solution is considered to be fit for adoption by the general public.


title: Solitons and peaked solitons for the general Degasperis-Procesi model

abstract: We consider the general Degasperis-Procesi model of shallow water out-flows. This six parametric family of conservation laws contains, in particular, KdV, Benjamin-Bona-Mahony, Camassa-Holm, and Degasperis-Procesi equations. The main result consists of criterions which guarantee the existence of solitary wave solutions: solitons and peakons ("peaked solitons").


title: A unified study of nonparametric inference for monotone functions

abstract: The problem of nonparametric inference on a monotone function has been extensively studied in many particular cases. Estimators considered have often been of so-called Grenander type, being representable as the left derivative of the greatest convex minorant or least concave majorant of an estimator of a primitive function. In this paper, we provide general conditions for consistency and pointwise convergence in distribution of a class of generalized Grenander-type estimators of a monotone function. This broad class allows the minorization or majoratization operation to be performed on a data-dependent transformation of the domain, possibly yielding benefits in practice. Additionally, we provide simpler conditions and more concrete distributional theory in the important case that the primitive estimator and data-dependent transformation function are asymptotically linear. We use our general results in the context of various well-studied problems, and show that we readily recover classical results established separately in each case. More importantly, we show that our results allow us to tackle more challenging problems involving parameters for which the use of flexible learning strategies appears necessary. In particular, we study inference on monotone density and hazard functions using informatively right-censored data, extending the classical work on independent censoring, and on a covariate-marginalized conditional mean function, extending the classical work on monotone regression functions. In addition to a theoretical study, we present numerical evidence supporting our large-sample results.


title: Exact solvability of two new 3D and 1D nonrelativistic potentials within the TRA framework

abstract: This work is concerned about introducing two new 1D and 3D confined potentials and present their solutions using the Tridiagonal Representation Approach (TRA). The wavefunction is written as a series in terms of square integrable basis functions which are expressed in terms of Jacobi polynomials. Moreover, the expansion coefficients are written in terms of new orthogonal polynomials that were introduced recently by Alhaidari, the analytical properties of these polynomials are yet to be derived. Moreover, we have computed the numerical eigen-energies for both potentials by considering specific choices of the potential parameters.


title: On Elo based prediction models for the FIFA Worldcup 2018

abstract: We propose an approach for the analysis and prediction of a football championship. It is based on Poisson regression models that include the Elo points of the teams as covariates and incorporates differences of team-specific effects. These models for the prediction of the FIFA World Cup 2018 are fitted on all football games on neutral ground of the participating teams since 2010. Based on the model estimates for single matches Monte-Carlo simulations are used to estimate probabilities for reaching the different stages in the FIFA World Cup 2018 for all teams. We propose two score functions for ordinal random variables that serve together with the rank probability score for the validation of our models with the results of the FIFA World Cups 2010 and 2014. All models favor Germany as the new FIFA World Champion. All possible courses of the tournament and their probabilities are visualized using a single Sankey diagram.


title: Order Preserving Interpolation for Summation-by-Parts Operators at Non-Conforming Grid Interfaces

abstract: We study non-conforming grid interfaces for summation-by-parts finite difference methods applied to partial differential equations with second derivatives in space. To maintain energy stability, previous efforts have been forced to accept a reduction of the global convergence rate by one order, due to large truncation errors at the non-conforming interface. We avoid the order reduction by generalizing the interface treatment and introducing order preserving interpolation operators. We prove that, given two diagonal-norm summation-by-parts schemes, order preserving interpolation operators with the necessary properties are guaranteed to exist, regardless of the grid-point distributions along the interface. The new methods retain the stability and global accuracy properties of the underlying schemes for conforming interfaces.


title: Experimental observation of symmetry protected bound state in the continuum in a chain of dielectric disks

abstract: Existence of bound states in the continuum (BIC) manifests a general wave phenomenon firstly predicted in quantum mechanics in 1929 by J. von Neumann and E. Wigner. Today it is being actively explored in photonics, radiophysics, acoustics, and hydrodynamics. In this paper, we report the first experimental observation of electromagnetic bound state in the radiation continuum in 1D array of dielectric particles. By measurement of the transmission spectra of the ceramic disk chain at GHz frequencies we demonstrate how a resonant state in the vicinity of the center of the Brillouin zone turns into a symmetry-protected BIC with increase the number of the disks. We estimate a number of the disks when the radiation losses become negligible in comparison to material absorption and, therefore, the chain could be considered practically as infinite. The presented analysis is supplemented by measurements of the near fields of the symmetry protected BIC. All measurements are in a good agreement with the results of numerical simulation and analytical model based on tight-binding approximation. The obtained results provide useful guidelines for practical implementations of structures with bound states in the continuum that opens up new horizons for the development of optical and radiofrequency metadevices.


title: Explainable Neural Networks based on Additive Index Models

abstract: Machine Learning algorithms are increasingly being used in recent years due to their flexibility in model fitting and increased predictive performance. However, the complexity of the models makes them hard for the data analyst to interpret the results and explain them without additional tools. This has led to much research in developing various approaches to understand the model behavior. In this paper, we present the Explainable Neural Network (xNN), a structured neural network designed especially to learn interpretable features. Unlike fully connected neural networks, the features engineered by the xNN can be extracted from the network in a relatively straightforward manner and the results displayed. With appropriate regularization, the xNN provides a parsimonious explanation of the relationship between the features and the output. We illustrate this interpretable feature--engineering property on simulated examples.


title: Global-in-time classical solutions and qualitative properties for the NNLIF neuron model with synaptic delay

abstract: The Nonlinear Noisy Leaky Integrate and Fire (NNLIF) model is widely used to describe the dynamics of neural networks after a diffusive approximation of the mean-field limit of a stochastic differential equation system. When the total activity of the network has an instantaneous effect on the network, in the average-excitatory case, a blow-up phenomenon occurs. This article is devoted to the theoretical study of the NNLIF model in the case where a delay in the effect of the total activity on the neurons is added. We first prove global-in-time existence and uniqueness of classical solutions, independently of the sign of the connectivity parameter, that is, for both cases: excitatory and inhibitory. Secondly, we prove some qualitative properties of solutions: asymptotic convergence to the stationary state for weak interconnections and a non-existence result for periodic solutions if the connectivity parameter is large enough. The proofs are mainly based on an appropriate change of variables to rewrite the NNLIF equation as a Stefan-like free boundary problem, constructions of universal super-solutions, the entropy dissipation method and Poincar\'e's inequality.


title: Exploring Feature Reuse in DenseNet Architectures

abstract: Densely Connected Convolutional Networks (DenseNets) have been shown to achieve state-of-the-art results on image classification tasks while using fewer parameters and computation than competing methods. Since each layer in this architecture has full access to the feature maps of all previous layers, the network is freed from the burden of having to relearn previously useful features, thus alleviating issues with vanishing gradients. In this work we explore the question: To what extent is it necessary to connect to all previous layers in order to reap the benefits of feature reuse? To this end, we introduce the notion of local dense connectivity and present evidence that less connectivity, allowing for increased growth rate at a fixed network capacity, can achieve a more efficient reuse of features and lead to higher accuracy in dense architectures.


title: Selection and Estimation Optimality in High Dimensions with the TWIN Penalty

abstract: We introduce a novel class of variable selection penalties called TWIN, which provides sensible data-adaptive penalization. Under a linear sparsity regime and random Gaussian designs we show that penalties in the TWIN class have a high probability of selecting the correct model and furthermore result in minimax optimal estimators. The general shape of penalty functions in the TWIN class is the key ingredient to its desirable properties and results in improved theoretical and empirical performance over existing penalties. In this work we introduce two examples of TWIN penalties that admit simple and efficient coordinate descent algorithms, making TWIN practical in large data settings. We demonstrate in challenging and realistic simulation settings with high correlations between active and inactive variables that TWIN has high power in variable selection while controlling the number of false discoveries, outperforming standard penalties.


title: Post-inflationary thermal histories and the refractive index of relic gravitons

abstract: We investigate the impact of the post-inflationary thermal histories on the cosmic graviton spectrum caused by the inflationary variation of their refractive index. Depending on the frequency band, the spectral energy distribution can be mildly red, blue or even violet. Wide portions of the parameter space lead to potentially relevant signals both in the audio range (probed by the advanced generation of terrestrial interferometers) and in the mHz band (where space-borne detectors could be operational within the incoming score year). The description of the refractive index in conformally related frames is clarified.


title: Fusion systems with Benson-Solomon components

abstract: The Benson-Solomon systems comprise a one-parameter family of simple exotic fusion systems at the prime $2$. The results we prove give significant additional evidence that these are the only simple exotic $2$-fusion systems, as conjectured by Solomon. We consider a saturated fusion system $\mathcal{F}$ having an involution centralizer with a component $\mathcal{C}$ isomorphic to a Benson-Solomon fusion system, and we show under rather general hypotheses that $\mathcal{F}$ cannot be simple. Furthermore, we prove that if $\mathcal{F}$ is almost simple with these properties, then $\mathcal{F}$ is isomorphic to the next larger Benson-Solomon system extended by a group of field automorphisms. Our results are situated within Aschbacher's program to provide a new proof of a major part of the classification of finite simple groups via fusion systems. One of the most important steps in this program is a proof of Walter's Theorem for fusion systems, and our first result is specifically tailored for use in the proof of that step. We then apply Walter's Theorem to treat the general Benson-Solomon component problem under the assumption that each component of an involution centralizer in $\mathcal{F}$ is on the list of currently known quasisimple $2$-fusion systems.


title: Stacks in Poisson Geometry

abstract: This thesis is divided into four chapters. The first chapter discusses the relationship between stacks on a site and groupoids internal to the site. It includes a rigorous proof of the folklore result that there is an equivalence between the bicategory of internal groupoids and the bicategory of geometric stacks. The second chapter discusses standard concepts in the theory of geometric stacks, including Morita equivalence, stack symmetries, and some Morita invariants. The third chapter introduces a new site of Dirac structures and provides a rigorous answer to the question: What is the stack associated to a symplectic groupoid? The last chapter discusses a remarkable class of Poisson manifolds, called b-symplectic manifolds, giving a classification of them up to Morita equivalence and computing their Picard group.


title: EIGEN: Ecologically-Inspired GENetic Approach for Neural Network Structure Searching from Scratch

abstract: Designing the structure of neural networks is considered one of the most challenging tasks in deep learning, especially when there is few prior knowledge about the task domain. In this paper, we propose an Ecologically-Inspired GENetic (EIGEN) approach that uses the concept of succession, extinction, mimicry, and gene duplication to search neural network structure from scratch with poorly initialized simple network and few constraints forced during the evolution, as we assume no prior knowledge about the task domain. Specifically, we first use primary succession to rapidly evolve a population of poorly initialized neural network structures into a more diverse population, followed by a secondary succession stage for fine-grained searching based on the networks from the primary succession. Extinction is applied in both stages to reduce computational cost. Mimicry is employed during the entire evolution process to help the inferior networks imitate the behavior of a superior network and gene duplication is utilized to duplicate the learned blocks of novel structures, both of which help to find better network structures. Experimental results show that our proposed approach can achieve similar or better performance compared to the existing genetic approaches with dramatically reduced computation cost. For example, the network discovered by our approach on CIFAR-100 dataset achieves 78.1% test accuracy under 120 GPU hours, compared to 77.0% test accuracy in more than 65, 536 GPU hours in [35].


title: The MESS of cosmological perturbations

abstract: We introduce two new effective quantities for the study of comoving curvature perturbations $\zeta$: the space dependent effective sound speed (SESS) and the momentum dependent effective sound speed (MESS) . We use the SESS and the MESS to derive a new set of equations which can be applied to any system described by an effective stress-energy-momentum tensor (EST), including multi-fields systems, supergravity and modified gravity theories. We show that this approach is completely equivalent to the standard one and it has the advantage of requiring to solve only one differential equation for $\zeta$ instead of a system, without the need of explicitly computing the evolution of entropy perturbations. The equations are valid for perturbations respect to any arbitrary flat spatially homogeneous background, including any inflationary and bounce model. As an application we derive the equation for $\zeta$ for multi-fields $KGB$ models and show that observed features of the primordial curvature perturbation spectrum are compatible with the effects of an appropriate local variation of the MESS in momentum space. The MESS is the natural quantity to parametrize in a model independent way the effects produced on curvature perturbations by multi-fields systems, particle production and modified gravity theories and could be conveniently used in the analysis of LSS observations, such as the ones from the upcoming EUCLID mission or CMB radiation measurements.


title: Solar wind dynamics around a comet: The paradigmatic inverse-square-law model

abstract: Aims. Observations of solar protons near comet 67P/Churyumov-Gerasimenko (67P) by the Rosetta spacecraft can be modelled by the planar motion in an effective magnetic field proportional to 1/r^2. We aim to provide a thorough study of such dynamics, with a clear description of the behaviour of an incoming flux of particles. We will be able, then, to calibrate the free parameters of the model to Rosetta observations. Methods. Basic tools of dynamical analysis are used. They lead to a definition of the relevant parameters for the system and a classification of the possible types of trajectories. Using the so-obtained formalism, the structures formed by a flux of particles coming from infinity can be studied. Results. All the trajectories are parametrised by two characteristic radii, rE and rC, derived from first integrals. There are three different types of motion possible divided by a separatrix corresponding to rE = rC. An analytical expression of the trajectories, defined by an integral, is developed. Using this formalism, the application to a flux of particles coming from infinity (modelling the incident solar wind) gives one free parameter only, the radius rE, which scales the problem. A circular cavity of radius 0.28 rE is created, as well as an overdensity curve (analogous to a caustic in optics). At each observation time, rE can be calibrated to Rosetta plasma measurements, giving a qualitative understanding of the solar particle dynamics (incoming direction, cavity and density map). We also deduce that, in order to properly capture the essence of the dynamics, numerical simulations of the solar wind around a comet must use simulation boxes much larger than rE and grids much finer than rE.


title: Cloning Algorithms: from Large Deviations to Population Dynamics

abstract: Population dynamics provides a numerical tool allowing for the study of rare events by means of simulating a large number of copies of the system, supplemented with a selection rule that favours the rare trajectories of interest. The cloning algorithm allows the estimation of a large deviation function (LDF) of additive observables in Markov processes. However, such algorithms are plagued by finite simulation time $t$ and finite population size $N_c$ effects that can render their use delicate. First, using a non-constant population approach, we analyze the small-$N_c$ effects in the initial transient regime. These effects play an important role in the numerical determination of LDF. We show how to overcome these effects by introducing a time delay in the evolution of populations, additional to the discarding of the initial regime of the population growth where these discreteness effects are strong. Then, the study of the finite-$t$ and finite-$N_c$ scalings in the LDF evaluation is done using two different versions of the algorithm, in discrete and continuous-time. We show that these scalings behave as $1/N_c$ and $1/t$ in the large-$N_c$ and large-$t$ asymptotics respectively. Moreover, we show that one can make use of this convergence speed in order to extract the asymptotic behavior in the infinite-$t$ and infinite-$N_c$ limits resulting in a better LDF estimation. These scalings are later generalized and evidence of a breakdown for large-size systems is presented.


title: Disrupted orbital order and the pseudo-gap in layered 1T-TaS$_2$

abstract: We present a state-of-the-art density functional theory (DFT) study which models crucial features of the partially disordered orbital order stacking in the prototypical layered transition metal dichalcogenide 1T-TaS2 . Our results not only show that DFT models with realistic assumptions about the orbital order perpendicular to the layers yield band structures which agree remarkably well with experiments. They also demonstrate that DFT correctly predicts the formation of an excitation pseudo-gap which is commonly attributed to Mott-Hubbard type electron-electron correlations. These results highlight the importance of interlayer interactions in layered transition metal dichalcogenides and serve as an intriguing example of how disorder within an electronic crystal can give rise to pseudo-gap features.


title: Spectrally and Power Efficient Optical Communication Systems

abstract: Increased traffic demands globally and in particular in short-reach links in data centers will require optical communication systems to continue scaling at an accelerated pace. Nevertheless, energy constraints start to limit the bit rate that can be practically transmitted over optical systems both at the shortest distances in data centers and at the longest distances in ultra-long submarine links. Short-reach links in data centers face strict constraints on power consumption, size, and cost, which will demand low-power solutions that scale to bit rates beyond 100 Gbit/s per wavelength, while accommodating increased losses due to longer fiber plant, multiplexing of more wavelengths, and possibly optical switching. At the longest distances, submarine optical cables longer than about 5,000 km face energy constraints due to power feed limits at the shores, which restricts the electrical power available to the undersea optical amplifiers, ultimately limiting the optical power and throughput per fiber. This dissertation presents strategies to mitigate these limitations. For short-reach links in data centers, we propose low-power coherent receiver architectures that completely avoid high-speed analog-to-digital converters and digital signal processors. For long-haul submarine links, we demonstrate how optimizing the channel power allocation under a constraint on the amplifier power maximizes the information-theoretic capacity per fiber.


title: Learning to Understand Goal Specifications by Modelling Reward

abstract: Recent work has shown that deep reinforcement-learning agents can learn to follow language-like instructions from infrequent environment rewards. However, this places on environment designers the onus of designing language-conditional reward functions which may not be easily or tractably implemented as the complexity of the environment and the language scales. To overcome this limitation, we present a framework within which instruction-conditional RL agents are trained using rewards obtained not from the environment, but from reward models which are jointly trained from expert examples. As reward models improve, they learn to accurately reward agents for completing tasks for environment configurations---and for instructions---not present amongst the expert data. This framework effectively separates the representation of what instructions require from how they can be executed. In a simple grid world, it enables an agent to learn a range of commands requiring interaction with blocks and understanding of spatial relations and underspecified abstract arrangements. We further show the method allows our agent to adapt to changes in the environment without requiring new expert examples.


title: A linear time method for the detection of point and collective anomalies

abstract: The challenge of efficiently identifying anomalies in data sequences is an important statistical problem that now arises in many applications. Whilst there has been substantial work aimed at making statistical analyses robust to outliers, or point anomalies, there has been much less work on detecting anomalous segments, or collective anomalies, particularly in those settings where point anomalies might also occur. In this article, we introduce Collective And Point Anomalies (CAPA), a computationally efficient approach that is suitable when collective anomalies are characterised by either a change in mean, variance, or both, and distinguishes them from point anomalies. Theoretical results establish the consistency of CAPA at detecting collective anomalies and, as a by-product, the consistency of a popular penalised cost based change in mean and variance detection method. Empirical results show that CAPA has close to linear computational cost as well as being more accurate at detecting and locating collective anomalies than other approaches. We demonstrate the utility of CAPA through its ability to detect exoplanets from light curve data from the Kepler telescope.


title: The estimation of impact ionization coefficients for $\beta $-Ga$_{2}$O$_{3}$

abstract: Impact ionization coefficients of anisotropic monoclinic $\beta $-Ga$_{2}$O$_{3}$ are estimated along four crystallographic directions and the plot for the $\left[ 010 \right]$ direction is shown. The approximation models were fitted to Baraff's universal plot for ionization rate in semiconductors and the values were obtained for $\beta $-Ga$_{2}$O$_{3}$. The phonon mean free path of $\beta $-Ga$_{2}$O$_{3}$ was estimated to be 5.2604 \si{\angstrom} using Gray medium approximation. The phonon group velocity takes the value of longitudinal acoustic phonons. The ionization rate has a maximum value of $ 3.98\times \,10^{6}\,cm^{-1}$ along the $\left[ 010 \right]$ direction over the applied electric field range $(1.43-4)\times10^{7}\,V\,\cdot\,cm^{-1}$. Contrary to expectations, the phonon mean free path along direction $\left[ \bar{2}01 \right]$ is the lowest since it has a lower thermal conductivity to phonon group velocity ratio. The plots were compared with GaN and 4H-SiC, which shows that as bandgap increases the field required for ionization increases. The critical electric field was estimated to be $ 0.921 \times10^{8}\,V\cdot\,cm^{-1}$ along $\left[010 \right]$ direction.


title: Reduced-Order Modeling through Machine Learning Approaches for Brittle Fracture Applications

abstract: In this paper, five different approaches for reduced-order modeling of brittle fracture in geomaterials, specifically concrete, are presented and compared. Four of the five methods rely on machine learning (ML) algorithms to approximate important aspects of the brittle fracture problem. In addition to the ML algorithms, each method incorporates different physics-based assumptions in order to reduce the computational complexity while maintaining the physics as much as possible. This work specifically focuses on using the ML approaches to model a 2D concrete sample under low strain rate pure tensile loading conditions with 20 preexisting cracks present. A high-fidelity finite element-discrete element model is used to both produce a training dataset of 150 simulations and an additional 35 simulations for validation. Results from the ML approaches are directly compared against the results from the high-fidelity model. Strengths and weaknesses of each approach are discussed and the most important conclusion is that a combination of physics-informed and data-driven features are necessary for emulating the physics of crack propagation, interaction and coalescence. All of the models presented here have runtimes that are orders of magnitude faster than the original high-fidelity model and pave the path for developing accurate reduced order models that could be used to inform larger length-scale models with important sub-scale physics that often cannot be accounted for due to computational cost.


title: Spin-unrestricted self-energy embedding theory

abstract: We present a new theoretical approach, unrestricted self-energy embedding theory (USEET) that is a Green's function embedding theory used to study problems in which an open, embedded system exchanges electrons with the environment. USEET has a high potential to be used in studies of strongly correlated systems with odd number of electrons and open shell systems such as transition metal complexes important in inorganic chemistry. In this paper, we show that USEET results agree very well with common quantum chemistry methods while avoiding typical bottlenecks present in these method.


title: A resolved warm/dense gas Schmidt-Kennicutt relationship in a binary HyLIRG at $z=2.41$

abstract: Hyperluminous infrared galaxies (HyLIRGs) provide crucial "long lever arm" constraints on galaxy evolution. H-ATLAS $J084933.4+021443$, a $z=2.41$ binary HyLIRG with at least two additional luminous companion galaxies, is thus an optimal test-ground for studies of star formation and galaxy evolution during "cosmic noon". We have used ALMA to obtain resolved imaging and kinematics of atomic and molecular emission lines, and rest-frame $340$ to $1160$GHz continuum emission, for the known luminous component galaxies in H-ATLAS $J084933.4+021443$: W, T, M, C. All four component galaxies are spatially ($\sim 0 .\!\!^{''} 3$ or $2.5$kpc) resolved in CO J:7-6, [C$_\mathrm{I}$] 2-1, H$_2$O and the millimetre (mm) to sub-mm continuum. Rotation-dominated gas kinematics is confirmed in W and T. The significant extension to component T, in gas and continuum, along its kinematic minor axis, is attributable to its lensing magnification. Spatially resolved sub-mm spectral energy distributions reveal that component W is well fit with greybody emission from dust at a single temperature over the full extent of the galaxy, despite it containing a powerful AGN, while component T requires an additional component of hotter nuclear dust and additional sources of emission in the mm. We confirm that [C$_\mathrm{I}$] 2-1 can be used as a rough tracer of warm/dense molecular gas in extreme systems, though the [C$_\mathrm{I}$] 2-1/CO luminosity ratio increases sub-linearly. We obtain an exquisite and unprecedented resolved ($2.5$-kpc-scale) "warm/dense molecular gas" Schmidt-Kennicutt (SK) relationship for components W and T. Gas exhaustion times for all apertures in W (T) are $1-4$Gyr ($0.5-2$Gyr). Both W and T follow a resolved "warm/dense gas" SK relationship with power law $n\sim1.7$, significantly steeper than the $n\sim1$ found previously via "cold" molecular gas in nearby "normal" star-forming galaxies.


title: Ultrastrongly dissipative quantum Rabi model

abstract: We discuss the equilibrium and out of equilibrium dynamics of cavity QED in presence of dissipation beyond the standard perturbative treatment of losses. Using the dynamical polaron \emph{ansatz} and Matrix Product State simulations, we discuss the case where both light-matter $g$-coupling and system-bath interaction are in the ultrastrong coupling regime. We provide a critical $g$ for the onset of Rabi oscillations. Besides, we demonstrate that the qubit is \emph{dressed} by the cavity and dissipation. That such dressing governs the dynamics and, thus, it can be measured. Finally, we sketch an implementation for our theoretical ideas within circuit QED technology.