title: Visualizing Invisible Dark Matter Annihilation with the CMB and Matter Power Spectrum

abstract: We study the cosmological signatures of Invisibly Annihilating Dark Matter (IAnDM), where DM annihilates into dark radiation particles that are decoupled from the Standard Model (SM). In the simple benchmark model we consider here, such invisible annihilation determines the relic abundance of DM via dark thermal freeze-out. We demonstrate that IAnDM may reveal itself through observable, novel signatures that are correlated: scale-dependent $\Delta N_\text{eff}$ (number of extra effective neutrinos) in the Cosmic Microwave Background (CMB) spectrum due to DM residual annihilation, while the phase of acoustic peaks shifts towards the opposite direction relative to that due to SM neutrinos, resembling the effect due to scattering (fluidlike) thermal dark radiation; in addition, IAnDM induces modifications to the matter power spectrum that resemble those due to warm dark matter. Current data are sensitive to IAnDM with masses up to $\sim200$ keV, while future observations will improve this reach, especially if the late-time DM annihilation cross section is enhanced relative to the standard thermal value, which can be realized in a variety of scenarios. This work also demonstrates a new possibility of realizing thermal sub-MeV DM with observable signals.


title: Particle Physics and the Cosmic Microwave Background

abstract: Temperature and polarization variations across the microwave sky include the fingerprints of quantum fluctuations in the early universe. They may soon reveal physics at unprecedented energy scales.


title: A nonlinear and time-dependent visco-elasto-plastic rheology model for studying shock-physics phenomena

abstract: We present a simple and efficient implementation of a viscous creep rheology based on diffusion creep, dislocation creep and the Peierls mechanism in conjunction with an elasto-plastic rheology model into a shock-physics code, the iSALE open-source impact code. Our approach is based on the calculation of an effective viscosity which is then used as a reference viscosity for any underlying viscoelastic (or even visco-elasto-plastic) model. Here we use a Maxwell-model which best describes stress relaxation and is therefore likely most important for the formation of large meteorite impact basins. While common viscoelastic behavior during mantle convection or other slow geodynamic or geological processes is mostly controlled by diffusion and dislocation creep, we showed that the Peierls mechanism dominates at the large strain rates that typically occur during meteorite impacts. Thus, the resulting visco-elasto-plastic rheology allows implementation of a more realistic mantle behavior in computer simulations, especially for those dealing with large meteorite impacts. The approach shown here opens the way to more faithful simulations of large impact basin formation, especially in elucidating the physics behind the formation of the external fault rings characteristic of large lunar basins.


title: Unsupervised Machine Learning Based on Non-Negative Tensor Factorization for Analyzing Reactive-Mixing

abstract: Analysis of reactive-diffusion simulations requires a large number of independent model runs. For each high-fidelity simulation, inputs are varied and the predicted mixing behavior is represented by changes in species concentration. It is then required to discern how the model inputs impact the mixing process. This task is challenging and typically involves interpretation of large model outputs. However, the task can be automated and substantially simplified by applying Machine Learning (ML) methods. In this paper, we present an application of an unsupervised ML method (called NTFk) using Non-negative Tensor Factorization (NTF) coupled with a custom clustering procedure based on k-means to reveal hidden features in product concentration. An attractive aspect of the proposed ML method is that it ensures the extracted features are non-negative, which are important to obtain a meaningful deconstruction of the mixing processes. The ML method is applied to a large set of high-resolution FEM simulations representing reaction-diffusion processes in perturbed vortex-based velocity fields. The applied FEM ensures that species concentration are always non-negative. The simulated reaction is a fast irreversible bimolecular reaction. The reactive-diffusion model input parameters that control mixing include properties of velocity field, anisotropic dispersion, and molecular diffusion. We demonstrate the applicability of the ML method to produce a meaningful deconstruction of model outputs to discriminate between different physical processes impacting the reactants, their mixing, and the spatial distribution of the product. The presented ML analysis allowed us to identify additive features that characterize mixing behavior.


title: Towards a Theory of Additive Eigenvectors

abstract: The standard approach in solving stochastic equations is eigenvector decomposition. Using separation ansatz $P(i,t)=u(i)e^{\mu t}$ one obtains standard equation for eigenvectors $Ku=\mu u$, where $K$ is the rate matrix of the master equation. While universally accepted, the standard approach is not the only possibility. Using additive separation ansatz $S(i,t)=W(i)-\nu t$ one arrives at additive eigenvectors. Here we suggest a theory of such eigenvectors. We argue that additive eigenvectors describe conditioned Markov processes and derive corresponding equations. The formalism is applied to one-dimensional stochastic process corresponding to the telegraph equation. We derive differential equations for additive eigenvectors and explore their properties. The proposed theory of additive eigenvectors provides a new description of stochastic processes with peculiar properties.


title: Weak shock propagation with accretion. I. Self-similar solutions and application to failed supernovae

abstract: We present solutions for the self-similar propagation of a shock wave in a hydrostatic, adiabatic medium with a point mass gravitational field. In contrast to the well-known, Sedov-Taylor blastwave, these solutions apply to the case when the shock Mach number is of order a few, and the energy of the shocked fluid is not conserved but self-consistently modified by the binding energy of the ambient medium that is swept up by the passage of the shock. Furthermore, we show that there is one solution (for a given ambient density profile) that smoothly passes through a sonic point in the post-shock flow and results in accretion onto the central object; in analogy with the Bondi problem, we propose that these solutions are the ones that are most relevant in astrophysical environments. We apply these accreting models to failed supernovae, in which neutron star formation does not unbind the envelope, but a weak shock is still generated in the outer layers of the star from neutrino-induced mass loss. We find excellent agreement between the predictions of our self-similar, shock propagation model and numerical simulations of the collapse of a yellow supergiant; the self-similar solutions reproduce the overall scaling of the shock speed, the time and space-dependent evolution of the velocity, density, and pressure behind the shock, and the accretion rate onto the black hole. Our results have important implications for the fallback and ejection of material in failed supernovae.


title: The Millimeter Continuum Size-Frequency Relationship in the UZ Tau E Disk

abstract: We present high spatial resolution observations of the continuum emission from the young multiple star system UZ Tau at frequencies from 6 to 340 GHz. To quantify the spatial variation of dust emission in the UZ Tau E circumbinary disk, the observed interferometric visibilities are modeled with a simple parametric prescription for the radial surface brightnesses at each frequency. We find evidence that the spectrum steepens with radius in the disk, manifested as a positive correlation between the observing frequency and the radius that encircles a fixed fraction of the emission ($R_{eff} \propto \nu^{0.34 \pm 0.08}$). The origins of this size--frequency relation are explored in the context of a theoretical framework for the growth and migration of disk solids. While that framework can reproduce a similar size--frequency relation, it predicts a steeper spectrum than is observed. Moreover, it comes closest to matching the data only on timescales much shorter ($\le 1$ Myr) than the putative UZ Tau age (~2-3 Myr). These discrepancies are the direct consequences of the rapid radial drift rates predicted by models of dust evolution in a smooth gas disk. One way to mitigate that efficiency problem is to invoke small-scale gas pressure modulations that locally concentrate drifting solids. If such particle traps reach high continuum optical depths at 30-340 GHz with a ~30-60% filling fraction in the inner disk ($r \lesssim20$ au), they can also explain the observed spatial gradient in the UZ Tau E disk spectrum.


title: Order in the chaos: eccentric black hole binary mergers in triples formed via strong binary-binary scatterings

abstract: Black hole (BH) triples represent one of the astrophysical pathways for BH mergers in the Universe detectable by LIGO and VIRGO. We study the formation of BH triples via binary-binary encounters in dense clusters, showing that one-third of triples formed through this channel are in a non-hierarchical, unstable configuration. We build a database of $32,000$ $N$-body simulations to investigate the evolution of these BH triples. Varying the mutual orbital inclination, the three BH masses and the inner and outer eccentricities, we show that retrograde, nearly planar configurations lead to a significant shrinkage of the inner binary. We find an universal trend of triple systems, namely that they tend to evolve toward prograde configurations and that the orbital flip, driven by the torque exerted on the inner BH binary (BHB) by the outer BH, leads in general to tighter inner orbits that, in some cases, merge within a Hubble time. We find that the mergers' mass distribution maps the original BH binary spectrum. A large fraction of merging BHBs with initial separation $1$ AU enter the $10^{-3}-10^{-1}$ Hz frequency band with large eccentricities, thus representing potential eccentric LISA sources. Mergers originating from initially tighter BHB ($a\sim 0.01$ AU), instead, often have eccentricities above 0.7 in the $1$ Hz band. This might imply that eccentric sources detected in high-frequency detectors are most likely connected with a high-velocity dispersion stellar environment, whereas eccentric sources detected in low-frequency detectors are likely to develop in low-density clusters. (Abridged)


title: On the orientation of Narrow Line Seyfert I

abstract: We study a sample of Narrow-line Seyfert 1 galaxies (NLS1) in their optical and radio features, in order to understand the differences between their radio silent, radio-loud and radio-quiet subclasses. We first show that the different redshift and mass distributions of radio-loud and -quiet NLS1s could be ascribed to observational biases. We then present a geometrical model according to which most of the different observational features of radio-loud and radio-quiet NLS1s are ascribed to the orientation of an intrinsically structured jet. We estimate the fraction of intrinsically jetted sources among NLS1s that justifies the observed radio-detected population. Noticeably, under the assumptions of the geometrical model, we derive a fraction of jetted sources significantly larger than in standard AGN.


title: Four-loop splitting functions at small $x$

