title: On Lie algebras responsible for zero-curvature representations and Backlund transformations of (1+1)-dimensional scalar evolution PDEs

abstract: Zero-curvature representations (ZCRs) are one of the main tools in the theory of integrable PDEs. In particular, Lax pairs for $(1+1)$-dimensional PDEs can be interpreted as ZCRs. In [arXiv:1303.3575], for any $(1+1)$-dimensional scalar evolution equation $E$, we defined a family of Lie algebras $F(E)$ which are responsible for all ZCRs of $E$ in the following sense. Representations of the algebras $F(E)$ classify all ZCRs of the equation $E$ up to local gauge transformations. Also, using these algebras, one obtains necessary conditions for existence of a Backlund transformation between two given equations. The algebras $F(E)$ are defined in [arXiv:1303.3575] in terms of generators and relations. In this approach, ZCRs may depend on partial derivatives of arbitrary order, which may be higher than the order of the equation $E$. The algebras $F(E)$ generalize Wahlquist-Estabrook prolongation algebras, which are responsible for a much smaller class of ZCRs. In this preprint we prove a number of results on $F(E)$ which were announced in [arXiv:1303.3575]. We present applications of $F(E)$ to the theory of Backlund transformations in more detail and describe the explicit structure (up to non-essential nilpotent ideals) of the algebras $F(E)$ for a number of equations of orders $3$ and $5$.


title: Analytic Two-Loop Higgs Amplitudes in Effective Field Theory and the Maximal Transcendentality Principle

abstract: We obtain for the first time the two-loop amplitudes for Higgs plus three gluons in Higgs effective field theory including dimension-seven operators. This provides the S-matrix elements for the top mass corrections for Higgs plus a jet production at LHC. The computation is based on the on-shell unitarity method combined with integration by parts reduction. We work in conventional dimensional regularization and obtain analytic expressions renormalized in the MS-bar scheme. The two-loop anomalous dimensions present operator mixing behavior. The infrared divergences agree with that predicted by Catani and the finite remainders take remarkably simple forms, where the maximally transcendental parts are identical to the corresponding results in N=4 super-Yang-Mills theory. The parts of lower transcendentality turn out to be also largely determined by the N=4 results.


title: A panoply of Schwinger-Keldysh transport

abstract: We classify all possible allowed constitutive relations of relativistic fluids in a statistical mechanical limit using the Schwinger-Keldysh effective action for hydrodynamics. We find that microscopic unitarity enforces genuinely new constraints on the allowed transport coefficients that are invisible in the classical hydrodynamic description; they are not implied by the second law or the Onsager relations. We term these conditions Schwinger-Keldysh positivity and provide explicit examples of the various allowed terms.


title: Eigenfunction distribution for the Rosenzweig-Porter model

abstract: The statistical distribution of eigenfunctions for the Rosenzweig-Porter model is derived for the region where eigenfunctions have fractal behaviour. The result is based on simple physical ideas and leads to transparent explicit formulas which agree very well with numerical calculations. It constitutes a rare case where a non-trivial eigenfunction distribution is obtained in a closed form.


title: 3D G-CNNs for Pulmonary Nodule Detection

abstract: Convolutional Neural Networks (CNNs) require a large amount of annotated data to learn from, which is often difficult to obtain in the medical domain. In this paper we show that the sample complexity of CNNs can be significantly improved by using 3D roto-translation group convolutions (G-Convs) instead of the more conventional translational convolutions. These 3D G-CNNs were applied to the problem of false positive reduction for pulmonary nodule detection, and proved to be substantially more effective in terms of performance, sensitivity to malignant nodules, and speed of convergence compared to a strong and comparable baseline architecture with regular convolutions, data augmentation and a similar number of parameters. For every dataset size tested, the G-CNN achieved a FROC score close to the CNN trained on ten times more data.


title: Galois Theory - a first course

abstract: These notes are a self-contained introduction to Galois theory, designed for the student who has done a first course in abstract algebra.


title: Francis Guthrie's approach to The Four Color Problem

abstract: The odd wheel is the only type of 4-critical graph in which one vertex always gets a unique color. This supports Frederic Guthrie's approach to the Four Color Problem.


title: Asynch-SGBDT: Asynchronous Parallel Stochastic Gradient Boosting Decision Tree based on Parameters Server

abstract: In AI research and industry, machine learning is the most widely used tool. One of the most important machine learning algorithms is Gradient Boosting Decision Tree, i.e. GBDT whose training process needs considerable computational resources and time. To shorten GBDT training time, many works tried to apply GBDT on Parameter Server. However, those GBDT algorithms are synchronous parallel algorithms which fail to make full use of Parameter Server. In this paper, we examine the possibility of using asynchronous parallel methods to train GBDT model and name this algorithm as asynch-SGBDT (asynchronous parallel stochastic gradient boosting decision tree). Our theoretical and experimental results indicate that the scalability of asynch-SGBDT is influenced by the sample diversity of datasets, sampling rate, step length and the setting of GBDT tree. Experimental results also show asynch-SGBDT training process reaches a linear speedup in asynchronous parallel manner when datasets and GBDT trees meet high scalability requirements.


title: Multiple Superadiabatic Transitions and Landau-Zener Formulas

abstract: We consider nonadiabatic systems in which the classical Born-Oppenheimer approximation breaks down. We present a general theory that accurately captures the full transmitted wavepacket after multiple transitions through either a single or distinct avoided crossings, including phase information and associated interference effects. Under suitable approximations, we recover both the celebrated Landau-Zener formula and standard surface-hopping algorithms. Our algorithm shows excellent agreement with the full quantum dynamics for a range of avoided crossing systems, and can also be applied to single full crossings with similar accuracy.

categories: physics.chem-ph math-ph math.MP quant-ph


title: Non-adiabatic transitions in multiple dimensions

abstract: We consider non-adiabatic transitions in multiple dimensions, which occur when the Born-Oppenheimer approximation breaks down. We present a general, multi-dimensional algorithm which can be used to accurately and efficiently compute the transmitted wavepacket at an avoided crossing. The algorithm requires only one-level Born-Oppenheimer dynamics and local knowledge of the potential surfaces. Crucially, in contrast to standard methods in the literature, we compute the whole wavepacket, including its phase, rather than simply the transition probability. We demonstrate the excellent agreement with full quantum dynamics for a range of examples in two dimensions. We also demonstrate surprisingly good agreement for a system with a full conical intersection.


title: Binary companions of evolved stars in APOGEE DR14: Search method and catalog of ~5,000 companions

abstract: Multi-epoch radial velocity measurements of stars can be used to identify stellar, sub-stellar, and planetary-mass companions. Even a small number of observation epochs can be informative about companions, though there can be multiple qualitatively different orbital solutions that fit the data. We have custom-built a Monte Carlo sampler (The Joker) that delivers reliable (and often highly multi-modal) posterior samplings for companion orbital parameters given sparse radial-velocity data. Here we use The Joker to perform a search for companions to 96,231 red-giant stars observed in the APOGEE survey (DR14) with $\geq 3$ spectroscopic epochs. We select stars with probable companions by making a cut on our posterior belief about the amplitude of the stellar radial-velocity variation induced by the orbit. We provide (1) a catalog of 320 companions for which the stellar companion properties can be confidently determined, (2) a catalog of 4,898 stars that likely have companions, but would require more observations to uniquely determine the orbital properties, and (3) posterior samplings for the full orbital parameters for all stars in the parent sample. We show the characteristics of systems with confidently determined companion properties and highlight interesting systems with candidate compact object companions.


title: The angular momentum-mass relation: a fundamental law from dwarf irregulars to massive spirals

abstract: In a $\Lambda$CDM Universe, the specific stellar angular momentum ($j_\ast$) and stellar mass ($M_\ast$) of a galaxy are correlated as a consequence of the scaling existing for dark matter haloes ($j_{\rm h}\propto M_{\rm h}^{2/3}$). The shape of this law is crucial to test galaxy formation models, which are currently discrepant especially at the lowest masses, allowing to constrain fundamental parameters, e.g. the retained fraction of angular momentum. In this study, we accurately determine the empirical $j_\ast-M_\ast$ relation (Fall relation) for 92 nearby spiral galaxies (from S0 to Irr) selected from the Spitzer Photometry and Accurate Rotation Curves (SPARC) sample in the unprecedented mass range $7 \lesssim \log M_\ast/M_\odot \lesssim 11.5$. We significantly improve all previous estimates of the Fall relation by determining $j_\ast$ profiles homogeneously for all galaxies, using extended HI rotation curves, and selecting only galaxies for which a robust $j_\ast$ could be measured (converged $j_\ast(<R)$ radial profile). We find the relation to be well described by a single, unbroken power-law $j_\ast\propto M_\ast^\alpha$ over the entire mass range, with $\alpha=0.55\pm 0.02$ and orthogonal intrinsic scatter of $0.17\pm 0.01$ dex. We finally discuss some implications for galaxy formation models of this fundamental scaling law and, in particular, the fact that it excludes models in which discs of all masses retain the same fraction of the halo angular momentum.


title: Modeling UV Radiation Feedback from Massive Stars: II. Dispersal of Star-Forming Giant Molecular Clouds by Photoionization and Radiation Pressure

abstract: UV radiation feedback from young massive stars plays a key role in the evolution of giant molecular clouds (GMCs) by photoevaporating and ejecting the surrounding gas. We conduct a suite of radiation hydrodynamic simulations of star cluster formation in marginally-bound, turbulent GMCs, focusing on the effects of photoionization and radiation pressure on regulating the net star formation efficiency (SFE) and cloud lifetime. We find that the net SFE depends primarily on the initial gas surface density, $\Sigma_0$, such that the SFE increases from 4% to 51% as $\Sigma_0$ increases from $13\,M_{\odot}\,{\rm pc}^{-2}$ to $1300\,M_{\odot}\,{\rm pc}^{-2}$. Cloud destruction occurs within $2$-$10\,{\rm Myr}$ after the onset of radiation feedback, or within $0.6$-$4.1$ freefall times (increasing with $\Sigma_0$). Photoevaporation dominates the mass loss in massive, low surface-density clouds, but because most photons are absorbed in an ionization-bounded Str\"{o}mgren volume the photoevaporated gas fraction is proportional to the square root of the SFE. The measured momentum injection due to thermal and radiation pressure forces is proportional to $\Sigma_0^{-0.74}$, and the ejection of neutrals substantially contributes to the disruption of low-mass and/or high-surface density clouds. We present semi-analytic models for cloud dispersal mediated by photoevaporation and by dynamical mass ejection, and show that the predicted net SFE and mass loss efficiencies are consistent with the results of our numerical simulations.


