title: New Anisotropic Sudden Singularities and Dimensional Reduction

abstract: We demonstrate the existence of sudden, finite-time, singularities, with constant scale factor, expansion rate and density, in expanding Bianchi type IX universes with free anisotropic pressures. A new type of non-simultaneous anisotropic sudden singularity arises because of the divergences of the pressures, which may be of barrel or pancake type. The effect of one or more directions of expansion hitting a sudden singularity is tantamount to dimensional reductions as the non-singular directions continue expanding and can see the sudden singularity in their past.


title: Factors of HOMFLY polynomials

abstract: We study factorizations of HOMFLY polynomials of certain knots and oriented links. We begin with a computer analysis of knots with at most 12 crossings, finding 17 non-trivial factorizations. Next, we give an irreducibility criterion for HOMFLY polynomials of oriented links associated to 2-connected plane graphs.


title: New Approximations and Hardness Results for Submodular Partitioning Problems

abstract: We consider the following class of submodular k-multiway partitioning problems: (Sub-$k$-MP) $\min \sum_{i=1}^k f(S_i): S_1 \uplus S_2 \uplus \cdots \uplus S_k = V \mbox{ and } S_i \neq \emptyset \mbox{ for all }i\in [k]$. Here $f$ is a non-negative submodular function, and $\uplus$ denotes the union of disjoint sets. Hence the goal is to partition $V$ into $k$ non-empty sets $S_1,S_2,\ldots,S_k$ such that $\sum_{i=1}^k f(S_i)$ is minimized. These problems were introduced by Zhao et al. partly motivated by applications to network reliability analysis, VLSI design, hypergraph cut, and other partitioning problems. In this work we revisit this class of problems and shed some light onto their hardness of approximation in the value oracle model. We provide new unconditional hardness results for Sub-$k$-MP in the special settings where the function $f$ is either monotone or symmetric. For symmetric functions we show that given any $\epsilon > 0$, any algorithm achieving a $(2 - \epsilon)$-approximation requires exponentially many queries in the value oracle model. For monotone objectives we show that given any $\epsilon > 0$, any algorithm achieving a $(4/3 - \epsilon)$-approximation requires exponentially many queries in the value oracle model. We then extend Sub-$k$-MP to a larger class of partitioning problems, where the functions $f_i(S_i)$ can be different, and there is a more general partitioning constraint $ S_1 \uplus S_2 \uplus \cdots \uplus S_k \in \mathcal{F}$ for some family $\mathcal{F} \subseteq 2^V$ of feasible sets. We provide a black box reduction that allows us to leverage several existing results from the literature; leading to new approximations for this class of problems.


title: Towards Entrepreneurial Ecosystem Indicators : Speed and Acceleration

abstract: We suggest the use of indicators to analyze entrepreneurial ecosystems, in a way similar to ecological indicators: simple, measurable, and actionable characteristics, used to convey relevant information to stakeholders and policymakers. We define 3 possible such indicators: Fundraising Speed, Acceleration and nth-year speed, all related to the ability of startups to develop more or less rapidly in a given ecosystem. Results based on these 3 indicators for 6 prominent ecosystems (Berlin, Israel, London, New York, Paris, Silicon Valley) exhibit markedly different situations and trajectories. Altogether, they contribute to confirm that such indicators can help shed new and interesting light on entrepreneurial ecosystems, to the benefit of potentially more grounded policy decisions, and all the more so in otherwise blurred and somewhat cacophonic environments.


title: Discovery of a near-infrared bar and a pseudobulge in the collisional ring galaxy Cartwheel

abstract: We report the discovery of a bar, a pseudobulge and unresolved point source in the archetype collisional ring galaxy Cartwheel using careful morphological analysis of a near-infrared (NIR) K$_s$ band image of excellent quality (seeing=0.42\arcsec) at the ESO archive. The bar is oval-shaped with a semi-major axis length of 3.23\arcsec\($\sim$2.09~kpc), with almost a flat light distribution along it. The bulge is almost round (ellipticity=0.21) with an effective radius of 1.62\arcsec\ ($\sim$1.05~kpc) and a Sersic index of 0.99, parameters typical of pseudobulges in late-type galaxies. The newly discovered bar is not recognisable as such in the optical images even with more than a factor of two higher spatial resolution of the Hubble Space Telescope, due to a combination of its red colour and the presence of dusty features. The observed bar and pseudobulge most likely belonged to the pre-collisional progenitor of the Cartwheel. The discovery of a bar in an archetype collisional ring galaxy Cartwheel is the first observational evidence to confirm the prediction that bars can survive a drop-through collision along with the morphological structures like a central bulge (pseudo).


title: An Error-limited NOMA-HARQ Approach using Short Packets

abstract: In this paper, we propose and analyze an error-constrained data transmission approach for uplink communication using non-orthogonal multiple access (NOMA). The results are presented for the cases with hybrid automatic repeat request (HARQ) protocols where, depending on the message decoding status in the previous rounds, the decoding scheme of the receiver as well as the data rates/transmission powers of the transmitters are adapted correspondingly. We perform finite block-length analysis of the proposed scheme and study the effect of the codewords length on the performance of NOMA-HARQ protocols. Particularly, considering repetition time diversity HARQ and codewords of finite length, we determine the expected error-limited achievable rate and the transmit power improvement of our proposed scheme, compared to the conventional NOMA-HARQ protocols. As we show, smart adaptation of the decoding scheme improves the spectral/energy efficiency of NOMA-HARQ protocols considerably with no additional receiver complexity.


title: Inferring median survival differences in general factorial designs via permutation tests

abstract: Factorial survival designs with right-censored observations are commonly inferred by Cox regression and explained by means of hazard ratios. However, in case of non-proportional hazards, their interpretation can become cumbersome; especially for clinicians. We therefore offer an alternative: median survival times are used to estimate treatment and interaction effects and null hypotheses are formulated in contrasts of their population versions. Permutation-based tests and confidence regions are proposed and shown to be asymptotically valid. Their type-1 error control and power behavior are investigated in extensive simulations, showing the new methods' wide applicability. The latter is complemented by an illustrative data analysis.


title: A Fast Finite Field Multiplier for SIKE

abstract: Various post-quantum cryptography algorithms have been recently proposed. Supersingluar isogeny Diffie-Hellman key exchange (SIKE) is one of the most promising candidates due to its small key size. However, the SIKE scheme requires numerous finite field multiplications for its isogeny computation, and hence suffers from slow encryption and decryption process. In this paper, we propose a fast finite field multiplier design that performs multiplications in GF(p) with high throughput and low latency. The design accelerates the computation by adopting deep pipelining, and achieves high hardware utilization through data interleaving. The proposed finite field multiplier demonstrates 4.48 times higher throughput than prior work based on the identical fast multiplication algorithm and 1.43 times higher throughput than the state-of-the-art fast finite field multiplier design aimed at SIKE.


title: A random walk on Area Restricted Search

abstract: These notes from a graduate class at the Unuversidad Autonoma de Madrid analyze a search behavior known as Area Resticted Search (ARS), widespread in the animal kingdom, and optimal when the resources that one is after are "patchy". In the first section we study the importance of the behavior in animal and its dependence on the dopamine as a indicator of reward. In the second section we put together a genetic algorithm to determine the optimality of ARS and its characteristics. Finally, we relate ARS to a type of random walks known as "Levy Walks", in which the probability of jumping at a distance d from the current location follows a power law distribution.


title: Deep Learning for Cornea Microscopy Blind Deblurring

abstract: The goal of this project is to build a deep-learning solution that deblurs cornea scans, used for medical examination. The spherical shape of the eye prevents ophtamologist from having completely sharp image. Provided with a stack of corneas from confocal images, our approach is to build a model that performs an upscaling of the images using an SR (Super Resolution) Network.


title: Analyzing Effect of Repeated Reading on Oral Fluency and Narrative Production for Computer-Assisted Language Learning

abstract: Repeated reading (RR) helps learners, who have little to no experience with reading fluently to gain confidence, speed and process words automatically. The benefits of repeated readings include helping all learners with fact recall, aiding identification of learners' main ideas and vocabulary, increasing comprehension, leading to faster reading as well as increasing word recognition accuracy, and assisting struggling learners as they transition from word-by-word reading to more meaningful phrasing. Thus, RR ultimately helps in improvements of learners' oral fluency and narrative production. However, there are no open audio datasets available on oral responses of learners based on their RR practices. Therefore, in this paper, we present our dataset, discuss its properties, and propose a method to assess oral fluency and narrative production for learners of English using acoustic, prosodic, lexical and syntactical characteristics. The results show that a CALL system can be developed for assessing the improvements in learners' oral fluency and narrative production.


title: Perfusion Quantification from Endoscopic Videos: Learning to Read Tumor Signatures

abstract: Intra-operative identification of malignant versus benign or healthy tissue is a major challenge in fluorescence guided cancer surgery. We propose a perfusion quantification method for computer-aided interpretation of subtle differences in dynamic perfusion patterns which can be used to distinguish between normal tissue and benign or malignant tumors intra-operatively in real-time by using multispectral endoscopic videos. The method exploits the fact that vasculature arising from cancer angiogenesis gives tumors differing perfusion patterns from the surrounding tissue, and defines a signature of tumor which could be used to differentiate tumors from normal tissues. Experimental evaluation of our method on a cohort of colorectal cancer surgery endoscopic videos suggests that the proposed tumor signature is able to successfully discriminate between healthy, cancerous and benign tissue with 95% accuracy.


title: Radiative cooling rates, ion fractions, molecule abundances and line emissivities including self-shielding and both local and metagalactic radiation fields