abstract: We consider the expansion of small-$x$ resummed DGLAP splitting functions at next-to-leading logarithmic (NLL) accuracy to four-loop order, namely next-to-next-to-next-to-leading order (N$^3$LO). From this, we extract the exact LL and NLL small-$x$ contributions to the yet unknown N$^3$LO splitting functions, both in the standard $\overline{MS}$ scheme and in the $Q_0 \overline{MS}$ scheme usually considered in small-$x$ literature. We show that the impact of unknown subleading logarithmic contributions (NNLL and beyond) at N$^3$LO is significant, thus motivating future work towards their computation. Our results will be also needed in future to match NLL resummation to N$^3$LO evolution. In turn, we propose an improved implementation of the small-$x$ resummation and therefore release a new version of the resummation code (HELL 3.0) which contains these changes.


title: Simulations of Jet Heating in Galaxy Clusters: Successes and Challenges

abstract: We study how jets driven by active galactic nuclei influence the cooling flow in Perseus-like galaxy cluster cores with idealised, non-relativistic, hydrodynamical simulations performed with the Eulerian code ATHENA using high-resolution Godunov methods with low numerical diffusion. We use novel analysis methods to measure the cooling rate, the heating rate associated to multiple mechanisms, and the power associated with adiabatic compression/expansion. A significant reduction of the cooling rate and cooling flow within 20 kpc from the centre can be achieved with kinetic jets. However, at larger scales and away from the jet axis, the system relaxes to a cooling flow configuration. Jet feedback is anisotropic and is mostly distributed along the jet axis, where the cooling rate is reduced and a significant fraction of the jet power is converted into kinetic power of heated outflowing gas. Away from the jet axis weak shock heating represents the dominant heating source. Turbulent heating is significant only near the cluster centre, but it becomes inefficient at 50 kpc scales where it only represents a few percent of the total heating rate. Several details of the simulations depend on the choice made for the hydro solver, a consequence of the difficulty of achieving proper numerical convergence for this problem: current physics implementations and resolutions do not properly capture multi-phase gas that develops as a consequence of thermal instability. These processes happen at the grid scale and leave numerical solutions sensitive to the properties of the chosen hydro solver.


title: Why is the HLR theory particle-hole symmetric?

abstract: Long wavelength descriptions of a half-filled lowest Landau level ($\nu = 1/2$) must be consistent with the experimental observation of particle-hole (PH) symmetry. The traditional description of the $\nu=1/2$ state pioneered by Halperin, Lee and Read (HLR) naively appears to break PH symmetry. However, recent studies have shown that the HLR theory with weak quenched disorder can exhibit an emergent PH symmetry. We find that such inhomogeneous configurations of the $\nu=1/2$ fluid, when described by HLR mean-field theory, are tuned to a topological phase transition between an integer quantum Hall state and an insulator of composite fermions with a dc Hall conductivity $\sigma_{xy}^{\rm (cf)} = - {1 \over 2} {e^2 \over h}$. Our observations help explain why the HLR theory exhibits PH symmetric dc response.


title: Holographic dual of hot Polchinski-Strassler quark-gluon plasma

abstract: We construct the supergravity dual of the hot quark-gluon plasma in the mass-deformed ${\cal N}=4$ Super-Yang-Mills theory (also known as ${\cal N}=1^*$). The full ten-dimensional type IIB holographic dual is described by 20 functions of two variables, which we determine numerically, and it contains a black hole with $S^5$ horizon topology. As we lower the temperature to around half of the mass of the chiral multiplets, we find evidence for (most likely a first-order) phase transition, which could lead either to one of the Polchinski-Strassler confining, screening, or oblique vacua with polarized branes, or to an intermediate phase corresponding to blackened polarized branes with an $S^2 \times S^3$ horizon topology. This phase transition is a feature that could in principle be seen by putting the theory on the lattice, and thus our result for the ratio of the chiral multiplet mass to the phase transition temperature, $m_c/T = 2.15467491205(6)$, constitutes the first prediction of string theory and AdS/CFT that could be independently checked via four-dimensional super-QCD lattice computation. We also construct the black-hole solution in certain five-dimensional gauged supergravity truncations and, without directly using uplift/reduction formulae, we find strong evidence that the five- and ten-dimensional solutions are the same. This indicates that five-dimensional gauged supergravity is powerful enough to capture the physics of the high-temperature deconfined phase of the Polchinski-Strassler quark-gluon plasma.


title: On the universality of late-time correlators in semi-classical 2d CFTs

abstract: In the framework of AdS$_3$/ CFT$_2$ correspondence, we present a systematic analysis of the late time thermalization of a two dimensional CFT state created by insertion of small number of heavy operators on the vacuum. We show that at late Lorentzian time, the universal features of this thermalization are solely captured by the eigenvalues of the monodromy matrix corresponding to the solutions of the uniformization equation. We discuss two different ways to extract the monodromy eigenvalues while bypassing the need for finding explicitly the full monodromy matrix - first, using a monodromy preserving diffeomorphism and second using Chen-Simons formulation of gravity in AdS$_3$. Both of the methods yield the same precise relation between the eigenvalues and the final black hole temperature at late Lorentzian time.


title: Color Unifed Dynamical Axion

abstract: We consider an enlarged color sector which solves the strong CP problem via new massless fermions. The spontaneous breaking of a unified color group into QCD and another confining group provides a source of naturally large axion mass $m_a$ due to small size instantons. This extra source of axion mass respects automatically the alignment of the vacuum, ensuring a low-energy CP-conserving vacuum. The mechanism does not appeal to a $Z_2$ 'mirror' copy of the SM, nor does it require any fine-tuning of the axion-related couplings at the unification scale. There is no light axion and uncharacteristically the lighter spectrum contains instead sterile fermions. The axion scale $f_a$ can be naturally brought down to a few TeV, with an exotic spectrum of colored pseudoscalars lighter than this scale, observable at colliders exclusively via strong interactions. The $\{m_a, f_a\}$ parameter space which allows a solution of the strong CP problem is thus enlarged well beyond that of invisible axion models.


title: Stellar-mass black holes in young massive and open stellar clusters and their role in gravitational-wave generation III: dissecting black hole dynamics

abstract: Stellar-remnant black holes (BH) in dense stellar clusters comprise a natural setup to trigger general-relativistic (GR) inspiral and merger of binary black holes (BBH), detectable by the LISA and the LIGO-Virgo, through dynamical encounters inside such environments. In this work, the intricacies of such dynamical interactions are probed utilizing realistic, self-consistent, post-Newtonian, direct N-body evolutionary models of young massive and open stellar clusters. Particularly, the configurations of the compact subsystems, that drive the in-cluster GR BBH coalescences, are tracked on the fly. Such an approach reveals that the GR coalescences within the open clusters take place primarily via chaotic interactions involving triple BH systems. Although less frequently, such mergers are found to happen also in higher-order subsystems such as quadruples and in subsystems involving non-BH members; the mergers can themselves be BH---non-BH, which events would leave electromagnetic signatures. Close, fly-by encounters inside the clusters can also make BBHs and other types of double-compact binaries temporarily post-Newtonian; such binaries would potentially contribute to the GW background for the LISA and the PTA. These calculations, furthermore, suggest that open clusters are potential hosts for not only detached BH---main-sequence binaries, as recently identified in the globular cluster NGC 3201, but also a wide variety of other types of remnant---non-remnant binaries, which are assembled via dynamical interactions inside the clusters and which have the prospects of being discovered in radial-velocity surveys.


title: Top Down Approach to 6D SCFTs

abstract: Six-dimensional superconformal field theories (6D SCFTs) occupy a central place in the study of quantum field theories encountered in high energy theory. This article reviews the top down construction and study of this rich class of quantum field theories, in particular, how they are realized by suitable backgrounds in string / M- / F-theory. We review the recent F-theoretic classification of 6D SCFTs, explain how to calculate physical quantities of interest such as the anomaly polynomial of 6D SCFTs, and also explain recent progress in understanding renormalization group flows for deformations of such theories. Additional topics covered by this review include some discussion on the (weighted and signed) counting of states in these theories via superconformal indices. We also include several previously unpublished results as well as a new variant on the swampland conjecture for general quantum field theories decoupled from gravity. The aim of the article is to provide a point of entry into this growing literature rather than an exhaustive overview.


title: Unscrambling Structured Chirality with Structured Light at Nanoscale Using Photo-induced Force

abstract: We introduce a microscopy technique that facilitates the prediction of spatial features of chirality of nanoscale samples by exploiting photo-induced optical force exerted on an achiral tip in the vicinity of the test specimen. The tip-sample interactive system is illuminated by structured light to probe both the transverse and longitudinal (with respect to the beam propagation direction) components of the sample magnetoelectric polarizability as the manifestation of its sense of handedness, i.e., chirality. We specifically prove that although circularly polarized waves are adequate to detect the transverse polarizability components of the sample, they are unable to probe the longitudinal component. To overcome this inadequacy, we propose a judiciously engineered combination of radially and azimuthally polarized beams, as optical vortices possessing pure longitudinal electric and magnetic field components along their vortex axis, respectively, hence probing longitudinal chirality. The proposed technique may benefit branches of science like stereochemistry, biomedicine, physical and material science, and pharmaceutics


title: Attempts at vacuum counterexamples to cosmic censorship in AdS

abstract: We consider vacuum solutions of four dimensional general relativity with $\Lambda < 0$. We numerically construct stationary solutions that asymptotically approach a boundary metric with differential rotation. Smooth solutions only exist up to a critical rotation. We thus argue that increasing the differential rotation by a finite amount will cause the curvature to grow without bound. This holds for both zero and nonzero temperature, and both compact and noncompact boundaries. However, the boundary metric always develops an ergoregion before reaching the critical rotation, which probably means that the energy is unbounded from below for these counterexamples to cosmic censorship.


title: Topological Band Engineering of Graphene Nanoribbons

