title: Dynamics of Equity, Efficiency, and Efficacy in Home Health Care with Patient and Caregiver Preferences

abstract: There are three main entities in an HHC setting: the firm, the caregivers, and the patients. Often times the interests of each group of entities conflict with the other two. Whether a non-profit or a for-profit HHC, the firm generally seeks efficiency in terms of maximum utilization of the resources and highest level of care provided to the patients. Caregivers, seek a fair utilization across the network in order to increase job satisfaction. Finally, the patients seek maximum satisfaction level and a fair share of the care provided by the firm. Therefore, looking at HHCs from either one of these lenses by themselves leads to decisions that may neglect the interests of the others. The goal of this paper is to integrate the objectives of these three groups into one place and study their dynamics.


title: Structure from Silence: Learning Scene Structure from Ambient Sound

abstract: From whirling ceiling fans to ticking clocks, the sounds that we hear subtly vary as we move through a scene. We ask whether these ambient sounds convey information about 3D scene structure and, if so, whether they provide a useful learning signal for multimodal models. To study this, we collect a dataset of paired audio and RGB-D recordings from a variety of quiet indoor scenes. We then train models that estimate the distance to nearby walls, given only audio as input. We also use these recordings to learn multimodal representations through self-supervision, by training a network to associate images with their corresponding sounds. These results suggest that ambient sound conveys a surprising amount of information about scene structure, and that it is a useful signal for learning multimodal features.


title: Impact of LIGO-Virgo black hole binaries on gravitational wave background searches

abstract: We study the impact of the black hole binary population currently probed by LIGO-Virgo on future searches for the primordial gravitational wave background. We estimate the foreground generated by the binaries using the observed event rate and a simple modeling of the black hole population. We subtract individually resolvable binaries from the foreground and utilize Fisher analysis to derive sensitivity curves for power-law signals including these astrophysical foregrounds. Even with optimistic assumptions, we find that the reach of future experiments will be severely reduced.


title: Chirality production during axion inflation

abstract: We study the generation of the chiral charge during axion inflation where the pseudoscalar inflaton field $\phi$ couples axially to the electromagnetic field through the term $(\beta/M_p)\phi\,\boldsymbol{E}\cdot\boldsymbol{B}$ with dimensionless coupling constant $\beta$. To describe the evolution of electromagnetic field and determine $\langle\boldsymbol{E}\cdot\boldsymbol{B}\rangle$ sourcing the chiral asymmetry during inflation due to the chiral anomaly, we employ the gradient expansion formalism. It operates with a set of vacuum expectation values of bilinear electromagnetic functions and allows us to take into account the backreaction of generated fields on the inflaton evolution as well as the Schwinger production of charged fermions. In addition, we include the chiral magnetic effect contribution to the electric current $\boldsymbol{j}_{\rm CME}=e^{2}/(2\pi^2)\mu_{5}\boldsymbol{B}$, where $\mu_5$ is the chiral chemical potential which quantifies the chiral charge production. Solving a set of equations for the inflaton field, scale factor, quadratic functions of the electromagnetic field, and the chiral charge density (chiral chemical potential), we find that the chirality production is quite efficient leading to the generation of a large chemical potential at the end of axion inflation.


title: Advances in Neural Rendering

abstract: Synthesizing photo-realistic images and videos is at the heart of computer graphics and has been the focus of decades of research. Traditionally, synthetic images of a scene are generated using rendering algorithms such as rasterization or ray tracing, which take specifically defined representations of geometry and material properties as input. Collectively, these inputs define the actual scene and what is rendered, and are referred to as the scene representation (where a scene consists of one or more objects). Example scene representations are triangle meshes with accompanied textures (e.g., created by an artist), point clouds (e.g., from a depth sensor), volumetric grids (e.g., from a CT scan), or implicit surface functions (e.g., truncated signed distance fields). The reconstruction of such a scene representation from observations using differentiable rendering losses is known as inverse graphics or inverse rendering. Neural rendering is closely related, and combines ideas from classical computer graphics and machine learning to create algorithms for synthesizing images from real-world observations. Neural rendering is a leap forward towards the goal of synthesizing photo-realistic image and video content. In recent years, we have seen immense progress in this field through hundreds of publications that show different ways to inject learnable components into the rendering pipeline. This state-of-the-art report on advances in neural rendering focuses on methods that combine classical rendering principles with learned 3D scene representations, often now referred to as neural scene representations. A key advantage of these methods is that they are 3D-consistent by design, enabling applications such as novel viewpoint synthesis of a captured scene. In addition to methods that handle static scenes, we cover neural scene representations for modeling non-rigidly deforming objects...


title: Towards Green Automated Machine Learning: Status Quo and Future Directions

abstract: Automated machine learning (AutoML) strives for the automatic configuration of machine learning algorithms and their composition into an overall (software) solution - a machine learning pipeline - tailored to the learning task (dataset) at hand. Over the last decade, AutoML has developed into an independent research field with hundreds of contributions. At the same time, AutoML is being criticised for its high resource consumption as many approaches rely on the (costly) evaluation of many machine learning pipelines, as well as the expensive large scale experiments across many datasets and approaches. In the spirit of recent work on Green AI, this paper proposes Green AutoML, a paradigm to make the whole AutoML process more environmentally friendly. Therefore, we first elaborate on how to quantify the environmental footprint of an AutoML tool. Afterward, different strategies on how to design and benchmark an AutoML tool wrt. their "greenness", i.e. sustainability, are summarized. Finally, we elaborate on how to be transparent about the environmental footprint and what kind of research incentives could direct the community into a more sustainable AutoML research direction. Additionally, we propose a sustainability checklist to be attached to every AutoML paper featuring all core aspects of Green AutoML.


title: Light propagation and magnon-photon coupling in optically dispersive magnetic media

abstract: Achieving strong coupling between light and matter excitations in hybrid systems is a benchmark for the implementation of quantum technologies. We recently proposed [arXiv:2110.02984] that strong single-particle coupling between magnons and light can be realized in a magnetized epsilon-near-zero (ENZ) medium, in which magneto-optical effects are enhanced. Here we present a detailed derivation of the magnon-photon coupling Hamiltonian in dispersive media both for degenerate and non-degenerate optical modes, and show the enhancement of the coupling near the ENZ frequency. Moreover, we show that the coupling of magnons to plane-wave non-degenerate Voigt modes vanishes at specific frequencies due to polarization selection rules tuned by dispersion. Finally, we present specific results using a Lorentz dispersion model. Our results pave the way for the design of dispersive optomagnonic systems, providing a general theoretical framework for describing engineering ENZ-based optomagnonic systems.


title: Light scalars in neutron star mergers

abstract: Due to their unique set of multimessenger signals, neutron star mergers have emerged as a novel environment for studies of new physics beyond the Standard Model (SM). As a case study, we consider the simplest extension of the SM scalar sector involving a light CP-even scalar singlet $S$ mixing with the SM Higgs boson. These $S$ particles can be produced abundantly in neutron star mergers via the nucleon bremsstrahlung process. We show that the $S$ particles may either be trapped in or stream freely out of the merger remnant, depending on the $S$ mass, its mixing with the SM Higgs boson, and the temperature and baryon density in the merger. In the free-streaming region, the scalar $S$ will provide an extra channel to cool down the merger remnant, with cooling timescales as small as ${\cal O}$(ms). On the other hand, in the trapped region, the Bose gas of $S$ particles could contribute a larger thermal conductivity than the trapped neutrinos in some parts of the parameter space, thus leading to faster thermal equilibration than expected. Therefore, future observations of the early postmerger phase of a neutron star merger could effectively probe a unique range of the $S$ parameter space, largely complementary to the existing and future laboratory and supernova limits. In view of these results, we hope the merger simulation community will be motivated to implement the effects of light CP-even scalars into their simulations in both the free-streaming and trapped regimes.


title: Complexity for Dynamical Anisotropic Sphere in f(G,T) Gravity

abstract: This paper is devoted to the formulation of a complexity factor for dynamical anisotropic sphere in the framework of $f(G,T)$ gravity, where $G$ is the Gauss-Bonnet invariant and $T$ is the trace of energy-momentum tensor. Inhomogeneous energy density, anisotropic pressure, heat dissipation and modified terms create complexity within the self-gravitating system. We evaluate the structure scalars by orthogonal splitting of the Riemann tensor to evaluate a complexity factor which incorporates all the fundamental properties of the system. Moreover, we examine the dynamics of the sphere by assuming homologous mode as the simplest pattern of evolution. We also discuss dissipative as well as non-dissipative scenarios corresponding to homologous and complexity free conditions. Finally, we establish a criterion under which the complexity free condition remains stable throughout the process of evolution. We conclude that the presence of dark source terms of $f(G,T)$ gravity increase the system's complexity.


title: Deep ReLU neural network approximation in Bochner spaces and applications to parametric PDEs

abstract: We investigate non-adaptive methods of deep ReLU neural network approximation in Bochner spaces $L_2({\mathbb U}^\infty, X, \mu)$ of functions on ${\mathbb U}^\infty$ taking values in a separable Hilbert space $X$, where ${\mathbb U}^\infty$ is either ${\mathbb R}^\infty$ equipped with the standard Gaussian probability measure, or ${\mathbb I}^\infty:= [-1,1]^\infty$ equipped with the Jacobi probability measure. Functions to be approximated are assumed to satisfy a certain weighted $\ell_2$-summability of the generalized chaos polynomial expansion coefficients with respect to the measure $\mu$. We prove the convergence rate of this approximation in terms of the size of approximating deep ReLU neural networks. These results then are applied to approximation of the solution to parametric elliptic PDEs with random inputs for the lognormal and affine cases.


title: A new radiomics feature: image frequency analysis

abstract: Radiomics is a promising technology that focuses on improvements of image analysis, using an automated high-throughput extraction of quantitative features. However, the character of lesion is affected by the surrounding tissue. A lesion on medical image should be characterized from the inter-relation between lesion and surrounding tissue as well as property of the lesion itself. The aim of this study is to introduce a new radiomics feature which quantitatively analyze the inter-relation between lesion and surrounding tissue focusing on the value change of rows and columns in a medical image.