title: First Subleading Power Resummation for Event Shapes

abstract: We derive and analytically solve renormalization group (RG) equations of gauge invariant non-local Wilson line operators which resum logarithms for event shape observables $\tau$ at subleading power in the $\tau\ll 1$ expansion. These equations involve a class of universal jet and soft functions arising through operator mixing, which we call $\theta$-jet and $\theta$-soft functions. An illustrative example involving these operators is introduced which captures the generic features of subleading power resummation, allowing us to derive the structure of the RG to all orders in $\alpha_s$, and provide field theory definitions of all ingredients. As a simple application, we use this to obtain an analytic leading logarithmic result for the subleading power resummed thrust spectrum for $H\to gg$ in pure glue QCD. This resummation determines the nature of the double logarithmic series at subleading power, which we find is still governed by the cusp anomalous dimension. We check our result by performing an analytic calculation up to ${\cal O}(\alpha_s^3)$. Consistency of the subleading power RG relates subleading power anomalous dimensions, constrains the form of the $\theta$-soft and $\theta$-jet functions, and implies an exponentiation of higher order loop corrections in the subleading power collinear limit. Our results provide a path for carrying out systematic resummation at subleading power for collider observables.


title: C IV BAL disappearance in a large SDSS QSO sample

abstract: Broad absorption lines (BALs) in the spectra of quasi-stellar objects (QSOs) originate from outflowing winds along our line of sight; winds are thought to originate from the inner regions of the QSO accretion disk, close to the central supermassive black hole (SMBH). Winds likely play a role in galaxy evolution and aid the accretion mechanism onto the SMBH. BAL equivalent widths can change on typical timescales from months to years; such variability is generally attributed to changes in the covering factor and/or in the ionization level of the gas. We investigate BAL variability, focusing on BAL disappearance. We analyze multi-epoch spectra of more than 1500 QSOs -the largest sample ever used for such a study- observed by different programs from the Sloan Digital Sky Survey-I/II/III (SDSS), and search for disappearing C IV BALs. The spectra rest-frame time baseline ranges from 0.28 to 4.9 yr; the source redshifts range from 1.68 to 4.27. We detect 73 disappearing BALs in the spectra of 67 sources. This corresponds to 3.9% of disappearing BALs, and 5.1% of our BAL QSOs exhibit at least one disappearing BAL. We estimate the average lifetime of a BAL along our line of sight (~ 80-100 yr), which appears consistent with the accretion disk orbital time at distances where winds are thought to originate. We inspect properties of the disappearing BALs and compare them to the properties of our main sample. We also investigate the existence of a correlation in the variability of multiple troughs in the same spectrum, and find it persistent at large velocity offsets between BAL pairs, suggesting that a mechanism extending on a global scale is necessary to explain the phenomenon. We select a more reliable sample of disappearing BALs following Filiz Ak et al. (2012), where a subset of our sample was analyzed, and compare the findings from the two works, obtaining generally consistent results.


title: NIHAO XV: The environmental impact of the host galaxy on galactic satellite and field dwarf galaxies

abstract: We study the impact of the host on dwarf galaxy properties using four new Milky Way-like, ultra high-resolution simulations, ($N_{\rm part} >10^7$) from the NIHAO project. We split our sample in satellite ($R<R_{200}$), nearby ($1<R/R_{200}<2.5$), and field ($R>2.5 R_{\rm 200}$) galaxies. Simulated galaxies from all three groups are in excellent agreement with Local Group dwarf galaxies in terms of: stellar mass-velocity dispersion, stellar mass-metallicity relation, star formation histories, and stellar mass functions. Satellites and nearby galaxies show lower velocity dispersions and gas fractions compared to field galaxies. While field galaxies follow global abundance matching relations, satellites and nearby galaxies deviate from them, showing lower dark matter masses for given stellar mass. The reason for this deficit in dark matter mass is substantial mass loss experienced by satellites and $\sim80$\% of the nearby galaxies, while orbiting inside $R_{200}$ at earlier times. However, both satellites and nearby objects fall back onto the relation for field galaxies if we use the maximum of their virial mass instead of the present-day value. This allows us to provide estimates for the peak masses of observed Local Group galaxies. Finally, using radial velocities, distances, and the velocity dispersion-stellar mass relation from our simulations, we derive a metric to distinguish between galaxies harassed by the central object and unaffected ones. Applying this metric to observed objects we find that even far away dwarf galaxies like Eri II ($D\approx$ 370 kpc) have a strong probability ($\approx 83$\%) of having been affected by the Milky Way in the past. This naturally explains the lack of gas and recent star formation seen in Eri II.


title: X-ray absorption towards high-redshift sources: probing the intergalactic medium with blazars

abstract: The role played by the intergalactic medium (IGM) in the X-ray absorption towards high-redshift sources has recently drawn more attention in spectral analysis studies. Here, we study the X-ray absorption towards 15 flat-spectrum radio quasars at $z>2$, relying on high counting statistic ($\gtrsim10\,000$ photons) provided by XMM-Newton, with additional NuSTAR (and simultaneous Swift-XRT) observations when available. Blazars can be confidently considered to have negligible X-ray absorption along the line of sight within the host galaxy, likely swept by the kpc-scale relativistic jet. This makes our sources ideal for testing the absorption component along the IGM. Our new approach is to revisit the origin of the soft X-ray spectral hardening observed in high-z blazars in terms of X-ray absorption occurring along the IGM, with the help of a low-z sample used as comparison. We verify that the presence of absorption in excess of the Galactic value is the preferred explanation to explain the observed hardening, while intrinsic energy breaks, predicted by blazars' emission models, can easily occur out of the observing energy band in most sources. First, we perform an indirect analysis comparing the inferred amount of absorption in excess of the Galactic value with a simulated IGM absorption contribution, that increases with redshift and includes both a minimum component from diffuse IGM metals, and the additional contribution of discrete denser intervening regions. Then, we directly investigate the warm-hot IGM with a spectral model on the best candidates of our sample, obtaining an average IGM density of $n_0=1.01^{+0.53}_{-0.72}\times10^{-7}$ cm$^{-3}$ and temperature of $\log(T/\text{K})=6.45^{+0.51}_{-2.12}$. A more dedicated study is currently beyond reach, but our results can be used as a stepping stone for future more accurate analysis, involving Athena.


title: Convex resource theory of non-Gaussianity

abstract: Continuous-variable systems realized in quantum optics play a major role in quantum information processing, and it is also one of the promising candidates for a scalable quantum computer. We introduce a resource theory for continuous-variable systems relevant to universal quantum computation. In our theory, easily implementable operations---Gaussian operations combined with feed-forward---are chosen to be the free operations, making the convex hull of the Gaussian states the natural free states. Since our free operations and free states cannot perform universal quantum computation, genuine non-Gaussian states---states not in the convex hull of Gaussian states---are the necessary resource states for universal quantum computation together with free operations. We introduce a monotone to quantify the genuine non-Gaussianity of resource states, in analogy to the stabilizer theory. A direct application of our resource theory is to bound the conversion rate between genuine non-Gaussian states. Finally, we give a protocol that probabilistically distills genuine non-Gaussianity---increases the genuine non-Gaussianity of resource states---only using free operations and postselection on Gaussian measurements, where our theory gives an upper bound for the distillation rate. In particular, the same protocol allows the distillation of cubic phase states, which enable universal quantum computation when combined with free operations.


title: General purpose ray-tracing and polarized radiative transfer in General Relativity

abstract: Ray-tracing is a central tool for constructing mock observations of compact object emission and for comparing physical emission models with observations. We present Arcmancer, a publicly available general ray-tracing and tensor algebra library, written in C++ and providing a Python interface. Arcmancer supports Riemannian and semi-Riemannian spaces of any dimension and metric, and has novel features such as support for multiple simultaneous coordinate charts, embedded geometric shapes, local coordinate systems and automatic parallel propagation. The Arcmancer interface is extensively documented and user-friendly. While these capabilities make the library well suited for a large variety of problems in numerical geometry, the main focus of this paper is in general relativistic polarized radiative transfer. The accuracy of the code is demonstrated in several code tests and in a comparison with GRTRANS, an existing ray-tracing code. We then use the library in several scenarios as a way to showcase the wide applicability of the code. We study a thin variable-geometry accretion disk model, and find that polarization carries information of the inner disk opening angle. Next, we study rotating neutron stars and determine that to obtain polarized light curves at better than $\sim1\%$ level of accuracy, the rotation needs to be taken into account both in the space-time metric as well as in the shape of the star. Finally, we investigate the observational signatures of an accreting black hole lensed by an orbiting black hole. We find that these systems exhibit a characteristic asymmetric twin-peak profile both in flux and polarization properties.


title: Ghost-free completion of an effective matter coupling in bimetric theory

abstract: We consider a particular set of ghost-free interactions for three spin-2 fields in which we freeze out the dynamics of the metric tensor that couples to the matter sector. Integrating out the non-dynamical degrees of freedom in vacuum results in a ghost-free bimetric theory. In the presence of the matter source, which we treat as a small perturbation, the equations for the non-dynamical field can be solved perturbatively for its vierbein. This results in ghost-free bimetric theory in vierbein formulation with a modified matter coupling. To lowest order in matter perturbations, we precisely obtain an effective matter coupling that has been suggested earlier in the literature. This coupling contains a linear combination of the two vierbeine whose corresponding metric fluctuation coincides with the massless spin-2 mode. In the past, bimetric theory with this symmetric coupling has been treated as an effective theory valid at low energies without a ghost-free completion. Our results demonstrate that the effective matter coupling can be rendered entirely ghost-free by including the higher-order corrections obtained from the trimetric setup.