abstract: We use the spectral synthesis code Cloudy to tabulate the properties of gas for an extensive range in redshift (z=0 to 9), temperature (log T [K] = 1 to 9.5), metallicity (log Z/$\mathrm{Z}_{\odot}$ = -4 to +0.5, Z=0), and density (log $n_{\mathrm{H}}$ [$\mathrm{cm}^{-3}$] = -8 to +6). This therefore includes gas with properties characteristic of the interstellar, circumgalactic and intergalactic media. The gas is exposed to a redshift-dependent UV/X-ray background, while for the self-shielded lower-temperature gas (i.e. ISM gas) an interstellar radiation field and cosmic rays are added. The radiation field is attenuated by a density- and temperature-dependent column of gas and dust. Motivated by the observed star formation law, this gas column density also determines the intensity of the interstellar radiation field and the cosmic ray density. The ionization balance, molecule fractions, cooling rates, line emissivities, and equilibrium temperatures are calculated self-consistently. We include dust, cosmic rays, and the interstellar radiation field step-by-step to study their relative impact. These publicly available tables are ideal for hydrodynamical simulations. They can be used stand alone or coupled to a non-equilibrium network for a subset of elements. The release includes a C routine to read in and interpolate the tables, as well as an easy to use python graphical user interface to explore the tables.


title: A room temperature optomechanical squeezer

abstract: One of the noise sources that currently limits gravitational wave (GW) detectors comes from the quantum nature of the light causing uncertain amplitude and phase. Phase uncertainty limits the precision of an interferometric measurement. This measurement is also subject to quantum back-action, caused by the radiation pressure force fluctuations produced by the amplitude uncertainty (QRPN). In order to lower this quantum noise, GW detectors plan to use squeezed light injection. In this thesis, I focus on using radiation-pressure-mediated optomechanical (OM) interaction to generate squeezed light. Creating squeezed states by using OM interaction enables wavelength-independent squeezed light sources that may also be more compact and robust than traditionally used non-linear crystals. We analyze the system with realistic imperfections (losses & classical noise), and use the concepts to design an experiment to obtain the most possible squeezing in a broad audio-frequency band at room temperature. This involves an optimization for the optical properties of the cavity and the mechanical properties of the oscillator. We then show its experimental implementation, and subsequent observation of QRPN as well as OM squeezing. These are the first ever direct observations of a room temperature oscillator's motion being overwhelmed by vacuum fluctuations. This is shown in the low frequency band, which is relevant to GW detectors, but poses its own technical challenges, and hence has not been done before. Being in the back-action dominated regime along with optimized optical properties has also enabled us to observe OM squeezing. That is the first direct observation of quantum noise suppression in a room temperature OM system. It is also the first direct evidence of quantum correlations in the audio frequency band, in a broad band at non-resonant frequencies.


title: Jet quenching parameter from a soft wall AdS/QCD model

abstract: We study the effect of chemical potential and nonconformality on the jet quenching parameter in a holographic QCD model with conformal invariance broken by a background dilaton. It turns out that the presence of chemical potential and nonconformality both increase the jet quenching parameter thus enhancing the energy loss, consistently with the findings of the drag force.


title: High-Dimensional Quadratic Discriminant Analysis under Spiked Covariance Model

abstract: Quadratic discriminant analysis (QDA) is a widely used classification technique that generalizes the linear discriminant analysis (LDA) classifier to the case of distinct covariance matrices among classes. For the QDA classifier to yield high classification performance, an accurate estimation of the covariance matrices is required. Such a task becomes all the more challenging in high dimensional settings, wherein the number of observations is comparable with the feature dimension. A popular way to enhance the performance of QDA classifier under these circumstances is to regularize the covariance matrix, giving the name regularized QDA (R-QDA) to the corresponding classifier. In this work, we consider the case in which the population covariance matrix has a spiked covariance structure, a model that is often assumed in several applications. Building on the classical QDA, we propose a novel quadratic classification technique, the parameters of which are chosen such that the fisher-discriminant ratio is maximized. Numerical simulations show that the proposed classifier not only outperforms the classical R-QDA for both synthetic and real data but also requires lower computational complexity, making it suitable to high dimensional settings.


title: Optimal control, contact dynamics and Herglotz variational problem

abstract: In this paper we combine two main topics in mechanics and optimal control theory: contact Hamiltonian systems and Pontryagin Maximum Principle. As an important result, among others, we develop a contact Pontryagin Maximum Principle that permits to deal with optimal control problems with dissipation. We also consider the Herglotz optimal control problem, which is simultaneously a generalization of the Herglotz variational principle and an optimal control problem. An application to the study of a thermodynamic system is provided.


title: Predicting hydration free energies of the FreeSolv database of druglike molecules with molecular density functional theory

abstract: We assess the performance of molecular densityfunctional theory (MDFT) to predict hydration freeenergies of the small drug-like molecules benchmark,FreeSolv. MDFT in the hyper-netted chain approx-imation (HNC) coupled with a pressure correctionpredicts experimental hydration free energies of theFreeSolv database within 1 kcal/mol with an averagecomputation time of two cpu.min per molecule. Thisis the same accuracy as for simulation based free en-ergy calculations that typically require hundreds ofcpu.h or tens of gpu.h per molecule.


title: Global polarization in heavy-ion collisions based on axial vortical effect

abstract: Global polarization of $\Lambda$ and $\bar{\Lambda}$ is calculated based on the axial vortical effect (AVE). Simulations are performed within the model of the three-fluid dynamics. Equations of state with the deconfinement transition result in a good agreement with STAR data for both $\Lambda$ and $\bar{\Lambda}$ polarization, in particular, with the $\Lambda$-$\bar{\Lambda}$ splitting. Suppression of the gravitational-anomaly contribution required for the data reproduction is in agreement with predictions of the QCD lattice simulations. Predictions for the global polarization in forthcoming experiments at lower collision energies are made. These forthcoming data will provide a critical test for the AVE and thermodynamic mechanisms of the polarization.


title: Differentially Private Health Tokens for Estimating COVID-19 Risk