title: Fokker-Planck equation for black holes in thermal potential

abstract: We construct a kind of thermal potential and then put the black hole thermodynamic system in it. In this regard, some thermodynamic properties of the black hole are related to the geometric characteristics of the thermal potential. Driven by the intrinsic thermodynamic fluctuations, the behavior of the black hole in the thermal potential is stochastic. With the help of solving the Fokker-Planck equation analytically, we obtain the discrete energy spectrum of Schwarzschild and Banados-Teitelboim-Zanelli (BTZ) black holes in the thermal potential. For Schwarzschild black hole, the energy spectrum is proportional to the temperature of the ensemble, which is an external parameter, and the ground state is non-zero. For BTZ black hole, the energy spectrum only depends on the AdS radius, which is the intrinsic parameter. Moreover, the ground state of BTZ black hole in thermal potential is zero. This also reflects the difference between three-dimensional gravity and four-dimensional gravity.


title: New Chandrasekhar transformation in Kerr spacetime

abstract: We construct a new type of Chandrasekhar transformation in Kerr spacetime using the different tortoise coordinate, which is useful for exact analysis to study Teukolsky equation with arbitrary frequency. We also give the interpretation of our transformation using the formalism of the quantum Seiberg-Witten geometry.


title: Fractional Quantum Hall edge dynamics from a Quantum Optics perspective

abstract: The goal of this master thesis is to give a different perspective on the problem of quasi-particle tunneling in Fractional Quantum Hall liquids by using a typical Quantum Optics tool such as the Truncated Wigner Approximation. Our novel approach provides a mapping of the quasi-particle tunneling problem into a stochastic classical field problem. This offers an intuitive physical picture which allows us to recover non-trivial features of FQH physics and predict novel phenomena.


title: PDMP Monte Carlo methods for piecewise-smooth densities

abstract: There has been substantial interest in developing Markov chain Monte Carlo algorithms based on piecewise-deterministic Markov processes. However existing algorithms can only be used if the target distribution of interest is differentiable everywhere. The key to adapting these algorithms so that they can sample from to densities with discontinuities is defining appropriate dynamics for the process when it hits a discontinuity. We present a simple condition for the transition of the process at a discontinuity which can be used to extend any existing sampler for smooth densities, and give specific choices for this transition which work with popular algorithms such as the Bouncy Particle Sampler, the Coordinate Sampler and the Zig-Zag Process. Our theoretical results extend and make rigorous arguments that have been presented previously, for instance constructing samplers for continuous densities restricted to a bounded domain, and we present a version of the Zig-Zag Process that can work in such a scenario. Our novel approach to deriving the invariant distribution of a piecewise-deterministic Markov process with boundaries may be of independent interest.


title: X-Ray Evidence Against the Hypothesis that the Hyper-Luminous z=6.3 Quasar J0100+2802 is Lensed

abstract: The $z=6.327$ quasar SDSS J010013.02+280225.8 (hereafter J0100+2802) is believed to be powered by a black hole more massive than $10^{10}\ {\rm M}_\odot$, making it the most massive black hole known in the first billion years of the Universe. However, recent high-resolution ALMA imaging shows four structures at the location of this quasar, potentially implying that it is lensed with a magnification of $\mu\sim450$ and thus its black hole is significantly less massive. Furthermore, for the underlying distribution of magnifications of $z\gtrsim6$ quasars to produce such an extreme value, theoretical models predict that a larger number of quasars in this epoch should be lensed, implying further overestimates of early black hole masses. To provide an independent constraint on the possibility that J0100+2802 is lensed, we re-analyzed archival XMM-Newton observations of the quasar and compared the expected ratios of X-ray luminosity to rest-frame UV and IR luminosities. For both cases, J0100+2802's X-ray flux is consistent with the no-lensing scenario; while this could be explained by J0100+2802 being X-ray faint, we find it does not have the X-ray or optical spectral features expected for an X-ray faint quasar. Finally, we compare the overall distribution of X-ray fluxes for known, typical $z\gtrsim6$ quasars. We find a $3\sigma$ tension between the observed and predicted X-ray-to-UV flux ratios when adopting the magnification probability distribution required to produce a $\mu=450$ quasar.


title: NRPyLaTeX: A LaTeX interface to computer algebra systems for general relativity

abstract: While each computer algebra system (CAS) contains its own unique syntax for inputting mathematical expressions, LaTeX is perhaps the most widespread language for typesetting mathematics. NRPyLaTeX (NL) enables direct LaTeX input of complex tensorial expressions (written in Einstein notation) relevant to general relativity and differential geometry into the SymPy CAS. As SymPy also supports output compatible with the Mathematica and Maple CASs, NL lowers the learning curve for inputting and manipulating tensorial expressions in three widely used CASs. LaTeX however is a typesetting language, and as such is not designed to resolve ambiguities in mathematical expressions. To address this, NL implements a convenient configuration interface that, e.g., defines variables with certain attributes. Configuration commands appear as LaTeX comments, so that entire NL workflows can fit seamlessly into the LaTeX source code of scientific papers without interfering with the rendered mathematical expressions. Further, NL adopts NRPy+'s rigid syntax for indexed symbols (e.g., tensors), which enables NL output to be directly converted into highly optimized C/C++-code kernels using NRPy+. Finally NL has robust and user-friendly error-handling, which catches common tensor indexing errors and reports unresolved ambiguities, further expediting the input and validation of LaTeX expressions into a CAS.


title: TRIDENT: A Rapid 3D Radiative Transfer Model for Exoplanet Transmission Spectra

abstract: Transmission spectroscopy is one of the premier methods used to probe the temperature, composition, and cloud properties of exoplanet atmospheres. Recent studies have demonstrated that the multidimensional nature of exoplanet atmospheres -- due to non-uniformities across the day-night transition and between the morning and evening terminators -- can strongly influence transmission spectra. However, the computational demands of 3D radiative transfer techniques have precluded their usage within atmospheric retrievals. Here we introduce TRIDENT, a new 3D radiative transfer model which rapidly computes transmission spectra of exoplanet atmospheres with day-night, morning-evening, and vertical variations in temperature, chemical abundances, and cloud properties. We also derive a general equation for transmission spectra, accounting for 3D atmospheres, refraction, multiple scattering, ingress/egress, grazing transits, stellar heterogeneities, and nightside thermal emission. After introducing TRIDENT's linear algebra-based approach to 3D radiative transfer, we propose new parametric prescriptions for 3D temperature and abundance profiles and 3D clouds. We show that multidimensional transmission spectra exhibit two significant observational signatures: (i) day-night composition gradients alter the relative amplitudes of absorption features; and (ii) morning-evening composition gradients distort the peak-to-wing contrast of absorption features. Finally, we demonstrate that these signatures of multidimensional atmospheres incur residuals > 100 ppm compared to 1D models, rendering them potentially detectable with JWST. TRIDENT's rapid radiative transfer, coupled with parametric multidimensional atmospheres, unlocks the final barrier to 3D atmospheric retrievals.


title: On the synthesis of heavy nuclei in protomagnetar outflows and implications for ultra-high energy cosmic rays

abstract: It has been suggested that strongly magnetised and rapidly rotating protoneutron stars (PNSs) may produce long duration gamma-ray bursts (GRBs) originating from stellar core collapse. We explore the steady-state properties and heavy element nucleosynthesis in neutrino-driven winds from such PNSs whose magnetic axis is generally misaligned with the axis of rotation. We consider a wide variety of central engine properties such as surface dipole field strength, initial rotation period and magnetic obliquity to show that heavy element nuclei can be synthesised in the radially expanding wind. This process is facilitated provided the outflow is Poynting-flux dominated such that its low entropy and fast expansion timescale enables heavy nuclei to form in a more efficient manner as compared to the equivalent thermal GRB outflows. We also examine the acceleration and survival of these heavy nuclei and show that they can reach sufficiently high energies $\gtrsim 10^{20}\ {\rm eV}$ within the same physical regions that are also responsible for powering gamma-ray emission, primarily through magnetic dissipation processes. Although these magnetised outflows generally fail to achieve the production of elements heavier than lanthanides for our explored electron fraction range 0.4-0.6, we show that they are more than capable of synthesizing nuclei near and beyond iron peak elements.


title: Utilizing cosmic-ray positron and electron observations to probe the averaged properties of Milky Way pulsars

abstract: Pulsars have long been studied in the electromagnetic spectrum. Their environments are rich in high-energy cosmic-ray electrons and positrons likely enriching the interstellar medium with such particles. In this work we use recent cosmic-ray observations from the AMS-02, CALET and DAMPE collaborations to study the averaged properties of the local Milky Way pulsar population. We perform simulations of the local Milky Way pulsar population, for interstellar medium assumptions in agreement with a range of cosmic-ray nuclei measurements. Each such simulation contains $\sim 10^{4}$ pulsars of unique age, location, initial spin-down power and cosmic-ray electron/positron spectra. We produce more than $7\times 10^{3}$ such Milky Way pulsar simulations. We account for and study i) the pulsars' birth rates and the stochastic nature of their birth, ii) their initial spin-down power distribution, iii) their time evolution in terms of their braking index and characteristic spin-down timescale, iv) the fraction of spin-down power going to cosmic-ray electrons and positrons and v) their propagation through the interstellar medium and the Heliosphere. We find that pulsars of ages $\sim 10^{5}-10^{7}$ yr, have a braking index that on average has to be 3 or larger. Given that electromagnetic spectrum observations of young pulsars find braking indices lower than 3, our work provides strong hints that pulsars' braking index increases on average as they age, allowing them to retain some of their rotational energy. Moreover, we find that pulsars have relatively uniform properties as sources of cosmic-ray electrons and positrons in terms of the spectra they produce and likely release O($10\%$) of their rotational energy to cosmic-rays in the ISM. Finally, we find at $\simeq$12 GeV positrons a spectral feature that suggests a new subpopulation of positron sources contributing at these energies.


title: MATBG as Topological Heavy Fermion: I. Exact Mapping and Correlated Insulators