title: Non-Hermitian Chern bands

abstract: The relation between chiral edge modes and bulk Chern numbers of quantum Hall insulators is a paradigmatic example of bulk-boundary correspondence. We show that the chiral edge modes are not strictly tied to the Chern numbers defined by a non-Hermitian Bloch Hamiltonian. This breakdown of conventional bulk-boundary correspondence stems from the non-Bloch-wave behavior of eigenstates (non-Hermitian skin effect), which generates pronounced deviations of phase diagrams from the Bloch theory. We introduce non-Bloch Chern numbers that faithfully predict the numbers of chiral edge modes. The theory is backed up by the open-boundary energy spectra, dynamics, and phase diagram of representative lattice models. Our results highlight a unique feature of non-Hermitian bands and suggest a non-Bloch framework to characterize their topology.

categories: cond-mat.mes-hall cond-mat.other cond-mat.quant-gas cond-mat.str-el physics.optics


title: ASASSN-18di: discovery of a $\Delta V \sim 10$ flare on a mid-M dwarf

abstract: We report and characterize a white-light superflare on a previously undiscovered M dwarf detected by the ASAS-SN survey. Employing various color-magnitude and color-spectral type relationships, we estimate several stellar parameters, including the quiescent V-band magnitude, from which we derive a flare amplitude of $\Delta V \sim 10$. We determine an r-band absolute magnitude of $M_{r} = 11.4$, consistent with a mid-M dwarf, and an approximate distance to the source of $2.2$ kpc. Using classical-flare models, we infer a flare energy of $E_{V} \simeq (4.1\pm 2.2)\times 10^{36}$ ergs, making this one of the strongest flares documented on an M dwarf.


title: Ballistic one-dimensional holes with strong g-factor anisotropy in germanium

abstract: We report experimental evidence of ballistic hole transport in one-dimensional quantum wires gate-defined in a strained SiGe/Ge/SiGe quantum well. At zero magnetic field, we observe conductance plateaus at integer multiples of 2e^2/h. At finite magnetic field, the splitting of these plateaus by Zeeman effect reveals largely anisotropic g-factors, with absolute values below 1 in the quantum-well plane, and exceeding 10 out of plane. This g-factor anisotropy is consistent with a heavy-hole character of the propagating valence-band states, in line with a predominant confinement in the growth direction. Remarkably, we observe quantized ballistic conductance in device channels up to 600 nm long. These findings mark an important step towards the realization of novel devices for applications in quantum spintronics.


title: An ALMA 3mm continuum census of Westerlund 1

abstract: Massive stars play an important role in both cluster and galactic evolution and the rate at which they lose mass is a key driver of both their own evolution and their interaction with the environment up to and including their SNe explosions. Young massive clusters provide an ideal opportunity to study a co-eval population of massive stars. We performed 3mm continuum observations with the Atacama Large Millimetre/submillimetre Array of the Galactic cluster Westerlund 1, to study the constituent massive stars and determine mass-loss rates for the diverse post-main sequence population. We detected emission from 50 stars in Westerlund 1, comprising all 21 Wolf-Rayets within the field of view, eight cool and 21 OB super-/hypergiants. Emission nebulae were associated with a number of the cool hypergiants while, unexpectedly, a number of hot stars also appear spatially resolved. We measured the mass-loss rates for a unique population of massive post-main sequence stars at every stage of evolution, confirming a significant increase as stars transition from OB supergiant to WR states. The range of spectral types exhibited provides a critical test of radiatively driven wind theory and the reality of the bi-stability jump. The extreme mass-loss rate inferred for the interacting binary Wd1-9 in comparison to other cluster members confirmed the key role binarity plays in massive stellar evolution. The presence of compact nebulae around a number of OB and WR stars is unexpected; by analogy to the cool super-/hypergiants we attribute this to confinement and sculpting of the stellar wind via interaction with the intra-cluster medium/wind. Given the morphology of core collapse SNe depend on the nature of the pre-explosion circumstellar environment, if this hypothesis is correct then the properties of the explosion depend not just on the progenitor, but also the environment in which it is located.


title: Topological unification of time-reversal and particle-hole symmetries in non-Hermitian physics

abstract: Topological phases are enriched in non-equilibrium open systems effectively described by non-Hermitian Hamiltonians. While several properties unique to non-Hermitian topological systems were uncovered, the fundamental role of symmetry in non-Hermitian physics has yet to be fully understood, and it has remained unclear how symmetry protects non-Hermitian topological phases. Here we show that two fundamental anti-unitary symmetries, time-reversal and particle-hole symmetries, are topologically equivalent in the complex energy plane and hence unified in non-Hermitian physics. A striking consequence of this symmetry unification is the emergence of unique non-equilibrium topological phases that have no counterparts in Hermitian systems. We illustrate this by presenting a non-Hermitian counterpart of the Majorana chain in an insulator with time-reversal symmetry and that of the quantum spin Hall insulator in a superconductor with particle-hole symmetry. Our work establishes a fundamental symmetry principle in non-Hermitian physics and paves the way towards a unified framework for non-equilibrium topological phases.

categories: quant-ph cond-mat.mes-hall cond-mat.quant-gas math-ph math.MP physics.optics


title: The Italian contributions to cosmic-ray physics from Bruno Rossi to the G-Stack. A new window into the inexhaustible wealth of nature

abstract: From the late 1920s to the early 1950s, cosmic rays were the main instrument to investigate what we now call "high-energy physics". In approximately 25 years, an intense experimental and theoretical work brought particle physics from its childhood to its maturity. The data collected at that time played a crucial role in order to outline the quantum theories of fundamental interactions - electromagnetic, weak and strong - and they were an excellent training field for many young scientists. The knowledge acquired on cosmic rays constituted the basis of the extraordinary progresses that were to be achieved in fundamental physics with the introduction of particles accelerators. Many papers on this subject have already been published, both by historians and by the scientists who personally worked at that time. However, a complete historical reconstruction of those 25 years of cosmic-ray researches has not been written yet. In particular, an overall description of the researches carried out with the contributions of Italian physicists is still missing. The present paper tries to outline such a description, taking into account the fact that in those very years, for scientific and non scientific reasons, physics researches were getting organised within international collaborations. Without pretending to be exhaustive at all, we would like this paper to be a stimulus for further researches on the subject.


title: Fast approaches for Bayesian estimation of size of hard-to-reach populations using Network Scale-up

abstract: The Network scale-up method is commonly used to overcome difficulties in estimating the size of hard-to-reach populations. The method uses indirect information based on social network of each participant taken from the general population, but in some applications a fast computational approach would be highly recommended. We propose a Gibbs sampling method and a Monte Carlo approach to sample from the random degree model. We applied the abovementioned analytical strategies to previous data on heavy drug users from Curitiba, Brazil.


title: On the theory of excitonic delocalization for robust vibronic dynamics in LH2

abstract: Nonlinear spectroscopy has revealed long-lasting oscillations in the optical response of a variety of photosynthetic complexes. Different theoretical models which involve the coherent coupling of electronic (excitonic) or electronic-vibrational (vibronic) degrees of freedom have been put forward to explain these observations. The ensuing debate concerning the relevance of either one or the other mechanism may have obscured their potential synergy. To illustrate this synergy, we quantify how the excitonic delocalization in the LH2 unit of Rhodopseudomonas Acidophila purple bacterium, leads to correlations of excitonic energy fluctuations, relevant coherent vibronic coupling and, importantly, a decrease in the excitonic dephasing rates. Combining these effects, we identify a feasible origin for the long-lasting oscillations observed in fluorescent traces from time-delayed two-pulse single molecule experiments performed on this photosynthetic complex.


title: Smoothness Conditions in Cohomogeneity manifolds

abstract: In this paper we discuss the smoothness conditions for metrics on a cohomogeneity one manifold, i.e. metrics invariant under a Lie group whose generic orbits are hypersurfaces. Along these hypersurfaces one describes the metrics in terms of a collection of functions defined along a geodesic normal to the hypersurfaces. In a neighborhood of a lower dimensional orbit the functions must satisfy certain smoothness conditions for the metric to extend smoothly. In general these can be quite complicated. We present a method that makes it straightforward to compute them, and illustrate it in several examples. This second version contains some improvements in exposition and a reformulation of Theorem B, with a proof added.


title: BPS Solutions of Six-Dimensional (1,0) Supergravity Coupled to Tensor Multiplets

abstract: We derive a general local form for supersymmetric solutions of six-dimensional (1,0) supergravity coupled to an arbitrary number of tensor multiplets. We consider some special cases in which the resulting equations can be solved explicitly. In particular we derive black string solutions and calculate their entropy. Upon reducing to five dimensions they yield spinning black hole solutions. We also discuss BPS pp-waves and black string solutions with traveling waves. Lastly, as an application, we study the attractor mechanism in this theory.


title: Concurrent formation of supermassive stars and globular clusters: implications for early self-enrichment