abstract: Topological insulators (TIs) are an emerging class of materials that host highly robust in-gap surface/interface states while maintaining an insulating bulk. While most notable scientific advancements in this field have been focused on TIs and related topological crystalline insulators in 2D and 3D, more recent theoretical work has predicted the existence of 1D symmetry-protected topological phases in graphene nanoribbons (GNRs). The topological phase of these laterally-confined, semiconducting strips of graphene is determined by their width, edge shape, and the terminating unit cell, and is characterized by a Z2 invariant (similar to 1D solitonic systems). Interfaces between topologically distinct GNRs characterized by different Z2 are predicted to support half-filled in-gap localized electronic states which can, in principle, be utilized as a tool for material engineering. Here we present the rational design and experimental realization of a topologically-engineered GNR superlattice that hosts a 1D array of such states, thus generating otherwise inaccessible electronic structure. This strategy also enables new end states to be engineered directly into the termini of the 1D GNR superlattice. Atomically-precise topological GNR superlattices were synthesized from molecular precursors on a Au(111) surface under ultra-high vacuum (UHV) conditions and characterized by low temperature scanning tunneling microscopy (STM) and spectroscopy (STS). Our experimental results and first-principles calculations reveal that the frontier band structure of these GNR superlattices is defined purely by the coupling between adjacent topological interface states. This novel manifestation of 1D topological phases presents an entirely new route to band engineering in 1D materials based on precise control of their electronic topology, and is a promising platform for future studies of 1D quantum spin physics.


title: Plasma Effects on Lasing of Uniform Ultralight Axion Condensate

abstract: Lasing of ultralight axion condensate into photons can be sensitive to the presence of a background plasma owing to its coupling to electromagnetism. Such a scenario is particularly relevant for superradiant axion condensate around stellar mass black holes since the axion mass can be within a few orders of magnitude of the plasma frequency of the surrounding medium. In this paper I discuss the properties of the plasma around a black hole and analyze its effects on the lasing of a uniform axion condensate of mass of the order of the plasma frequency.


title: Survival probability in Generalized Rosenzweig-Porter random matrix ensemble

abstract: We study analytically and numerically the dynamics of the generalized Rosenzweig-Porter model, which is known to possess three distinct phases: ergodic, multifractal and localized phases. Our focus is on the survival probability $R(t)$, the probability of finding the initial state after time $t$. In particular, if the system is initially prepared in a highly-excited non-stationary state (wave packet) confined in space and containing a fixed fraction of all eigenstates, we show that $R(t)$ can be used as a dynamical indicator to distinguish these three phases. Three main aspects are identified in different phases. The ergodic phase is characterized by the standard power-law decay of $R(t)$ with periodic oscillations in time, surviving in the thermodynamic limit, with frequency equals to the energy bandwidth of the wave packet. In multifractal extended phase the survival probability shows an exponential decay but the decay rate vanishes in the thermodynamic limit in a non-trivial manner determined by the fractal dimension of wave functions. Localized phase is characterized by the saturation value of $R(t\to\infty)=k$, finite in the thermodynamic limit $N\rightarrow\infty$, which approaches $k=R(t\to 0)$ in this limit.


title: Nine tiny star clusters in Gaia DR1, PS1 and DES

abstract: We present the results of a systematic Milky Way satellite search performed across an array of publicly available wide-area photometric surveys. Our aim is to complement previous searches by widening the parameter space covered. Specifically, we focus on objects smaller than $1'$ and include old, young, metal poor and metal rich stellar population masks. As a result we find 9 new likely genuine stellar systems in data from GAIA, DES, and Pan-STARRS, which were picked from the candidate list because of conspicuous counterparts in the cut-out images. The presented systems are all very compact ($r_h<1'$) and faint ($M_V\gtrsim-3$), and are associated either with the Galactic disk, or the Magellanic Clouds. While most of the stellar systems look like Open Clusters, their exact classification is, as of today, unclear. With these discoveries, we extend the parameter space occupied by star clusters to sizes and luminosities previously unexplored and demonstrate that rather than two distinct classes of Globular and Open clusters, there appears to be a continuity of objects, unmarked by a clear decision boundary.


title: Understanding Federation: An Analytical Framework for the Interoperability of Social Networking Sites

abstract: Although social networking has become a remarkable feature in the Web, full interoperability has not arrived. This work explores the main 5 paradigms of interoperability across social networking sites, corresponding to the layers in which we an find interoperability. Building on those, a novel analytical framework for SNS interoperability is introduced. Seven representative interoperability SNS technologies are compared using the proposed framework. The analysis exposes an overwhelming disparity and fragmentation in the solutions for tackling the same problems. Although there are a few solutions where consensus is reached and are widely adopted (e.g. in object IDs), there are multiple central issues that are still far from being widely standarized (e.g. in profile representation). In addition, several areas have been identified where there is clear room for improvement, such as privacy controls or data synchronization.


title: Impact of Type II Spicules in the Corona: Simulations and Synthetic Observables

abstract: The role of type II spicules in the corona has been a much debated topic in recent years. This paper aims to shed light on the impact of type II spicules in the corona using novel 2.5D radiative MHD simulations including ion-neutral interaction effects with the Bifrost code. We find that the formation of simulated type II spicules, driven by the release of magnetic tension, impacts the corona in various manners. Associated with the formation of spicules, the corona exhibits 1) magneto-acoustic shocks and flows which supply mass to coronal loops, and 2) transversal magnetic waves and electric currents that propagate at Alfv\'en speeds. The transversal waves and electric currents, generated by the spicule's driver and lasting for many minutes, are dissipated and heat the associated loop. These complex interactions in the corona can be connected with blue shifted secondary components in coronal spectral lines (Red-Blue asymmetries) observed with Hinode/EIS and SOHO/SUMER, as well as the EUV counterpart of type II spicules and propagating coronal disturbances (PCDs) observed with the 171~\AA\ and 193~\AA\ SDO/AIA channels.


title: Control of Second Harmonic Generation in Doubly Resonant Aluminum Nitride Microrings to Address Rubidium Two-Photon Clock Transition

abstract: Nonlinear optical effects have been studied extensively in microresonators as more photonics applications transition to integrated on-chip platforms. Due to low optical losses and small mode volumes, microresonators are demonstrably the state-of-the-art platform for second harmonic generation (SHG). However, the working bandwidth of such microresonator-based devices are relatively small, presenting a challenge for applications where a specifically targeted wavelength needs to be addressed. In this work, we analyzed the phase-matching window and resonance wavelength with respect to varying microring width, radius and temperature. A chip with precise design parameters was fabricated with phase-matching realized at the exact wavelength of two-photon transition of 85-Rubidium. This procedure can be generalized to any target pump wavelength in the telecom-band with picometer precision.


title: The role of molecular gas in the nuclear regions of IRAS 00183-7111. ALMA and X-ray investigations of an ultraluminous infrared galaxy

abstract: We present a multi-frequency study of the ultraluminous infrared galaxy (ULIRG) IRAS 00183-7111 (z=0.327), selected from the Spoon diagnostic diagram as a highly obscured active galactic nucleus (AGN) candidate. ALMA Cycle 0 millimetre and X-ray observations are used; the main aim is to verify at what level the molecular gas, traced by the CO, may be responsible for the obscuration observed at X-ray energies. Theory and observations both suggest that galaxy-scale absorption may play a role in the AGN obscuration at intermediate (i.e. Compton-thin) column densities. We derived a molecular gas column density of $(8.0\pm0.9)\times10^{21}$ cm$^{-2}$ from the ALMA CO$_{(1-0)}$ detection, while the best-fit column density of cold gas obtained from X-ray spectral fitting is $6.8^{+2.1}_{-1.5}\times10^{22}$ cm$^{-2}$. The two quantities suggest that the molecular gas may contribute only a fraction of the AGN obscuration; however, the link between them is not straightforward. The nuclear regions of IRAS 00183-7111 are likely stratified into different layers of matter: one inner and highly ionized by the strong radiation field of the AGN (as inferred from the high-ionization iron line found in the X-ray spectra), and one outer and colder, extending more than 5~kpc from the nucleus (as traced by the molecular gas observed with ALMA). The molecular gas regions also give rise to a vigorous starburst with SFR$\sim260\pm28$ M$_{\odot}$ yr$^{-1}$. The complexity of this nuclear environment makes it difficult to identify the origin of the AGN obscuration given the quality of the data currently available. Higher resolution observations in the millimetre regime are needed to deeply investigate this issue.


title: Semi-parametric Bayesian change-point model based on the Dirichlet process

abstract: In this work we introduce a semi-parametric Bayesian change-point model, defining its time dynamic as a latent Markov process based on the Dirichlet process. We treat the number of change point as a random variable and we estimate it during model fitting. Posterior inference is carried out using a Markov chain Monte Carlo algorithm based on a marginalized version of the proposed model. The model is illustrated using simulated examples and two real datasets, namely the coal- mining disasters, that is a widely used dataset for illustrative purpose, and a dataset of indoor radon recordings. With the simulated examples we show that the model is able to recover the parameters and number of change points, and we compare our results with the ones of the-state- of-the-art models, showing a clear improvement in terms of change points identification. The results obtained on the coal-mining disasters and radon data are coherent with previous literature.


title: GalWeight: A New and Effective Weighting Technique for Determining Galaxy Cluster and Group Membership

abstract: We introduce GalWeight, a new technique for assigning galaxy cluster membership. This technique is specifically designed to simultaneously maximize the number of bona fide cluster members while minimizing the number of contaminating interlopers. The GalWeight technique can be applied to both massive galaxy clusters and poor galaxy groups. Moreover, it is effective in identifying members in both the virial and infall regions with high efficiency. We apply the GalWeight technique to MDPL2 \& Bolshoi N-body simulations, and find that it is $> 98\%$ accurate in correctly assigning cluster membership. We show that GalWeight compares very favorably against four well-known existing cluster membership techniques (shifting gapper, den Hartog, caustic, SIM). We also apply the GalWeight technique to a sample of twelve Abell clusters (including the Coma cluster) using observations from the Sloan Digital Sky Survey. We end by discussing GalWeight's potential for other astrophysical applications.


title: In-situ soil parametrization from multi-layer moisture data