abstract: Magic-angle ($\theta=1.05^\circ$) twisted bilayer graphene (MATBG) has shown two seemingly contradictory characters: the localization and quantum-dot-like behavior in STM experiments, and delocalization in transport experiments. We construct a model, which naturally captures the two aspects, from the Bistritzer-MacDonald (BM) model in a first principle spirit. A set of local flat-band orbitals ($f$) centered at the AA-stacking regions are responsible to the localization. A set of extended topological conduction bands ($c$), which are at small energetic separation from the local orbitals, are responsible to the delocalization and transport. The topological flat bands of the BM model appear as a result of the hybridization of $f$- and $c$-electrons. This model then provides a new perspective for the strong correlation physics, which is now described as strongly correlated $f$-electrons coupled to nearly free topological semimetallic $c$-electrons - we hence name our model as the topological heavy fermion model. Using this model, we obtain the U(4) and U(4)$\times$U(4) symmetries as well as the correlated insulator phases and their energies. Simple rules for the ground states and their Chern numbers are derived. Moreover, features such as the large dispersion of the charge $\pm1$ excitations and the minima of the charge gap at the $\Gamma_M$ point can now, for the first time, be understood both qualitatively and quantitatively in a simple physical picture. Our mapping opens the prospect of using heavy-fermion physics machinery to the superconducting physics of MATBG.


title: Structure and kinematics of tidally limited satellite galaxies in LCDM

abstract: We use N-body simulations to model the tidal evolution of dark matter-dominated dwarf spheroidal galaxies embedded in cuspy Navarro-Frenk-White subhalos. Tides gradually peel off stars and dark matter from a subhalo, trimming it down according to their initial binding energy. This process strips preferentially particles with long orbital times, and comes to an end when the remaining bound particles have crossing times shorter than a fraction of the orbital time at pericentre. The properties of the final stellar remnant thus depend on the energy distribution of stars in the progenitor subhalo, which in turn depends on the initial density profile and radial segregation of the initial stellar component. The stellar component may actually be completely dispersed if its energy distribution does not extend all the way to the subhalo potential minimum, although a bound dark remnant may remain. These results imply that 'tidally-limited' galaxies, defined as systems whose stellar components have undergone substantial tidal mass loss, neither converge to a unique structure nor follow a single tidal track, as claimed in earlier work. On the other hand, tidally limited dwarfs do have characteristic sizes and velocity dispersions that trace directly the characteristic radius ($r_{max}$) and circular velocity ($V_{max}$) of the subhalo remnant. This result places strong upper limits on the size of satellites whose unusually low velocity dispersions are often ascribed to tidal effects. In particular, the large size of kinematically-cold 'feeble giant' satellites like Crater 2 or Antlia 2 cannot be explained as due to tidal effects alone in the Lambda Cold Dark Matter scenario.


title: Pulsar-timing measurement of the circular polarization of the stochastic gravitational-wave background

abstract: Pulsar-timing arrays (PTAs) are in the near future expected to detect a stochastic gravitational-wave background (SGWB) produced by a population of inspiralling super-massive black hole binaries. In this work, we consider a background that can be anisotropic and circularly polarized. We use the expansion of the intensity and the circular polarization in terms of spherical harmonics and the overlap reduction functions for each term in this expansion. We propose an unbiased real-space estimator that can separate the intensity and circular-polarization contributions of the SGWB to pulsar-timing-residual correlations and then validate the estimator on simulated data. We compute the signal-to-noise ratio of a circular-polarization component that has a dipole pattern under different assumptions about the PTA. We find that a nearly-maximal circular-polarization dipole may be detectable, which can aid in determining whether or not the background is dominated by a handful of bright sources.


title: Probing Neutrino-Portal Dark Matter at the Forward Physics Facility

abstract: The Forward Physics Facility (FPF), planned to operate near the ATLAS interaction point at the LHC, offers exciting new terrain to explore neutrino properties at TeV energy scales. It will reach an unprecedented regime for terrestrial neutrino experiments and provide the opportunity to reveal new physics of neutrinos at higher energy scales. We demonstrate that future detectors at the FPF have the potential to discover new mediators that couple predominantly to neutrinos, with masses between 0.3 and 20 GeV and small couplings not yet probed by existing searches. Such a neutrinophilic mediator is well motivated for addressing the origin of several neutrino-portal dark matter candidates, including thermal freeze-out and sterile-neutrino dark matter scenarios. Experimentally, the corresponding signatures include neutrino charged-current scattering events associated with large missing transverse momentum, and excessive apparent tau-neutrino events. We discuss the FPF detector capabilities needed for this search, most importantly the hadronic energy resolution.


title: Efficient estimation method for time evolution of proto-neutron star mass and radius from supernova neutrino signal

abstract: In this paper we present a novel method to estimate the time evolution of proto-neutron star (PNS) structure from the neutrino signal in core-collapse supernovae (CCSN). Employing recent results of multi-dimensional CCSN simulations, we delve into a relation between total emitted neutrino energy (TONE) and PNS mass/radius, and we find that they are strongly correlated with each other. We fit the relation by simple polynomial functions connecting TONE to PNS mass and radius as a function of time. By combining another fitting function representing the correlation between TONE and cumulative number of event at each neutrino observatory, PNS mass and radius can be retrieved from purely observed neutrino data. We demonstrate retrievals of PNS mass and radius from mock data of neutrino signal, and we assess the capability of our proposed method. While underlining the limitations of the method, we also discuss the importance of the joint analysis with gravitational wave signal. This would reduce uncertainties of parameter estimations in our method, and may narrow down the possible neutrino oscillation model. The proposed method is a very easy and inexpensive computation, which will be useful in real data analysis of CCSN neutrino signal.

categories: astro-ph.HE astro-ph.IM gr-qc hep-ex physics.data-an


title: Interplay between topological protected valley and quantum Hall edge transport

abstract: An established way of realizing topologically protected states in a two-dimensional electron gas is by applying a perpendicular magnetic field thus creating quantum Hall edge channels. In electrostatically gapped bilayer graphene intriguingly, even in the absence of a magnetic field topologically protected electronic states can emerge at naturally occurring stacking domain walls. While individually both types of topologically protected states have been investigated, their intriguing interplay remains poorly understood. Here, we focus on the interplay between topological domain wall states and quantum Hall edge transport within the eight-fold degenerate zeroth Landau level of high-quality suspended bilayer graphene. We find that the two-terminal conductance remains approximately constant for low magnetic fields throughout the distinct quantum Hall states since the conduction channels are traded between domain wall and device edges. For high magnetic fields, however, we observe evidence of transport suppression at the domain wall, which can be attributed to the emergence of spectral minigaps. This indicates that stacking domain walls do potentially do not correspond to a topological domain wall in the order parameter.


title: VLT/MUSE observations of SDSS J1029+2623: towards a high-precision strong lensing model

abstract: We present a strong lensing analysis of the galaxy cluster SDSS J1029+2623 at $z=0.588$, one of the few currently known lens clusters with multiple images of a background ($z=2.1992$) quasar with a measured time delay. We use archival Hubble Space Telescope multi-band imaging and new Multi Unit Spectroscopic Explorer follow-up spectroscopy to build an accurate lens mass model, a crucial step towards future cosmological applications. The spectroscopic data enable the secure identification of 57 cluster members and of two nearby perturbers along the line-of-sight. We estimate the inner kinematics of a sub-set of 20 cluster galaxies to calibrate the scaling relations parametrizing the sub-halo mass component. We also reliably determine the redshift of 4 multiply imaged sources, provide a tentative measurement for one system, and report the discovery of a new four-image system. The final catalog comprises 26 multiple images from 7 background sources, spanning a wide redshift range, from 1.02 to 5.06. We present two parametric lens models, with slightly different cluster mass parametrizations. The observed positions of the multiple images are accurately reproduced within approximately $0''.2$, the three image positions of the quasar within only $\sim0''.1$. We estimate a cluster projected total mass of $M(<300~ {\rm kpc}) \sim 2.1 \times 10^{14}~ M_{\odot}$, with a statistical uncertainty of a few percent. Both models, that include a small galaxy close to one of the quasar images, predict magnitude differences and time delays between the quasar images that are consistent with the observations.


title: Multi-gap topology and non-Abelian braiding of phonons from first principles

abstract: Non-Abelian states of matter, in which the final state depends on the order of the interchanges of two quasiparticles, can encode information immune from environmental noise with the potential to provide a robust platform for topological quantum computation. We demonstrate that phonons can carry non-Abelian frame charges at the band crossing points of their frequency spectrum, and that external stimuli can drive their braiding. We present a general framework to understand the topological configurations of phonons from first principles calculations using a topological invariant called Euler class, and provide a complete analysis of phonon braiding by combining different topological configurations. Taking a well-known dielectric material, Al$_2$O$_3$, as a representative example, we demonstrate that electrostatic doping gives rise to phonon band inversions that can induce redistribution of the frame charges, leading to non-Abelian braiding of phonons. Our work provides a new quasiparticle platform for realizable non-Abelian braiding in reciprocal space, and expands the toolset for studying braiding processes.

categories: cond-mat.mes-hall cond-mat.mtrl-sci physics.comp-ph quant-ph


title: Asymptotically flat black hole solution in modified gravity

abstract: From time to time, different observations suggest that Einstein's theory of general relativity (GR) may not be the ultimate theory of gravity. Various researchers have suggested that the $f(R)$ theory of gravity is the best alternative to replace GR. Using $f(R)$ gravity, one can elucidate the various unexplained physics of compact objects, such as black holes, neutron stars, and white dwarfs. Researchers have already put effort into finding the vacuum solution around a black hole in $f(R)$ gravity. However, for a long time, they could not find an asymptotically flat vacuum solution. In this article, we show that the asymptotically flat vacuum solution of $f(R)$ gravity is possible and thereby use it to explain the spherical accretion flow around the black hole.


title: A Hierarchy for Replica Quantum Advantage