abstract: We present a model for the concurrent formation of globular clusters (GCs) and supermassive stars (SMSs, $>10^3\,{\rm M}_\odot$) to address the origin of the HeCNONaMgAl abundance anomalies in GCs. GCs form in converging gas flows and accumulate low-angular momentum gas, which accretes onto protostars. This leads to an adiabatic contraction of the cluster and an increase of the stellar collision rate. A SMS can form via runaway collisions if the cluster reaches sufficiently high density before two-body relaxation halts the contraction. This condition is met if the number of stars $\gtrsim10^6$ and the gas accretion rate $\gtrsim10^5\,{\rm M}_\odot$/Myr, reminiscent of GC formation in high gas-density environments, such as -- but not restricted to -- the early Universe. The strong SMS wind mixes with the inflowing pristine gas, such that the protostars accrete diluted hot-hydrogen burning yields of the SMS. Because of continuous rejuvenation, the amount of processed material liberated by the SMS can be an order of magnitude higher than its maximum mass. This `conveyor-belt' production of hot-hydrogen burning products provides a solution to the mass budget problem that plagues other scenarios. Additionally, the liberated material is mildly enriched in helium and relatively rich in other hot-hydrogen burning products, in agreement with abundances of GCs today. Finally, we find a super-linear scaling between the amount of processed material and cluster mass, providing an explanation for the observed increase of the fraction of processed material with GC mass. We discuss open questions of this new GC enrichment scenario and propose observational tests.


title: Bounds on the largest Kronecker and induced multiplicities of finite groups

abstract: We give new bounds and asymptotic estimates on the largest Kronecker and induced multiplicities of finite groups. The results apply to large simple groups of Lie type and other groups with few conjugacy classes.


title: Greybody factors for a minimally coupled scalar field in three-dimensional Einstein-power-Maxwell black hole background

abstract: In the present work we study the propagation of a probe minimally coupled scalar field in Einstein-power-Maxwell charged black hole background in (1+2) dimensions. We find analytical expressions for the reflection coefficient as well as for the absorption cross-section in the low energy regime, and we show graphically their behaviour as functions of the frequency for several values of the free parameters of the theory.


title: Timing Recollision in Nonsequential Double Ionization by Intense Elliptically Polarized Laser Pulses

abstract: We examine correlated electron and doubly charged ion momentum spectra from strong field double ionization of Neon employing intense elliptically polarized laser pulses. An ellipticity-dependent asymmetry of correlated electron and ion momentum distributions has been observed. Using a 3D semiclassical model, we demonstrate that our observations reflect the sub-cycle dynamics of the recollision process. Our work reveals a general physical picture for recollision-impact double ionization with elliptical polarization, and demonstrates the possibility of ultrafast control of the recollision dynamics.


title: Separation Of Edge Effects In Highly Non-Dipole Radiation Spectra

abstract: Edge effects in spectra of non-dipole radiation from ultra-relativistic electrons passing through various finite targets are analyzed from a unified point of view. Examples include radiation in a finite magnet, bremsstrahlung at double scattering, and bremsstrahlung in an amorphous plate. A generalization of the encountered spectral decomposition property is proposed. The relevance of phenomena of jet-interjet radiation interference and electron time delay in the target is emphasized. Features produced by them in the radiation spectra are discussed.


title: Cross-Domain Visual Recognition via Domain Adaptive Dictionary Learning

abstract: In real-world visual recognition problems, the assumption that the training data (source domain) and test data (target domain) are sampled from the same distribution is often violated. This is known as the domain adaptation problem. In this work, we propose a novel domain-adaptive dictionary learning framework for cross-domain visual recognition. Our method generates a set of intermediate domains. These intermediate domains form a smooth path and bridge the gap between the source and target domains. Specifically, we not only learn a common dictionary to encode the domain-shared features, but also learn a set of domain-specific dictionaries to model the domain shift. The separation of the common and domain-specific dictionaries enables us to learn more compact and reconstructive dictionaries for domain adaptation. These dictionaries are learned by alternating between domain-adaptive sparse coding and dictionary updating steps. Meanwhile, our approach gradually recovers the feature representations of both source and target data along the domain path. By aligning all the recovered domain data, we derive the final domain-adaptive features for cross-domain visual recognition. Extensive experiments on three public datasets demonstrates that our approach outperforms most state-of-the-art methods.


title: Topological data analysis and diagnostics of compressible MHD turbulence

abstract: The predictions of mean-field electrodynamics can now be probed using direct numerical simulations of random flows and magnetic fields. When modelling astrophysical MHD, it is important to verify that such simulations are in agreement with observations. One of the main challenges in this area is to identify robust \it{quantitative} measures to compare structures found in simulations with those inferred from astrophysical observations. A similar challenge is to compare quantitatively results from different simulations. Topological data analysis offers a range of techniques, including the Betti numbers and persistence diagrams, that can be used to facilitate such a comparison. After describing these tools, we first apply them to synthetic random fields and demonstrate that, when the data are standardized in a straightforward manner, some topological measures are insensitive to either large-scale trends or the resolution of the data. Focusing upon one particular astrophysical example, we apply topological data analysis to HI observations of the turbulent interstellar medium (ISM) in the Milky Way and to recent MHD simulations of the random, strongly compressible ISM. We stress that these topological techniques are generic and could be applied to any complex, multi-dimensional random field.


title: The VMC survey XXX. Stellar proper motions in the central parts of the Small Magellanic Cloud

abstract: We present the first spatially resolved map of stellar proper motions within the central ($\sim$3.1 $\times$ 2.4 kpc) regions of the Small Magellanic Cloud (SMC). The data used for this study encompasses four tiles from the ongoing near-infrared VISTA survey of the Magellanic Clouds system and covers a total contiguous area on the sky of $\sim$6.81 deg$^2$. Proper motions have been calculated independently in two dimensions from the spatial offsets in the $K_s$ filter over time baselines between 22 and 27 months. The reflex motions of $\sim$33~000 background galaxies are used to calibrate the stellar motions to an absolute scale. The resulting catalog is composed of more than 690 000 stars which have been selected based on their position in the $(J-K_s, K_s)$ color-magnitude diagram. For the median absolute proper motion of the SMC, we find ($\mu_{\alpha}\mathrm{cos}(\delta)$, $\mu_{\delta}$) = (1.087 $\pm$ 0.192 (sys.) $\pm$ 0.003 (stat.), $-$1.187 $\pm$ 0.008 (sys.) $\pm$ 0.003 (stat.)) mas yr$^{-1}$, consistent with previous studies. Mapping the proper motions as a function of position within the SMC reveals a non uniform velocity pattern indicative of a tidal feature behind the main body of the SMC and a flow of stars in the South-East moving predominantly along the line-of-sight.


title: Aesthetical Attributes for Segmenting Arabic Word

abstract: The connected allograph representing calligraphic Arabic word does not appear individually in any calligraphic resource but in association with other letters all adapted to each other. The graphic segmentation of the word by respecting aesthetical attributes indicating the grapheme of every letter is far from being an obvious task. The question consists in discovering every letter constituting the word, points of cutting which separate its grapheme from other constituents of word's shape. The obtained segment must be a complete drawing of the represented letter. This segmentation according to contextual graphic and qualitative criteria connecting the attached allograph will have to satisfy typographic constraints varying in conformity with the possibilities offered by the wanted technology. In this paper, we develop an approach for segmenting Arabic word from which the purpose is to extract graphemes respecting the design of Arabic letters such as it is in the calligraphic literature. The procedure bases itself on the principle that the Arabic connected letters have a common part included in the cursive area, which must not be lost during the process of cutting.


title: Problem of Multiple Diacritics Design for Arabic Script

abstract: This study focuses on the design of multiple Arabic diacritical marks and to developing a model that generates the stacking of multiples Arabic diacritics in order to integrate it into a system of Arabic composition. The problem concerns the presence of multiple diacritics on a single basic letter. This model is based on the layering composition. The combination of diacritics with letters requires a basic layering to combine any diacritics in the word with their base letter, without having to deal individually and separately each pair of base letter and diacritics.


title: Determining the Absolute Magnitudes of Galactic-Bulge Red Clump Giants in the Z and Y Filters of the Vista Sky Surveys and the IRAC Filters of the Spitzer Sky Surveys

abstract: The properties of red clump giants in the central regions of the Galactic bulge are investigated in the photometric $Z$ and $Y$ bands of the infrared VVV (VISTA/ESO) survey and the [3.6],[4.5], [5.8] and [8.0] bands of the GLIMPSE (Spitzer/IRAC) Galactic plane survey. The absolute magnitudes for objects of this class have been determined in these bands for the first time: $M_{Z}= -0.20\pm0.04$ and $M_{Y}=-0.470\pm0.045$, $M_{[3.6]}= -1.70\pm0.03$, $M_{[4.5]}= -1.60\pm0.03$, $M_{[5.8]}= -1.67\pm0.03$ and $M_{[8.0]}= -1.70\pm0.03$. A comparison of the measured magnitudes with the predictions of theoretical models for the spectra of the objects under study has demonstrated good mutual agreement and has allowed some important constraints to be obtained for the properties of bulge red clump giants. In particular, a comparison with evolutionary tracks has shown that we are dealing predominantly with the high-metallicity subgroup of bulge red clump giants. Their metallicity is slightly higher than has been thought previously, $[M/H]\simeq0.40$ (${\bf Z} \simeq 0.038$) with an error of $[M/H]\simeq 0.1$ dex, while the effective temperature is $4250\pm150$ K. Stars with an age of 9-10 Gyr are shown to dominate among the red clump giants, although some number of younger objects with an age of $\sim8$ can also be present. In addition, the distances to several Galactic bulge regions have been measured, as D = 8200-8500 pc, and the extinction law in these directions is shown to differ noticeably from the standard one.


title: On the largest Kronecker and Littlewood--Richardson coefficients

abstract: We give new bounds and asymptotic estimates for Kronecker and Littlewood--Richardson coefficients. Notably, we resolve Stanley's questions on the shape of partitions attaining the largest Kronecker and Littlewood--Richardson coefficients. We apply the results to asymptotics of the number of standard Young tableaux of skew shapes.


title: A Variational U-Net for Conditional Appearance and Shape Generation

abstract: Deep generative models have demonstrated great performance in image synthesis. However, results deteriorate in case of spatial deformations, since they generate images of objects directly, rather than modeling the intricate interplay of their inherent shape and appearance. We present a conditional U-Net for shape-guided image generation, conditioned on the output of a variational autoencoder for appearance. The approach is trained end-to-end on images, without requiring samples of the same object with varying pose or appearance. Experiments show that the model enables conditional image generation and transfer. Therefore, either shape or appearance can be retained from a query image, while freely altering the other. Moreover, appearance can be sampled due to its stochastic latent representation, while preserving shape. In quantitative and qualitative experiments on COCO, DeepFashion, shoes, Market-1501 and handbags, the approach demonstrates significant improvements over the state-of-the-art.


title: Quantum corrections to conductivity of disordered electrons due to inelastic scattering off magnetic impurities

abstract: We study the quantum corrections to the conductivity of the two-dimensional disordered interacting electron system in the diffusive regime due to inelastic scattering off rare magnetic impurities. We focus on the case of very different g-factors for electrons and magnetic impurities. Within the Born approximation for the inelastic scattering off magnetic impurities we find additional temperature-dependent corrections to the conductivity of the Altshuler-Aronov type. Our results demonstrate that the low temperature transport in interacting disordered electron systems with rare magnetic impurities is more interesting than it was commonly believed on the basis of treatment of magnetic impurity spins as classical ones.


title: Efficient Model Identification for Tensegrity Locomotion

abstract: This paper aims to identify in a practical manner unknown physical parameters, such as mechanical models of actuated robot links, which are critical in dynamical robotic tasks. Key features include the use of an off-the-shelf physics engine and the Bayesian optimization framework. The task being considered is locomotion with a high-dimensional, compliant Tensegrity robot. A key insight, in this case, is the need to project the model identification challenge into an appropriate lower dimensional space for efficiency. Comparisons with alternatives indicate that the proposed method can identify the parameters more accurately within the given time budget, which also results in more precise locomotion control.


title: Two-color interpolation of absorption response for quantitative Acousto-Optic imaging

abstract: Diffuse Optical Tomography (DOT) is a reliable and widespread technique for monitoring qualitative changes in absorption inside highly scattering media. It has been shown, however, that Acousto-Optic (AO) imaging can provide significantly more qualitative information, without the need for inversion algorithms, due to the spatial resolution afforded by ultrasound probing. In this article, we show how, by using multiple-wavelength AO imaging, it is also possible to perform quantitative measurements of absorber concentration inside scattering media.


title: Design and fabrication of low cost Disdrometer and turbulence sensor

abstract: Disdrometers detects size and speed of falling rain droplets in a rainfall. They detect droplets either by impact or optical based methods. Optical based ones, use light source and an optical detector. Droplet crossing the light creates shadow on the detector. This change in the intensity and duration of change provides the size and speed of the droplet. Whereas the impact based method utilizes peizo elements to sense the droplet by physical momentum transfer at the time of droplets impact with the sensor. Apart from the technical differences in them, both of these disdrometers fall in the similar cost range. Wind turbulence sensors works on detection of gradients in wind speeds by employing various working principles like thermal anemometer, laser Doppler anemometer, etc. In thermo anemometers, a heating element is maintained at a constant temperature. As the wind flows, it carries some heat energy away from the element. This causes a change in the resistance in the heating element providing further change in voltages across it. Laser Doppler anemometers work on change in wavelength of light scattered by a moving particle which is caught in the wind flow. These techniques are highly accurate in measuring the turbulence in the air but at elevated costs. Aimed at cost efficacy and acceptable performance of the device, a novel concept was implemented in a single device to study its feasibility to work as a disdrometer and a turbulence sensor. It employs the principles of electromagnetic induction and response of a cantilever to impulses. Based on these principles, a device was designed and fabricated to detect rain droplet size and wind turbulence. It was fabricated from readily available materials. The testing of the device in laboratory conditions yielded satisfactory results as compare to the standard commercial instruments.


title: Stein kernels and moment maps

abstract: We describe a construction of Stein kernels using moment maps, which are solutions to a variant of the Monge-Amp\`ere equation. As a consequence, we show how regularity bounds on these maps control the rate of convergence in the classical central limit theorem, and derive new rates in Kantorovitch-Wasserstein distance in the log-concave situation, with explicit polynomial dependence on the dimension.


title: An investigation of the non-trivial zeros of the Riemann zeta function

abstract: While many zeros of the Riemann zeta function are located on the critical line $\Re(s)=1/2$, the non-existence of zeros in the remaining part of the critical strip $\Re(s) \in \, ]0, 1[$ is the main scope to be proven for the Riemann hypothesis. The Riemann zeta functional leads to a relation between the zeros on either sides of the critical line. Given $s$ a complex number and $\bar{s}$ its complex conjugate, if $s$ is a zero of the Riemann zeta function in the critical strip $\Re(s) \in \, ]0, 1[$, then $\zeta(s) = \zeta(1-\bar{s})$, as a key proposition to prove the Riemann hypothesis.


title: Reputation in M2M Economy

abstract: Triggered by modern technologies, our possibilities may now expand beyond the unthinkable. Cars externally may look similar to decades ago, but a dramatic revolution happened inside the cabin as a result of their computation, communications, and storage capabilities. With the advent of Electric Autonomous Vehicles (EAVs), Artificial Intelligence and ecological technologies found the best synergy. Several transportation problems may be solved (accidents, emissions, and congestion among others), and the foundation of Machine-to-Machine (M2M) economy could be established, in addition to value-added services such as infotainment (information and entertainment). In the world where intelligent technologies are pervading everyday life, software and algorithms play a major role. Software has been lately introduced in virtually every technological product available on the market, from phones to television sets to cars and even housing. Artificial Intelligence is one of the consequences of this pervasive presence of algorithms. The role of software is becoming dominant and technology is, at times pervasive, of our existence. Concerns, such as privacy and security, demand high attention and have been already explored to some level of detail. However, intelligent agents and actors are often considered as perfect entities that will overcome human error-prone nature. This may not always be the case and we advocate that the notion of reputation is also applicable to intelligent artificial agents, in particular to EAVs.


title: Representations of monomial matrices and restriction from $GL_n$ to $S_n$

abstract: We outline an approach to understanding restrictions of polynomial representations of $GL_n(\mathbb{C})$ to $S_n$ by first restricting to $T \rtimes S_n$, the subgroup of $n \times n$ monomial matrices. Using this approach we give a combinatorial interpretation for the decomposition of a tensor product of symmetric powers of the defining representation.


title: Rost nilpotence and free theories

abstract: We introduce coherent cohomology theories h_* and prove that if such a theory is moreover generically constant then the Rost nilpotence principle holds for projective homogeneous varieties in the category of h_*-motives. Examples of such theories are algebraic cobordism and its descendants the free theories.


title: Generating Sequences and Semigroups of Valuations on 2-Dimensional Normal Local Rings

abstract: In this paper we develop a method for constructing generating sequences for valuations dominating the ring of a two dimensional quotient singularity. Suppose that $K$ is an algebraically closed field of characteristic zero, $K[X,Y]$ is a polynomial ring over $K$ and $\nu$ is a rational rank 1 valuation of the field $K(X,Y)$ which dominates $K[X,Y]_{(X,Y)}$. Given a finite Abelian group $H$ acting diagonally on $K[X,Y]$, and a generating sequence of $\nu$ in $K[X,Y]$ whose members are eigenfunctions for the action of $H$, we compute a generating sequence for the invariant ring $K[X,Y]^H$. We use this to compute the semigroup $S^{K[X,Y]^H}$ of values of elements of $K[X,Y]^H$. We further determine when $S^{K[X,Y]} (\nu)$ is a finitely generated $S^{K[X,Y]^H} (\nu)$ -module.


title: On surfaces endowed with canonical principal direction in Euclidean spaces

abstract: In this paper, we introduce canonical principal direction (CPD) submanifolds with higher codimension in Euclidean spaces. We obtain the complete classification of surfaces endowed with CPD in the Euclidean 4-space.


title: Towards a complete study of central exclusive production of $K^{+}K^{-}$ pairs in proton-proton collisions within the tensor Pomeron approach

abstract: We present a study of the central exclusive production of the $K^{+} K^{-}$ pairs in proton-proton collisions at high energies. We consider diffractive mechanisms including the $K^{+} K^{-}$ continuum, the dominant scalar $f_{0}(980)$, $f_{0}(1500)$, $f_{0}(1710)$ and tensor $f_{2}(1270)$, $f'_{2}(1525)$ resonances decaying into the $K^{+} K^{-}$ pairs. We include also photoproduction mechanisms for the non-resonant (Drell-S\"{o}ding) and the $\phi(1020)$ resonance contributions. The theoretical results are calculated within the tensor-pomeron approach including both pomeron and reggeon exchanges. Predictions for planned or current experiments at RHIC and LHC are presented. We discuss the influence of the experimental cuts on the integrated cross section and on various differential distributions for outgoing particles. The distributions in two-kaon invariant mass, in a special `glueball filter variable', as well as examples of angular distributions in the $K^{+}K^{-}$ rest frame are presented. We compare the $\phi(1020)$ and continuum photoproduction contributions to the $f_{0}(980)$ and continuum diffractive contributions and discuss whether the $\phi(1020)$ resonance could be extracted experimentally. For the determination of some model parameters we also include a discussion of $K$-nucleon scattering, in particular total cross sections, and of $\phi(1020)$ photoproduction.