abstract: In the fight against Covid-19, many governments and businesses are in the process of evaluating, trialling and even implementing so-called immunity passports. Also known as antibody or health certificates, there is a clear demand for any technology that could allow people to return to work and other crowded places without placing others at risk. One of the major criticisms of such systems is that they could be misused to unfairly discriminate against those without immunity, allowing the formation of an `immuno-privileged' class of people. In this work we are motivated to explore an alternative technical solution that is non-discriminatory by design. In particular we propose health tokens -- randomised health certificates which, using methods from differential privacy, allow individual test results to be randomised whilst still allowing useful aggregate risk estimates to be calculated. We show that health tokens could mitigate immunity-based discrimination whilst still presenting a viable mechanism for estimating the collective transmission risk posed by small groups of users. We evaluate the viability of our approach in the context of identity-free and identity-binding use cases and then consider a number of possible attacks. Our experimental results show that for groups of size 500 or more, the error associated with our method can be as low as 0.03 on average and thus the aggregated results can be useful in a number of identity-free contexts. Finally, we present the results of our open-source prototype which demonstrates the practicality of our solution.


title: Time-varying Graph Representation Learning via Higher-Order Skip-Gram with Negative Sampling

abstract: Representation learning models for graphs are a successful family of techniques that project nodes into feature spaces that can be exploited by other machine learning algorithms. Since many real-world networks are inherently dynamic, with interactions among nodes changing over time, these techniques can be defined both for static and for time-varying graphs. Here, we build upon the fact that the skip-gram embedding approach implicitly performs a matrix factorization, and we extend it to perform implicit tensor factorization on different tensor representations of time-varying graphs. We show that higher-order skip-gram with negative sampling (HOSGNS) is able to disentangle the role of nodes and time, with a small fraction of the number of parameters needed by other approaches. We empirically evaluate our approach using time-resolved face-to-face proximity data, showing that the learned time-varying graph representations outperform state-of-the-art methods when used to solve downstream tasks such as network reconstruction, and to predict the outcome of dynamical processes such as disease spreading. The source code and data are publicly available at https://github.com/simonepiaggesi/hosgns.


title: A Theoretical Framework for Target Propagation

abstract: The success of deep learning, a brain-inspired form of AI, has sparked interest in understanding how the brain could similarly learn across multiple layers of neurons. However, the majority of biologically-plausible learning algorithms have not yet reached the performance of backpropagation (BP), nor are they built on strong theoretical foundations. Here, we analyze target propagation (TP), a popular but not yet fully understood alternative to BP, from the standpoint of mathematical optimization. Our theory shows that TP is closely related to Gauss-Newton optimization and thus substantially differs from BP. Furthermore, our analysis reveals a fundamental limitation of difference target propagation (DTP), a well-known variant of TP, in the realistic scenario of non-invertible neural networks. We provide a first solution to this problem through a novel reconstruction loss that improves feedback weight training, while simultaneously introducing architectural flexibility by allowing for direct feedback connections from the output to each hidden layer. Our theory is corroborated by experimental results that show significant improvements in performance and in the alignment of forward weight updates with loss gradients, compared to DTP.


title: Dynamical Systems, Representation of Particle Beams

abstract: An overview of dynamical systems in accelerator physics is presented with a suggestion of a few issues to be addressed. Also mentioned are a few possible developments in the future. Technical details supporting the views are not presented.


title: Data characterization in dynamical inverse problem for the 1d wave equation with matrix potential

abstract: The dynamical system under consideration is \begin{align*} & u_{tt}-u_{xx}+Vu=0,\qquad x>0,\,\,\,t>0;\\ & u|_{t=0}=u_t|_{t=0}=0,\,\,x\geqslant 0;\quad u|_{x=0}=f,\,\,t\geqslant 0, \end{align*} where $V=V(x)$ is a matrix-valued function ({\it potential}); $f=f(t)$ is an $\mathbb R^N$-valued function of time ({\it boundary control}); $u=u^f(x,t)$ is a {\it trajectory} (an $\mathbb R^N$-valued function of $x$ and $t$). The input/output map of the system is a {\it response operator} $R:f\mapsto u^f_x(0,\cdot),\,\,\,t\geqslant0$. The {\it inverse problem} is to determine $V$ from given $R$. To characterize its data is to provide the necessary and sufficient conditions on $R$ that ensure its solvability. The procedure that solves this problem has long been known and the characterization has been announced (Avdonin and Belishev, 1996). However, the proof was not provided and, moreover, it turned out that the formulation must be corrected. Our paper fills this gap.


title: Wavevector analysis of plasmon-assisted distributed nonlinear photoluminescence along Au nanowire antennas

abstract: We report a quantitative analysis of the wavevector diagram emitted by nonlinear photoluminescence generated by a tightly focused pulsed laser beam and distributed along Au nanowire via the mediation of surface plasmon polaritions. The nonlinear photoluminescence is locally excited at key locations along the nanowire in order to understand the different contributions constituting the emission pattern measured in a conjugate Fourier plane of the microscope. Polarization-resolved measurements reveal that the nanowire preferentially emits nonlinear photoluminescence polarized transverse to the long axis at close to the detection limit wavevectors with a small azimuthal spread in comparison to the signal polarized along the long axis. We utilize finite element method to simulate the observed directional scattering by using localized incoherent sources placed on the nanowire. Simulation results faithfully mimic the directional emission of the nonlinear signal emitted by the different portions of the nanowire.


title: Thermal radiation and inclusive production in the running coupling $k_T$ -- factorization approach

abstract: The characteristics of the thermal radiation are investigated using a two - component model, with the hard component being described by the Color Glass Condensate formalism. The inclusive transverse momentum spectra of charged hadrons produced in proton - proton and proton - nucleus collisions at LHC energies and large - $p_T$ are estimated using the running coupling $k_T$ - factorization formula and the solution of the Balitsky - Kovchegov equation. Our results indicate that the thermal term is necessary to describe the experimental data and that the effective thermal temperature has an energy dependence similar to the saturation scale. We demonstrate that the enhancement of the thermal temperature in $pPb$ collisions is consistent with that predicted by the saturation scale.


title: Multi-wavelength Photometry and Progenitor Analysis of the Nova V906 Car

abstract: We present optical and infrared photometry of the classical nova V906 Car, also known as Nova Car 2018 and ASASSN-18fv, discovered by ASASS-SN survey on 16.32 March 2018 UT (MJD 58193.0). The nova reached its maximum on MJD 58222.56 at $V_{\rm{max}} = 5.84 \pm 0.09$ mag and had decline times of $t_{2,V} = 26.2 $ d and $t_{3,V} = 33.0 $ d. The data from Evryscope shows that the nova had already brightened to $g'\simeq 13$\,mag five days before discovery, as compared to its quiescent magnitude of $g=$20.13$\pm$0.03. The extinction towards the nova, as derived from high resolution spectroscopy, shows an estimate consistent with foreground extinction to the Carina Nebula of $A_V = 1.11_{-0.39}^{+0.54}$. The light curve resembles a rare C (cusp) class nova with a steep decline slope of $\alpha=-3.94$ post cusp flare. From the lightcurve decline rate, we estimate the mass of white dwarf to be $M_{WD}$ = $ < 0.8$M\textsubscript{\(\odot\)}, consistent with $M_{WD}=0.71^{+0.23}_{-0.19}$ derived from modelling the accretion disk of the system in quiescence. The donor star is likely a K-M dwarf of 0.23-0.43\,\Msun, which is being heated by its companion.


title: Secure quantum key distribution with a subset of malicious devices

abstract: The malicious manipulation of quantum key distribution (QKD) hardware is a serious threat to its security, as, typically, neither end users nor QKD manufacturers can validate the integrity of every component of their QKD system in practice. One possible approach to re-establish the security of QKD is to use a redundant number of devices. Following this idea, we introduce an efficient distributed QKD post-processing protocol and prove its security in a variety of corruption models of the possibly malicious devices. We find that, compared to the most conservative model of active and collaborative corrupted devices, natural assumptions lead to a significant enhancement of the secret key rate and considerably simpler QKD setups. Furthermore, we show that, for most practical situations, the resulting finite-size secret key rate is similar to that of the standard scenario assuming trusted devices.


title: Classical and quantum stochastic thermodynamics

abstract: The stochastic thermodynamics provides a framework for the description of systems that are out of thermodynamic equilibrium. It is based on the assumption that the elementary constituents are acted by random forces that generate a stochastic dynamics, which is here represented by a Fokker-Planck-Kramers equation. We emphasize the role of the irreversible probability current, the vanishing of which characterizes the thermodynamic equilibrium and yields a special relation between fluctuation and dissipation. The connection to thermodynamics is obtained by the definition of the energy function and the entropy as well as the rate at which entropy is generated. The extension to quantum systems is provided by a quantum evolution equation which is a canonical quantization of the Fokker-Planck-Kramers equation. An example of an irreversible systems is presented which shows a nonequilibrium stationary state with an unceasing production of entropy. A relationship between the fluxes and the path integral is also presented.


title: In-plane magnetic domains and N\'eel-like domain walls in thin flakes of the room temperature CrTe$_2$ van der Waals ferromagnet

abstract: The recent discovery of magnetic van der Waals materials has triggered a wealth of investigations in materials science, and now offers genuinely new prospects for both fundamental and applied research. Although the catalogue of van der Waals ferromagnets is rapidly expanding, most of them have a Curie temperature below 300 K, a notable disadvantage for potential applications. Combining element-selective x-ray magnetic imaging and magnetic force microscopy, we resolve at room temperature the magnetic domains and domains walls in micron-sized flakes of the CrTe$_2$ van der Waals ferromagnet. Flux-closure magnetic patterns suggesting in-plane six-fold symmetry are observed. Upon annealing the material above its Curie point (315 K), the magnetic domains disappear. By cooling back down the sample, a different magnetic domain distribution is obtained, indicating material stability and lack of magnetic memory upon thermal cycling. The domain walls presumably have N\'eel texture, are preferentially oriented along directions separated by 120 degrees, and have a width of several tens of nanometers. Besides microscopic mapping of magnetic domains and domain walls, the coercivity of the material is found to be of a few mT only, showing that the CrTe$_2$ compound is magnetically soft. The coercivity is found to increase as the volume of the material decreases.


title: Linear Optics Computations

abstract: In this Chapter we briefly recall and summarize the main linear optics concepts of the accelerators beam dynamics theory. In doing so we put our emphasis on the related computational aspects: the reader will be provided with the basic elements to write a linear optics code. To this aim, we complement the text with few practical examples and code listings.


title: Solving frustrated Ising models using tensor networks

abstract: Motivated by the recent success of tensor networks to calculate the residual entropy of spin ice and kagome Ising models, we develop a general framework to study frustrated Ising models in terms of infinite tensor networks %, i.e. tensor networks that can be contracted using standard algorithms for infinite systems. This is achieved by reformulating the problem as local rules for configurations on overlapping clusters chosen in such a way that they relieve the frustration, i.e. that the energy can be minimized independently on each cluster. We show that optimizing the choice of clusters, including the weight on shared bonds, is crucial for the contractibility of the tensor networks, and we derive some basic rules and a linear program to implement them. We illustrate the power of the method by computing the residual entropy of a frustrated Ising spin system on the kagome lattice with next-next-nearest neighbour interactions, vastly outperforming Monte Carlo methods in speed and accuracy. The extension to finite-temperature is briefly discussed.


title: Hecke eigensystems of automorphic forms (mod $p$) of Hodge type and algebraic modular forms

abstract: We show that the systems of prime-to-$p$ Hecke eigenvalues arising from automorphic forms$\pmod p$ for a good prime $p$ associated to an algebraic group $G/\mathbb Q$ of Hodge type are the same as those arising from algebraic modular forms$\pmod p$ associated to an inner form of $G$.


title: Spatio-temporal Inversion using the Selection Kalman Model

abstract: Data assimilation in models representing spatio-temporal phenomena poses a challenge, particularly if the spatial histogram of the variable appears with multiple modes. The traditional Kalman model is based on a Gaussian initial distribution and Gauss-linear dynamic and observation models. This model is contained in the class of Gaussian distribution and is therefore analytically tractable. It is however unsuitable for representing multimodality. We define the selection Kalman model that is based on a selection-Gaussian initial distribution and Gauss-linear dynamic and observation models. The selection-Gaussian distribution can be seen as a generalization of the Gaussian distribution and may represent multimodality, skewness and peakedness. This selection Kalman model is contained in the class of selection-Gaussian distributions and therefore it is analytically tractable. An efficient recursive algorithm for assessing the selection Kalman model is specified. The synthetic case study of spatio-temporal inversion of an initial state, inspired by pollution monitoring, containing an extreme event suggests that the use of the selection Kalman model offers significant improvements compared to the traditional Kalman model when reconstructing discontinuous initial states.


title: Thickness gradient promotes the performance of Si-based anode material for lithium-ion battery

abstract: The large volume change of the silicon (Si) during lithiation and delithiation process has long been a problem impeding its application as one of the most promising anode materials for LIBs. In this paper, we proposed a conceptually new idea to address this problem simply by varying the thickness of the electrode material film. The resulting thickness-gradient electrode exhibits considerable enhancement in the electrochemical performances including capacity, capacity retention, Coulombic efficiency, and rate capability in comparison to the traditional counterparts with uniform thickness. Such enhancement in the electrochemical performance can be attributed to the lessening of the stress concentration on the interface between the electrode film and the current collector upon the volume change of Si taking place in the lithiation and delithiation process. To make the best use of this strategy, optimal design of the gradient thickness is proposed based on the theory of stress homogenization, followed by the experimental verification. The results of this paper provide a facile, cost-effective, and scalable way for enhancing the performance of Si-based anodes for LIBs. This strategy can be further extended to the other anode materials suffering from the similar lithiation-induced volume change problem.


title: Collaborative Boundary-aware Context Encoding Networks for Error Map Prediction

abstract: Medical image segmentation is usually regarded as one of the most important intermediate steps in clinical situations and medical imaging research. Thus, accurately assessing the segmentation quality of the automatically generated predictions is essential for guaranteeing the reliability of the results of the computer-assisted diagnosis (CAD). Many researchers apply neural networks to train segmentation quality regression models to estimate the segmentation quality of a new data cohort without labeled ground truth. Recently, a novel idea is proposed that transforming the segmentation quality assessment (SQA) problem intothe pixel-wise error map prediction task in the form of segmentation. However, the simple application of vanilla segmentation structures in medical image fails to detect some small and thin error regions of the auto-generated masks with complex anatomical structures. In this paper, we propose collaborative boundaryaware context encoding networks called AEP-Net for error prediction task. Specifically, we propose a collaborative feature transformation branch for better feature fusion between images and masks, and precise localization of error regions. Further, we propose a context encoding module to utilize the global predictor from the error map to enhance the feature representation and regularize the networks. We perform experiments on IBSR v2.0 dataset and ACDC dataset. The AEP-Net achieves an average DSC of 0.8358, 0.8164 for error prediction task,and shows a high Pearson correlation coefficient of 0.9873 between the actual segmentation accuracy and the predicted accuracy inferred from the predicted error map on IBSR v2.0 dataset, which verifies the efficacy of our AEP-Net.


title: Fast General Distributed Transactions with Opacity using Global Time

abstract: Transactions can simplify distributed applications by hiding data distribution, concurrency, and failures from the application developer. Ideally the developer would see the abstraction of a single large machine that runs transactions sequentially and never fails. This requires the transactional subsystem to provide opacity (strict serializability for both committed and aborted transactions), as well as transparent fault tolerance with high availability. As even the best abstractions are unlikely to be used if they perform poorly, the system must also provide high performance. Existing distributed transactional designs either weaken this abstraction or are not designed for the best performance within a data center. This paper extends the design of FaRM - which provides strict serializability only for committed transactions - to provide opacity while maintaining FaRM's high throughput, low latency, and high availability within a modern data center. It uses timestamp ordering based on real time with clocks synchronized to within tens of microseconds across a cluster, and a failover protocol to ensure correctness across clock master failures. FaRM with opacity can commit 5.4 million neworder transactions per second when running the TPC-C transaction mix on 90 machines with 3-way replication.


title: Epoch-evolving Gaussian Process Guided Learning

abstract: In this paper, we propose a novel learning scheme called epoch-evolving Gaussian Process Guided Learning (GPGL), which aims at characterizing the correlation information between the batch-level distribution and the global data distribution. Such correlation information is encoded as context labels and needs renewal every epoch. With the guidance of the context label and ground truth label, GPGL scheme provides a more efficient optimization through updating the model parameters with a triangle consistency loss. Furthermore, our GPGL scheme can be further generalized and naturally applied to the current deep models, outperforming the existing batch-based state-of-the-art models on mainstream datasets (CIFAR-10, CIFAR-100, and Tiny-ImageNet) remarkably.


title: Audeo: Audio Generation for a Silent Performance Video

abstract: We present a novel system that gets as an input video frames of a musician playing the piano and generates the music for that video. Generation of music from visual cues is a challenging problem and it is not clear whether it is an attainable goal at all. Our main aim in this work is to explore the plausibility of such a transformation and to identify cues and components able to carry the association of sounds with visual events. To achieve the transformation we built a full pipeline named `\textit{Audeo}' containing three components. We first translate the video frames of the keyboard and the musician hand movements into raw mechanical musical symbolic representation Piano-Roll (Roll) for each video frame which represents the keys pressed at each time step. We then adapt the Roll to be amenable for audio synthesis by including temporal correlations. This step turns out to be critical for meaningful audio generation. As a last step, we implement Midi synthesizers to generate realistic music. \textit{Audeo} converts video to audio smoothly and clearly with only a few setup constraints. We evaluate \textit{Audeo} on `in the wild' piano performance videos and obtain that their generated music is of reasonable audio quality and can be successfully recognized with high precision by popular music identification software.