abstract: We prove that given the ability to make entangled measurements on at most $k$ replicas of an $n$-qubit state $\rho$ simultaneously, there is a property of $\rho$ which requires at least order $2^n$ measurements to learn. However, the same property only requires one measurement to learn if we can make an entangled measurement over a number of replicas polynomial in $k, n$. Because the above holds for each positive integer $k$, we obtain a hierarchy of tasks necessitating progressively more replicas to be performed efficiently. We introduce a powerful proof technique to establish our results, and also use this to provide new bounds for testing the mixedness of a quantum state.


title: Higgs Branches of rank-0 5d theories from M-theory on $(A_j,A_l)$ and $(A_k,D_n)$ singularities

abstract: We study the dynamics of M-theory on isolated non-toric Calabi-Yau threefold singularities of type $(A_j,A_l)$ and $(A_k,D_n)$, engineering five-dimensional rank-zero SCFTs. Our approach consists in mapping these backgrounds to type IIA string theory with D6 branes at angles and O6$^{-}$ planes, computing the five-dimensional open string modes located at the brane intersections. This permits us to characterize the Higgs Branches of these theories as algebraic varieties, determine the flavour and gauge group and inspect subtleties related to T-branes. Our methods apply for all the considered singularities: we give a closed formula for the $(A_j,A_l)$ Higgs Branches, and tables for the Higgs Branches of the $(A_k,D_n)$ series.


title: Effective limits on single scalar extensions in the light of recent LHC data

abstract: In this paper, we work with 16 different single scalar particle extensions of the Standard Model. We present the sets of dimension-6 effective operators and the associated Wilson coefficients as functions of model parameters after integrating out the heavy scalars up to 1-loop, including the heavy-light mixing, for each such scenario. Using the correspondence between the effective operators and the observables at electroweak scale, and employing Bayesian statistics, we compute the allowed ranges of new physics parameters that are further translated and depicted in 2-dimensional Wilson coefficient space in the light of the latest CMS and ATLAS data up to $137 \text{ fb}^{-1}$ and $139\text{ fb}^{-1}$, respectively. We also adjudge the status of those new physics extensions that offer similar sets of relevant effective operators. In addition, we provide a model-independent fit of $23$ Standard Model effective field theory Wilson coefficients using electroweak precision observables, single and di-Higgs data as well as kinematic distributions of di-boson production.


title: Spectral shift technique for strongly correlated lattice fermions

abstract: We present a development of strong-coupling diagrammatic techniques which relies on integrating out mean-field-like paths prior to conducting the expansion. This makes it possible to expand around a state with a quasiparticle spectrum that takes into account all kinetic effects that do not explicitly depend on nonlocal correlations. These paths contribute most of the kinetic energy in correlated systems, and so this protocol provides a starting point for the expansion that more closely resembles the full theory. Comparisons to existing exact results clearly confirm this.


title: Natural evidence for fuzzy sphere noncommutative geometry: super-Chandrasekhar white dwarfs

abstract: Noncommutative geometry is one of the quantum gravity theories, which various researchers have been using to describe different physical and astrophysical systems. However, so far, no direct observations can justify its existence, and this theory remains a hypothesis. On the other hand, over the past two decades, more than a dozen over-luminous type Ia supernovae have been observed, which indirectly predict that they originate from white dwarfs with super-Chandrasekhar masses $2.1-2.8 \rm\,M_\odot$. In this article, we discuss that considering white dwarfs as squashed fuzzy spheres, a class of noncommutative geometry, helps in accumulating more mass than the Chandrasekhar mass-limit. The length-scale beyond which the effect of noncommutativity becomes prominent is an emergent phenomenon, which depends only on the inter-electron separations in the white dwarf.


title: Insights on Entanglement Entropy in $1+1$ Dimensional Causal Sets

abstract: Entanglement entropy in causal sets offers a fundamentally covariant characterisation of quantum field degrees of freedom. A known result in this context is that the degrees of freedom consist of a number of contributions that have continuum-like analogues, in addition to a number of contributions that do not. The latter exhibit features below the discreteness scale and are excluded from the entanglement entropy using a "truncation scheme". This truncation is necessary to recover the standard spatial area law of entanglement entropy. In this paper we build on previous work on the entanglement entropy of a massless scalar field on a causal set approximated by a 1+1D causal diamond in Minkowski spacetime. We present new insights into the truncated contributions, including evidence that they behave as fluctuations and encode features specific to a particular causal set sprinkling. We extend previous results in the massless theory to include R\'enyi entropies and include new results for the massive theory. We also discuss the implications of our work for the treatment of entanglement entropy in causal sets in more general settings.


title: Generalized Tevelev degrees of $\mathbb{P}^1$

abstract: Let $(C,p_1,\ldots,p_n)$ be a general curve. We consider the problem of enumerating covers of the projective line by $C$ subject to incidence conditions at the marked points. These counts have been obtained by the first named author with Pandharipande and Schmitt via intersection theory on Hurwitz spaces and by the second named author with Farkas via limit linear series. In this paper, we build on these two approaches to generalize these counts to the situation where the covers are constrained to have arbitrary ramification profiles: that is, additional ramification conditions are imposed at the marked points, and some collections of marked points are constrained to have equal image.


title: Exponential separations between learning with and without quantum memory

abstract: We study the power of quantum memory for learning properties of quantum systems and dynamics, which is of great importance in physics and chemistry. Many state-of-the-art learning algorithms require access to an additional external quantum memory. While such a quantum memory is not required a priori, in many cases, algorithms that do not utilize quantum memory require much more data than those which do. We show that this trade-off is inherent in a wide range of learning problems. Our results include the following: (1) We show that to perform shadow tomography on an $n$-qubit state rho with $M$ observables, any algorithm without quantum memory requires $\Omega(\min(M, 2^n))$ samples of rho in the worst case. Up to logarithmic factors, this matches the upper bound of [HKP20] and completely resolves an open question in [Aar18, AR19]. (2) We establish exponential separations between algorithms with and without quantum memory for purity testing, distinguishing scrambling and depolarizing evolutions, as well as uncovering symmetry in physical dynamics. Our separations improve and generalize prior work of [ACQ21] by allowing for a broader class of algorithms without quantum memory. (3) We give the first tradeoff between quantum memory and sample complexity. We prove that to estimate absolute values of all $n$-qubit Pauli observables, algorithms with $k < n$ qubits of quantum memory require at least $\Omega(2^{(n-k)/3})$ samples, but there is an algorithm using $n$-qubit quantum memory which only requires $O(n)$ samples. The separations we show are sufficiently large and could already be evident, for instance, with tens of qubits. This provides a concrete path towards demonstrating real-world advantage for learning algorithms with quantum memory.


title: A Histopathology Study Comparing Contrastive Semi-Supervised and Fully Supervised Learning

abstract: Data labeling is often the most challenging task when developing computational pathology models. Pathologist participation is necessary to generate accurate labels, and the limitations on pathologist time and demand for large, labeled datasets has led to research in areas including weakly supervised learning using patient-level labels, machine assisted annotation and active learning. In this paper we explore self-supervised learning to reduce labeling burdens in computational pathology. We explore this in the context of classification of breast cancer tissue using the Barlow Twins approach, and we compare self-supervision with alternatives like pre-trained networks in low-data scenarios. For the task explored in this paper, we find that ImageNet pre-trained networks largely outperform the self-supervised representations obtained using Barlow Twins.


title: Towards the use of conservative thermodynamic variables in data assimilation: a case study using ground-based microwave radiometer measurements

abstract: This study aims at introducing two conservative thermodynamic variables (moist-air entropy potential temperature and total water content) into a one-dimensional variational data assimilation system (1D-Var) to demonstrate the benefit for future operational assimilation schemes. This system is assessed using microwave brightness temperatures from a ground-based radiometer installed during the field campaign SOFOG3D dedicated to fog forecast improvement. An underlying objective is to ease the specification of background error covariance matrices that are highly dependent on weather conditions when using classical variables, making difficult the optimal retrievals of cloud and thermodynamic properties during fog conditions. Background error covariance matrices for these new conservative variables have thus been computed by an ensemble approach based on the French convective scale model AROME, for both all-weather and fog conditions. A first result shows that the use of these matrices for the new variables reduces some dependencies to the meteorological conditions (diurnal cycle, presence or not of clouds) compared to usual variables (temperature, specific humidity). Then, two 1D-Var experiments (classical vs. conservative variables) are evaluated over a full diurnal cycle characterized by a stratus-evolving radiative fog situation, using hourly brightness temperatures. Results show, as expected, that analysed brightness temperatures by the 1D-Var are much closer to the observed ones than background values for both variable choices. This is especially the case for channels sensitive to water vapour and liquid water. On the other hand, analysis increments in model space (water vapour, liquid water) show significant differences between the two sets of variables.


title: Search in Imperfect Information Games

abstract: From the very dawn of the field, search with value functions was a fundamental concept of computer games research. Turing's chess algorithm from 1950 was able to think two moves ahead, and Shannon's work on chess from $1950$ includes an extensive section on evaluation functions to be used within a search. Samuel's checkers program from 1959 already combines search and value functions that are learned through self-play and bootstrapping. TD-Gammon improves upon those ideas and uses neural networks to learn those complex value functions -- only to be again used within search. The combination of decision-time search and value functions has been present in the remarkable milestones where computers bested their human counterparts in long standing challenging games -- DeepBlue for Chess and AlphaGo for Go. Until recently, this powerful framework of search aided with (learned) value functions has been limited to perfect information games. As many interesting problems do not provide the agent perfect information of the environment, this was an unfortunate limitation. This thesis introduces the reader to sound search for imperfect information games.


title: Predicting Lattice Phonon Vibrational Frequencies Using Deep Graph Neural Networks

abstract: Lattice vibration frequencies are related to many important materials properties such as thermal and electrical conductivity as well as superconductivity. However, computational calculation of vibration frequencies using density functional theory (DFT) methods is too computationally demanding for a large number of samples in materials screening. Here we propose a deep graph neural network-based algorithm for predicting crystal vibration frequencies from crystal structures with high accuracy. Our algorithm addresses the variable dimension of vibration frequency spectrum using the zero padding scheme. Benchmark studies on two data sets with 15,000 and 35,552 samples show that the aggregated $R^2$ scores of the prediction reaches 0.554 and 0.724 respectively. Our work demonstrates the capability of deep graph neural networks to learn to predict phonon spectrum properties of crystal structures in addition to phonon density of states (DOS) and electronic DOS in which the output dimension is constant.