title: Unveiling the chemistry of interstellar CH: Spectroscopy of the 2 THz $N=2\leftarrow 1$ ground state line

abstract: The methylidyne radical CH is commonly used as a proxy for H$_2$ in the cold, neutral phase of the interstellar medium. The optical spectroscopy of CH is limited by interstellar extinction, whereas far-infrared observations provide an integral view through the Galaxy. While the HF ground state absorption, another H$_2$ proxy in diffuse gas, frequently suffers from saturation, CH remains transparent both in spiral-arm crossings and high-mass star forming regions, turning this light hydride into a universal surrogate for H$_2$. However, in slow shocks and in regions dissipating turbulence its abundance is expected to be enhanced by an endothermic production path, and the idea of a "canonical" CH abundance needs to be addressed. The $N=2\leftarrow 1$ ground state transition of CH at $\lambda 149\,\mu$m has become accessible to high-resolution spectroscopy thanks to GREAT aboard SOFIA. Its unsaturated absorption and the absence of emission makes it an ideal candidate for the determination of column densities with a minimum of assumptions. Here we present an analysis of four sightlines towards distant, far-infrared bright Galactic star forming regions. If combined with the sub-millimeter line of CH at $\lambda 560\,\mu$m, environments forming massive stars can be analyzed. For this we present a case study on the "proto-Trapezium" cluster W3 IRS5, and demonstrate that the sub-millimeter/far-infrared lines of CH reliably trace not only diffuse but also dense, molecular gas. While we confirm the global correlation between the column densities of HF and those of CH, clear signposts of an over-abundance of CH are observed towards lower densities. A quiescent ion-neutral chemistry alone cannot account for this over-abundance. Vortices forming in turbulent, diffuse gas may be the setting for an enhanced production path.