abstract: Inversion methodology has been used to obtain, from multi-layer soil probes records, a complete soil parametrisation, namely water retention curve, unsaturated conductivity curve and bulk density at 4 depths. The approach integrates water dynamics, hysteresis and the effect of bulk density on conductivity to extract soil parameters required from most simulation models. The method is applied to sub-sets of data collection, allowing to understand that not every data-sets contains the information required for method convergence. A comparison with experimental bulk-density values show that inversion could give information even with a better adherence to model, as it considers the effect of roots and skeleton. The method may be applied to any type of multi-layer water content probes giving the opportunity to enrich soil parameter availability and reliability.


title: A new focal-plane 3D imaging method based on temporal ghost imaging

abstract: A new focal-plane three-dimensional (3D) imaging method based on temporal ghost imaging is proposed and demonstrated. By exploiting the advantages of temporal ghost imaging, this method enables slow integrating cameras have an ability of 3D surface imaging in the framework of sequential flood-illumination and focal-plane detection. The depth information of 3D objects is easily lost when imaging with traditional cameras, but it can be reconstructed with high-resolution by temporal correlation between received signals and reference signals. Combining with a two-dimensional (2D) projection image obtained by one single shot, a 3D image of the object can be achieved. The feasibility and performance of this focal-plane 3D imaging method have been verified through theoretical analysis and numerical experiments in this paper.


title: Two-loop corrections to false vacuum decay in scalar field theory

abstract: We consider radiative corrections to false vacuum decay in a four-dimensional scalar field theory with cubic and quartic potential. Using planar thin wall approximation we were able to get analytical expression for the decay rate up to two loop order. The results obtained employ dimensional regularization and $\overline{MS}$ renormalization scheme.


title: A new convexity-based inequality, characterization of probability distributions and some free-of-distribution tests

abstract: A new inequality between some functional of probability distribution functions is given. The inequality is based on strict convexity of a function used in functional definition. Equality sign in the inequality gives a characteristic property of some probability distributions. This fact together with special character of functional is used to construct free-of-distribution two sample tests. Key words: convex functions; probability distances; characterization of distributions; Cram\'{e}r - von Mises distance; statistical tests.


title: Klyshko's Advanced-Wave Picture in Stimulated Parametric Down-Conversion with a Spatially Structured Pump Beam

abstract: The advanced-wave picture is "... an intuitive treatment of two-photon correlation with the help of the concept of an effective field acting upon one of the two detectors and formed by parametric conversion of the advanced wave emitted by the second detector ..." [A. V. Belinskii and D. N. Klyshko, JETP 78, 259 (1994)]. This quote from Belinskii and Klyshko nicely describes the concept of the advanced-wave picture; an intuitive tool for designing and predicting results from coincidence-based two-photon experiments. Up to now, the advanced-wave picture has been considered primarily for the case of an ideal plane-wave pump beam and only for design purposes. Here we study the advanced wave picture for a structured pump beam and in the context of stimulated emission provoked by an auxiliary input laser beam. This suggests stimulated parametric down-conversion as a useful experimental tool for testing the experimental sets designed with the advanced-wave picture. We present experimental results demonstrating the strategy of designing the experiment with advanced-wave picture and testing with stimulated emission.


title: QuaterNet: A Quaternion-based Recurrent Model for Human Motion

abstract: Deep learning for predicting or generating 3D human pose sequences is an active research area. Previous work regresses either joint rotations or joint positions. The former strategy is prone to error accumulation along the kinematic chain, as well as discontinuities when using Euler angle or exponential map parameterizations. The latter requires re-projection onto skeleton constraints to avoid bone stretching and invalid configurations. This work addresses both limitations. Our recurrent network, QuaterNet, represents rotations with quaternions and our loss function performs forward kinematics on a skeleton to penalize absolute position errors instead of angle errors. On short-term predictions, QuaterNet improves the state-of-the-art quantitatively. For long-term generation, our approach is qualitatively judged as realistic as recent neural strategies from the graphics literature.


title: The origin and early evolution of a bipolar magnetic region in the solar photosphere

abstract: Finding the formation mechanisms for bipolar configurations of strong local magnetic field under control of the relatively weak global magnetic field of the Sun is a key problem of the physics of solar activity. This study is aimed at discriminating whether the magnetic field or fluid motion plays a primary, active role in this process. The very origin and early development stage of Active Region 12548 are investigated based on SDO/HMI observations of 2016 May 20--25. Full-vector magnetic and velocity fields are analyzed in parallel. The leading and trailing magnetic polarities are found to grow asymmetrically in terms of their amplitude, magnetic flux, and the time variation of these quantities. The leading-polarity magnetic element originates as a compact feature against the background of a distributed trailing-polarity field, with an already existing trailing-polarity magnetic element. No signs of strong horizontal magnetic fields are detected between the two magnetic poles. No predominant upflow between their future locations precedes the origin of this bipolar magnetic region (BMR); instead, upflows and downflows are mixed, with some prevalence of downflows. Any signs of a large-scale horizontal divergent flow from the area where the BMR develops are missing; in contrast, a normal supergranulation and mesogranulation pattern is preserved. This scenario of early BMR evolution is in strong contradiction with the expectations based on the model of a rising $\Omega$-shaped loop of a flux tube of strong magnetic field, and an \emph{in situ} mechanism of magnetic-field amplification and structuring should operate in this case.


title: Evolution of macromolecular structure: a 'double tale' of biological accretion

abstract: The evolution of structure in biology is driven by accretion and change. Accretion brings together disparate parts to form bigger wholes. Change provides opportunities for growth and innovation. Here we review patterns and processes that are responsible for a 'double tale' of evolutionary accretion at various levels of complexity, from proteins and nucleic acids to high-rise building structures in cities. Parts are at first weakly linked and associate variously. As they diversify, they compete with each other and are selected for performance. The emerging interactions constrain their structure and associations. This causes parts to self-organize into modules with tight linkage. In a second phase, variants of the modules evolve and become new parts for a new generative cycle of higher-level organization. Evolutionary genomics and network biology support the 'double tale' of structural module creation and validate an evolutionary principle of maximum abundance that drives the gain and loss of modules.


title: Valid and Approximately Valid Confidence Intervals for Current Status Data

abstract: We introduce a new framework for creating point-wise confidence intervals for the distribution of event times for current status data. Existing methods are based on asymptotics. Our framework is based on binomial properties and motivates confidence intervals that are very simple to apply and are valid, i.e., guarantee nominal coverage. Although these confidence intervals are necessarily conservative for small sample sizes, asymptotically their coverage rate approaches the nominal one. This binomial framework also motivates approximately valid confidence intervals, and simulations show that these approximate intervals generally have coverage rates closer to the nominal level with shorter length than existing intervals, including the likelihood ratio-based confidence interval. Unlike previous asymptotic methods that require different asymptotic distributions for continuous or grid-based assessment, the binomial framework can be applied to either type of assessment distribution.


title: Deterministic transformations of coherent states under incoherent operations

abstract: It is well known that the majorization condition is the necessary and sufficient condition for the deterministic transformations of both pure bipartite entangled states by local operations and coherent states under incoherent operations. In this paper, we present two explicit protocols for these transformations. We first present a permutation-based protocol which provides a method for the single-step transformation of $d$-dimensional coherent states. We also obtain generalized solutions of this protocol for some special cases of $d$-level systems. Then, we present an alternative protocol where we use $d'$-level ($d'$ $<$ $d$) subspace solutions of the permutation-based protocol to achieve the complete transformation as a sequence of coherent-state transformations. We show that these two protocols also provide solutions for deterministic transformations of pure bipartite entangled states.


title: Octet and decuplet contribution to the proton self energy

abstract: Within the hadronic language of Chiral Perturbation Theory we present the full leading-order octet-baryon$-$meson and decuplet-baryon$-$meson contribution to the proton self energy and thus to its wave function renormalization factor $Z$. By Fock-expanding the physical proton state into its bare and hadron-cloud part, we show how each individual baryon-meson probability depend on the average momenta of the particles in the fluctuation. We present how the results depend on the choice of the form factor involved in the regularization (Gaussian or Besselian) and how they depend on the cut-off parameter. We also show how the results vary with respect to a variation of the decuplet coupling constant $h_A$. The momentum distributions of the fluctuations are given and the fluctuations' relative probabilities are presented. We show that for reasonable values of the cut-off parameter, the Delta-pion fluctuation is of the same strength as the nucleon-pion fluctuation.


title: Deterministic chaos in an ytterbium-doped mode-locked fiber laser

abstract: We experimentally study the nonlinear dynamics of a femtosecond ytterbium doped mode-locked fiber laser. With the laser operating in the pulsed regime a route to chaos is presented, starting from stable mode-locking, period two, period four, chaos and period three regimes. Return maps and bifurcation diagrams were extracted from time series for each regime. The analysis of the time series with the laser operating in the quasi mode-locked regime presents deterministic chaos described by an unidimensional Rossler map. A positive Lyapunov exponent $\lambda = 0.14$ confirms the deterministic chaos of the system. We suggest an explanation about the observed map by relating gain saturation and intra-cavity loss.


title: Biorthogonal Bulk-Boundary Correspondence in Non-Hermitian Systems

abstract: Non-Hermitian systems exhibit striking exceptions from the paradigmatic bulk-boundary correspondence, including the failure of bulk Bloch band invariants in predicting boundary states and the (dis)appearance of boundary states at parameter values far from those corresponding to gap closings in periodic systems without boundaries. Here, we provide a comprehensive framework to unravel this disparity based on the notion of biorthogonal quantum mechanics: While the properties of the left and right eigenstates corresponding to boundary modes are individually decoupled from the bulk physics in non-Hermitian systems, their combined biorthogonal density penetrates the bulk precisely when phase transitions occur. This leads to generalized bulk-boundary correspondence and a quantized biorthogonal polarization that is formulated directly in systems with open boundaries. We illustrate our general insights by deriving the phase diagram for several microscopic open boundary models, including exactly solvable non-Hermitian extensions of the Su-Schrieffer-Heeger model and Chern insulators.