title: Isotropy statistics of CMB hot and cold spots

abstract: Statistical Isotropy of the Cosmic Microwave Background (CMB) radiation has been studied and debated extensively in recent years. Under this assumption, the hot spots and cold spots of the CMB are expected to be uniformly distributed over a 2-sphere. We use the orientation matrix, first proposed by Watson (1965) and Scheidegger (1965) and associated shape and strength parameters (Woodcock, 1977) to analyse whether the hot and cold spots of the observed CMB temperature anisotropy field are uniformly placed. We demonstrate the usefulness of our estimators by using simulated toy models containing non-uniform data. We apply our method on several foreground minimized CMB maps observed by WMAP and Planck over large angular scales. The shape and strength parameters constrain geometric features of possible deviations from uniformity (isotropy) and the power of the anomalous signal. We find that distributions of hot or cold spots in cleaned maps show no unusual signature of clustering or girdling. Instead, we notice a strikingly uniform distribution of hot spots over the full sky. The signal remains robust with respect to the four cleaned maps used and presence or absence of the non-Gaussian cold spot (NGCS). On the partial sky with WMAP KQ75 and Planck U73 masks we find anomalously low strength of non-uniformity for cold spots which is found to be robust with respect to various cleaning methods, masks applied, instruments, frequencies, and the presence or absence of the NGCS. Interestingly we find that the signal of anomalously weak non-uniformity could be due to contributions from the quadrupole and octupole and may be related with the low CMB temperature variance anomaly.


title: The Covariance of Squeezed Bispectrum Configurations

abstract: We measure the halo bispectrum covariance in a large set of N-body simulations and compare it with theoretical expectations. We find a large correlation among (even mildly) squeezed halo bispectrum configurations. A similarly large correlation can be found between squeezed triangles and the long-wavelength halo power spectrum. This shows that the diagonal Gaussian contribution fails to describe, even approximately, the full covariance in these cases. We compare our numerical estimate with a model that includes, in addition to the Gaussian one, only the non-Gaussian terms that are large for squeezed configurations. We find that accounting for these large terms in the modeling greatly improves the agreement of the full covariance with simulations. We apply these results to a simple Fisher matrix forecast, and find that constraints on primordial non-Gaussianity are degraded by a factor of $\sim 2$ when a non-Gaussian covariance is assumed instead of the diagonal, Gaussian approximation.


title: Normalizers of chains of discrete $p$-toral subgroups in compact Lie groups

abstract: In this paper we study the normalizer decomposition of a compact Lie group $G$ using the information of the fusion system $\mathcal{F}$ of $G$ on a maximal discrete $p$-toral subgroup. We prove that there is an injective map from the set of conjugacy classes of chains of $\mathcal{F}$-centric, $\mathcal{F}$-radical discrete $p$-toral subgroups to the set of conjugacy classes of chains of $p$-centric, $p$-stubborn continuous $p$-toral subgroups. The map is a bijection when $\pi_0(G)$ is a finite $p$-group. We also prove that the classifying space of the normalizer of a chain of discrete $p$-toral subgroups of $G$ is mod $p$ equivalent to the classifying space of the normalizer of the corresponding chain of $p$-toral subgroups.


title: Very-High Dynamic Range, 10,000 frames/second Pixel Array Detector for Electron Microscopy

abstract: Precision and accuracy of quantitative scanning transmission electron microscopy (STEM) methods such as ptychography, and the mapping of electric, magnetic and strain fields depend on the dose. Reasonable acquisition time requires high beam current and the ability to quantitatively detect both large and minute changes in signal. A new hybrid pixel array detector (PAD), the second-generation Electron Microscope Pixel Array Detector (EMPAD-G2), addresses this challenge by advancing the technology of a previous generation PAD, the EMPAD. The EMPAD-G2 images continuously at a frame-rates up to 10 kHz with a dynamic range that spans from low-noise detection of single electrons to electron beam currents exceeding 180 pA per pixel, even at electron energies of 300 keV. The EMPAD-G2 enables rapid collection of high-quality STEM data that simultaneously contain full diffraction information from unsaturated bright field disks to usable Kikuchi bands and higher-order Laue zones. Test results from 80 to 300 keV are presented, as are first experimental results demonstrating ptychographic reconstructions, strain and polarization maps. We introduce a new information metric, the Maximum Usable Imaging Speed (MUIS), to identify when a detector becomes electron-starved, saturated or its pixel count is mismatched with the beam current.


title: Multimodal End-to-End Group Emotion Recognition using Cross-Modal Attention

abstract: Classifying group-level emotions is a challenging task due to complexity of video, in which not only visual, but also audio information should be taken into consideration. Existing works on multimodal emotion recognition are using bulky approach, where pretrained neural networks are used as a feature extractors and then extracted features are being fused. However, this approach does not consider attributes of multimodal data and feature extractors cannot be fine-tuned for specific task which can be disadvantageous for overall model accuracy. To this end, our impact is twofold: (i) we train model end-to-end, which allows early layers of neural network to be adapted with taking into account later, fusion layers, of two modalities; (ii) all layers of our model was fine-tuned for downstream task of emotion recognition, so there were no need to train neural networks from scratch. Our model achieves best validation accuracy of 60.37% which is approximately 8.5% higher, than VGAF dataset baseline and is competitive with existing works, audio and video modalities.


title: A Portable and Passive Gravity Compensation Arm Support for Drone Teleoperation

abstract: Gesture-based interfaces are often used to achieve a more natural and intuitive teleoperation of robots. Yet, sometimes, gesture control requires postures or movements that cause significant fatigue to the user. In a previous user study, we demonstrated that na\"ive users can control a fixed-wing drone with torso movements while their arms are spread out. However, this posture induced significant arm fatigue. In this work, we present a passive arm support that compensates the arm weight with a mean torque error smaller than 0.005 N/kg for more than 97% of the range of motion used by subjects to fly, therefore reducing muscular fatigue in the shoulder of on average 58%. In addition, this arm support is designed to fit users from the body dimension of the 1st percentile female to the 99th percentile male. The performance analysis of the arm support is described with a mechanical model and its implementation is validated with both a mechanical characterization and a user study, which measures the flight performance, the shoulder muscle activity and the user acceptance.


title: Concurrent Donsker-Varadhan and hydrodynamical large deviations

abstract: We consider the weakly asymmetric exclusion process on the $d$-dimensional torus. We prove a large deviations principle for the time averaged empirical density and current in the joint limit in which both the time interval and the number of particles diverge. This result is obtained both by analyzing the variational convergence, as the number of particles diverges, of the Donsker-Varadhan functional for the empirical process and by considering the large time behavior of the hydrodynamical rate function. The large deviations asymptotic of the time averaged current is then deduced by contraction principle. The structure of the minimizers of this variational problem corresponds to the possible occurrence of dynamical phase transitions.


title: Asymptotic resemblance relations on groups

abstract: In this paper, we study properties of asymptotic resemblance relations induced by compatible coarse structures on groups. We generalize the notion of asymptotic dimensiongrad for groups with compatible coarse structures and show this notion is coarse invariant. We end by defining the notion of set theoretic coupling for groups with compatible coarse structures and showing this notion is the generalization of the notion of topological coupling for finitely generated groups. We show if two groups with compatible coarse structures admit a set theoretic coupling then they are asymptotic equivalent.


title: Graph Neural Network Training with Data Tiering

abstract: Graph Neural Networks (GNNs) have shown success in learning from graph-structured data, with applications to fraud detection, recommendation, and knowledge graph reasoning. However, training GNN efficiently is challenging because: 1) GPU memory capacity is limited and can be insufficient for large datasets, and 2) the graph-based data structure causes irregular data access patterns. In this work, we provide a method to statistical analyze and identify more frequently accessed data ahead of GNN training. Our data tiering method not only utilizes the structure of input graph, but also an insight gained from actual GNN training process to achieve a higher prediction result. With our data tiering method, we additionally provide a new data placement and access strategy to further minimize the CPU-GPU communication overhead. We also take into account of multi-GPU GNN training as well and we demonstrate the effectiveness of our strategy in a multi-GPU system. The evaluation results show that our work reduces CPU-GPU traffic by 87-95% and improves the training speed of GNN over the existing solutions by 1.6-2.1x on graphs with hundreds of millions of nodes and billions of edges.


title: A Generic Deep Learning Based Cough Analysis System from Clinically Validated Samples for Point-of-Need Covid-19 Test and Severity Levels

abstract: We seek to evaluate the detection performance of a rapid primary screening tool of Covid-19 solely based on the cough sound from 8,380 clinically validated samples with laboratory molecular-test (2,339 Covid-19 positives and 6,041 Covid-19 negatives). Samples were clinically labeled according to the results and severity based on quantitative RT-PCR (qRT-PCR) analysis, cycle threshold, and lymphocytes count from the patients. Our proposed generic method is an algorithm based on Empirical Mode Decomposition (EMD) with subsequent classification based on a tensor of audio features and a deep artificial neural network classifier with convolutional layers called DeepCough'. Two different versions of DeepCough based on the number of tensor dimensions, i.e. DeepCough2D and DeepCough3D, have been investigated. These methods have been deployed in a multi-platform proof-of-concept Web App CoughDetect to administer this test anonymously. Covid-19 recognition results rates achieved a promising AUC (Area Under Curve) of 98.800.83%, sensitivity of 96.431.85%, and specificity of 96.201.74%, and 81.08%5.05% AUC for the recognition of three severity levels. Our proposed web tool and underpinning algorithm for the robust, fast, point-of-need identification of Covid-19 facilitates the rapid detection of the infection. We believe that it has the potential to significantly hamper the Covid-19 pandemic across the world.


title: Generating faster algorithms for d-Path Vertex Cover