title: RFI Flagging Implications for Short-Duration Transients

abstract: With their wide fields of view and often relatively long coverage of any position in the sky in imaging survey mode, modern radio telescopes provide a data stream that is naturally suited to searching for rare transients. However, Radio Frequency Interference (RFI) can show up in the data stream in similar ways to such transients, and thus the normal pre-treatment of filtering RFI (flagging) may also remove astrophysical transients from the data stream before imaging. In this paper we investigate how standard flagging affects the detectability of such transients by examining the case of transient detection in an observing mode used for Low Frequency Array (LOFAR; \citep{LOFAR}) surveys. We quantify the fluence range of transients that would be detected, and the reduction of their SNR due to partial flagging. We find that transients with a duration close to the integration sampling time, as well as bright transients with durations on the order of tens of seconds, are completely flagged. For longer transients on the order of several tens of seconds to minutes, the flagging effects are not as severe, although part of the signal is lost. For these transients, we present a modified flagging strategy which mitigates the effect of flagging on transient signals. We also present a script which uses the differences between the two strategies, and known differences between transient RFI and astrophysical transients, to notify the observer when a potential transient is in the data stream.


title: Rotation of the asymptotic giant branch star R Doradus

abstract: High resolution observations of the extended atmospheres of asymptotic giant branch (AGB) stars can now directly confront the theories that describe stellar mass loss. Using Atacama Large Millimeter/submillimeter Array (ALMA) high angular resolution ($30\times42$~mas) observations we have, for the first time, resolved stellar rotation of an AGB star, R~Dor. We measure an angular rotation velocity of $\omega_R\sin{i}=(3.5\pm0.3)\times10^{-9}$~rad~s$^{-1}$ which indicates a rotational velocity of $|\upsilon_{\rm rot}\sin{i}|=1.0\pm0.1$~km~s$^{-1}$ at the stellar surface ($R_*=31.2$~mas at $214$~GHz). The rotation axis projected on the plane of the sky has a position angle $\Phi=7\pm6^\circ$. We find that the rotation of R Dor is two orders of magnitude faster than expected for a solitary AGB star that will have lost most of its angular momentum. Its rotational velocity is consistent with angular momentum transfer from a close companion. As a companion has not been directly detected we thus suggest R~Dor has a low-mass, close-in, companion. The rotational velocity approaches the critical velocity, set by the local sound speed in the extended envelope, and is thus expected to affect the mass loss characteristics of R~Dor.


title: Fully Distributed Secondary Voltage Control in Inverter-Based Microgrids

abstract: Centralized secondary voltage control in a power system has been replaced by the distributed controller in the recent literature due to its high dependency on extensive communication messages. Although in the new method each distributed generator only communicate with its neighbors to control the voltage, yet the messages are circulating among the whole system. In this paper, we have utilized distributed controller locally so that it will work as a fully distributed control system. This controller has been justified by being studied within a case study including 6 distributed generators.


title: High-Temperature Majorana Corner States

abstract: Majorana bound states often occur at the end of 1D topological superconductor or at the $\pi$ Josephson junction mediated by a helical edge state. Validated by a new bulk invariant and an intuitive edge argument, we show the emergence of one Majorana Kramers pair at each corner of a square-shaped 2D topological insulator proximitized by an $s_\pm$-wave (e.g., Fe-based) superconductor. We obtain a phase diagram that emphasizes the roles of bulk parameters and edge orientations. We propose several experimental realizations in lattice-matched candidate materials. Our scheme offers a high-temperature platform for exploring higher-order non-Abelian quasiparticles.

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


title: Emergent three-dimensional sperm motility: Coupling calcium dynamics and preferred curvature in a Kirchhoff rod model

abstract: Changes in calcium concentration along the sperm flagellum regulate sperm motility and hyperactivation, characterized by an increased flagellar bend amplitude and beat asymmetry, enabling the sperm to reach and penetrate the ovum (egg). The signaling pathways by which calcium increases within the flagellum are well established. However, the exact mechanisms of how calcium regulates flagellar bending are still under investigation. We extend our previous model of planar flagellar bending by developing a fluid-structure interaction model that couples the three-dimensional motion of the flagellum in a viscous, Newtonian fluid with the evolving calcium concentration. The flagellum is modeled as a Kirchhoff rod: an elastic rod with preferred curvature and twist. The calcium dynamics are represented as a one-dimensional reaction-diffusion model on a moving domain, the centerline of the flagellum. The two models are coupled assuming that the preferred curvature and twist of the sperm flagellum depend on the local calcium concentration. To investigate the effect of calcium on sperm motility, we compare model results of flagellar bend amplitude and swimming speed for three cases: planar, helical (spiral with equal amplitude in both directions), and quasi-planar (spiral with small amplitude in one direction). We observe that for the same parameters, the planar swimmer is faster and a turning motion is more clearly observed when calcium coupling is accounted for in the model. In the case of flagellar bending coupled to the calcium concentration, we observe emergent trajectories that can be characterized as a hypotrochoid for both quasi-planar and helical bending.


title: New method of bandlimited extrapolation

abstract: The paper deals with numerical solution of the Fredholm integral equation associated with the classical problem of extrapolating bandlimited functions known on $(-1,1)$ to the entire real line. The approach presented can be characterized as the degenerate kernel method using the spherical Bessel functions as basis functions where the Tikhonov regularization is applied at the discrete level in order to deal with the ill-posedness of the problem.


title: Generalised Sobolev Stable Flux Reconstruction

abstract: A new set of symmetric correction functions is presented for high-order flux reconstruction, that expands upon, while incorporating, all previous correction function sets and opens the possibility for improved performance. By considering FR applied to the linear advection equation, and through modification to the Sobolev norm, criteria are presented for a wider set of correction functions. Legendre polynomials are then used to fulfil these criterion and realise functions for third to fifth order FR. The sufficient conditions for the existence of the modified norms are also explored, before Fourier and Von Neumann analysis are applied to analytically find temporal stability limits for various Runge-Kutta temporal integration schemes. For all cases, correction functions are found that extend the temporal stability of FR. Two application-inspired investigations are performed that aim to explore the effect of aliasing and non-linear equations. In both cases unique correction functions could be found that give good performance, compared to previous FR schemes, while also improving upon the temporal stability limit.


title: Sound Event Detection and Time-Frequency Segmentation from Weakly Labelled Data

abstract: Sound event detection (SED) aims to detect when and recognize what sound events happen in an audio clip. Many supervised SED algorithms rely on strongly labelled data which contains the onset and offset annotations of sound events. However, many audio tagging datasets are weakly labelled, that is, only the presence of the sound events is known, without knowing their onset and offset annotations. In this paper, we propose a time-frequency (T-F) segmentation framework trained on weakly labelled data to tackle the sound event detection and separation problem. In training, a segmentation mapping is applied on a T-F representation, such as log mel spectrogram of an audio clip to obtain T-F segmentation masks of sound events. The T-F segmentation masks can be used for separating the sound events from the background scenes in the time-frequency domain. Then a classification mapping is applied on the T-F segmentation masks to estimate the presence probabilities of the sound events. We model the segmentation mapping using a convolutional neural network and the classification mapping using a global weighted rank pooling (GWRP). In SED, predicted onset and offset times can be obtained from the T-F segmentation masks. As a byproduct, separated waveforms of sound events can be obtained from the T-F segmentation masks. We remixed the DCASE 2018 Task 1 acoustic scene data with the DCASE 2018 Task 2 sound events data. When mixing under 0 dB, the proposed method achieved F1 scores of 0.534, 0.398 and 0.167 in audio tagging, frame-wise SED and event-wise SED, outperforming the fully connected deep neural network baseline of 0.331, 0.237 and 0.120, respectively. In T-F segmentation, we achieved an F1 score of 0.218, where previous methods were not able to do T-F segmentation.


title: The pre-symplectic geometry of opers and the holonomy map

abstract: In this paper, we construct the moduli space of marked oper structures on a closed, oriented smooth surface of negative Euler characteristic as a holomorphic fiber bundle over Teichm\"{u}ller space. We prove that the holonomy map from the space of marked oper structures to the moduli space of reductive flat bundles is a holomorphic immersion, generalizing the known results for the moduli space of marked complex projective structures. Finally, we prove that the symplectic structure on the moduli space of marked complex projective structures extends to a pre-symplectic structure on the moduli space of marked opers whose reduced phase space is the space of marked complex projective structures.


title: Arithmetic models for Shimura varieties

abstract: We describe recent work on the construction of well-behaved arithmetic models for large classes of Shimura varieties and report on progress in the study of these models.


title: On inverse problem with phase retrieval for an inclined line in the parabolic approximation

abstract: The inverse problem of amplitude reconstruction on an inclined line based on the values of amplitude or its module as recorded on semi-infinite line orthogonal to the beam propagation direction is considered within the framework of 2D parabolic equation. It is demonstrated that this inverse problem, in case when the complex image plane amplitude is known, can be reduced to a singular Cauchy type integral equation. The existence of its solutions requires that certain conditions be met but if a solution exists it is necessary unique. The obtained integral equation is then approximated piece-wisely and the resulting linear algebraic system is solved numerically while applying necessary regularization procedures to enhance the stability of its solutions. Finally, an iterative method of phase retrieval is developed and a set of numerical experiments is performed.


title: Target detection in synthetic aperture radar imagery: a state-of-the-art survey