title: Flavor Symmetry and Ferroelectric Nematics in Transition Metal Dichalcogenides

abstract: Recent magneto-transport experiments have provided compelling evidence for the presence of an energetically isolated threefold Q-valley degeneracy in few-layer transition metal dichalcogenides. We study the flavor SU(3) symmetry breaking when each Landau level triplet is one-third filled or empty and predict that a pure flavor nematic phase and a flavorless charge-density-wave phase will occur respectively below and above a critical magnetic field. Surprisingly, electrons carry flavor-dependent electric dipole moments even at zero magnetic field, rendering the nematics ferroelectric, allowing electric-field manipulation of the flavors, and leading to the concept of flavortronics.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.str-el cond-mat.supr-con


title: Stripping a supermultiplet of all but one scalar: a microscopic view

abstract: We explain how a single scalar degree of freedom can be obtained from projecting out all the other components of a chiral superfield using simple operators in the microscopic theory. We use the Fayet-Iliopoulos model as an example of the origin for the necessary supersymmetry breaking. We also comment on one peculiar aspect where non-linear realisation can be useful: the freeze-in scenario for gravitino dark matter.


title: Geometric Phase using Diagonal Coherent State Representation

abstract: Glauber-Sudarshan diagonal coherent state P-representation has been used to determine geometric phase for non-classical states of light. For a given density operator $\hat{\rho_1}$ of two mode optical beam, we evolve it in complex projective ray space $(\cal R)$ to $\hat{\rho_2}$ and to $\hat{\rho_3}$ by changing its state of polarisation using unitary operator $\mathrm{\hat{U}_p}(\theta)$. The diagonal coherent state basis has been utilized to represent the density operators instead of fock state basis as in the fock state basis the state vector in present work evolve under unitary operator produe infinitely numerous terms which make the density operators very messy to handle. This cumbersome situation can be easily avoided by using the approach proposed above. The trace of the product of $\hat{\rho_1}$, $\hat{\rho_2}$ and $\hat{\rho_3}$ is taken to get three vertex Bargmann invariant. Argument of this gives the geometric phase which is represented in terms of phase space variables containing a notion of symplectic area.


title: Herschel Planetary Nebula Survey (HerPlaNS): Hydrogen Recombination Laser Lines in Mz 3

abstract: The bipolar nebula Menzel 3 (Mz 3) was observed as part of the \textit{Herschel} Planetary Nebula Survey (\textit{HerPlaNS}), which used the PACS and SPIRE instruments aboard the \textit{Herschel Space Observatory} to study a sample of planetary nebulae (PNe). In this paper, one of the series describing \textit{HerPlaNS} results, we report the detection of H I recombination lines (HRLs) in the spectrum of Mz 3. Inspection of the spectrum reveals the presence of 12 HRLs in the 55 to 680 $\mu$m range covered by the PACS and SPIRE instruments (H11$\alpha$ to H21$\alpha$ and H14$\beta$). The presence of HRLs in this range is unusual for PNe and has not been reported in Mz 3 before. Our analysis indicates that the HRLs we observed are enhanced by laser effect occurring in the core of Mz 3. Our arguments for this are: (i) the available Mz 3 optical to submillimetre HRL $\alpha$ line intensity ratios are not well reproduced by the spontaneous emission of optically thin ionized gas, as would be typical for nebular gas in PNe; (ii) the compact core of Mz 3 is responsible for a large fraction of the Herschel HRLs emission; (iii) the line intensity ratios for Mz 3 are very similar to those in the core emission of the well known star MWC 349A, where laser effect is responsible for the enhancement of HRLs in the Herschel wavelength range; (iv) the physical characteristics relevant to cause laser effect in the core of MWC 349A are very similar to those in the core of Mz 3.


title: Deconvolution of dust mixtures by latent Dirichlet allocation in forensic science

abstract: Dust particles recovered from the soles of shoes may be indicative of the sites recently visited by an individual, and, in particular, of the presence of an individual at a particular site of interest, e.g., the scene of a crime. By describing the dust profile of a given site by a multinomial distribution over a fixed number of dust particle types, we can define the probability distribution of the mixture of dust recovered from the sole of a shoe via Latent Dirichlet Allocation. We use Variational Bayesian Inference to study the parameters of the model, and use their resulting posterior distributions to make inference on (a) the contributions of sites of interest to a dust mixture, and (b) the particle profiles associated with these sites.


title: Utility maximization with proportional transaction costs under model uncertainty

abstract: We consider a discrete time financial market with proportional transaction costs under model uncertainty, and study a num\'eraire-based semi-static utility maximization problem with an exponential utility preference. The randomization techniques recently developed in \cite{BDT17} allow us to transform the original problem into a frictionless counterpart on an enlarged space. By suggesting a different dynamic programming argument than in \cite{bartl2016exponential}, we are able to prove the existence of the optimal strategy and the convex duality theorem in our context with transaction costs. In the frictionless framework, this alternative dynamic programming argument also allows us to generalize the main results in \cite{bartl2016exponential} to a weaker market condition. Moreover, as an application of the duality representation, some basic features of utility indifference prices are investigated in our robust setting with transaction costs.


title: QoE-Aware Beamforming Design for Massive MIMO Heterogeneous Networks

abstract: One of the main goals of the future wireless networks is improving the users quality of experience (QoE). In this paper, we consider the problem of QoE-based resource allocation in the downlink of a massive multiple-input multiple-output (MIMO) heterogeneous network (HetNet). The network consists of a macro cell with a number of small cells embedded in it. The small cells base stations (BSs) are equipped with a few antennas, while the macro BS is equipped with a massive number of antennas. We consider the two services Video and Web Browsing and design the beamforming vectors at the BSs. The objective is to maximize the aggregated Mean Opinion Score (MOS) of the users under constraints on the BSs powers and the required quality of service (QoS) of the users. We also consider extra constraints on the QoE of users to more strongly enforce the QoE in the beamforming design. To reduce the complexity of the optimization problem, we suggest suboptimal and computationally efficient solutions. Our results illustrate that increasing the number of antennas at the BSs and also increasing the number of small cells antennas in the network leads to a higher user satisfaction.


title: Fluidization of epithelial sheets by active cell rearrangements

abstract: We theoretically explore fluidization of epithelial tissues by active T1 neighbor exchanges. We show that the geometry of cell-cell junctions encodes important information about the local features of the energy landscape, which we support by an elastic theory of T1 transformations. Using a 3D vertex model, we show that the degree of active noise driving forced cell rearrangements governs the stress-relaxation time-scale of the tissue. We study tissue response to in-plane shear at different time scales. At short time, the tissue behaves as a solid, whereas its long-time fluid behavior can be associated with an effective viscosity which scales with the rate of active T1 transformations. Furthermore, we develop a coarse-grained theory, where we treat the tissue as an active fluid and confirm the results of the vertex model. The impact of cell rearrangements on tissue shape is illustrated by studying axial compression of an epithelial tube.


title: Evolution of shell structure in exotic nuclei

abstract: The atomic nucleus is a quantum many-body system whose constituent nucleons (protons and neutrons) are subject to complex nucleon-nucleon interactions that include spin- and isospin-dependent components. For stable nuclei, already several decades ago, emerging seemingly regular patterns in some observables could be described successfully within a shell-model picture that results in particularly stable nuclei at certain magic fillings of the shells with protons and/or neutrons: N,Z = 8, 20, 28, 50, 82, 126. However, in short-lived, so-called exotic nuclei or rare isotopes, characterized by a large N/Z asymmetry and located far away from the valley of beta stability on the nuclear chart, these magic numbers, viewed through observables, were shown to change. These changes in the regime of exotic nuclei offer an unprecedented view at the roles of the various components of the nuclear force when theoretical descriptions are confronted with experimental data on exotic nuclei where certain effects are enhanced. This article reviews the driving forces behind shell evolution from a theoretical point of view and connects this to experimental signatures.


title: First Experiments with Neural Translation of Informal to Formal Mathematics

abstract: We report on our experiments to train deep neural networks that automatically translate informalized LaTeX-written Mizar texts into the formal Mizar language. To the best of our knowledge, this is the first time when neural networks have been adopted in the formalization of mathematics. Using Luong et al.'s neural machine translation model (NMT), we tested our aligned informal-formal corpora against various hyperparameters and evaluated their results. Our experiments show that our best performing model configurations are able to generate correct Mizar statements on 65.73\% of the inference data, with the union of all models covering 79.17\%. These results indicate that formalization through artificial neural network is a promising approach for automated formalization of mathematics. We present several case studies to illustrate our results.


title: Weight Initialization in Neural Language Models

abstract: Semantic Similarity is an important application which finds its use in many downstream NLP applications. Though the task is mathematically defined, semantic similarity's essence is to capture the notions of similarity impregnated in humans. Machines use some heuristics to calculate the similarity between words, but these are typically corpus dependent or are useful for specific domains. The difference between Semantic Similarity and Semantic Relatedness motivates the development of new algorithms. For a human, the word car and road are probably as related as car and bus. But this may not be the case for computational methods. Ontological methods are good at encoding Semantic Similarity and Vector Space models are better at encoding Semantic Relatedness. There is a dearth of methods which leverage ontologies to create better vector representations. The aim of this proposal is to explore in the direction of a hybrid method which combines statistical/vector space methods like Word2Vec and Ontological methods like WordNet to leverage the advantages provided by both.


title: Analogical Reasoning on Chinese Morphological and Semantic Relations

abstract: Analogical reasoning is effective in capturing linguistic regularities. This paper proposes an analogical reasoning task on Chinese. After delving into Chinese lexical knowledge, we sketch 68 implicit morphological relations and 28 explicit semantic relations. A big and balanced dataset CA8 is then built for this task, including 17813 questions. Furthermore, we systematically explore the influences of vector representations, context features, and corpora on analogical reasoning. With the experiments, CA8 is proved to be a reliable benchmark for evaluating Chinese word embeddings.