abstract: Many algorithms which exactly solve hard problems require branching on more or less complex structures in order to do their job. Those who design such algorithms often find themselves doing a meticulous analysis of numerous different cases in order to identify these structures and design suitable branching rules, all done by hand. This process tends to be error prone and often the resulting algorithm may be difficult to implement in practice. In this work, we aim to automate a part of this process and focus on simplicity of the resulting implementation. We showcase our approach on the following problem. For a constant $d$, the $d$-Path Vertex Cover problem ($d$-PVC) is as follows: Given an undirected graph and an integer $k$, find a subset of at most $k$ vertices of the graph, such that their deletion results in a graph not containing a path on $d$ vertices as a subgraph. We develop a fully automated framework to generate parameterized branching algorithms for the problem and obtain algorithms outperforming those previously known for $3 \le d \le 8$. E.g., we show that $5$-PVC can be solved in $O(2.7^k\cdot n^{O(1)})$ time.


title: Persia: An Open, Hybrid System Scaling Deep Learning-based Recommenders up to 100 Trillion Parameters

abstract: Deep learning based models have dominated the current landscape of production recommender systems. Furthermore, recent years have witnessed an exponential growth of the model scale--from Google's 2016 model with 1 billion parameters to the latest Facebook's model with 12 trillion parameters. Significant quality boost has come with each jump of the model capacity, which makes us believe the era of 100 trillion parameters is around the corner. However, the training of such models is challenging even within industrial scale data centers. This difficulty is inherited from the staggering heterogeneity of the training computation--the model's embedding layer could include more than 99.99% of the total model size, which is extremely memory-intensive; while the rest neural network is increasingly computation-intensive. To support the training of such huge models, an efficient distributed training system is in urgent need. In this paper, we resolve this challenge by careful co-design of both the optimization algorithm and the distributed system architecture. Specifically, in order to ensure both the training efficiency and the training accuracy, we design a novel hybrid training algorithm, where the embedding layer and the dense neural network are handled by different synchronization mechanisms; then we build a system called Persia (short for parallel recommendation training system with hybrid acceleration) to support this hybrid training algorithm. Both theoretical demonstration and empirical study up to 100 trillion parameters have conducted to justified the system design and implementation of Persia. We make Persia publicly available (at https://github.com/PersiaML/Persia) so that anyone would be able to easily train a recommender model at the scale of 100 trillion parameters.


title: Beyond Importance Scores: Interpreting Tabular ML by Visualizing Feature Semantics

abstract: Interpretability is becoming an active research topic as machine learning (ML) models are more widely used to make critical decisions. Tabular data is one of the most commonly used modes of data in diverse applications such as healthcare and finance. Much of the existing interpretability methods used for tabular data only report feature-importance scores -- either locally (per example) or globally (per model) -- but they do not provide interpretation or visualization of how the features interact. We address this limitation by introducing Feature Vectors, a new global interpretability method designed for tabular datasets. In addition to providing feature-importance, Feature Vectors discovers the inherent semantic relationship among features via an intuitive feature visualization technique. Our systematic experiments demonstrate the empirical utility of this new method by applying it to several real-world datasets. We further provide an easy-to-use Python package for Feature Vectors.


title: On power integral bases of certain pure number fields defined by $X^{60}-m$

abstract: Let $K$ be a pure number field generated by a complex root of a monic irreducible polynomial $F(x)=x^{60}-m\in \mathbb{Z}[x]$, with $m\neq \pm1$ a square free integer. In this paper, we study the monogeneity of $K$. We prove that if $m\not\equiv 1\md{4}$, $m\not\equiv \mp 1 \md{9} $ and $\overline{m}\not\in\{\mp 1,\mp 7\} \md{25}$, then $K$ is monogenic. But if $m\equiv 1\md{4}$, $m\equiv \mp1 \md{9}$, or $m\equiv \mp 1\md{25}$, then $K$ is not monogenic. Our results are illustrated by examples.


title: Local Normal Forms of Noncommutative Functions

abstract: This article describes local normal forms of functions in noncommuting variables, up to equivalence generated by isomorphism of noncommutative Jacobi algebras, extending singularity theory in the style of Arnold's commutative local normal forms into the noncommutative realm. This generalisation unveils many new phenomena, including an ADE classification when the Jacobi ring has dimension zero and, by suitably taking limits, a further ADE classification in dimension one. These are natural generalisations of the simple singularities and those with infinite multiplicity in Arnold's classification. We obtain normal forms away from some exceptional Type E cases. Remarkably these normal forms have no moduli, and the key new feature is that the noncommutative world affords larger families. This theory has a range of immediate consequences to the birational geometry of 3-folds. The normal forms of dimension zero are the analytic classification of smooth 3-fold flops, and one outcome of NC singularity theory is the first complete list of all Type D flopping germs, with variants covering Type E. The normal forms of dimension one have further ramifications to divisorial contractions to a curve. In addition, the general techniques also give strong evidence towards new contractibility criteria for rational curves.


title: An Extensive Study of User Identification via Eye Movements across Multiple Datasets

abstract: Several studies have reported that biometric identification based on eye movement characteristics can be used for authentication. This paper provides an extensive study of user identification via eye movements across multiple datasets based on an improved version of method originally proposed by George and Routray. We analyzed our method with respect to several factors that affect the identification accuracy, such as the type of stimulus, the IVT parameters (used for segmenting the trajectories into fixation and saccades), adding new features such as higher-order derivatives of eye movements, the inclusion of blink information, template aging, age and gender.We find that three methods namely selecting optimal IVT parameters, adding higher-order derivatives features and including an additional blink classifier have a positive impact on the identification accuracy. The improvements range from a few percentage points, up to an impressive 9 % increase on one of the datasets.


title: Heteroclinic cycling and extinction in May-Leonard models with demographic stochasticity

abstract: May and Leonard (SIAM J. Appl. Math 1975) introduced a three-species Lotka-Volterra type population model that exhibits heteroclinic cycling. Rather than producing a periodic limit cycle, the trajectory takes longer and longer to complete each "cycle", passing closer and closer to unstable fixed points in which one population dominates and the others approach zero. Aperiodic heteroclinic dynamics have subsequently been studied in ecological systems (side-blotched lizards; colicinogenic E. coli), in the immune system, in neural information processing models ("winnerless competition"), and in models of neural central pattern generators. Yet as May and Leonard observed "Biologically, the behavior (produced by the model) is nonsense. Once it is conceded that the variables represent animals, and therefore cannot fall below unity, it is clear that the system will, after a few cycles, converge on some single population, extinguishing the other two." Here, we explore different ways of introducing discrete stochastic dynamics based on May and Leonard's ODE model, with application to ecological population dynamics, and to a neuromotor central pattern generator system. We study examples of several quantitatively distinct asymptotic behaviors, including total extinction of all species, extinction to a single species, and persistent cyclic dominance with finite mean cycle length.


title: Particle Acceleration by Pickup Process Upstream of Relativistic Shocks

abstract: Particle acceleration at magnetized purely perpendicular relativistic shocks in electron-ion plasmas are studied by means of two-dimensional particle-in-cell simulations. Magnetized shocks with the upstream bulk Lorentz factor $\gamma_1 \gg 1$ are known to emit intense electromagnetic waves from the shock front, which induce electrostatic plasma waves (wakefield) and transverse filamentary structures in the upstream region via the stimulated/induced Raman scattering and the filamentation instability, respectively. The wakefield and filaments inject a fraction of incoming particles into a particle acceleration process, in which particles are once decoupled from the upstream bulk flow by the wakefield, and are piked up again by the flow. The picked-up particles are accelerated by the motional electric field. The maximum attainable Lorentz factor is estimated as $\gamma_{max,e} \sim \alpha\gamma_1^3$ for electrons and $\gamma_{max,i} \sim (1+m_e\gamma_1/m_i)\gamma_1^2$ for ions, where $\alpha \sim 10$ is determined from our simulation results. $\alpha$ can increase up to $\gamma_1$ for weakly magnetized shock if $\gamma_1$ is sufficiently large. This result indicates that highly relativistic astrophysical shocks such as external shocks of gamma-ray bursts can be an efficient particle accelerator.


title: An Optimal Engagement Zone Avoidance Scenario in 2-D

abstract: In this paper, an optimal control problem is considered where a target vehicle aims to reach a desired location in minimum time while avoiding a dynamic engagement zone. Using simple motion, four potential approaches are considered. First, the min-time strategy which ignores the engagement zone is posed and solved. Second, the min-time strategy which avoids the engagement zone entirely is considered. Third, the min-time strategy which allows for some time in the engagement zone; but, still strives to stay away from the center of the engagement zone is posed. Lastly, a fixed final-time strategy is considered, wherein the target tries to avoid the engagement zone; but, is required to arrive at the desired location at a specific time. Using a nonlinear program solver, the optimal strategies are numerically solved. From the results of the numeric solutions, the optimal strategies are discussed and comparisons are drawn.


title: Model independent search for transient multimessenger events with AMON using outlier detection methods

abstract: The Astrophysical Multimessenger Observatory Network (AMON) receives subthreshold data from multiple observatories in order to look for coincidences. Combining more than two datasets at the same time is challenging because of the range of possible signals (time windows, energies, number of events...). However, outlier detection methods can circumvent this issue by identifying any signal divergent from the background (e.g. scrambled data). We propose to use these methods to make a model independent combination of the subthreshold data of neutrino and gamma ray experiments. Using the python outlier detection (PyOD) package, it allows us to test several methods from a simple "k-nearest neighbours" algorithm to a more sophisticated Generative Adversarial Active Learning neural networks which generates data points to better discriminate inliers from outliers.


title: Optimal control of a heroin epidemic mathematical model

abstract: A heroin epidemic mathematical model with prevention information and treatment, as control interventions, is analyzed, assuming that an individual's behavioral response depends on the spreading of information about the effects of heroin. Such information creates awareness, which helps individuals to participate in preventive education and self-protective schemes with additional efforts. We prove that the basic reproduction number is the threshold of local stability of a drug-free and endemic equilibrium. Then, we formulate an optimal control problem to minimize the total number of drug users and the cost associated with prevention education measures and treatment. We prove existence of an optimal control and derive its characterization through Pontryagin's maximum principle. The resulting optimality system is solved numerically. We observe that among all possible strategies, the most effective and cost-less is to implement both control policies.