abstract: Target detection is the front-end stage in any automatic target recognition system for synthetic aperture radar (SAR) imagery (SAR-ATR). The efficacy of the detector directly impacts the succeeding stages in the SAR-ATR processing chain. There are numerous methods reported in the literature for implementing the detector. We offer an umbrella under which the various research activities in the field are broadly probed and taxonomized. First, a taxonomy for the various detection methods is proposed. Second, the underlying assumptions for different implementation strategies are overviewed. Third, a tabular comparison between careful selections of representative examples is introduced. Finally, a novel discussion is presented, wherein the issues covered include suitability of SAR data models, understanding the multiplicative SAR data models, and two unique perspectives on constant false alarm rate (CFAR) detection: signal processing and pattern recognition. From a signal processing perspective, CFAR is shown to be a finite impulse response band-pass filter. From a statistical pattern recognition perspective, CFAR is shown to be a suboptimal one-class classifier: a Euclidian distance classifier and a quadratic discriminant with a missing term for one-parameter and two-parameter CFAR, respectively. We make a contribution toward enabling an objective design and implementation for target detection in SAR imagery.


title: Fast Prefix Search in Little Space, with Applications

abstract: It has been shown in the indexing literature that there is an essential difference between prefix/range searches on the one hand, and predecessor/rank searches on the other hand, in that the former provably allows faster query resolution. Traditionally, prefix search is solved by data structures that are also dictionaries---they actually contain the strings in $S$. For very large collections stored in slow-access memory, we propose much more compact data structures that support \emph{weak} prefix searches---they return the ranks of matching strings provided that \emph{some} string in $S$ starts with the given prefix. In fact, we show that our most space-efficient data structure is asymptotically space-optimal. Previously, data structures such as String B-trees (and more complicated cache-oblivious string data structures) have implicitly supported weak prefix queries, but they all have query time that grows logarithmically with the size of the string collection. In contrast, our data structures are simple, naturally cache-efficient, and have query time that depends only on the length of the prefix, all the way down to constant query time for strings that fit in one machine word. We give several applications of weak prefix searches, including exact prefix counting and approximate counting of tuples matching conjunctive prefix conditions.


title: Econophysics Beyond General Equilibrium: the Business Cycle Model

abstract: Current business cycle theory is an application of the general equilibrium theory. This paper presents the business cycle model without using general equilibrium framework. We treat agents risk assessments as their coordinates x on economic space and establish distribution of all economic agents by their risk coordinates. We suggest aggregation of agents and their variables by scales large to compare with risk scales of single agents and small to compare with economic domain on economic space. Such model is alike to transition from kinetic description of multi-particle system to hydrodynamic approximation. Aggregates of agents extensive variables with risk coordinate x determine macro variables as functions of x alike to hydrodynamic variables. Economic and financial transactions between agents define evolution of their variables. Aggregation of transactions between agents with risk coordinates x and y determine macro transactions as functions of x and y and define evolution of macro variables at points x and y. We describe evolution and interactions between macro transactions by hydrodynamic-like system of economic equations. We show that business cycles are described as consequence of the system of economic equations on macro transactions. As example we describe Credit transactions CL(tax,y) that provide Loans from Creditors at point x to Borrowers at point y and Loan-Repayment transactions LR(t,x,y) that describe repayments from Borrowers at point y to Creditors at point x. We use hydrodynamic-like economic equations and derive from them the system of ordinary differential equations that describe business cycle fluctuations of macro Credits C(t) and macro Loan-Repayments LR(t) of the entire economics. The nature of business cycle fluctuations is explained as oscillations of "mean risk" of economic variables on bounded economic domain of economic space.


title: Testing the Weak Equivalence Principle using Optical and Near-Infrared Crab Pulses

abstract: The Weak Equivalence Principle states that the geodesics of a test particle in a gravitational field are independent of the particle's constitution. To constrain violations of the Weak Equivalence Principle, we use the one-meter telescope at Table Mountain Observatory near Los Angeles to monitor the relative arrival times of pulses from the Crab Pulsar in the optical ($\lambda \approx 585$ nm) and near-infrared ($\lambda \approx 814$ nm) using an instrument which detects single photons with nanosecond-timing resolution in those two bands. The infrared pulse arrives slightly before the visible pulse. Our three analysis methods give delays with statistical errors of $\Delta t_{obs} = 7.41 \pm 0.58$, $0.4 \pm 3.6$, and $7.35 \pm 4.48$ microseconds (at most 1/4000 of the pulsar period). We attribute this discrepancy to systematic error from the fact that the visible and infrared pulses have slightly different shapes. Whether this delay emerges from the pulsar, is caused by passing through wavelength-dependent media, or is caused by a violation of the equivalence principle, unless there is a fine-tuned cancellation among these, we set the first upper limit on the differential post-Newtonian parameter at these wavelengths of $\Delta \gamma < 1.07 \times 10^{-10}~(3\sigma)$. This result falls in an unexplored region of parameter space and complements existing limits on equivalence-principle violation from fast radio bursts, gamma ray bursts, as well as previous limits from the Crab.


title: Lower semicontinuity of the ADM mass in dimensions two through seven

abstract: The semicontinuity phenomenon of the ADM mass under pointed (i.e., local) convergence of asymptotically flat metrics is of interest because of its connections to nonnegative scalar curvature, the positive mass theorem, and Bartnik's mass-minimization problem in general relativity. In this paper, we extend a previously known semicontinuity result in dimension three for $C^2$ pointed convergence to higher dimensions, up through seven, using recent work of S. McCormick and P. Miao (which itself builds on the Riemannian Penrose inequality of H. Bray and D. Lee). For a technical reason, we restrict to the case in which the limit space is asymptotically Schwarzschild. In a separate result, we show that semicontinuity holds under weighted, rather than pointed, $C^2$ convergence, in all dimensions $n \geq 3$, with a simpler proof independent of the positive mass theorem. Finally, we also address the two-dimensional case for pointed convergence, in which the asymptotic cone angle assumes the role of the ADM mass.


title: Escher squares and lattice links

abstract: We give a shorter and simpler proof of the result of [2], which gives a necessary and sufficient condition for when a lattice diagram is the projection of a lattice link.


title: Network-based protein structural classification

abstract: Experimental determination of protein function is resource-consuming. As an alternative, computational prediction of protein function has received attention. In this context, protein structural classification (PSC) can help, by allowing for determining structural classes of currently unclassified proteins based on their features, and then relying on the fact that proteins with similar structures have similar functions. Existing PSC approaches rely on sequence-based or direct 3-dimensional (3D) structure-based protein features. In contrast, we first model 3D structures of proteins as protein structure networks (PSNs). Then, we use network-based features for PSC. We propose the use of graphlets, state-of-the-art features in many research areas of network science, in the task of PSC. Moreover, because graphlets can deal only with unweighted PSNs, and because accounting for edge weights when constructing PSNs could improve PSC accuracy, we also propose a deep learning framework that automatically learns network features from weighted PSNs. When evaluated on a large set of ~9,400 CATH and ~12,800 SCOP protein domains (spanning 36 PSN sets), our proposed approaches are superior to existing PSC approaches in terms of accuracy, with comparable running time.


title: Compression stiffening in biological tissues: on the possibility of classic elasticity origins

abstract: Compression stiffening, or an increase in shear modulus with increasing compressive strain, has been observed in recent rheometry experiments on brain, liver, and fat tissues. Here, we extend the known types of biomaterials exhibiting this phenomenon to include agarose gel and fruit flesh. Further, we show that two different results from classic elasticity theory can account for the phenomenon of linear compression stiffening. One result is due to Barron and Klein, extended here to the relevant geometry and pre-stresses; the other is due to Birch. For incompressible materials, there are no adjustable parameters in either theory. Which one applies to a given situation is a matter of reference state, suggesting that the reference state is determined by the tendency of the material to develop, or not develop, axial stress (in excess of the applied pre-stress) when subjected to torsion at constant axial strain. Our experiments and analysis also strengthen the notion that seemingly distinct animal and plant tissues can have mechanically similar behavior under certain conditions.


title: Spaceflight from Super-Earths is difficult

abstract: Many rocky exoplanets are heavier and larger than the Earth, and have higher surface gravity. This makes space-flight on these worlds very challenging, because the required fuel mass for a given payload is an exponential function of planetary surface gravity. We find that chemical rockets still allow for escape velocities on Super-Earths up to 10x Earth mass. More massive rocky worlds, if they exist, would require other means to leave the planet, such as nuclear propulsion.


title: Two-mode squeezing operator in circuit QED

abstract: We theoretically investigate the implementation of the two-mode squeezing operator in circuit quantum electrodynamics. Inspired by a previous scheme for optical cavities [Phys. Rev. A $\textbf{73}$, 043803(2006)], we employ a superconducting qubit coupled to two nondegenerate quantum modes and use a driving field on the qubit to adequately control the resonator-qubit interaction. Based on the generation of two-mode squeezed vacuum states, firstly we analyze the validity of our model in the ideal situation and then we investigate the influence of the dissipation mechanisms on the generation of the two-mode squeezing operation, namely the qubit and resonator mode decays and qubit dephasing. We show that our scheme allows the generation of highly squeezed states even with the state-of-the-art parameters, leading to a theoretical prediction of more than 10 dB of two-mode squeezing. Furthermore, our protocol is able to squeeze an arbitrary initial state of the resonators, which makes our scheme attractive for future applications in continuous-variable quantum information processing and quantum metrology in the realm of circuit quantum electrodynamics.


title: Jet Lag Recovery: Synchronization of Circadian Oscillators as a Mean Field Game

abstract: The Suprachiasmatic Nucleus (SCN) is a region in the brain that is responsible for controlling circadian rhythms. The SCN contains on the order of 10^4 neuronal oscillators which have a preferred period slightly longer than 24 hours. The oscillators try to synchronize with each other as well as responding to external stimuli such as sunlight exposure. A mean field game model for these neuronal oscillators is formulated with two goals in mind: 1) to understand the long time behavior of the oscillators when an individual remains in the same time zone, and 2) to understand how the oscillators recover from jet lag when the individual has traveled across time zones. In particular, we would like to study the claim that jet lag is worse after traveling east than west. Finite difference schemes are used to find numerical approximations to the mean field game solutions. Numerical results are presented and conjectures are posed. The numerics suggest the time to recover from jet lag is about the same for east versus west trips, but the cost the oscillators accrue while recovering is larger for eastward trips.