title: Realization of Third Order Exceptional Singularities in a Three level non-Hermitian System: Towards Cascaded State Conversion

abstract: The appearance of topological singularities, namely exceptional points (EPs) is an intriguing feature of parameter-dependent open quantum or wave systems. EPs are the special type of nonHermitian degeneracies where two (or more) eigenstates of the underlying system coalesce. In this paper, we present a three level non-Hermitian Hamiltonian which hosts three interacting eigenstates. The matrix elements are optimized in such a way that the intermediate eigenstate interacts with both the other states and the underlying system hosts at least two different second order EPs. The impact of quasi-static parameter variation along a cyclic contour around the embedded EPs on the dynamics of interacting eigenvalues is well investigated in the context of cascaded state conversion. Such dynamics of the eigenvalues shows a clear signature of the third order EP with a combined effect of both the second order EPs. Moreover, we examine the accumulation of phases around the identified EPs and study the hallmark of phase exchange during cascaded state conversions accompanied by the parametric encirclement of the third order EP.


title: Emergence of Benford's Law in Music

abstract: We analyzed a large selection of classical musical pieces composed by Bach, Beethoven, Mozart, Schubert and Tchaikovsky, and found a surprising connection with mathematics. For each composer, we extracted the time intervals each note was played in each piece and found that the corresponding data sets are Benford distributed. Remarkably, the logarithmic distribution is present not only for the leading digits, but also for all digits.


title: Nowhere-differentiability of the solution map of 2D Euler equations on bounded spatial domain

abstract: We consider the incompressible 2D Euler equations on bounded spatial domain $S$, and study the solution map on the Sobolev spaces $H^k(S)$ ($k > 2$). Through an elaborate geometric construction, we show that for any $T >0$, the time $T$ solution map $u_0 \mapsto u(T)$ is nowhere locally uniformly continuous and nowhere Fr\'echet differentiable.

categories: math.AP math-ph math.DS math.MP nlin.CD physics.flu-dyn


title: Sex-Disaggregated Systematics in Canadian Time Allocation Committee Telescope Proposal Reviews