title: The concentration inequality for a discrete height function model perturbed by random potential

abstract: We present a proof of the concentration inequality for a discrete random surface model, where the underlying potential is perturbed by an additive random potential. The proof is based on annealing the random potential, and follows the method of [CEP96] and other works. Our result demonstrates the robustness of this method.


title: Weak degeneracy of graphs

abstract: Motivated by the study of greedy algorithms for graph coloring, we introduce a new graph parameter, which we call weak degeneracy. By definition, every $d$-degenerate graph is also weakly $d$-degenerate. On the other hand, if $G$ is weakly $d$-degenerate, then $\chi(G) \leq d + 1$ (and, moreover, the same bound holds for the list-chromatic and even the DP-chromatic number of $G$). It turns out that several upper bounds in graph coloring theory can be phrased in terms of weak degeneracy. For example, we show that planar graphs are weakly $4$-degenerate, which implies Thomassen's famous theorem that planar graphs are $5$-list-colorable. We also prove a version of Brooks's theorem for weak degeneracy: a connected graph $G$ of maximum degree $d \geq 3$ is weakly $(d-1)$-degenerate unless $G \cong K_{d + 1}$. (By contrast, all $d$-regular graphs have degeneracy $d$.) We actually prove an even stronger result, namely that for every $d \geq 3$, there is $\epsilon > 0$ such that if $G$ is a graph of weak degeneracy at least $d$, then either $G$ contains a $(d+1)$-clique or the maximum average degree of $G$ is at least $d + \epsilon$. Finally, we show that graphs of maximum degree $d$ and either of girth at least $5$ or of bounded chromatic number are weakly $(d - \Omega(\sqrt{d}))$-degenerate, which is best possible up to the value of the implied constant.


title: Rise and fall of Yu-Shiba-Rusinov bound-states in charge conserving $s$-wave one-dimensional superconductors

abstract: We re-examine the problem of a magnetic impurity coupled to a superconductor focusing on the role of quantum fluctuations. We study in detail, a system that consists of a one-dimensional charge conserving spin-singlet superconductor coupled to a boundary magnetic impurity. Our main finding is that quantum fluctuations lead to the destruction of Yu-Shiba-Rusinov (YSR) intra-gap bound-states in all but a narrow region of the phase diagram. We carry out our analysis in three stages, increasing the role of the quantum fluctuations at each stage. First we consider the limit of a classical impurity and study the bulk semiclassically, finding YSR states throughout the phase diagram, a situation similar to conventional BCS superconductors. In the second stage, we reintroduce quantum fluctuations in the bulk and find that the YSR state is suppressed over half of the phase diagram, existing only around the quantum critical point separating the unscreened and the partially screened phases. In the final stage we solve exactly the full interacting model with arbitrary coupling constants using Bethe Ansatz. We find that including both the quantum fluctuating bulk and quantum impurity destabilizes the YSR state over most of the phase diagram allowing it to exist only in a small region, the YSR regime, between a Kondo-screened and an unscreened regime. Within the YSR regime a first order phase transition occurs between a spin singlet and doublet ground state. We also find that for large enough impurity spin exchange interaction a renormalized Kondo-screened regime is established. In this regime, not found for BCS superconductors, there is no YSR state and a renormalized Kondo temperature scale is generated.


title: Combining Electron Spin Resonance Spectroscopy with Scanning Tunneling Microscopy at High Magnetic Fields

abstract: Magnetic media remain a key in information storage and processing. The continuous increase of storage densities and the desire for quantum memories and computers pushes the limits of magnetic characterisation techniques. Ultimately, a tool which is capable of coherently manipulating and detecting individual quantum spins is needed. The scanning tunnelling microscope (STM) is the only technique which unites the prerequisites of high spatial and energy resolution, low temperature and high magnetic fields to achieve this goal. Limitations in the available frequency range for electron spin resonance STM (ESR-STM) mean that many instruments operate in the thermal noise regime. We resolve challenges in signal delivery to extend the operational frequency range of ESR-STM by more than a factor of two and up to 100GHz, making the Zeeman energy the dominant energy scale at achievable cryogenic temperatures of a few hundred millikelvin. We present a general method for augmenting existing instruments into ESR-STMs to investigate spin dynamics in the high-field limit. We demonstrate the performance of the instrument by analysing inelastic tunnelling in a junction driven by a microwave signal and provide proof of principle measurements for ESR-STM.


title: Activity cycles in RS CVn type stars

abstract: We compile a list of 121 RS CVn type stars from the bibliography in southern hemisphere, to search for activity cycles, covering a large range of luminosities and rotation periods. For each system of the list, we obtain photometric data from the All Sky Automated Survey (ASAS), and we complement it with our own photometry, obtained with the Optical Robotic Observatory (ORO). We analyze this data with the Generalized Lomb-Scargle periodogram to determine the rotation period and possible activity cycles for each system. We found rotation periods for 102 systems and activity cycles for 91 systems. From the statistical analysis of the results, we found that giant stars behave differently than subgiants and main-sequence stars, and that there is a good correlation between rotation periods and luminosities.


title: Functional Connectivity of the Brain Across Rodents and Humans

abstract: Resting-state functional magnetic resonance imaging (rs-fMRI), which measures the spontaneous fluctuations in the blood oxygen level-dependent (BOLD) signal, is increasingly utilized for the investigation of the brain's physiological and pathological functional activity. Rodents, as a typical animal model in neuroscience, play an important role in the studies that examine the neuronal processes that underpin the spontaneous fluctuations in the BOLD signal and the functional connectivity that results. Translating this knowledge from rodents to humans requires a basic knowledge of the similarities and differences across species in terms of both the BOLD signal fluctuations and the resulting functional connectivity. This review begins by examining similarities and differences in anatomical features, acquisition parameters, and preprocessing techniques, as factors that contribute to functional connectivity. Homologous functional networks are compared across species, and aspects of the BOLD fluctuations such as the topography of the global signal and the relationship between structural and functional connectivity are examined. Time-varying features of functional connectivity, obtained by sliding windowed approaches, quasi-periodic patterns, and coactivation patterns, are compared across species. Applications demonstrating the use of rs-fMRI as a translational tool for cross-species analysis are discussed, with an emphasis on neurological and psychiatric disorders. Finally, open questions are presented to encapsulate the future direction of the field.


title: An Agmon-Allegretto-Piepenbrink principle for Schroedinger operators

abstract: We prove that each Borel function $V : \Omega \to [-\infty, +\infty]$ defined on an open subset $\Omega \subset \mathbb{R}^{N}$ induces a decomposition $\Omega = S \cup \bigcup_{i} D_{i}$ such that every function in $W^{1,2}_{0}(\Omega) \cap L^{2}(\Omega; V^{+} dx)$ is zero almost everywhere on $S$ and existence of nonnegative supersolutions of $-\Delta + V$ on each component $D_{i}$ yields nonnegativity of the associated quadratic form $\int_{D_{i}} (|\nabla \xi|^2+V\xi^2)$.


title: Monte Carlo Event Generator updates, for tau pair events at Belle II energies

abstract: The Monte Carlo for lepton pair production and tau decays consist of KKMC for lepton pair production, tauola for tau lepton decays and photos for radiative corrections in decays. An effort for adaptation of the system for precision data to be collected at Belle II experiment lead to extension of phase space generation modules both in photos and tauola to enable decays and/or radiative corrections of additional light lepton pairs. The phase-space and matrix element parts of the programs are separated, that is why extension to processes where lepton pair is produced through narrow resonances, like dark photon, was straight forward. In the present version of tauola, list of tau decay channels is enriched with multitude of exotic ones, useful for searches of new physics. The hadronic currents parametrizations of main decay channels is prepared for basic simulation in the experiment. The basis for future work on precise fits of hadronic currents including Machine Learning is retained, but development of necessary software solutions is left for the forthcoming years. Presented programs versions are available in stand-alone format from gitlab or through the basf2 system of Belle II software. Official distribution web pages, documenting programs tests, are retained but not necessarily up to the date.


title: The Central 300 pc of the Galaxy Probed by Infrared Spectra of H$_3^+$ and CO: III. Locations of Sgr B2 and Star $\iota$

abstract: Using a simple relation between the radial expansion velocity of diffuse gas in the Central Molecular Zone (CMZ) of the Galaxy and its distance from Sgr A$^\ast$ we estimate the physical depths within the CMZ of Star $\iota$ (2MASS $J$17470898-2829561) and two Sgr B2 far-infrared continuum sources with respect to the location of Sgr A$^\ast$. To do this we use velocity profiles of mid-infrared absorption lines of H$_3^+$ and of far-infrared absorption lines of H$_2$O$^+$, OH$^+$, and $^{13}$CH$^+$. The distances to Star $\iota$ and to the Sgr B2 sources are found to be $\sim$90 pc greater than the distance to Sgr A$^\ast$. Our conclusion that Sgr B2 lies toward the rear of the CMZ is contrary to many previous models in which it has been placed shallower than Sgr A$^\ast$.


title: Dance In the Wild: Monocular Human Animation with Neural Dynamic Appearance Synthesis

abstract: Synthesizing dynamic appearances of humans in motion plays a central role in applications such as AR/VR and video editing. While many recent methods have been proposed to tackle this problem, handling loose garments with complex textures and high dynamic motion still remains challenging. In this paper, we propose a video based appearance synthesis method that tackles such challenges and demonstrates high quality results for in-the-wild videos that have not been shown before. Specifically, we adopt a StyleGAN based architecture to the task of person specific video based motion retargeting. We introduce a novel motion signature that is used to modulate the generator weights to capture dynamic appearance changes as well as regularizing the single frame based pose estimates to improve temporal coherency. We evaluate our method on a set of challenging videos and show that our approach achieves state-of-the art performance both qualitatively and quantitatively.