title: Secure certification of mixed quantum states with application to two-party randomness generation

abstract: We investigate sampling procedures that certify that an arbitrary quantum state on $n$ subsystems is close to an ideal mixed state $\varphi^{\otimes n}$ for a given reference state $\varphi$, up to errors on a few positions. This task makes no sense classically: it would correspond to certifying that a given bitstring was generated according to some desired probability distribution. However, in the quantum case, this is possible if one has access to a prover who can supply a purification of the mixed state. In this work, we introduce the concept of mixed-state certification, and we show that a natural sampling protocol offers secure certification in the presence of a possibly dishonest prover: if the verifier accepts then he can be almost certain that the state in question has been correctly prepared, up to a small number of errors. We then apply this result to two-party quantum coin-tossing. Given that strong coin tossing is impossible, it is natural to ask "how close can we get". This question has been well studied and is nowadays well understood from the perspective of the bias of individual coin tosses. We approach and answer this question from a different---and somewhat orthogonal---perspective, where we do not look at individual coin tosses but at the global entropy instead. We show how two distrusting parties can produce a common high-entropy source, where the entropy is an arbitrarily small fraction below the maximum (except with negligible probability).


title: MPSM: Multi-prospective PaaS Security Model

abstract: Cloud computing has brought a revolution in the field of information technology and improving the efficiency of computational resources. It offers computing as a service enabling huge cost and resource efficiency. Despite its advantages, certain security issues still hinder organizations and enterprises from it being adopted. This study mainly focused on the security of Platform-as-a-Service (PaaS) as well as the most critical security issues that were documented regarding PaaS infrastructure. The prime outcome of this study was a security model proposed to mitigate security vulnerabilities of PaaS. This security model consists of a number of tools, techniques and guidelines to mitigate and neutralize security issues of PaaS. The security vulnerabilities along with mitigation strategies were discussed to offer a deep insight into PaaS security for both vendor and client that may facilitate future design to implement secure PaaS platforms.


title: Multimodal Unsupervised Image-to-Image Translation

abstract: Unsupervised image-to-image translation is an important and challenging problem in computer vision. Given an image in the source domain, the goal is to learn the conditional distribution of corresponding images in the target domain, without seeing any pairs of corresponding images. While this conditional distribution is inherently multimodal, existing approaches make an overly simplified assumption, modeling it as a deterministic one-to-one mapping. As a result, they fail to generate diverse outputs from a given source domain image. To address this limitation, we propose a Multimodal Unsupervised Image-to-image Translation (MUNIT) framework. We assume that the image representation can be decomposed into a content code that is domain-invariant, and a style code that captures domain-specific properties. To translate an image to another domain, we recombine its content code with a random style code sampled from the style space of the target domain. We analyze the proposed framework and establish several theoretical results. Extensive experiments with comparisons to the state-of-the-art approaches further demonstrates the advantage of the proposed framework. Moreover, our framework allows users to control the style of translation outputs by providing an example style image. Code and pretrained models are available at https://github.com/nvlabs/MUNIT


title: Lattices With Internal Resonator Defects

abstract: We consider a variety of settings involving chains with one or more defects stemming from the introduction of nodes bearing internal resonators. Motivated by experimental results in woodpile elastic lattices with one or two defects, we consider a variety of different theoretical scenarios. These include multi-defect chains and their ability to transmit, reflect, and especially trap energy; they also include settings with linear vs. nonlinear defects of variable interaction exponent. Moreover, they involve defects which are spatially separated and either statically, or more effectively dynamically, enable the confinement of energy between the separated defects. Wherever possible, comparisons of the experiments with numerical simulations, as well as with theoretical intuition are also offered, to provide a justification for the observed findings.


title: Large deviations for fast transport stochastic RDEs with applications to the exit problem

abstract: We study reaction diffusion equations with a deterministic reaction term as well as two random reaction terms, one that acts on the interior of the domain, and another that acts only on the boundary of the domain. We are interested in the regime where the relative sizes of the diffusion and reaction terms are different. Specifically, we consider the case where the diffusion rate is much larger than the rate of reaction, and the deterministic rate of reaction is much larger than either of the random rate of reactions.


title: Holographic pump probe spectroscopy

abstract: We study the non-linear response of a 2+1 dimensional holographic model with weak momentum relaxation and finite charge density to an oscillatory electric field pump pulse. Following the time evolution of one point functions after the pumping has ended, we find that deviations from thermality are well captured within the linear response theory. For electric pulses with a negligible zero frequency component the response approaches the instantaneously thermalizing form typical of holographic Vaidya models. We link this to the suppression of the amplitude of the quasinormal mode that governs the approach to equilibrium. In the large frequency limit, we are also able to show analytically that the holographic geometry takes the Vaidya form. A simple toy model captures these features of our holographic setup. Computing the out-of-equilibrium probe optical conductivity after the pump pulse, we similarly find that for high-frequency pulses the optical conductivity reaches its final equilibrium value effectively instantaneously. Pulses with significant DC components show exponential relaxation governed by twice the frequency of the vector quasinormal mode that governs the approach to equilibrium for the background solution. We explain this numerical factor in terms of a simple symmetry argument.

categories: hep-th cond-mat.str-el cond-mat.supr-con quant-ph


title: Adaptive Ensemble Biomolecular Simulations at Scale

abstract: Recent advances in both theory and methods have created opportunities to simulate biomolecular processes more efficiently using adaptive ensemble simulations. Ensemble-based simulations are used widely to compute a number of individual simulation trajectories and analyze statistics across them. Adaptive ensemble simulations offer a further level of sophistication and flexibility by enabling high-level algorithms to control simulations based on intermediate results. Novel high-level algorithms require sophisticated approaches to utilize the intermediate data during runtime. Thus, there is a need for scalable software systems to support adaptive ensemble-based applications. We describe the operations in executing adaptive workflows, classify different types of adaptations, and describe challenges in implementing them in software tools. We enhance Ensemble Toolkit (EnTK) -- an ensemble execution system -- to support the scalable execution of adaptive workflows on HPC systems, and characterize the adaptation overhead in EnTK. We implement two high-level adaptive ensemble algorithms -- expanded ensemble and Markov state modeling, and execute upto $2^{12}$ ensemble members, on thousands of cores on three distinct HPC platforms. We highlight scientific advantages enabled by the novel capabilities of our approach. To the best of our knowledge, this is the first attempt at describing and implementing multiple adaptive ensemble workflows using a common conceptual and implementation framework.


title: A Scalable Shared-Memory Parallel Simplex for Large-Scale Linear Programming

abstract: The Simplex tableau has been broadly used and investigated in the industry and academia. With the advent of the big data era, ever larger problems are posed to be solved in ever larger machines whose architecture type did not exist in the conception of this algorithm. In this paper, we present a shared-memory parallel implementation of the Simplex tableau algorithm for dense large-scale Linear Programming (LP) problems for use in modern multi-core architectures. We present the general scheme and explain the strategies taken to parallelize each step of the standard simplex algorithm, emphasizing the solutions found to solve performance bottlenecks. We analyzed the speedup and the parallel efficiency for the proposed implementation relative to the standard Simplex algorithm using a shared-memory system with 64 processing cores. The experiments were performed for several different problems, with up to 8192 variables and constraints, in their primal and dual formulations. The results show that the performance is mostly much better when we use the formulation with more variables than inequality constraints. Also, they show that the parallelization strategies applied to avoid bottlenecks lead the implementation to scale well with the problem size and the core count up to a certain limit of problem size. Further analysis showed that this scaling limit was an effect of resource limitation. Even though, our implementation was able to reach speedups in the order of 19x.


title: An iterative reconstruction of cosmological initial density fields

abstract: We present an iterative method to reconstruct the linear-theory initial conditions from the late-time cosmological matter density field, with the intent of improving the recovery of the cosmic distance scale from the baryon acoustic oscillations (BAOs). We present tests using the dark matter density field in both real and redshift space generated from an $N$-body simulation. In redshift space at $z = 0.5$, we find that the reconstructed displacement field using our iterative method are more than 80\% correlated with the true displacement field of the dark matter particles on scales $k < 0.10h\ {\rm Mpc}^{-1}$. Furthermore, we show that the two-point correlation function of our reconstructed density field matches that of the initial density field substantially better, especially on small scales ($< 40h^{-1}\ {\rm Mpc}$). Our redshift-space results are improved if we use an anisotropic smoothing so as to account for the reduced small-scale information along the line of sight in redshift space.


title: Efficient algorithms for tensor scaling, quantum marginals and moment polytopes

abstract: We present a polynomial time algorithm to approximately scale tensors of any format to arbitrary prescribed marginals (whenever possible). This unifies and generalizes a sequence of past works on matrix, operator and tensor scaling. Our algorithm provides an efficient weak membership oracle for the associated moment polytopes, an important family of implicitly-defined convex polytopes with exponentially many facets and a wide range of applications. These include the entanglement polytopes from quantum information theory (in particular, we obtain an efficient solution to the notorious one-body quantum marginal problem) and the Kronecker polytopes from representation theory (which capture the asymptotic support of Kronecker coefficients). Our algorithm can be applied to succinct descriptions of the input tensor whenever the marginals can be efficiently computed, as in the important case of matrix product states or tensor-train decompositions, widely used in computational physics and numerical mathematics. We strengthen and generalize the alternating minimization approach of previous papers by introducing the theory of highest weight vectors from representation theory into the numerical optimization framework. We show that highest weight vectors are natural potential functions for scaling algorithms and prove new bounds on their evaluations to obtain polynomial-time convergence. Our techniques are general and we believe that they will be instrumental to obtain efficient algorithms for moment polytopes beyond the ones consider here, and more broadly, for other optimization problems possessing natural symmetries.


title: Persistent Homology and Euler Integral Transforms

abstract: The Euler calculus -- an integral calculus based on Euler characteristic as a valuation on constructible functions -- is shown to be an incisive tool for answering questions about injectivity and invertibility of recent transforms based on persistent homology for shape characterization.