title: Extended Phase Space Thermodynamics for Black Holes in a Cavity

abstract: In this paper, we extend the phase space of black holes enclosed by a spherical cavity of radius $r_{B}$ to include $V\equiv4\pi r_{B}^{3}/3$ as a thermodynamic volume. The thermodynamic behavior of Schwarzschild and Reissner-Nordstrom (RN) black holes is then investigated in the extended phase space. In a canonical ensemble at constant pressure, we find that the aforementioned thermodynamic behavior is remarkably similar to that of the anti-de Sitter (AdS) counterparts with the cosmological constant being interpreted as a pressure. Specifically, a first-order Hawking-Page-like phase transition occurs for a Schwarzschild black hole in a cavity. The phase structure of a RN black hole in a cavity shows a strong resemblance to that of the van der Waals fluid. Our results may provide a new perspective for the extended thermodynamics of AdS black holes by analogy with black holes in a cavity.


title: Data-dependent Pruning to find the Winning Lottery Ticket

abstract: The Lottery Ticket Hypothesis postulates that a freshly initialized neural network contains a small subnetwork that can be trained in isolation to achieve similar performance as the full network. Our paper examines several alternatives to search for such subnetworks. We conclude that incorporating a data dependent component into the pruning criterion in the form of the gradient of the training loss -- as done in the SNIP method -- consistently improves the performance of existing pruning algorithms.


title: Digraphs whose m-step competition graphs are trees

abstract: In this paper, we completely characterize the digraphs of order $n$ whose $m$-step competition graphs are star graphs for positive integers $2\leq m < n$. This result in matrix version identifies the solution set to the matrix equation $X^m(X^T)^m= \Lambda_n+I_n$ for positive integers $2\leq m < n$ where $I_n$ is the identity matrix of order $n$ and $\Lambda_n$ is a $(0,1)$ Boolean matrix such that the first row and the first column consist of $1$'s except $(1,1)$-entry and the remaining entries are $0$, which is the adjacency matrix of a star graph of order $n$. We also derive meaningful properties of the digraphs whose $m$-step competition graphs are trees. In the process, we extend a result of Helleloid~[Connected triangle-free $m$-step competition graphs, Discrete Appl.\ Math.\ 145 (2005) 376--383] by showing that for all positive integers $m \geq 2$ and $n$, the connected triangle-free $m$-step competition graph on $n$ vertices is a tree.


title: HARMer: Cyber-attacks Automation and Evaluation