abstract: Recent studies have shown that the proposal peer review processes employed by a variety of organizations to allocate astronomical telescope time produce outcomes that are systematically biased depending on whether proposal's principal investigator (PI) is a man or a woman. Using Canada-France-Hawaii Telescope (CFHT) and Gemini Observatory proposal statistics from Canada over 10 recent proposal cycles, we assess whether or not the mean proposal scores assigned by the National Research Council's (NRC's) Canadian Time Allocation Committee (CanTAC) also correlate significantly with PI sex. Classical t-tests, bootstrap and jackknife replications show that proposals submitted by women were rated significantly worse than those submitted by men. We subdivide the data in order to investigate sex-disaggregated statistics in relation to PI career stage (faculty vs. non-faculty), telescope requested, scientific review panel, observing semester, and the PhD year of faculty PIs. Consistent with the bivariate results, a multivariate regression analysis controlling for other covariates confirmed that PI sex is the only significant predictor of proposal rating scores for the sample as a whole, although differences emerge for proposals submitted by faculty and non-faculty PIs. While further research is needed to explain our results, it is possible that implicit social cognition is at work. NRC and CanTAC have taken steps to mitigate this possibility by altering proposal author lists in order to conceal the PI's identity among co-investigators. We recommend that the impact of this measure on mitigating bias in future observing semesters be quantitatively assessed using statistical techniques such as those employed here.


title: Density Response from Kinetic Theory and Time Dependent Density Functional Theory for Matter Under Extreme Conditions

abstract: The density linear response function for an inhomogeneous system of electrons in equilibrium with an array of fixed ions is considered. Two routes to its evaluation for extreme conditions (e.g., warm dense matter) are considered. The first is from a recently developed short-time kinetic equation; the second is from time-dependent density functional theory (tdDFT). The result from the latter approach agrees with that from kinetic theory in the "adiabatic approximation", providing support and context for each. Both provide a connection to the phenomenological Kubo-Greenwood method for calculating transport properties. A brief proof of the van Leeuwen theorem (an essential underpinning of tdDFT) extended to the mixed states of equilibrium ensembles is given.


title: Facebook Reaction-Based Emotion Classifier as Cue for Sarcasm Detection

abstract: Online social media users react to content in them based on context. Emotions or mood play a significant part of these reactions, which has filled these platforms with opinionated content. Different approaches and applications to make better use of this data are continuously being developed. However, due to the nature of the data, the variety of platforms, and dynamic online user behavior, there are still many issues to be dealt with. It remains a challenge to properly obtain a reliable emotional status from a user prior to posting a comment. This work introduces a methodology that explores semi-supervised multilingual emotion detection based on the overlap of Facebook reactions and textual data. With the resulting emotion detection system we evaluate the possibility of using emotions and user behavior features for the task of sarcasm detection. More than 1 million English and Chinese comments from over 62,000 public Facebook pages posts have been collected and processed, conducted experiments show acceptable performance metrics.


title: Improving End-of-turn Detection in Spoken Dialogues by Detecting Speaker Intentions as a Secondary Task

abstract: This work focuses on the use of acoustic cues for modeling turn-taking in dyadic spoken dialogues. Previous work has shown that speaker intentions (e.g., asking a question, uttering a backchannel, etc.) can influence turn-taking behavior and are good predictors of turn-transitions in spoken dialogues. However, speaker intentions are not readily available for use by automated systems at run-time; making it difficult to use this information to anticipate a turn-transition. To this end, we propose a multi-task neural approach for predicting turn- transitions and speaker intentions simultaneously. Our results show that adding the auxiliary task of speaker intention prediction improves the performance of turn-transition prediction in spoken dialogues, without relying on additional input features during run-time.


title: The Broken Stick Project

abstract: The broken stick problem is the following classical question. You have a segment $[0,1]$. You choose two points on this segment at random. They divide the segment into three smaller segments. Show that the probability that the three segments form a triangle is $1/4$. The MIT PRIMES program, together with Art of Problem Solving, organized a high school research project where participants worked on several variations of this problem. Participants were generally high school students who posted ideas and progress to the Art of Problem Solving forums over the course of an entire year, under the supervision of PRIMES mentors. This report summarizes the findings of this CrowdMath project.


title: Constraints on primordial black hole dark matter from Galactic center X-ray observations

abstract: Surprisingly high masses of the black holes inferred from the LIGO & Virgo gravitational wave measurements have lead to speculations that the observed mergers might be due to ${\cal O}(10) M_\odot$ primordial black holes (PBHs). Furthermore, it has been suggested that the whole amount of dark matter (DM) might be in that exotic form. We investigate constraints on the PBH DM using NuSTAR Galactic center (GC) X-ray data. We used a robust Monte Carlo approach in conjunction with a radiatively inefficient PBH accretion model with commonly accepted model parameters. Compared to previous studies we allowed for multiple forms of DM density profiles. Most importantly, our study includes treatment of the gas turbulence, which significantly modifies the relative velocity between PBHs and gas. We show that inclusion of the effects of gas turbulence and the uncertainties related to the DM density profile reduces significantly the gas accretion onto PBHs compared to the claimed values in previous papers. It is highly improbable to obtain accreting PBHs brighter than the NuSTAR point source limit using observationally determined gas velocities. As such, one can safely conclude that GC X-ray observations cannot rule out ${\cal O}(10) M_\odot$ PBH DM.


title: Affine Matsuki correspondence for sheaves

abstract: We lift the affine Matsuki correspondence between real and symmetric loop group orbits in affine Grassmannians to an equivalence of derived categories of sheaves. In analogy with the finite-dimensional setting, our arguments depend upon the Morse theory of energy functions obtained from symmetrizations of coadjoint orbits. The additional fusion structures of the affine setting lead to further equivalences with Schubert constructible derived categories of sheaves on real affine Grassmannians.


title: Remote Source Coding under Gaussian Noise : Dueling Roles of Power and Entropy Power

abstract: The distributed remote source coding (so-called CEO) problem is studied in the case where the underlying source, not necessarily Gaussian, has finite differential entropy and the observation noise is Gaussian. The main result is a new lower bound for the sum-rate-distortion function under arbitrary distortion measures. When specialized to the case of mean-squared error, it is shown that the bound exactly mirrors a corresponding upper bound, except that the upper bound has the source power (variance) whereas the lower bound has the source entropy power. Bounds exhibiting this pleasing duality of power and entropy power have been well known for direct and centralized source coding since Shannon's work. While the bounds hold generally, their value is most pronounced when interpreted as a function of the number of agents in the CEO problem.


title: Dependence upon initial conditions

abstract: This article discusses dependence on initial conditions in natural and social sciences with focus on physical science. The main focus is on the newly discovered rough dependence on initial data.

categories: nlin.CD math-ph math.AP math.MP physics.ao-ph physics.flu-dyn


title: Using three-partite GHZ states for partial quantum error-detection in entanglement-based protocols

abstract: The problem of noise incidence on qubits taking part of bipartite entanglement-based protocols is addressed. It is shown that the use of a three-partite GHZ state and measurements instead of their EPR counterparts allows the experimenter to detect $2/3$ of the times whenever one of the qubits involved in the measurement is affected by bit-flip noise through the mere observation of unexpected outcomes in the teleportation and superdense coding protocols when compared to the ideal case. It is shown that the use of post-selection after the detection of noise leads to an enhancement in the efficiency of the protocols. The idea is extended to any protocol using entangled states and measurements. Furthermore it is provided a generalization in which GHZ states and measurements with an arbitrary amount of qubits are used instead of EPR pairs, and remarkably, it is concluded that the optimal number of qubits is only three.


title: On displacement of viscous liquid in a system of parallel tubes

abstract: An explicitly solvable quasi 1D model of oil displacement is studied. The problem of recovering of the reservoir geometry is solved by means of a fixed point algorithm. The stability of solution is studied in various functional classes.


title: Hypergeometric expansions of the general Heun function governed by two-term recurrence relations

abstract: We show that there exist infinitely many particular choices of parameters for which the three-term recurrence relations governing the expansions of the solutions of the general Heun equation in terms of the Gauss hypergeometric functions become two-term. In these cases the coefficients are explicitly written in terms of the gamma functions.


title: $L^p$ norms of the lattice point discrepancy

abstract: We estimate the $L^{p}$ norms of the discrepancy between the volume and the number of integer points in $r\Omega-x$, a dilated by a factor $r$ and translated by a vector $x$ of a convex body $\Omega$ in $\mathbb{R}^{d}$ with smooth boundary with strictly positive curvature, \[ \left\{ {\displaystyle\int_{\mathbb R}}{\displaystyle\int_{\mathbb{T}^{d}}}\left\vert \sum_{k\in\mathbb{Z}^{d}}\chi _{r\Omega-x}(k)-r^{d}\left\vert \Omega\right\vert \right\vert ^{p}dxd\mu(r-R) \right\} ^{1/p}, \] where $\mu$ is a Borel measure compactly supported on the positive real axis and $R\to+\infty$.


title: Composite Semantic Relation Classification

abstract: Different semantic interpretation tasks such as text entailment and question answering require the classification of semantic relations between terms or entities within text. However, in most cases it is not possible to assign a direct semantic relation between entities/terms. This paper proposes an approach for composite semantic relation classification, extending the traditional semantic relation classification task. Different from existing approaches, which use machine learning models built over lexical and distributional word vector features, the proposed model uses the combination of a large commonsense knowledge base of binary relations, a distributional navigational algorithm and sequence classification to provide a solution for the composite semantic relation classification problem.


title: Semantic Relatedness for All (Languages): A Comparative Analysis of Multilingual Semantic Relatedness Using Machine Translation

abstract: This paper provides a comparative analysis of the performance of four state-of-the-art distributional semantic models (DSMs) over 11 languages, contrasting the native language-specific models with the use of machine translation over English-based DSMs. The experimental results show that there is a significant improvement (average of 16.7% for the Spearman correlation) by using state-of-the-art machine translation approaches. The results also show that the benefit of using the most informative corpus outweighs the possible errors introduced by the machine translation. For all languages, the combination of machine translation over the Word2Vec English distributional model provided the best results consistently (average Spearman correlation of 0.68).


title: End-to-end Learning of a Convolutional Neural Network via Deep Tensor Decomposition

abstract: In this paper we study the problem of learning the weights of a deep convolutional neural network. We consider a network where convolutions are carried out over non-overlapping patches with a single kernel in each layer. We develop an algorithm for simultaneously learning all the kernels from the training data. Our approach dubbed Deep Tensor Decomposition (DeepTD) is based on a rank-1 tensor decomposition. We theoretically investigate DeepTD under a realizable model for the training data where the inputs are chosen i.i.d. from a Gaussian distribution and the labels are generated according to planted convolutional kernels. We show that DeepTD is data-efficient and provably works as soon as the sample size exceeds the total number of convolutional weights in the network. We carry out a variety of numerical experiments to investigate the effectiveness of DeepTD and verify our theoretical findings.


title: Hybrid Adaptive Fuzzy Extreme Learning Machine for text classification

abstract: In traditional ELM and its improved versions suffer from the problems of outliers or noises due to overfitting and imbalance due to distribution. We propose a novel hybrid adaptive fuzzy ELM(HA-FELM), which introduces a fuzzy membership function to the traditional ELM method to deal with the above problems. We define the fuzzy membership function not only basing on the distance between each sample and the center of the class but also the density among samples which based on the quantum harmonic oscillator model. The proposed fuzzy membership function overcomes the shortcoming of the traditional fuzzy membership function and could make itself adjusted according to the specific distribution of different samples adaptively. Experiments show the proposed HA-FELM can produce better performance than SVM, ELM, and RELM in text classification.


title: Text classification based on ensemble extreme learning machine

abstract: In this paper, we propose a novel approach based on cost-sensitive ensemble weighted extreme learning machine; we call this approach AE1-WELM. We apply this approach to text classification. AE1-WELM is an algorithm including balanced and imbalanced multiclassification for text classification. Weighted ELM assigning the different weights to the different samples improves the classification accuracy to a certain extent, but weighted ELM considers the differences between samples in the different categories only and ignores the differences between samples within the same categories. We measure the importance of the documents by the sample information entropy, and generate cost-sensitive matrix and factor based on the document importance, then embed the cost-sensitive weighted ELM into the AdaBoost.M1 framework seamlessly. Vector space model(VSM) text representation produces the high dimensions and sparse features which increase the burden of ELM. To overcome this problem, we develop a text classification framework combining the word vector and AE1-WELM. The experimental results show that our method provides an accurate, reliable and effective solution for text classification.


title: Sun flux variation due to the effects of orbiting planets. Case of study of a non-compact planetary system

abstract: We study the phase curves for the planets of our Solar System; which, is considered as a non-compact planetary system. We focus on modeling the small variations of the light curve, based on the three photometric effects: reflection, ellipsoidal, and Doppler beaming. Theoretical predictions for these photometric variations are proposed, as if a hypothetical external observer would measure them. In contrast to similar studies of multi-planetary systems, the physical and geometrical parameters for each planet of the Solar System are well-known. Therefore, we can evaluate with accuracy the mathematical relations that shape the planetary light curves for an external fictitious observer. Our results suggest that in all the planets of study the ellipsoidal effect is very weak, while the Doppler beaming effect is in general dominant. In fact, the latter effect seems to be confirmed as the principal cause of variations of the light curves for the planets. This affirmation could not be definitive in Mercury or Venus where the Doppler beaming and the reflection effects have similar amplitudes. The obtained phase curves for the Solar System planets show interesting new features that have not been presented before, so the results presented here are relevant in their application to other non-compact systems, since they allow us to have an idea of what it is expected to find in their light curves.


title: Inequalities and separation for covariant Schr\"odinger operators

abstract: We consider a differential expression $L^{\nabla}_{V}=\nabla^{\dagger}\nabla+V$, where $\nabla$ is a metric covariant derivative on a Hermitian bundle $E$ over a geodesically complete Riemannian manifold $(M,g)$ with metric $g$, and $V$ is a linear self-adjoint bundle map on $E$. In the language of Everitt and Giertz, the differential expression $L^{\nabla}_{V}$ is said to be separated in $L^p(E)$ if for all $u\in L^p(E)$ such that $L^{\nabla}_{V}u\in L^p(E)$, we have $Vu\in L^p(E)$. We give sufficient conditions for $L^{\nabla}_{V}$ to be separated in $L^2(E)$. We then study the problem of separation of $L^{\nabla}_{V}$ in the more general $L^p$-spaces, and give sufficient conditions for $L^{\nabla}_{V}$ to be separated in $L^p(E)$, when $1<p<\infty$.


title: Pulsating Fronts for a Bistable Lotka-Volterra Competition System with Advection in a Periodic Habitat

abstract: This paper is concerned with the following Lotka-Volterra competition system with advection in a periodic habitat \begin{equation*} \begin{cases} \frac{\partial u_1}{\partial t} =d_1(x)\frac{\partial^2 u_1}{\partial x^2}-a_1(x)\frac{\partial u_1}{\partial x}+u_1\left(b_1(x)-a_{11}(x)u_1-a_{12}(x)u_2\right),\\ \frac{\partial u_2}{\partial t} =d_2(x)\frac{\partial^2 u_2}{\partial x^2}-a_2(x)\frac{\partial u_2}{\partial x}+u_2\left(b_2(x)-a_{21}(x)u_1-a_{22}(x)u_2\right), \end{cases} t>0,~x\in\Bbb R, \end{equation*} where $d_i(\cdot)$, $a_i(\cdot)$, $b_i(\cdot)$, $a_{ij}(\cdot)$ $(i,j=1,2)$ are $L$-periodic functions in $C^\nu(\Bbb{R})$ with some $\nu\in(0,1)$. Under certain assumptions, the system admits two periodic locally stable steady states $(u_1^*(x),0)$ and $(0,u_2^*(x))$. In this work, we first establish the existence of the pulsating front $U(x,x+ct)=(U_1(x,x+ct),U_2(x,x+ct))$ connecting two periodic solutions $(0,u_2^*(x))$ and $(u_1^*(x),0)$ at infinities. By using a dynamical method, we confirm further that the pulsating front is asymptotically stable for front-like initial values. As a consequence of the global asymptotically stability, we finally show that the pulsating front is unique up to translation.


title: Slicing Virtualized EPC-based 5G Core Network for Content Delivery

abstract: Traditional Content Delivery Networks (CDNs) built with traditional Internet technology are less and less able to cope with today's tremendous growth of content. Information Centric Networks (ICN), a proposed future Internet technology, may aid in remedying the situation. Unlike the current Internet, it decouples information from its sources and provides in-network storage. We expect traditional CDN and ICN-based CDN to co-exist in the foreseeable future, especially as it is now known that it might be possible to evolve traditional CDNs to gain the benefits promised by ICN. 5G providers must therefore aim to offer core network slices on which both ICN-based CDNs and traditional CDNs can be built. These slices could of course also be offered to providers of other applications with requirements similar to those of content delivery. This paper tackles the problem of slicing 5G for content delivery over ICN-based CDNs and traditional CDNs. Only virtualized Evolved Packet Core (EPC)-based 5G is considered. The problem is defined as a resource allocation problem which aims at minimizing the cost of slice assignment, while meeting QoS requirements. An Integer linear programming (ILP) formulation is provided and evaluated in a small-scale scenario.


title: Improving the Gaussian Mechanism for Differential Privacy: Analytical Calibration and Optimal Denoising

abstract: The Gaussian mechanism is an essential building block used in multitude of differentially private data analysis algorithms. In this paper we revisit the Gaussian mechanism and show that the original analysis has several important limitations. Our analysis reveals that the variance formula for the original mechanism is far from tight in the high privacy regime ($\varepsilon \to 0$) and it cannot be extended to the low privacy regime ($\varepsilon \to \infty$). We address these limitations by developing an optimal Gaussian mechanism whose variance is calibrated directly using the Gaussian cumulative density function instead of a tail bound approximation. We also propose to equip the Gaussian mechanism with a post-processing step based on adaptive estimation techniques by leveraging that the distribution of the perturbation is known. Our experiments show that analytical calibration removes at least a third of the variance of the noise compared to the classical Gaussian mechanism, and that denoising dramatically improves the accuracy of the Gaussian mechanism in the high-dimensional regime.


title: Continuum models of collective cell migration

abstract: Collective cell migration plays a central role in tissue development, morphogenesis, wound repair and cancer progression. With the growing realization that physical forces mediate cell motility in development and physiology, a key biological question is how cells integrate molecular activities for force generation on multicellular scales. In this review we discuss recent advances in modeling collective cell migration using quantitative tools and approaches rooted in soft matter physics. We focus on theoretical models of cell aggregates as continuous active media, where the feedback between mechanical forces and regulatory biochemistry gives rise to rich collective dynamical behavior. This class of models provides a powerful predictive framework for the physiological dynamics that underlies many developmental processes, where cells need to collectively migrate like a viscous fluid to reach a target region, and then stiffen to support mechanical stresses and maintain tissue cohesion.

categories: physics.bio-ph cond-mat.soft q-bio.CB q-bio.TO


title: Beyond 5G with UAVs: Foundations of a 3D Wireless Cellular Network

abstract: In this paper, a novel concept of three-dimensional (3D) cellular networks, that integrate drone base stations (drone-BS) and cellular-connected drone users (drone-UEs), is introduced. For this new 3D cellular architecture, a novel framework for network planning for drone-BSs as well as latency-minimal cell association for drone-UEs is proposed. For network planning, a tractable method for drone-BSs' deployment based on the notion of truncated octahedron shapes is proposed that ensures full coverage for a given space with minimum number of drone-BSs. In addition, to characterize frequency planning in such 3D wireless networks, an analytical expression for the feasible integer frequency reuse factors is derived. Subsequently, an optimal 3D cell association scheme is developed for which the drone-UEs' latency, considering transmission, computation, and backhaul delays, is minimized. To this end, first, the spatial distribution of the drone-UEs is estimated using a kernel density estimation method, and the parameters of the estimator are obtained using a cross-validation method. Then, according to the spatial distribution of drone-UEs and the locations of drone-BSs, the latency-minimal 3D cell association for drone-UEs is derived by exploiting tools from optimal transport theory. Simulation results show that the proposed approach reduces the latency of drone-UEs compared to the classical cell association approach that uses a signal-to-interference-plus-noise ratio (SINR) criterion. In particular, the proposed approach yields a reduction of up to 46% in the average latency compared to the SINR-based association. The results also show that the proposed latency-optimal cell association improves the spectral efficiency of a 3D wireless cellular network of drones.


title: Modeling Naive Psychology of Characters in Simple Commonsense Stories

abstract: Understanding a narrative requires reading between the lines and reasoning about the unspoken but obvious implications about events and people's mental states - a capability that is trivial for humans but remarkably hard for machines. To facilitate research addressing this challenge, we introduce a new annotation framework to explain naive psychology of story characters as fully-specified chains of mental states with respect to motivations and emotional reactions. Our work presents a new large-scale dataset with rich low-level annotations and establishes baseline performance on several new tasks, suggesting avenues for future research.


title: Career Transitions and Trajectories: A Case Study in Computing

abstract: From artificial intelligence to network security to hardware design, it is well-known that computing research drives many important technological and societal advancements. However, less is known about the long-term career paths of the people behind these innovations. What do their careers reveal about the evolution of computing research? Which institutions were and are the most important in this field, and for what reasons? Can insights into computing career trajectories help predict employer retention? In this paper we analyze several decades of post-PhD computing careers using a large new dataset rich with professional information, and propose a versatile career network model, R^3, that captures temporal career dynamics. With R^3 we track important organizations in computing research history, analyze career movement between industry, academia, and government, and build a powerful predictive model for individual career transitions. Our study, the first of its kind, is a starting point for understanding computing research careers, and may inform employer recruitment and retention mechanisms at a time when the demand for specialized computational expertise far exceeds supply.


title: Stabilization rates for the damped wave equation with H\"older-regular damping

abstract: We study the decay rate of the energy of solutions to the damped wave equation in a setup where the geometric control condition is violated. We consider damping coefficients which are $0$ on a strip and vanish like polynomials, $x^{\beta}$. We prove that the semigroup cannot be stable at rate faster than $1/t^{(\beta+2)/(\beta+3)}$ by producing quasimodes of the associated stationary damped wave equation. We also prove that the semigroup is stable at rate at least as fast as $1/t^{(\beta+2)/(\beta+4)}$. These two results establish an explicit relation between the rate of vanishing of the damping and rate of decay of solutions. Our result partially generalizes a decay result of Nonnemacher in which the damping is an indicator function on a strip.


title: Are BLEU and Meaning Representation in Opposition?

abstract: One of possible ways of obtaining continuous-space sentence representations is by training neural machine translation (NMT) systems. The recent attention mechanism however removes the single point in the neural network from which the source sentence representation can be extracted. We propose several variations of the attentive NMT architecture bringing this meeting point back. Empirical evaluation suggests that the better the translation quality, the worse the learned sentence representations serve in a wide range of classification and similarity tasks.


title: Direct transcription methods based on fractional integral approximation formulas for solving nonlinear fractional optimal control problems

abstract: This paper presents three direct methods based on Gr\"{u}nwald-Letnikov, trapezoidal and Simpson fractional integral formulas to solve fractional optimal control problems (FOCPs). At first, the fractional integral form of FOCP is considered, then the fractional integral is approximated by Gr\"{u}nwald-Letnikov, trapezoidal and Simpson formulas in a matrix approach. Thereafter, the performance index is approximated either by trapezoidal or Simpson quadrature. As a result, FOCP are reduced to nonlinear programming problems, which can be solved by many well-developed algorithms. To improve the efficiency of the presented method, the gradient of the objective function and the Jacobian of constraints are prepared in closed forms. It is pointed out that the implementation of the methods is simple and, due to the fact that there is no need to derive necessary conditions, the methods can be simply and quickly used to solve a wide class of FOCPs. The efficiency and reliability of the presented methods are assessed by ample numerical tests involving a free final time with path constraint FOCP, a bang-bang FOCP and an optimal control of a fractional-order HIV-immune system.


title: The Photophobic ALP

abstract: We explore models and phenomenology of a photophobic axion-like particle (ALP), an axion whose coupling to photons is maximally suppressed without fine-tuning of the underlying parameters. We demonstrate that photophobia can be a natural UV property of ALP models and determine the irreducible coupling of photophobic ALPs to photons induced by violations of the axion shift symmetry. The signatures of photophobic axions are radically different from those of typical ALPs and are of particular interest for collider-based experiments, for which Standard Model triboson measurements provide a significant probe at higher masses. A variety of terrestrial and astrophysical measurements constrain the parameter space of photophobic ALPs, though bounds are typically much weaker compared to typical ALPs. We discuss implications for particle production relaxion models based on photophobic ALPs, finding that they are in mild tension with existing experimental constraints.


title: Beyond Structural Causal Models: Causal Constraints Models

abstract: Structural Causal Models (SCMs) provide a popular causal modeling framework. In this work, we show that SCMs are not flexible enough to give a complete causal representation of dynamical systems at equilibrium. Instead, we propose a generalization of the notion of an SCM, that we call Causal Constraints Model (CCM), and prove that CCMs do capture the causal semantics of such systems. We show how CCMs can be constructed from differential equations and initial conditions and we illustrate our ideas further on a simple but ubiquitous (bio)chemical reaction. Our framework also allows to model functional laws, such as the ideal gas law, in a sensible and intuitive way.


title: An exact particle splitting algorithm for PIC codes

abstract: Splitting and merging are long standing issues in PIC codes. I propose a novel algorithm devoted to exact splitting for Particle-In-Cell (PIC) codes relying on Adaptive Mesh Refinement (AMR) grids. AMR grids have - by definition - a constant refinement factor between levels, this strong property combined with B-spline shape functions makes possible exact splitting of macro-particles in space. The proposed exact splitting method is applicable for any integer refinement factor and B-spline shape functions of arbitrary order. Besides, this exact algorithm removes by construction the usual numerical noise increase resulting from approximate particle splitting strategies.


title: Towards Malware Detection via CPU Power Consumption: Data Collection Design and Analytics (Extended Version)

abstract: This paper presents an experimental design and data analytics approach aimed at power-based malware detection on general-purpose computers. Leveraging the fact that malware executions must consume power, we explore the postulate that malware can be accurately detected via power data analytics. Our experimental design and implementation allow for programmatic collection of CPU power profiles for fixed tasks during uninfected and infected states using five different rootkits. To characterize the power consumption profiles, we use both simple statistical and novel, sophisticated features. We test a one-class anomaly detection ensemble (that baselines non-infected power profiles) and several kernel-based SVM classifiers (that train on both uninfected and infected profiles) in detecting previously unseen malware and clean profiles. The anomaly detection system exhibits perfect detection when using all features and tasks, with smaller false detection rate than the supervised classifiers. The primary contribution is the proof of concept that baselining power of fixed tasks can provide accurate detection of rootkits. Moreover, our treatment presents engineering hurdles needed for experimentation and allows analysis of each statistical feature individually. This work appears to be the first step towards a viable power-based detection capability for general-purpose computers, and presents next steps toward this goal.