title: Recognition of Patient Groups with Sleep Related Disorders using Bio-signal Processing and Deep Learning

abstract: Accurately diagnosing sleep disorders is essential for clinical assessments and treatments. Polysomnography (PSG) has long been used for detection of various sleep disorders. In this research, electrocardiography (ECG) and electromayography (EMG) have been used for recognition of breathing and movement-related sleep disorders. Bio-signal processing has been performed by extracting EMG features exploiting entropy and statistical moments, in addition to developing an iterative pulse peak detection algorithm using synchrosqueezed wavelet transform (SSWT) for reliable extraction of heart rate and breathing-related features from ECG. A deep learning framework has been designed to incorporate EMG and ECG features. The framework has been used to classify four groups: healthy subjects, patients with obstructive sleep apnea (OSA), patients with restless leg syndrome (RLS) and patients with both OSA and RLS. The proposed deep learning framework produced a mean accuracy of 72% and weighted F1 score of 0.57 across subjects for our formulated four-class problem.


title: Unbounded Algebraic Derivators

abstract: We show that the unbounded derived category of a Grothendieck category with enough projective objects is the base category of a derivator whose category of diagrams is the full 2-category of small categories. With this structure, we give a description of the localization functor associated to a specialization closed subset of the spectrum of a commutative noetherian ring. In addition, using the derivator of modules, we prove some basic theorems of group cohomology for complexes of representations over an arbitrary base ring.


title: Prediction of the excitation energies of the 2$^+_1$ states for superheavy nuclei based on the microscopically derived Grodzins relation

abstract: As the result of synthesis of nuclei with large proton numbers a new region of investigations of the structure of nuclei has been discovered. Due to the recent significant increase in the yield of superheavy nuclei their gamma-spectroscopic studies became possible. The purpose of paper is to predict the excitation energies of the $2^+_1$ states of nuclei with Z$\ge 100$ using the microscopic variant of the Grodzins relation derived based on the geometrical collective model. The excitation energies of the $2^+_1$ states of the even-even nuclei from $^{256}$Fm to $^{296}_{120}$X which differ from each other in the number of $\alpha$-particles are predicted. It is shown that at the beginning of the chain of the studied nuclei the excitation energies of the $2^+_1$ states don't exceed 100 keV. Then $E(2^+_1)$ sharply increases with $A$ and reaches maximum value of $400-500$ keV in $^{284}$Fl or $^{292}$Og.


title: Influence of the Magnetic Sub-Lattices in the Double Perovskite Compound LaCaNiReO$_6$

abstract: The magnetism of double perovskites is a complex phenomenon, determined from intra- or interatomic magnetic moment interactions, and strongly influenced by geometry. We take advantage of the complementary length and time scales of the muon spin rotation, relaxation and resonance ($\mu^+$SR) microscopic technique and bulk AC/DC magnetic susceptibility measurements to study the magnetic phases of the LaCaNiReO$_6$ double perovskite. As a result we are able to discern and report a newly found dynamic phase transition and the formation of magnetic domains below and above the known magnetic transition of this compound at T$_N$ = 103 K. $\mu^+$SR, serving as a local probe at crystallographic interstitial sites, reveals a transition from a metastable ferrimagnetic ordering below T = 103 K to a stable one below T = 30 K. The fast and slow collective dynamic state of this system are investigated. Between 103 K < T < 230 K, the following two magnetic environments appear, a dense spin region and a static-dilute spin region. The paramagnetic state is obtained only above T > 270 K. An evolution of the interaction between Ni and Re magnetic sublattices in this geometrically frustrated fcc perovskite structure, is revealed as a function of temperature and magnetic field, through the critical behaviour and thermal evolution of microscopic and macroscopic physical quantities.


title: The many faces (phases) of strong correlations

abstract: There has been considerable recent progress in discovering and understanding quantum phases and fluctuations produced by strong correlations. Heavy fermion systems are an ideal platform for systematic studies because low and competing energy scales make them highly tunable. As such the phases (faces) of strong correlations transform continuously into one another.


title: Exact-WKB analysis for SUSY and quantum deformed potentials: Quantum mechanics with Grassmann fields and Wess-Zumino terms

abstract: Quantum deformed potentials arise naturally in quantum mechanical systems of one bosonic coordinate coupled to $N_f$ Grassmann valued fermionic coordinates, or to a topological Wess-Zumino term. These systems decompose into sectors with a classical potential plus a quantum deformation. Using exact WKB, we derive exact quantization condition and its median resummation. The solution of median resummed form gives physical Borel-Ecalle resummed results, as we show explicitly in quantum deformed double- and triple- well potentials. Despite the fact that instantons are finite action, for generic quantum deformation, they do not contribute to the energy spectrum at leading order in semi-classics. For certain quantized quantum deformations, where the alignment of levels to all order in perturbation theory occurs, instantons contribute to the spectrum. If deformation parameter is not properly quantized, their effect disappears, but higher order effects in semi-classics survive. In this sense, we classify saddle contributions as fading and robust. Finally, for quantum deformed triple-well potential, we demonstrate the P-NP relation, by computing period integrals and Mellin transform.


title: The Decidability and Complexity of Interleaved Bidirected Dyck Reachability

abstract: Dyck reachability is the standard formulation of a large domain of static analyses, as it achieves the sweet spot between precision and efficiency, and has thus been studied extensively. Interleaved Dyck reachability (denoted $D_k\odot D_k$) uses two Dyck languages for increased precision (e.g., context and field sensitivity) but is well-known to be undecidable. As many static analyses yield a certain type of bidirected graphs, they give rise to interleaved bidirected Dyck reachability problems. Although these problems have seen numerous applications, their decidability and complexity has largely remained open. In a recent work, Li et al. made the first steps in this direction, showing that (i) $D_1\odot D_1$ reachability (i.e., when both Dyck languages are over a single parenthesis and act as counters) is computable in $O(n^7)$ time, while (ii) $D_k\odot D_k$ reachability is NP-hard. In this work we address the decidability and complexity of all variants of interleaved bidirected Dyck reachability. First, we show that $D_1\odot D_1$ reachability can be computed in $O(n^3\cdot \alpha(n))$ time, significantly improving over the existing $O(n^7)$ bound. Second, we show that $D_k\odot D_1$ reachability (i.e., when one language acts as a counter) is decidable, in contrast to the non-bidirected case where decidability is open. We further consider $D_k\odot D_1$ reachability where the counter remains linearly bounded. Our third result shows that this bounded variant can be solved in $O(n^2\cdot \alpha(n))$ time, while our fourth result shows that the problem has a (conditional) quadratic lower bound, and thus our upper bound is essentially optimal. Fifth, we show that full $D_k\odot D_k$ reachability is undecidable. This improves the recent NP-hardness lower-bound, and shows that the problem is equivalent to the non-bidirected case.


title: Numerical modeling of anisotropic ferroelectric materials with hybridizable discontinuous Galerkin methods

abstract: We investigate a gradient flow structure of the Ginzburg--Landau--Devonshire (GLD) model for anisotropic ferroelectric materials by reconstructing its energy form. We show that the modified energy form admits at least one minimizer. Under some regularity assumptions for the electric charge distribution and the initial polarization field, we prove that the $L^2$ gradient flow structure has a unique solution. To simulate the GLD model numerically, we propose an energy-stable semi-implicit time-stepping scheme and a hybridizable discontinuous Galerkin method for space discretization. Some numerical tests are provided to verify the stability and convergence of the proposed numerical scheme as well as some properties of ferroelectric materials.


title: Generalizing the Quantum Information Model for Dynamic Diffraction

abstract: The development of novel neutron optics devices that rely on perfect crystals and nano-scale features are ushering a new generation of neutron science experiments, from fundamental physics to material characterization of emerging quantum materials. However, the standard theory of dynamical diffraction (DD) that analyzes neutron propagation through perfect crystals does not consider complex geometries, deformations, and/or imperfections which are now becoming a relevant systematic effect in high precision interferometric experiments. In this work, we expand upon a quantum information (QI) model of DD that is based on propagating a particle through a lattice of unitary quantum gates. We show that the model output is mathematically equivalent to the spherical wave solution of the Takagi-Taupin equations when in the appropriate limit, and that the model can be extended to the Bragg as well as the Laue-Bragg geometry where it is consistent with experimental data. The presented results demonstrate the universality of the QI model and its potential for modeling scenarios that are beyond the scope of the standard theory of DD.


title: Quantifying Microstructural Evolution via Time-Dependent Reduced-Dimension Metrics Based on Hierarchical $n$-Point Polytope Functions

abstract: We devise reduced-dimension metrics for effectively measuring the distance between two points (i.e., microstructures) in the microstructure space and quantifying the pathway associated with microstructural evolution, based on a recently introduced set of hierarchical $n$-point polytope functions $P_n$. The $P_n$ functions provide the probability of finding particular $n$-point configurations associated with regular $n$-polytopes in the material system, and a special sub-set of the standard $n$-point correlation functions $S_n$ that effectively decomposes the structural features in the systems into regular polyhedral basis with different symmetry. The $n$-th order metric $\Omega_n$ is defined as the $\mathbb{L}_1$ norm associated with the $P_n$ functions of two distinct microstructures. By choosing a reference initial state (i.e., a microstructure associated with $t_0 = 0$), the $\Omega_n(t)$ set quantifies the evolution of distinct polyhedral symmetries and can in principle capture emerging polyhedral symmetries that are not apparent in the initial state. To demonstrate their utility, we apply the $\Omega_n$ metrics to a 2D binary system undergoing spinodal decomposition to extract the phase separation dynamics via the temporal scaling behavior of the corresponding $\Omega_n(t)$, which reveals mechanisms governing the evolution. Moreover, we employ $\Omega_n(t)$ to analyze pattern evolution during vapor-deposition of phase-separating alloy films with different surface contact angles, which exhibit rich evolution dynamics including both unstable and oscillating patterns. The $\Omega_n$ metrics have potential applications in establishing quantitative processing-structure-property relationships, as well as real-time processing control and optimization of complex heterogeneous material systems.