abstract: With the increasing growth of cyber-attack incidences, it is important to develop innovative and effective techniques to assess and defend networked systems against cyber attacks. One of the well-known techniques for this is performing penetration testing which is carried by a group of security professionals (i.e, red team). Penetration testing is also known to be effective to find existing and new vulnerabilities, however, the quality of security assessment can be depending on the quality of the red team members and their time and devotion to the penetration testing. In this paper, we propose a novel automation framework for cyber-attacks generation named `HARMer' to address the challenges with respect to manual attack execution by the red team. Our novel proposed framework, design, and implementation is based on a scalable graphical security model called Hierarchical Attack Representation Model (HARM). (1) We propose the requirements and the key phases for the automation framework. (2) We propose security metrics-based attack planning strategies along with their algorithms. (3) We conduct experiments in a real enterprise network and Amazon Web Services. The results show how the different phases of the framework interact to model the attackers' operations. This framework will allow security administrators to automatically assess the impact of various threats and attacks in an automated manner.


title: Dynamics of Causal Fermion Systems. Field Equations and Correction Terms for a New Unified Physical Theory

abstract: The theory of causal fermion systems is a new physical theory which aims to describe a fundamental level of physical reality. Its mathematical core is the causal action principle. In this thesis, we develop a formalism which connects the causal action principle to a suitable notion of fields on space-time. We derive field equations from the causal action principle and find that the dynamics induced by the field equations conserve a symplectic form which gives rise to an Hamiltonian time evolution if the causal fermion system admits a notion of time. In this way, we establish the dynamics of causal fermion systems. Remarkably, the causal action principle implies that there are correction terms to the field equations, which we subsequently derive and study. In particular, we prove that there is a stochastic and a non-linear correction term and investigate how they relate to the Hamiltonian time evolution. Furthermore, we give theorems which generalize the connection between symmetries and conservation laws in Noether's theorems to the theory of causal fermion systems. The appearance of the particular correction terms is reminiscent of dynamical collapse models in quantum theory.


title: Velocity-dependent phase shift in a light-pulse atom interferometer

abstract: Atom interferometry relies on the separation and recombination of atom wavepackets. When the two paths overlap perfectly at the end of the interferometer, the phase is insensitive to the atomic velocity distribution. Here, we show that, when the separation and recombination is performed using a Raman transition there is a displacement of the atomic wavepacket due to a phase shift during light pulses. Because of the variation of the laser intensity seen by the atoms, there is an imperfect cancellation of these displacements. The observation of a velocity-dependent phase shift on the interferometer is the signature of this effect, which has been modeled. Thanks to the signature we have identified, we are able to compensate for this effect by applying laser power ramps during the interferometer to mitigate intensity variations.


title: Controlling collective synchrony in oscillatory ensembles by precisely timed pulses

abstract: We present an efficient technique for control of synchrony in a globally coupled ensemble by pulsatile action. We assume that we can observe the collective oscillation and can stimulate all elements of the ensemble simultaneously. We pay special attention to the minimization of intervention into the system. The key idea is to stimulate only at the most sensitive phase. To find this phase we implement an adaptive feedback control. Estimating the instantaneous phase of the collective mode on the fly, we achieve efficient suppression using a few pulses per oscillatory cycle. We discuss the possible relevance of the results for neuroscience, namely for the development of advanced algorithms for deep brain stimulation, a medical technique used to treat Parkinson's disease.


title: Coherent excitation of the highly forbidden electric octupole transition in ${}^{172}$Yb$^+$

abstract: We report on the first coherent excitation of the highly forbidden $^2S_{1/2}\rightarrow{}^2F_{7/2}$ electric octupole (E3) transition in a single trapped ${}^{172}$Yb$^+$ ion, an isotope without nuclear spin. Using the transition in ${}^{171}$Yb$^+$ as a reference, we determine the transition frequency to be $642\,116\,784\,950\,887.6(2.4)\,$Hz. We map out the magnetic field environment using the forbidden $^2S_{1/2} \rightarrow{}^2D_{5/2}$ electric quadrupole (E2) transition and determine its frequency to be $729\,476\,867\,027\,206.8(4.4)\,$Hz. Our results are a factor of $1\times10^5$ ($3\times10^{5}$) more accurate for the E2 (E3) transition compared to previous measurements. The results open up the way to search for new physics via precise isotope shift measurements and improved tests of local Lorentz invariance using the metastable $^2F_{7/2}$ state of Yb$^+$.


title: Hierarchical stock assessment methods improve management performance in multi-species, data-limited fisheries

abstract: Management performance of five alternative stock assessment methods was evaluated by using them to set harvest levels targeting multi-species maximum yield in a multi-species flatfish fishery, including single-species and hierarchical multi-species models, and methods that pooled data across species and spatial strata, with catch outcomes of each method under three data scenarios compared to catch under an omniscient manager simulation. Operating models included technical interactions between species intended to produce choke effects often observed in output controlled multi-species fisheries. Hierarchical multi-species models outperformed all other methods under data-poor and data-moderate scenarios, and outperformed single-species models under the data-rich scenario. Hierarchical models were least sensitive to prior precision, sometimes improving in performance when prior precision was reduced. Choke effects were found to both positive and negative effects, sometimes leading to underfishing of non-choke species, but at other times preventing overfishing of non-choke species. We highlight the importance of including technical interactions in multi-species assessment models and management objectives, how choke species can indicate mismatches between management objectives and system dynamics, and recommend hierarchical multi-species models for multi-species fishery management systems.


title: An empirical study of availability and reliability properties of the Bitcoin Lightning Network

abstract: The Bitcoin Lightning network is a mechanism to enable fast and inexpensive off-chain Bitcoin transactions using peer-to-peer (P2P) channels between nodes that can also be composed into a routing path. Although the resulting possible channel graphs are well-studied, there is no empirical data on the network's reliability in terms of being able to successfully route payments at a given moment in time. In this paper we address this gap and investigate two forms of availability that are a necessary ingredient to achieve such reliability. We first study the Lightning network's ability to route payments of various sizes to nearly every participating node, over most available channels. We establish an inverse relationship between payment volume and success rate and show that at best only about a third of destination nodes can be successfully reached. The routing is hampered by a number of possible errors, both transient and permanent. We then study the availability of nodes in the network longitudinally and determine how long-lived they are. Churn in the network is actually low, and a considerable number of nodes are hosted on cloud providers. By testing node liveness, we find that the propagated network information is relatively often stale, however, both for IP addresses and Tor onion addresses. We provide recommendations how the Lightning network can be improved, including considerations which trade-offs between privacy and decentralization on the one hand and reliability on the other hand should at least be reconsidered by the community developing the Lightning network.


title: Reconstruction of signal amplitudes in the CMS electromagnetic calorimeter in the presence of overlapping proton-proton interactions

abstract: A template fitting technique for reconstructing the amplitude of signals produced by the lead tungstate crystals of the CMS electromagnetic calorimeter is described. This novel approach is designed to suppress the increased out-of-time pileup contribution to the signal following the reduction of the accelerator bunch spacing from 50 to 25 ns at the start of Run 2 of the LHC. Execution of the algorithm is sufficiently fast for it to be employed in the CMS high-level trigger. It is also used in the offline event reconstruction. Results obtained from simulations and from Run 2 collision data (2015-2018) demonstrate a substantial improvement in the energy resolution of the calorimeter over a range of energies extending from a few GeV to several tens of GeV.


title: Stability Enhanced Privacy and Applications in Private Stochastic Gradient Descent

abstract: Private machine learning involves addition of noise while training, resulting in lower accuracy. Intuitively, greater stability can imply greater privacy and improve this privacy-utility tradeoff. We study this role of stability in private empirical risk minimization, where differential privacy is achieved by output perturbation, and establish a corresponding theoretical result showing that for strongly-convex loss functions, an algorithm with uniform stability of $\beta$ implies a bound of $O(\sqrt{\beta})$ on the scale of noise required for differential privacy. The result applies to both explicit regularization and to implicitly stabilized ERM, such as adaptations of Stochastic Gradient Descent that are known to be stable. Thus, it generalizes recent results that improve privacy through modifications to SGD, and establishes stability as the unifying perspective. It implies new privacy guarantees for optimizations with uniform stability guarantees, where a corresponding differential privacy guarantee was previously not known. Experimental results validate the utility of stability enhanced privacy in several problems, including application of elastic nets and feature selection.


title: An updated version of "Leader-following consensus for linear multi-agent systems via asynchronous sampled-data control," IEEE Transactions on Automatic Control, DOI:10.1109/TAC.2019.2948256

abstract: In this article, we update the reference [14] in two aspects. First, we note that in order for the control law (12) in [14] to be equivalent to the control law (3) in [14], we need to assume that the samplings for all subsystems must be synchronous, i.e., we need to assume that $T_{i}=T$ for all $i=1,\cdots,N$. Second, we extend our results from periodic sampling to aperiodic sampling.


title: Dissociation in strong field: a quantum analysis of the relation between angular momentum and angular distribution of fragments

abstract: An ab initio simulation of strong-field photodissociation of diatomic molecules was developed, inspired by recent dissociation experiments of F2-. The transition between electronic states was modeled, including the laser pulse and transition dipole, and the angle between them. The initial conditions of the system were set to be thermal and to include different rovibrational states. Carefully designed absorbing boundary conditions were applied to describe the boundary conditions of the experiment. We studied the influence of field intensity on the direction of the outcoming fragments and laboratory-fixed axis, defined by the field polarization. At high intensities, the angular distribution became more peaked with a marginal influence on kinetic energy release.


title: SOAC: The Soft Option Actor-Critic Architecture

abstract: The option framework has shown great promise by automatically extracting temporally-extended sub-tasks from a long-horizon task. Methods have been proposed for concurrently learning low-level intra-option policies and high-level option selection policy. However, existing methods typically suffer from two major challenges: ineffective exploration and unstable updates. In this paper, we present a novel and stable off-policy approach that builds on the maximum entropy model to address these challenges. Our approach introduces an information-theoretical intrinsic reward for encouraging the identification of diverse and effective options. Meanwhile, we utilize a probability inference model to simplify the optimization problem as fitting optimal trajectories. Experimental results demonstrate that our approach significantly outperforms prior on-policy and off-policy methods in a range of Mujoco benchmark tasks while still providing benefits for transfer learning. In these tasks, our approach learns a diverse set of options, each of whose state-action space has strong coherence.


title: Finite-Sample Analysis of Proximal Gradient TD Algorithms

abstract: In this paper, we analyze the convergence rate of the gradient temporal difference learning (GTD) family of algorithms. Previous analyses of this class of algorithms use ODE techniques to prove asymptotic convergence, and to the best of our knowledge, no finite-sample analysis has been done. Moreover, there has been not much work on finite-sample analysis for convergent off-policy reinforcement learning algorithms. In this paper, we formulate GTD methods as stochastic gradient algorithms w.r.t.~a primal-dual saddle-point objective function, and then conduct a saddle-point error analysis to obtain finite-sample bounds on their performance. Two revised algorithms are also proposed, namely projected GTD2 and GTD2-MP, which offer improved convergence guarantees and acceleration, respectively. The results of our theoretical analysis show that the GTD family of algorithms are indeed comparable to the existing LSTD methods in off-policy learning scenarios.


title: Gaussian approximation potentials for body-centered-cubic transition metals

abstract: We develop a set of machine-learning interatomic potentials for elemental V, Nb, Mo, Ta, and W using the Gaussian approximation potential framework. The potentials show good accuracy and transferability for elastic, thermal, liquid, defect, and surface properties. All potentials are augmented with accurate repulsive potentials, making them applicable to radiation damage simulations involving high-energy collisions. We study melting and liquid properties in detail and use the potentials to provide melting curves up to 400 GPa for all five elements.


title: On the intermediate dimensions of Bedford-McMullen carpets

abstract: The intermediate dimensions of a set $\Lambda$, elsewhere denoted by $\dim_{\theta}\Lambda$, interpolates between its Hausdorff and box dimensions using the parameter $\theta\in[0,1]$. Determining a precise formula for $\dim_{\theta}\Lambda$ is particularly challenging when $\Lambda$ is a Bedford-McMullen carpet with distinct Hausdorff and box dimension. In this direction, answering a question of Fraser, we show that $\dim_{\theta}\Lambda$ is strictly less than the box dimension of $\Lambda$ for every $\theta<1$, moreover, the derivative of the upper bound is strictly positive at $\theta=1$. We also improve on the lower bound obtained by Falconer, Fraser and Kempton.


title: Fluid Flow on Vegetated Hillslope

abstract: In this paper, we present a deduction of swallow water equations in the presence of vegetation based on spatial averaging techniques starting from the general principles of conservation of mass and momentum. For this purpose, we worked in the hydrostatic approximation of the pressure field and we considered certain hypotheses of kinematic and topographical nature and assumptions on the structure of the vegetation. Some elements of differential geometry necessary to facilitate the reading of the paper can be found in the Appendix.


title: Magnetized cosmological model with variable deceleration parameter

abstract: In this paper, we have derived the field equations in an extended theory of gravity in an anisotropic space time background and in the presence of magnetic field. The physical and geometrical parameters of the models are determined with respect to the Hubble parameter using some algebraic approaches. A time varying scale factor has been introduced to analyze the behavior of the model. From some diagnostic approach, we found that the model behaves as $\Lambda_{CDM}$ model at late time of cosmic evolution.


title: Shadowing for codimension one sectional-Anosov flows

abstract: In hyperbolic dynamics, a well-known result is that every hyperbolic attracting set, have a finite pseudo-orbit tracing property (FPOTP). It's natural to wonder if this result is maintained in the sectional-hyperbolic dynamics; Komuro in [Lorenz attractors do not have the pseudo-orbit tracing property], provides a negative answer for this question, by proving that the geometric Lorenz Attractor doesn't have a FPOTP. In this paper, we generalized the result of Komuro, we prove that every codimension one sectional-hyperbolic attractor set with a unique singularity Lorenz-like, which is of boundary-type, does not have FPOTP.


title: Harmonic Standing-Wave Excitations of Simply-Supported Thick-Walled Hollow Elastic Circular Cylinders: Exact 3D Linear Elastodynamic Response

abstract: The forced-vibration response of a simply-supported isotropic thick-walled hollow elastic circular cylinder subjected to two-dimensional harmonic standing-wave excitations on its curved surfaces is studied within the framework of linear elastodynamics. Exact semi-analytical solutions for the steady-state displacement field of the cylinder are constructed using recently-published parametric solutions to the Navier-Lam\'{e} equation. Formal application of the standing-wave boundary conditions generates three parameter-dependent $6\times6$ linear systems, each of which can be numerically solved in order to determine the parametric response of the cylinder's displacement field under various conditions. The method of solution is direct and demonstrates a general approach that can be applied to solve many other elastodynamic forced-response problems involving isotropic elastic cylinders. As an application, and considering several examples, the obtained solution is used to compute the steady-state frequency response in a few specific low-order excitation cases. In each case, the solution generates a series of resonances that are in exact correspondence with a unique subset of the natural frequencies of the simply-supported cylinder. The considered problem is of general theoretical interest in structural mechanics and acoustics and more practically serves as a benchmark forced-vibration problem involving a thick-walled hollow elastic cylinder.


title: Accelerating Training in Artificial Neural Networks with Dynamic Mode Decomposition

abstract: Training of deep neural networks (DNNs) frequently involves optimizing several millions or even billions of parameters. Even with modern computing architectures, the computational expense of DNN training can inhibit, for instance, network architecture design optimization, hyper-parameter studies, and integration into scientific research cycles. The key factor limiting performance is that both the feed-forward evaluation and the back-propagation rule are needed for each weight during optimization in the update rule. In this work, we propose a method to decouple the evaluation of the update rule at each weight. At first, Proper Orthogonal Decomposition (POD) is used to identify a current estimate of the principal directions of evolution of weights per layer during training based on the evolution observed with a few backpropagation steps. Then, Dynamic Mode Decomposition (DMD) is used to learn the dynamics of the evolution of the weights in each layer according to these principal directions. The DMD model is used to evaluate an approximate converged state when training the ANN. Afterward, some number of backpropagation steps are performed, starting from the DMD estimates, leading to an update to the principal directions and DMD model. This iterative process is repeated until convergence. By fine-tuning the number of backpropagation steps used for each DMD model estimation, a significant reduction in the number of operations required to train the neural networks can be achieved. In this paper, the DMD acceleration method will be explained in detail, along with the theoretical justification for the acceleration provided by DMD. This method is illustrated using a regression problem of key interest for the scientific machine learning community: the prediction of a pollutant concentration field in a diffusion, advection, reaction problem.


title: Solving Differential Equations Using Neural Network Solution Bundles

abstract: The time evolution of dynamical systems is frequently described by ordinary differential equations (ODEs), which must be solved for given initial conditions. Most standard approaches numerically integrate ODEs producing a single solution whose values are computed at discrete times. When many varied solutions with different initial conditions to the ODE are required, the computational cost can become significant. We propose that a neural network be used as a solution bundle, a collection of solutions to an ODE for various initial states and system parameters. The neural network solution bundle is trained with an unsupervised loss that does not require any prior knowledge of the sought solutions, and the resulting object is differentiable in initial conditions and system parameters. The solution bundle exhibits fast, parallelizable evaluation of the system state, facilitating the use of Bayesian inference for parameter estimation in real dynamical systems.


title: Study of lockdown/testing mitigation strategies on stochastic SIR model and its comparison with South Korea, Germany and New York data

abstract: We are currently facing a highly critical case of a world-wide pandemic. The novel coronavirus (SARS-CoV-2, a.k.a. COVID-19) has proved to be extremely contagious and the original outbreak from Asia has now spread to all continents. This situation will fruitfully profit from the study in regards of the spread of the virus, assessing effective countermeasures to weight the impact of the adopted strategies. The standard Susceptible-Infectious-Recovered (SIR) model is a very successful and widely used mathematical model for predicting the spread of an epidemic. We adopt the SIR model on a random network and extend the model to include control strategies {\em lockdown} and {\em testing} -- two often employed mitigation strategies. The ability of these strategies in controlling the pandemic spread is investigated by varying the effectiveness with which they are implemented. The possibility of a second outbreak is evaluated in detail after the mitigation strategies are withdrawn. We notice that, in any case, a sudden interruption of such mitigation strategies will likely induce a resurgence of a second outbreak, whose peak will be correlated to the number of susceptible individuals. In fact, we find that a population will remain vulnerable to the infection until the herd immunity is achieved. We also test our model with real statistics and information on the epidemic spread in South Korea, Germany, and New York and find a remarkable agreement with the simulation data.


title: Discontinuous and Smooth Depth Completion with Binary Anisotropic Diffusion Tensor

abstract: We propose an unsupervised real-time dense depth completion from a sparse depth map guided by a single image. Our method generates a smooth depth map while preserving discontinuity between different objects. Our key idea is a Binary Anisotropic Diffusion Tensor (B-ADT) which can completely eliminate smoothness constraint at intended positions and directions by applying it to variational regularization. We also propose an Image-guided Nearest Neighbor Search (IGNNS) to derive a piecewise constant depth map which is used for B-ADT derivation and in the data term of the variational energy. Our experiments show that our method can outperform previous unsupervised and semi-supervised depth completion methods in terms of accuracy. Moreover, since our resulting depth map preserves the discontinuity between objects, the result can be converted to a visually plausible point cloud. This is remarkable since previous methods generate unnatural surface-like artifacts between discontinuous objects.


title: Asymmetric Stokes flow induced by a transverse point-force acting near a finite-sized elastic membrane

abstract: A deep understanding of the physical interactions between nanoparticles and target cell membranes is important in designing efficient nanocarrier systems for drug delivery applications. Here, we present a theoretical framework to describe the hydrodynamic flow field induced by a point-force singularity (Stokeslet) directed parallel to a finite-sized elastic membrane endowed with shear and bending rigidities. We formulate the elastohydrodynamic problem as a mixed-boundary-value problem, which we then reduce into a well-behaved system of integro-differential equations. It follows that shear and bending linearly decouple so that the solution of the overall flow problem can be obtained by linear superposition of the contributions arising from these modes of deformation. Additionally, we probe the effect of the membrane on the hydrodynamic drag acting on a nearby particle, finding that, in a certain range of parameters, translational motion near an elastic membrane with only energetic resistance toward shear can, surprisingly, be sped up compared to bulk fluid. Our results may find applications in microrheological characterizations of colloidal systems near elastic confinements.


title: Automatic Tuning of Stochastic Gradient Descent with Bayesian Optimisation

abstract: Many machine learning models require a training procedure based on running stochastic gradient descent. A key element for the efficiency of those algorithms is the choice of the learning rate schedule. While finding good learning rates schedules using Bayesian optimisation has been tackled by several authors, adapting it dynamically in a data-driven way is an open question. This is of high practical importance to users that need to train a single, expensive model. To tackle this problem, we introduce an original probabilistic model for traces of optimisers, based on latent Gaussian processes and an auto-/regressive formulation, that flexibly adjusts to abrupt changes of behaviours induced by new learning rate values. As illustrated, this model is well-suited to tackle a set of problems: first, for the on-line adaptation of the learning rate for a cold-started run; then, for tuning the schedule for a set of similar tasks (in a classical BO setup), as well as warm-starting it for a new task.


title: Spectral structure of the Neumann--Poincar\'e operator on thin domains in two dimensions

abstract: We consider the spectral structure of the Neumann--Poincar\'e operators defined on the boundaries of thin domains of rectangle shape in two dimensions. We prove that as the aspect ratio of the domains tends to $\infty$, or equivalently, as the domains get thinner, the spectra of the Neumann--Poincar\'e operators are densely distributed in the interval $[-1/2,1/2]$.


title: A Canonical Transform for Strengthening the Local $L^p$-Type Universal Approximation Property

abstract: Most $L^p$-type universal approximation theorems guarantee that a given machine learning model class $\mathscr{F}\subseteq C(\mathbb{R}^d,\mathbb{R}^D)$ is dense in $L^p_{\mu}(\mathbb{R}^d,\mathbb{R}^D)$ for any suitable finite Borel measure $\mu$ on $\mathbb{R}^d$. Unfortunately, this means that the model's approximation quality can rapidly degenerate outside some compact subset of $\mathbb{R}^d$, as any such measure is largely concentrated on some bounded subset of $\mathbb{R}^d$. This paper proposes a generic solution to this approximation theoretic problem by introducing a canonical transformation which "upgrades $\mathscr{F}$'s approximation property" in the following sense. The transformed model class, denoted by $\mathscr{F}\text{-tope}$, is shown to be dense in $L^p_{\mu,\text{strict}}(\mathbb{R}^d,\mathbb{R}^D)$ which is a topological space whose elements are locally $p$-integrable functions and whose topology is much finer than usual norm topology on $L^p_{\mu}(\mathbb{R}^d,\mathbb{R}^D)$; here $\mu$ is any suitable $\sigma$-finite Borel measure $\mu$ on $\mathbb{R}^d$. Next, we show that if $\mathscr{F}$ is any family of analytic functions then there is always a strict "gap" between $\mathscr{F}\text{-tope}$'s expressibility and that of $\mathscr{F}$, since we find that $\mathscr{F}$ can never dense in $L^p_{\mu,\text{strict}}(\mathbb{R}^d,\mathbb{R}^D)$. In the general case, where $\mathscr{F}$ may contain non-analytic functions, we provide an abstract form of these results guaranteeing that there always exists some function space in which $\mathscr{F}\text{-tope}$ is dense but $\mathscr{F}$ is not, while, the converse is never possible. Applications to feedforward networks, convolutional neural networks, and polynomial bases are explored.


title: $L^{p} \rightarrow L^{q}$ estimates for maximal functions associated with nonisotropic dilations of hypersurfaces in $\mathbb{R}^3$

abstract: The goal of this article is to establish $L^{p} \rightarrow L^{q}$ estimates for maximal functions associated with nonisotropic dilations $\delta_t(x)=(t^{a_1}x_1,t^{a_2}x_2,t^{a_3}x_3)$ of hypersurfaces $(x_{1}, x_{2},\Phi(x_1,x_2))$ in $\mathbb{R}^3$, where the Gaussian curvatures of the hypersurfaces are allowed to vanish. When $2 \alpha_{2} = \alpha_{3}$, this problem is reduced to study of the $L^{p} \rightarrow L^{q}$ estimates for maximal functions along the curve $\gamma(x)=(x,x^2(1+\phi(x)))$ and associated dilations $\delta_t(x)=(tx_1,t^2x_2)$. The corresponding maximal function shows features related to the Bourgain circular maximal function, whose $L^{p} \rightarrow L^{q}$ estimate has been considered by [Schlag, JAMS, 1997], [Schlag-Sogge, MRL, 1997] and [Lee, PAMS, 2003]. However, in the study of the maximal function related to the mentioned curve $\gamma(x)$ and associated dilations, we get the $L^{p} \rightarrow L^{q}$ regularity properties for a family of corresponding Fourier integral operators which fail to satisfy the "cinematic curvature condition" uniformly, which means that classical local smoothing estimates could not be directly applied to our problem. What's more, the $L^{p} \rightarrow L^{q}$ estimates are also new for maximal functions associated with isotropic dilations of hypersurfaces $(x_{1}, x_{2},\Phi(x_1,x_2))$ mentioned before.


title: Deep Convolutional GANs for Car Image Generation

abstract: In this paper, we investigate the application of deep convolutional GANs on car image generation. We improve upon the commonly used DCGAN architecture by implementing Wasserstein loss to decrease mode collapse and introducing dropout at the end of the discrimiantor to introduce stochasticity. Furthermore, we introduce convolutional layers at the end of the generator to improve expressiveness and smooth noise. All of these improvements upon the DCGAN architecture comprise our proposal of the novel BoolGAN architecture, which is able to decrease the FID from 195.922 (baseline) to 165.966.


title: Modeling and measuring the non-ideal characteristics of transmission lines

abstract: We describe a simple method to experimentally determine the frequency dependencies of the per-unit-length resistance and conductance of transmission lines. The experiment is intended as a supplement to the classic measurement of the transient response of a transmission line to a voltage step or pulse. In the transient experiment, an ideal (lossless) model of the transmission line is used to determine the characteristic impedance and signal propagation speed. In our experiment, the insertion losses of various coaxial cables are measured as a function of frequency from 1 to 2000 MHz. A full distributed circuit model of the transmission line that includes both conductor and dielectric losses is needed to fit the frequency dependence of the measured insertion losses. Our model assumes physically-sensible frequency dependencies for the per-unit-length resistance and conductance that are determined by the geometry of the coaxial transmission lines used in the measurements.


title: Coordination of OLTC and Smart Inverters for Optimal Voltage Regulation of Unbalanced Distribution Networks

abstract: Photovoltaic (PV) smart inverters can improve the voltage profile of distribution networks. A multi-objective optimization framework for coordination of reactive power injection of smart inverters and tap operations of on-load tap changers (OLTCs) for multi-phase unbalanced distribution systems is proposed. The optimization objective is to minimize voltage deviations and the number of tap operations simultaneously. A novel linearization method is proposed to linearize power flow equations and to convexify the problem, which guarantees convergence of the optimization and less computation costs. The optimization is modeled and solved using mixed-integer linear programming (MILP). The proposed method is validated against conventional rule-based autonomous voltage regulation (AVR) on the highly-unbalanced modified IEEE 37 bus test system and a large California utility feeder. Simulation results show that the proposed method accurately estimates feeder voltage, significantly reduces voltage deviations, mitigates over-voltage problems, and reduces voltage unbalance while eliminating unnecessary tap operations. The robustness of the method is validated against various levels of forecast error. The computational efficiency and scalability of the proposed approach are also demonstrated through the simulations on the large utility feeder.


title: Fractional derivatives and the fundamental theorem of Fractional Calculus

abstract: In this paper, we address the one-parameter families of the fractional integrals and derivatives defined on a finite interval. First we remind the reader of the known fact that under some reasonable conditions, there exists precisely one unique family of the fractional integrals, namely, the well-known Riemann-Liouville fractional integrals. As to the fractional derivatives, their natural definition follows from the fundamental theorem of the Fractional Calculus, i.e., they are introduced as the left-inverse operators to the Riemann-Liouville fractional integrals. Until now, three families of such derivatives were suggested in the literature: the Riemann-Liouville fractional derivatives, the Caputo fractional derivatives, and the Hilfer fractional derivatives. We clarify the interconnections between these derivatives on different spaces of functions and provide some of their properties including the formulas for their projectors and the Laplace transforms. However, it turns out that there exist infinitely many other families of the fractional derivatives that are the left-inverse operators to the Riemann-Liouville fractional integrals. In this paper, we focus on an important class of these fractional derivatives and discuss some of their properties.


title: Dual-Free Stochastic Decentralized Optimization with Variance Reduction

abstract: We consider the problem of training machine learning models on distributed data in a decentralized way. For finite-sum problems, fast single-machine algorithms for large datasets rely on stochastic updates combined with variance reduction. Yet, existing decentralized stochastic algorithms either do not obtain the full speedup allowed by stochastic updates, or require oracles that are more expensive than regular gradients. In this work, we introduce a Decentralized stochastic algorithm with Variance Reduction called DVR. DVR only requires computing stochastic gradients of the local functions, and is computationally as fast as a standard stochastic variance-reduced algorithms run on a $1/n$ fraction of the dataset, where $n$ is the number of nodes. To derive DVR, we use Bregman coordinate descent on a well-chosen dual problem, and obtain a dual-free algorithm using a specific Bregman divergence. We give an accelerated version of DVR based on the Catalyst framework, and illustrate its effectiveness with simulations on real data.


title: Computationally Efficient Attitude Estimation with Extended $\mathcal{H}_2$ Filtering

abstract: Accurate state estimation using low-cost MEMS (Micro Electro- Mechanical Systems) sensors present on Commercial-off-the-shelf (COTS) drones is a challenging problem. Most UAV systems use a combination of a gyroscope, an accelerometer, and a magnetometer to obtain measurements and estimate attitude. Under this paradigm of sensor fusion, the Extended Kalman Filter (EKF) is the most popular algorithm for attitude estimation in UAVs. In this work, we propose a novel estimation technique called extended H2 filter that can overcome the limitations of the EKF, specifically with respect to computational speed, memory usage, and root mean squared error. We formulate our attitude-estimation algorithm using unit quaternions. The H2 optimal filter gain is designed offline about a nominal operating point by solving a convex optimization problem, and the filter dynamics is implemented using the nonlinear system dynamics. This implementation of this H2 optimal estimator is referred as the extended H2 estimator. The proposed technique is tested on four cases corresponding to long time-scale motion, fast time-scale motion, transition from hover to forward flight for VTOL aircrafts, and an entire flight cycle (from take-off to landing). Its results are compared against that of the EKF in terms of the aforementioned performance metrics.


title: Periodic String Complexes over String Algebras

abstract: In this paper we develop combinatorial techniques for the case of string algebras with the aim to give a characterization of string complexes with infinite minimal projective resolution. These complexes will be called \textit{periodic string complexes}. As a consequence of this characterization, we give two important applications. The first one, is a sufficient condition for a string algebra to have infinite global dimension. In the second one, we exhibit a class of indecomposable objects in the derived category for a special case of string algebras. Every construction, concept and consequence in this paper is followed by some illustrative examples.


title: Normalized ground states for semilinear elliptic systems with critical and subcritical nonlinearities

abstract: In the present paper, we study the normalized solutions with least energy to the following system: $$\begin{cases} -\Delta u+\lambda_1u=\mu_1 |u|^{p-2}u+\beta r_1|u|^{r_1-2}|v|^{r_2}u\quad &\hbox{in}\;\mathbb R^N,\\ -\Delta v+\lambda_2v=\mu_2 |v|^{q-2}v+\beta r_2|u|^{r_1}|v|^{r_2-2}v\quad&\hbox{in}\;\mathbb R^N,\\ \int_{\mathbb R^N}u^2=a_1^2\quad\hbox{and}\;\int_{\mathbb R^N}v^2=a_2^2, \end{cases}$$ where $p,q,r_1+r_2$ can be Sobolev critical. To this purpose, we study the geometry of the Pohozaev manifold and the associated minimizition problem. Under some assumption on $a_1,a_2$ and $\beta$, we obtain the existence of the positive normalized ground state solution to the above system. We have solved some unsolved open problems in this area.


title: A study of $c c\bar{c}\bar{c}$ tetraquark decays in 4 muons and in $D^{(*)} \bar{D}^{(*)}$ at LHC

abstract: We perform a quantitative analysis of the decays of $cc\bar c\bar c$ tetraquarks with $J^{PC}=0^{++}, 2^{++}$ into 4 muons and into hidden- and open-charm mesons and estimate, for the first time, the fully charmed tetraquark decay width. The calculated cross section upper limit is $\sim 40$ fb for the 4 muons channel, and $\sim 28$ nb for the $D^{(*)} \bar D^{(*)} \to e\mu$ channel. On the basis of our results, with the present sensitivity LHCb should detect both signals, for $0^{++}$ and $2^{++}$ fully-charmed tetraquarks.


title: Reinforcement Learning for Non-Stationary Markov Decision Processes: The Blessing of (More) Optimism

abstract: We consider un-discounted reinforcement learning (RL) in Markov decision processes (MDPs) under drifting non-stationarity, i.e., both the reward and state transition distributions are allowed to evolve over time, as long as their respective total variations, quantified by suitable metrics, do not exceed certain variation budgets. We first develop the Sliding Window Upper-Confidence bound for Reinforcement Learning with Confidence Widening (SWUCRL2-CW) algorithm, and establish its dynamic regret bound when the variation budgets are known. In addition, we propose the Bandit-over-Reinforcement Learning (BORL) algorithm to adaptively tune the SWUCRL2-CW algorithm to achieve the same dynamic regret bound, but in a parameter-free manner, i.e., without knowing the variation budgets. Notably, learning non-stationary MDPs via the conventional optimistic exploration technique presents a unique challenge absent in existing (non-stationary) bandit learning settings. We overcome the challenge by a novel confidence widening technique that incorporates additional optimism.


title: A New Modal Autoencoder for Functionally Independent Feature Extraction

abstract: Autoencoders have been widely used for dimensional reduction and feature extraction. Various types of autoencoders have been proposed by introducing regularization terms. Most of these regularizations improve representation learning by constraining the weights in the encoder part, which maps input into hidden nodes and affects the generation of features. In this study, we show that a constraint to the decoder can also significantly improve its performance because the decoder determines how the latent variables contribute to the reconstruction of input. Inspired by the structural modal analysis method in mechanical engineering, a new modal autoencoder (MAE) is proposed by othogonalising the columns of the readout weight matrix. The new regularization helps to disentangle explanatory factors of variation and forces the MAE to extract fundamental modes in data. The learned representations are functionally independent in the reconstruction of input and perform better in consecutive classification tasks. The results were validated on the MNIST variations and USPS classification benchmark suite. Comparative experiments clearly show that the new algorithm has a surprising advantage. The new MAE introduces a very simple training principle for autoencoders and could be promising for the pre-training of deep neural networks.


title: Ladder operators and a second--order difference equation for general discrete Sobolev orthogonal polynomials

abstract: We consider a general discrete Sobolev inner product involving the Hahn difference operator, so this includes the well--known difference operators $\mathscr{D}_{q}$ and $\Delta$ and, as a limit case, the derivative operator. The objective is twofold. On the one hand, we construct the ladder operators for the corresponding nonstandard orthogonal polynomials and we obtain the second--order difference equation satisfied by these polynomials. On the other hand, we generalise some related results appeared in the literature as we are working in a more general framework. Moreover, we will show that all the functions involved in these constructions can be computed explicitly.


title: Spectral analysis of the multi-dimensional diffusion operator with random jumps from the boundary

abstract: We develop a Hilbert-space approach to the diffusion process of the Brownian motion in a bounded domain with random jumps from the boundary introduced by Ben-Ari and Pinsky in 2007. The generator of the process is introduced by a diffusion elliptic differential operator in the space of square-integrable functions, subject to non-self-adjoint and non-local boundary conditions expressed through a probability measure on the domain. We obtain an expression for the difference between the resolvent of the operator and that of its Dirichlet realization. We prove that the numerical range is the whole complex plane, despite the fact that the spectrum is purely discrete and is contained in a half-plane. Furthermore, for the class of absolutely continuous probability measures with square-integrable densities we characterise the adjoint operator and prove that the system of root vectors is complete. Finally, under certain assumptions on the densities, we obtain enclosures for the non-real spectrum and find a sufficient condition for the non-zero eigenvalue with the smallest real part to be real. The latter supports the conjecture of Ben-Ari and Pinsky that this eigenvalue is always real.


title: The GAPS Programme at TNG -- XXIV. An eccentric Neptune-mass planet near the inner edge of the BD-11 4672 habitable zone

abstract: With the growth of comparative exoplanetology, it is increasingly clear that investigating the relationships between inner and outer planets plays a key role in discriminating between competing formation and evolution models. To do so, it is important to probe the inner region of systems hosting long-period giants in search for undetected lower-mass planetary companions. In this work we present our results on the K-dwarf star BD-11~4672, already known to host a long-period giant planet, as the first output of a subsample of the GAPS programme specifically aimed at assessing the impact of inefficient migration of planets formed beyond the snowline by searching for Neptune-mass and super-Earths planetary companions of known longer-period giants. The high-precision HARPS-N observations of BD-11~4672 are used in conjunction with literature time series in order to search for additional inner planetary signals to be fitted using differential evolution Markov chain Monte Carlo. The long-term stability of the new orbital solutions is tested by N-body dynamical simulations. We report the detection of BD-11~4672~c, a new Neptune-mass planet with an orbital period of $74.20_{-0.08}^{+0.06}\,\rm{d}$, eccentricity $0.40_{-0.15}^{+0.13}$, semimajor axis $0.30\pm0.01\,\rm{au}$ and minimum mass $15.37_{-2.81}^{+2.97}\,M_\oplus$ orbiting slightly outside the inner edge of the optimistic circumstellar habitable zone. In order to assess its impact on the dynamical stability of the habitable zone we compute the angular momentum deficit of the system, showing that planet c has a severe negative impact on the stability of possible additional lower-mass temperate planets. The BD-11~4672 system is notable for its architecture, hosting both a long-period giant planet and an inner lower-mass planet, the latter being also among the most eccentric Neptune-mass planets known at similar periods.


title: Transition magnon modes in thin ferromagnetic nanogratings

abstract: This work presents micromagnetic simulations in ferromagnetic nanogratings for the full range of directions of an applied in-plane external magnetic field. We focus on the modification of the magnon mode characteristics when the magnetic field orientation is gradually changed between the classical Damon-Eshbach (DE) and backward-volume (BV) geometries. We found that in a specific range of field directions, the magnon mode parameters differ significantly from the parameters in the classical cases, namely, the modes are characterized by complex spatial distributions and have low group velocities. The center of this range corresponds to the direction of the external magnetic field, which gives the maximal nonuniform distribution of the static magnetization in the nanogratings.

categories: cond-mat.mes-hall cond-mat.other physics.app-ph physics.comp-ph


title: Training Variational Networks with Multi-Domain Simulations: Speed-of-Sound Image Reconstruction

abstract: Speed-of-sound has been shown as a potential biomarker for breast cancer imaging, successfully differentiating malignant tumors from benign ones. Speed-of-sound images can be reconstructed from time-of-flight measurements from ultrasound images acquired using conventional handheld ultrasound transducers. Variational Networks (VN) have recently been shown to be a potential learning-based approach for optimizing inverse problems in image reconstruction. Despite earlier promising results, these methods however do not generalize well from simulated to acquired data, due to the domain shift. In this work, we present for the first time a VN solution for a pulse-echo SoS image reconstruction problem using diverging waves with conventional transducers and single-sided tissue access. This is made possible by incorporating simulations with varying complexity into training. We use loop unrolling of gradient descent with momentum, with an exponentially weighted loss of outputs at each unrolled iteration in order to regularize training. We learn norms as activation functions regularized to have smooth forms for robustness to input distribution variations. We evaluate reconstruction quality on ray-based and full-wave simulations as well as on tissue-mimicking phantom data, in comparison to a classical iterative (L-BFGS) optimization of this image reconstruction problem. We show that the proposed regularization techniques combined with multi-source domain training yield substantial improvements in the domain adaptation capabilities of VN, reducing median RMSE by 54% on a wave-based simulation dataset compared to the baseline VN. We also show that on data acquired from a tissue-mimicking breast phantom the proposed VN provides improved reconstruction in 12 milliseconds.


title: Q-NET: A Network for Low-Dimensional Integrals of Neural Proxies

abstract: Many applications require the calculation of integrals of multidimensional functions. A general and popular procedure is to estimate integrals by averaging multiple evaluations of the function. Often, each evaluation of the function entails costly computations. The use of a \emph{proxy} or surrogate for the true function is useful if repeated evaluations are necessary. The proxy is even more useful if its integral is known analytically and can be calculated practically. We propose the use of a versatile yet simple class of artificial neural networks -- sigmoidal universal approximators -- as a proxy for functions whose integrals need to be estimated. We design a family of fixed networks, which we call Q-NETs, that operate on parameters of a trained proxy to calculate exact integrals over \emph{any subset of dimensions} of the input domain. We identify transformations to the input space for which integrals may be recalculated without resampling the integrand or retraining the proxy. We highlight the benefits of this scheme for a few applications such as inverse rendering, generation of procedural noise, visualization and simulation. The proposed proxy is appealing in the following contexts: the dimensionality is low ($<10$D); the estimation of integrals needs to be decoupled from the sampling strategy; sparse, adaptive sampling is used; marginal functions need to be known in functional form; or when powerful Single Instruction Multiple Data/Thread (SIMD/SIMT) pipelines are available for computation.


title: Approximate controllability of the semilinear reaction-diffusion equation governed by a multiplicative control

abstract: In this paper we are concerned with the approximate controllability of a multidimensional semilinear reaction-diffusion equation governed by a multiplicative control, which is locally distributed in the reaction term. For a given initial state we provide sufficient conditions on the desirable state to be approximately reached within an arbitrarily small time interval. Moreover, in the case of a globally supported control, we prove the approximate controllability within any time-interval given in advance which does not depend on the initial and target states. Our approaches are based on linear semigroup theory and some results on uniform approximation with smooth functions.


title: Periodic waves in the fractional modified Korteweg--de Vries equation

abstract: Periodic waves in the modified Korteweg-de Vries (mKdV) equation are revisited in the setting of the fractional Laplacian. Two families of solutions in the local case are given by the sign-definite dnoidal and sign-indefinite cnoidal solutions. Both solutions can be characterized in the general fractional case as global minimizers of the quadratic part of the energy functional subject to the fixed $L^4$ norm: the sign-definite (sign-indefinite) solutions are obtained in the subspace of even (odd) functions. Morse index is computed for both solutions and the spectral stability criterion is derived. We show numerically that the family of sign-definite solutions has a generic fold bifurcation for the fractional Laplacian of lower regularity and the family of sign-indefinite solutions has a generic symmetry-breaking bifurcation both in the fractional and local cases.

categories: math.AP math-ph math.DS math.MP nlin.PS nlin.SI


title: Modifications on parameters of $Z(4430)$ in a dense medium

abstract: The charmonium-like resonance $ Z_c(3900) $ and its excited state $ Z(4430) $ are among the particles that are serious candidates for double heavy tetraquarks. Calculations of different parameters associated with these states both in the vacuum and the medium with finite density are of great importance. Such investigations help us clarify their nature, internal quark-gluon organization and quantum numbers. In this accordance, we extend our previous analyses on the ground state $ Z_c(3900) $ to investigate the medium modifications on different parameters of the excited $ Z(4430) $ state. In particular, we calculate the mass, vector self-energy and current coupling of $ Z(4430) $ in terms of density, up to a density comparable to the density of the cores of massive neutron stars. The obtained results may help experimental groups aiming to study the behavior of exotic states at higher densities.