title: Ergodicity and invariant measures for a diffusing passive scalar advected by a random channel shear flow and the connection between the Kraichnan-Majda model and Taylor-Aris Dispersion

abstract: We study the long time behavior of an advection-diffusion equation with a random shear flow which depends on a stationary Ornstein-Uhlenbeck (OU) process in parallel-plate channels enforcing the no-flux boundary conditions. We derive a closed form formula for the long time asymptotics of the arbitrary $N$-point correlator using the ground state eigenvalue perturbation approach proposed in \cite{bronski1997scalar}. In turn, appealing to the conclusion of the Hausdorff moment problem \cite{shohat1943problem}, we discover a diffusion equation with a random drift and deterministic enhanced diffusion possessing the exact same probability distribution function at long times. Such equations enjoy many ergodic properties which immediately translate to ergodicity results for the original problem. In particular, we establish that the first two Aris moments using a single realization of the random field can be used to explicitly construct all ensemble averaged moments. Also, the first two ensemble averaged moments explicitly predict any long time centered Aris moment. Our formulae quantitatively depict the dependence of the deterministic effective diffusion on the interaction between spatial structure of flow and random temporal fluctuation. This noteably contrasts the white noise case from the OU case. Further, this approximation provides many identities regarding the stationary OU process dependent time integral. We derive explicit formulae for the decaying passive scalar's long time limiting probability distribution function (PDF) for different types of initial conditions (e.g. deterministic and random). All results are verified by Monte-Carlo simulations.


title: A Stringent Upper Limit on Dark Matter Self-Interaction Cross Section from Cluster Strong Lensing

abstract: We analyze strongly lensed images in 8 galaxy clusters to measure their dark matter density profiles in the radial region between 10 kpc and 150 kpc, and use this to constrain the self-interaction cross section of dark matter (DM) particles. We infer the mass profiles of the central DM halos, bright central galaxies, key member galaxies, and DM subhalos for the member galaxies for all 8 clusters using the Qlens code. The inferred DM halo surface densities are fit to a self-interacting dark matter (SIDM) model, which allows us to constrain the self-interaction cross section over mass $\sigma$/m. When our full method is applied to mock data generated from two clusters in the Illustris-TNG simulation, we find results consistent with no dark matter self-interactions as expected. For the eight observed clusters with average relative velocities of $1458_{-81}^{+80}$ km/s, we infer $\sigma$/m = $0.082_{-0.021}^{+0.027}$ cm$^2$/g and $\sigma$/m < 0.13 cm$^2$/g at the 95% confidence level.


title: Electronic excited states of benzene in interaction with water clusters : influence of structure and size

abstract: This work is dedicated to the theoretical investigation of the influence of water clusters' organisation and size on the electronic spectrum of an interacting benzene (Bz) molecule using both TD-DFT and CASPT2 approaches. Two series of geometries, namely $Geo_{IEI}$ and $Geo_{IED}$ were extracted from two Bz-hexagonal ice configurations leading to maximum and minimum ionization energies respectively. An appropriate basis set containing atomic diffuse and polarisation orbitals and describing the Rydberg states of Bz was determined. The TD-DFT approach was carefully benchmarked against CASPT2 results for the smallest systems.Despite some discrepancies, the trends were found to be similar at both levels of theory: the positions and intensities of the main $\pi \rightarrow \pi^{\star}$ transitions were found slightly split due to symmetry breaking. For the smallest systems, our results clearly show the dependence of the electronic transitions on the clusters' structures. Of particular interest, low energy transitions of non negligible oscillator strength from a Bz $\pi$ orbital to a virtual orbital of Rydberg character, also involving atomic diffuse functions and partially expanded on the water cluster, were found for the $Geo_{IED}$ series. The energies of such transitions were determined to be more than 2\,eV below the ionization potential of Bz. When the cluster's size increases, similar transitions were found for all structures, the virtual orbitals becoming mainly developed on the H atoms of the water molecules at the edge of the cluster. Given their nature and energy, such transitions could play a role in the photochemistry of aromatic species in interaction with water clusters or ice, such processes being of astrophysical interest.


title: Beyond Scaling: Calculable Error Bounds of the Power-of-Two-Choices Mean-Field Model in Heavy-Traffic

abstract: This paper provides a recipe for deriving calculable approximation errors of mean-field models in heavy-traffic with the focus on the well-known load balancing algorithm -- power-of-two-choices (Po2). The recipe combines Stein's method for linearized mean-field models and State Space Concentration (SSC) based on geometric tail bounds. In particular, we divide the state space into two regions, a neighborhood near the mean-field equilibrium and the complement of that. We first use a tail bound to show that the steady-state probability being outside the neighborhood is small. Then, we use a linearized mean-field model and Stein's method to characterize the generator difference, which provides the dominant term of the approximation error. From the dominant term, we are able to obtain an asymptotically-tight bound and a nonasymptotic upper bound, both are calculable bounds, not order-wise scaling results like most results in the literature. Finally, we compared the theoretical bounds with numerical evaluations to show the effectiveness of our results. We note that the simulation results show that both bounds are valid even for small size systems such as a system with only ten servers.


title: Improving fretting corrosion resistance of CoCrMo alloy with TiSiN and ZrN coatings for orthopedic applications

abstract: Total hip replacement is the most effective treatment for late stage osteoarthritis. However, adverse local tissue reactions (ALTRs) associated with fretting corrosion have been observed in patients with modular total hip implants. The purpose of this study is to increase the fretting corrosion resistance of the CoCrMo alloy and the associated metal ion release by applying hard coatings to the surface. Cathodic arc evaporation technique (arc-PVD) was used to deposit TiSiN and ZrN hard coatings on to CoCrMo substrates. The morphology, chemical composition, crystal structures and residual stress of the coatings were characterized by scanning electron microscopy, energy dispersive x-ray spectroscopy, and X-ray diffractometry. Hardness, elastic modulus, and adhesion of the coatings were measured by nano-indentation, nano-scratch test, and the Rockwell C test. Fretting corrosion resistance tests of coated and uncoated CoCrMo discs against Ti6Al4V spheres were conducted on a four-station fretting testing machine in simulated body fluid at 1Hz for 1 million cycles. Post-fretting samples were analyzed for morphological changes, volume loss and metal ion release. Our analyses showed better surface finish and lower residual stress for ZrN coating, but higher hardness and better scratch resistance for TiSiN coating. Fretting results demonstrated substantial improvement in fretting corrosion resistance of CoCrMo with both coatings. ZrN and TiSiN decreased fretting volume loss by more than 10 times and 1000 times, respectively. Both coatings showed close to 90% decrease of Co ion release during fretting corrosion tests. Our results suggest that hard coating deposition on CoCrMo alloy can significantly improve its fretting corrosion resistance and could thus potentially alleviate ALTRs in metal hip implants.


title: Probabilistic Conditional System Invariant Generation with Bayesian Inference

abstract: Invariants are a set of properties over program attributes that are expected to be true during the execution of a program. Since developing those invariants manually can be costly and challenging, there are a myriad of approaches that support automated mining of likely invariants from sources such as program traces. Existing approaches, however, are not equipped to capture the rich states that condition the behavior of autonomous mobile robots, or to manage the uncertainty associated with many variables in these systems. This means that valuable invariants that appear only under specific states remain uncovered. In this work we introduce an approach to infer conditional probabilistic invariants to assist in the characterization of the behavior of such rich stateful, stochastic systems. These probabilistic invariants can encode a family of conditional patterns, are generated using Bayesian inference to leverage observed trace data against priors gleaned from previous experience and expert knowledge, and are ranked based on their surprise value and information content. Our studies on two semi-autonomous mobile robotic systems show how the proposed approach is able to generate valuable and previously hidden stateful invariants.


title: Experimental QCD summary (ICHEP 2020)

abstract: This writeup summarizes the main experimental studies of the strong interaction, theoretically described by quantum chromodynamics (QCD), that were presented during the ICHEP-2020 conference. The latest results, measured mostly in p-p collisions at the LHC, are categorized in seven broad topics: (i) Extractions of the strong coupling constant $\alpha_s(m_{\rm Z})$; (ii) Comparison of data to fixed-order (N$^{\rm n}$LO) perturbative QCD calculations; (iii) Determinations of parton distribution functions (PDFs); (iv) Comparison of data to resummed (N$^{\rm n}$LL) pQCD calculations; (v) Parton showering and jet substructure analyses; (vi) Semihard (double parton interactions, multiparton interactions, hard diffraction), and soft (elastic and diffractive) QCD scatterings; and (vii) Studies of parton hadronization in $\rm e^+e^-$ and p-p collisions.


title: Pullback and uniform attractors for nonautonomous reaction-diffusion equation in Dumbbell domains

abstract: This work is devoted to the study of the asymptotic behavior of nonautonomous reaction-diffusion equations in Dumbbell domains $\Omega_{\varepsilon} \subset \mathbb{R}^{N}$. Each $\Omega_{\varepsilon}$ is the union of a fixed open set $\Omega$ and a channel $R_{\varepsilon}$ that collapses to a line segment $R_0$ as $\varepsilon \rightarrow 0^{+}$. We first establish the global existence of solution for each problem by using two properties of the parabolic equation considered, which are the positivity of the solutions and comparison results for them. We prove the existence of pullback and uniform attractors and we obtain uniform bounds (in $\varepsilon$) for them.


title: Slowly rotating black holes in Quasi-topological gravity

abstract: While cubic Quasi-topological gravity is unique, there is a family of quartic Quasi-topological gravities in five dimensions. These theories are defined by leading to a first order equation on spherically symmetric spacetimes, resembling the structure of the equations of Lovelock theories in higher-dimensions, and are also ghost free around AdS. Here we construct slowly rotating black holes in these theories, and show that the equations for the off-diagonal components of the metric in the cubic theory are automatically of second order, while imposing this as a restriction on the quartic theories allows to partially remove the degeneracy of these theories, leading to a three-parameter family of Lagrangians of order four in the Riemann tensor. This shows that the parallel with Lovelock theory observed on spherical symmetry, extends to the realm of slowly rotating solutions. In the quartic case, the equations for the slowly rotating black hole are obtained from a consistent, reduced action principle. These functions admit a simple integration in terms of quadratures. Finally, in order to go beyond the slowly rotating regime, we explore the consistency of the Kerr-Schild ansatz in cubic Quasi-topological gravity. Requiring the spacetime to asymptotically match with the rotating black hole in GR, for equal oblateness parameters, the Kerr-Schild deformation of an AdS vacuum, does not lead to a solution for generic values of the coupling. This result suggests that in order to have solutions with finite rotation in Quasi-topological gravity, one must go beyond the Kerr-Schild ansatz.


title: A systematic DECam search for RR Lyrae in the outer halo of the Milky Way

abstract: The discovery of very distant stars in the halo of the Milky Way provides valuable tracers on the Milky Way mass and its formation. Beyond 100 kpc from the Galactic center, most of the stars are likely to be in faint dwarf galaxies or tidal debris from recently accreted dwarfs, making the outer reaches of the Galaxy important for understanding the Milky Way's accretion history. However, distant stars in the halo are scarce. In that context, RR Lyrae are ideal probes of the distant halo as they are intrinsically bright and thus can be seen at large distances, follow well-known period-luminosity relations that enable precise distance measurements, and are easily identifiable in time-series data. Therefore, a detailed study of RR Lyrae will help us understand the accreted outskirts of the Milky Way. In this contribution, we present the current state of our systematic search for distant RR Lyrae stars in the halo using the DECam imager at the 4m telescope on Cerro Tololo (Chile). The total surveyed area consists of more than 110 DECam fields (~ 350 sq. deg) and includes two recent independent campaigns carried out in 2017 and 2018 with which we have detected > 650 candidate RR Lyrae stars. Here we describe the methodology followed to analyze the two latest campaigns. Our catalog contains a considerable number of candidate RR Lyrae beyond 100 kpc, and reaches out up to ~ 250 kpc. The number of distant RR Lyrae found is consistent with recent studies of the outer halo. These stars provide a set of important probes of the mass of the Milky Way, the nature of the halo, and the accretion history of the Galactic outskirts.


title: The Michigan Infrared Test Thermal ELT N-band (MITTEN) Cryostat

abstract: We introduce the Michigan Infrared Test Thermal ELT N-band (MITTEN) Cryostat, a new facility for testing infrared detectors with a focus on mid-infrared (MIR) wavelengths (8-13 microns). New generations of large format, deep well, fast readout MIR detectors are now becoming available to the astronomical community. As one example, Teledyne Imaging Sensors (TIS) has introduced a long-wave Mercury-Cadmium-Telluride (MCT) array, GeoSnap, with high quantum efficiency (> 65 %) and improved noise properties compared to previous generation Si:As blocked impurity band (BIB) detectors. GeoSnap promises improved sensitivities, and efficiencies, for future background-limited MIR instruments, in particular with future extremely large telescopes (ELTs). We describe our new test facility suitable for measuring characteristics of these detectors, such as read noise, dark current, linearity, gain, pixel operability, quantum efficiency, and point source imaging performance relative to a background scene, as well as multiple point sources of differing contrast. MITTEN has an internal light source, and soon an accompanying filter wheel and aperture plate, reimaged onto the detector using an Offner relay. The baseline temperature of the cryostat interior is maintained < 40 K and the optical bench maintains a temperature of 16 K using a two-stage pulse-tube cryocooler package from Cryomech. No measurable background radiation from the cryostat interior has yet been detected.


title: The 2017 May 20$^{\rm th}$ stellar occultation by the elongated centaur (95626) 2002 GZ$_{32}$

abstract: We predicted a stellar occultation of the bright star Gaia DR1 4332852996360346368 (UCAC4 385-75921) (m$_{\rm V}$= 14.0 mag) by the centaur 2002 GZ$_{32}$ for 2017 May 20$^{\rm th}$. Our latest shadow path prediction was favourable to a large region in Europe. Observations were arranged in a broad region inside the nominal shadow path. Series of images were obtained with 29 telescopes throughout Europe and from six of them (five in Spain and one in Greece) we detected the occultation. This is the fourth centaur, besides Chariklo, Chiron and Bienor, for which a multi-chord stellar occultation is reported. By means of an elliptical fit to the occultation chords we obtained the limb of 2002 GZ$_{32}$ during the occultation, resulting in an ellipse with axes of 305 $\pm$ 17 km $\times$ 146 $\pm$ 8 km. From this limb, thanks to a rotational light curve obtained shortly after the occultation, we derived the geometric albedo of 2002 GZ$_{32}$ ($p_{\rm V}$ = 0.043 $\pm$ 0.007) and a 3-D ellipsoidal shape with axes 366 km $\times$ 306 km $\times$ 120 km. This shape is not fully consistent with a homogeneous body in hydrostatic equilibrium for the known rotation period of 2002 GZ$_{32}$. The size (albedo) obtained from the occultation is respectively smaller (greater) than that derived from the radiometric technique but compatible within error bars. No rings or debris around 2002 GZ$_{32}$ were detected from the occultation, but narrow and thin rings cannot be discarded.


title: Phase transition in Modified Newtonian Dynamics (MONDian) self-gravitating systems

abstract: We study the statistical mechanics of binary systems under gravitational interaction of the Modified Newtonian Dynamics (MOND) in three-dimensional space. Considering the binary systems, in the microcanonical and canonical ensembles, we show that in the microcanonical systems, unlike the Newtonian gravity, there is a sharp phase transition, with a high-temperature homogeneous phase and a low temperature clumped binary one. Defining an order parameter in the canonical systems, we find a smoother phase transition and identify the corresponding critical temperature in terms of physical parameters of the binary system.


title: DPG approximation of eigenvalue problems

abstract: In this paper, the discontinuous Petrov--Galerkin approximation of the Laplace eigenvalue problem is discussed. We consider in particular the primal and ultra weak formulations of the problem and prove the convergence together with a priori error estimates. Moreover, we propose two possible error estimators and perform the corresponding a posteriori error analysis. The theoretical results are confirmed numerically and it is shown that the error estimators can be used to design an optimally convergent adaptive scheme.


title: Boron Content and the Superconducting Critical Temperature of Carbon-Based Materials

abstract: In this paper, we present results on magnetization properties of boron nitride-carbon (BN-C) and boron carbide-carbon (B4C-C) granular mixtures. The temperature-dependent magnetization for field-cooled during cooling and field-cooled during warming shows a kind of thermal hysteresis that is always seen around a metamagnetic phase transition from an antiferromagnetic martensite to a ferromagnetic austenite phase. The low-temperature magnetization has an upward turn that can be attributed to superparamagnetism, diamagnetic shielding, and trapped flux characteristic to high-temperature superconducting materials. After subtracting the diamagnetic background, the field-dependent magnetization loops M(B) are ferromagnetic-like, more significant for the BN-C than for the B4C-C mixture. In addition, the magnetization loops show the kink feature characteristic to granular superconductivity. The irreversibility temperature for a B4C-C mixture having 37.5 wt% B is Tc = 76 K. Combining our data with previous results on B-doped diamond and Q-carbon, we find that Tc increases linearly with the B concentration.


title: Reconstruction of Pairwise Interactions using Energy-Based Models

abstract: Pairwise models like the Ising model or the generalized Potts model have found many successful applications in fields like physics, biology, and economics. Closely connected is the problem of inverse statistical mechanics, where the goal is to infer the parameters of such models given observed data. An open problem in this field is the question of how to train these models in the case where the data contain additional higher-order interactions that are not present in the pairwise model. In this work, we propose an approach based on Energy-Based Models and pseudolikelihood maximization to address these complications: we show that hybrid models, which combine a pairwise model and a neural network, can lead to significant improvements in the reconstruction of pairwise interactions. We show these improvements to hold consistently when compared to a standard approach using only the pairwise model and to an approach using only a neural network. This is in line with the general idea that simple interpretable models and complex black-box models are not necessarily a dichotomy: interpolating these two classes of models can allow to keep some advantages of both.


title: On the symmetric lamination convex and quasiconvex hull for the coplanar n-well problem in two dimensions

abstract: We study some particular cases of the $n$-well problem in two-dimensional linear elasticity. Assuming that every well in $\mathcal{U}\subset\mathbb{R}^{2\times 2}_\text{sym}$ belong to the same two-dimensional affine subspace, we characterize the symmetric lamination convex hull $L^e(\mathcal{U})$ for any number of wells in terms of the symmetric lamination convex hull of all three-well subsets contained in $\mathcal{U}$. For a family of four-well sets where two pairs of wells are rank-one compatible, we show that the symmetric lamination convex and quasiconvex hulls coincide, but are strictly contained in its convex hull $C(\mathcal{U})$. We extend this result to some particular configurations of $n$ wells. Most of the proofs are constructive, and we also present explicit examples.


title: Cohomological $\chi$-independence for moduli of one-dimensional sheaves and moduli of Higgs bundles

abstract: We prove that the intersection cohomology (together with the perverse and the Hodge filtrations) for the moduli space of one-dimensional semistable sheaves supported in an ample curve class on a toric del Pezzo surface is independent of the Euler characteristic of the sheaves. We also prove an analogous result for the moduli space of semistable Higgs bundles with respect to an effective divisor $D$ of degree $\mathrm{deg}(D)>2g-2$. Our results confirm the cohomological $\chi$-independence conjecture by Bousseau for $\mathbb{P}^2$, and verify Toda's conjecture for Gopakumar-Vafa invariants for certain local curves and local surfaces. For the proof, we combine a generalized version of Ng\^o's support theorem, a dimension estimate for the stacky Hilbert-Chow morphism, and a splitting theorem for the morphism from the moduli stack to the good GIT quotient.


title: Sat2Vid: Street-view Panoramic Video Synthesis from a Single Satellite Image

abstract: We present a novel method for synthesizing both temporally and geometrically consistent street-view panoramic video from a single satellite image and camera trajectory. Existing cross-view synthesis approaches focus on images, while video synthesis in such a case has not yet received enough attention. For geometrical and temporal consistency, our approach explicitly creates a 3D point cloud representation of the scene and maintains dense 3D-2D correspondences across frames that reflect the geometric scene configuration inferred from the satellite view. As for synthesis in the 3D space, we implement a cascaded network architecture with two hourglass modules to generate point-wise coarse and fine features from semantics and per-class latent vectors, followed by projection to frames and an upsampling module to obtain the final realistic video. By leveraging computed correspondences, the produced street-view video frames adhere to the 3D geometric scene structure and maintain temporal consistency. Qualitative and quantitative experiments demonstrate superior results compared to other state-of-the-art synthesis approaches that either lack temporal consistency or realistic appearance. To the best of our knowledge, our work is the first one to synthesize cross-view images to video.


title: A Comparison of Proton Stopping Power Measured with Proton CT and x-ray CT in Fresh Post-Mortem Porcine Structures

abstract: Purpose: Currently, calculations of proton range in proton therapy patients are based on a conversion of CT Hounsfield Units of patient tissues into proton relative stopping power. Uncertainties in this conversion necessitate larger proximal and distal planned target volume margins. Proton CT can potentially reduce these uncertainties by directly measuring proton stopping power. We aim to demonstrate proton CT imaging with complex porcine samples, to analyze in detail three-dimensional regions of interest, and to compare proton stopping powers directly measured by proton CT to those determined from x-ray CT scans. Methods: We have used a prototype proton imaging system with single proton tracking to acquire proton radiography and proton CT images of a sample of porcine pectoral girdle and ribs, and a pig's head. We also acquired close in time x-ray CT scans of the same samples, and compared proton stopping power measurements from the two modalities. In the case of the pig's head, we obtained x-ray CT scans from two different scanners, and compared results from high-dose and low-dose settings. Results: Comparing our reconstructed proton CT images with images derived from x-ray CT scans, we find agreement within 1% to 2% for soft tissues, and discrepancies of up to 6% for compact bone. We also observed large discrepancies, up to 40%, for cavitated regions with mixed content of air, soft tissue, and bone, such as sinus cavities or tympanic bullae. Conclusions: Our images and findings from a clinically realistic proton CT scanner demonstrate the potential for proton CT to be used for low-dose treatment planning with reduced margins.


title: Universal collective modes from strong electronic correlations: Modified $1/\mathcal{N}_f$ theory with application to high-$T_c$ cuprates

abstract: A nonzero-temperature technique for strongly correlated electron lattice systems, combining elements of both variational wave function (VWF) approach and expansion in the inverse number of fermionic flavors ($1/\mathcal{N}_f$), is developed. The departure point, VWF method, goes beyond the renormalized mean-field theory and provides semi-quantitative description of principal equilibrium properties of high-$T_c$ superconducting cuprates. The developed here scheme of VWF+$1/\mathcal{N}_f$, in the leading order provides dynamical spin and charge responses around the VWF solution, generalizing the weak-coupling spin-fluctuation theory to the regime of strong correlations. Thermodynamic corrections to the correlated saddle-point state arise systematically at consecutive orders. Explicitly, VWF+$1/\mathcal{N}_f$ is applied to evaluate dynamical response functions for the hole-doped Hubbard model and compared with available determinant quantum-Monte-Carlo data, yielding a good overall agreement in the regime of coherent collective-mode dynamics. The emergence of well-defined spin and charge excitations from the incoherent continua is explicitly demonstrated and a non-monotonic dependence of the charge-excitation energy on the interaction magnitude is found. The charge-mode energy saturates slowly when approaching the strong-coupling limit, which calls for a reevaluation of the $t$-$J$-model approach to the charge dynamics in favor of more general $t$-$J$-$U$ and $t$-$J$-$U$-$V$ models. The results are also related to recent inelastic resonant $X$-ray and neutron scattering experiments for the high-$T_c$ cuprates.


title: Entangled symmetric states and copositive matrices

abstract: Entanglement in symmetric quantum states and the theory of copositive matrices are intimately related concepts. For the simplest symmetric states, i.e., the diagonal symmetric (DS) states, it has been shown that there exists a correspondence between exceptional (non-exceptional) copositive matrices and non-decomposable (decomposable) Entanglement Witnesses (EWs). Here we show that EWs of symmetric, but not DS, states can also be constructed from extended copositive matrices, providing new examples of bound entangled symmetric states, together with their corresponding EWs, in arbitrary odd dimensions.


title: Simple model for sequential multiphoton ionization by ultra-intense x-rays

abstract: A simple model for sequential multiphoton ionization by ultra-intense x-rays is presented. The derived scaling of the ion yield with pulse energy quantitatively reproduces the experimental data, which shows that the ion yield increases according to the "power law" behavior typical of multiphoton ionization, followed by saturation at high pulse energies. The calculated average time interval between ionizations for producing ions at a certain charge state is found to be proportional to the pulse duration and independent of all other x-ray pulse parameters. This agrees with previous studies where the kinetic energy of fragment ions with a given charge state produced by intense x-ray ionization of molecules was found to be independent of the pulse energy, but to increase with smaller pulse duration due to the smaller time interval between ionizations.


title: Biodiversity of marine microbes is safeguarded by phenotypic variability in ecological traits

abstract: Why, contrary to theoretical predictions, do marine microbe communities harbor tremendous phenotypic heterogeneity? How can so many marine microbe species competing in the same niche coexist? We discovered a unifying explanation for both phenomena by investigating a non-cooperative game that interpolates between individual-level competitions and species-level outcomes. We identified all equilibrium strategies of the game. These strategies are characterized by maximal phenotypic heterogeneity. They are also neutral towards each other in the sense that an unlimited number of species can co-exist while competing according to the equilibrium strategies. Whereas prior theory predicts that natural selection would minimize trait variation around an optimum value, here we obtained a rigorous mathematical proof that species with maximally variable traits are those that endure. This discrepancy may reflect a disparity between predictions from models developed for larger organisms in contrast to our microbe-centric model. Rigorous mathematics proves that phenotypic heterogeneity is itself a mechanistic underpinning of microbial diversity. This discovery has fundamental ramifications for microbial ecology and may represent an adaptive reservoir sheltering biodiversity in changing environmental conditions.


title: Axion search with BabyIAXO in view of IAXO

abstract: Axions are a natural consequence of the Peccei-Quinn mechanism, the most compelling solution to the strong-CP problem. Similar axion-like particles (ALPs) also appear in a number of possible extensions of the Standard Model, notably in string theories. Both axions and ALPs are very well motivated candidates for Dark Matter, and in addition, they would be copiously produced at the sun's core. A relevant effort during the last decade has been the CAST experiment at CERN, the most sensitive axion helioscope to-date. The International Axion Observatory (IAXO) is a large-scale 4th generation helioscope. As its primary physics goal, IAXO will look for solar axions or ALPs with a signal to background ratio of about 5 orders of magnitude higher than CAST. Recently the IAXO collaboration has proposed and intermediate experimental stage, BabyIAXO, conceived to test all IAXO subsystems (magnet, optics, detectors and sun-tracking systems) at a relevant scale for the final system and thus serve as pathfinder for IAXO but at the same time as a fully-fledged helioscope with record and relevant physics reach in itself with potential for discovery. BabyIAXO was endorsed by the Physics Review committee of DESY last May 2019. Here we will review the status and prospects of BabyIAXO and its potential to probe the most physics motivated regions of the axion & ALPs parameter space.


title: DataVault: A Data Storage Infrastructure for the Einstein Toolkit

abstract: Data sharing is essential in the numerical simulations research. We introduce a data repository, DataVault, that is designed for data sharing, search and analysis. A comparative study of existing repositories is performed to analyze features that are critical to a data repository. We describe the architecture, workflow, and deployment of DataVault, and provide three use-case scenarios for different communities to facilitate the use and application of DataVault. Potential features are proposed and we outline the future development for these features.


title: Smashed products of topological left quasigroups

abstract: In this article smashed products of topological left quasigroups are scrutinized. Quotients of topological fan quasigroups are investigated. Skew smashed products of fan quasigroups and their structure also are studied.


title: Low energy implications of cosmological data in $U(1)_X$ Higgs inflation

abstract: A scalar field having the Coleman-Weinberg type effective potential arises in various contexts of particle physics and serves as a useful framework for discussing cosmic inflation. According to recent studies based on the Markov chain Monte Carlo analysis, the coefficients of such an effective potential are severely constrained by the cosmological data. We investigate the impact of this observation on the physics beyond the Standard Model, focusing on an inflationary model based on the $U(1)_X$-extended Standard Model as a well-motivated example. We examine the parameter region that is not excluded by the Large Hadron Collider (LHC) Run-2 at 139 fb${}^{-1}$ integrated luminosity, and show that the model parameters can be further constrained by the High-Luminosity LHC experiments in the near future. We also comment on the possible reheating mechanism and the dark matter candidates of this scenario.


title: Enhanced high-dispersion coronagraphy with KPIC phase II: design, assembly and status of sub-modules

abstract: The Keck Planet Imager and Characterizer (KPIC) is a purpose-built instrument for high-dispersion coronagraphy in the K and L bands on Keck. This instrument will provide the first high resolution (R$>$30,000) spectra of known directly imaged exoplanets and low-mass brown dwarf companions visible in the northern hemisphere. KPIC is developed in phases. Phase I is currently at Keck in the early operations stage, and the phase II upgrade will deploy in late 2021. The goal of phase II is to maximize the throughput for planet light and minimize the stellar leakage, hence reducing the exposure time needed to acquire spectra with a given signal-to-noise ratio. To achieve this, KPIC phase II exploits several innovative technologies that have not been combined this way before. These include a 1000-element deformable mirror for wavefront correction and speckle control, a set of lossless beam shaping optics to maximize coupling into the fiber, a pupil apodizer to suppress unwanted starlight, a pupil plane vortex mask to enable the acquisition of spectra at and within the diffraction limit, and an atmospheric dispersion compensator. These modules, when combined with the active fiber injection unit present in phase I, will make for a highly efficient exoplanet characterization platform. In this paper, we will present the final design of the optics and opto-mechanics and highlight some innovative solutions we implemented to facilitate all the new capabilities. We will provide an overview of the assembly and laboratory testing of the sub-modules and some of the results. Finally, we will outline the deployment timeline.


title: Octonion random functions and integration of stochastic PDEs

abstract: In the article random functions in modules over the octonion algebra and Cayley-Dickson algebras are investigated. For their study transition measures with values in the octonion algebra and Cayley-Dickson algebras are used. Stochastic integrals over these algebras are studied. They are applied to integration of stochastic PDEs. This approach permits subsequently to analyze and integrate PDEs of orders higher than two of different types including parabolic, elliptic and hyperbolic.


title: Strong and tunable spin-orbit interaction in a single crystalline InSb nanosheet

abstract: A dual-gate InSb nanosheet field-effect device is realized and is used to investigate the physical origin and the controllability of the spin-orbit interaction in a narrow bandgap semiconductor InSb nanosheet. We demonstrate that by applying a voltage over the dual gate, efficiently tuning of the spin-orbit interaction in the InSb nanosheet can be achieved. We also find the presence of an intrinsic spin-orbit interaction in the InSb nanosheet at zero dual-gate voltage and identify its physical origin as a build-in asymmetry in the device layer structure. Having a strong and controllable spin-orbit interaction in an InSb nanosheet could simplify the design and realization of spintronic deceives, spin-based quantum devices and topological quantum devices.


title: Boundary Layer Circumplanetary Accretion: How Fast Could an Unmagnetized Planet Spin Up Through Its Disk?

abstract: Gas giant planets are expected to accrete most of their mass via a circumplanetary disk. If the planet is unmagnetized and initially slowly rotating, it will accrete gas via a radially narrow boundary layer and rapidly spin up. Radial broadening of the boundary layer as the planet spins up reduces the specific angular momentum of accreted gas, allowing the planet to find a terminal rotation rate short of the breakup rate. Here, we use axisymmetric viscous hydrodynamic simulations to quantify the terminal rotation rate of planets accreting from their circumplanetary disks. For an isothermal planet-disk system with a disk scale height $h/r =0.1$ near the planetary surface, spin up switches to spin down at between 70\% and 80\% of the planet's breakup angular velocity. In a qualitative difference from vertically-averaged models -- where spin down can co-exist with mass accretion -- we observe \emph{decretion} accompanying solutions where angular momentum is being lost. The critical spin rate depends upon the disk thickness near the planet. For an isothermal system with a disk scale height of $h/r = 0.15$ near the planet, the critical spin rate drops to between 60\% and 70\% of the planet's breakup angular velocity. In the disk outside the boundary layer, we identify meridional circulation flows, which are unsteady and instantaneously asymmetric across the mid-plane. The simulated flows are strong enough to vertically redistribute solid material in early-stage satellite formation. We discuss how extrasolar planetary rotation measurements, when combined with spectroscopic and variability studies of protoplanets with circumplanetary disks, could determine the role of magnetic and non-magnetic processes in setting giant planet spins.


title: Trefftz Functions for Nonlocal Electrostatics

abstract: Electrostatic interactions in solvents play a major role in biophysical systems. There is a consensus in the literature that the dielectric response of aqueous solutions is nonlocal: polarization depends on the electric field not only at a given point but in the vicinity of that point as well. This is typically modeled via a convolution of the electric field with an appropriate integral kernel. A primary problem with nonlocal models is high computational cost. A secondary problem is restriction of convolution integrals to the solvent, as opposed to their evaluation over the whole space. The paper develops a computational tool alleviating the "curse of nonlocality" and helping to handle the integration correctly. This tool is Trefftz approximations, which tend to furnish much higher accuracy than traditional polynomial ones. In the paper, Trefftz approximations are developed for problems of nonlocal electrostatics, with the goal of numerically "localizing" the original nonlocal problem. This approach can be extended to nonlocal problems in other areas of computational mathematics, physics and engineering.


title: On the radial abundance gradients of nitrogen and oxygen in the inner Galactic disc

abstract: We present optical spectra of nine Galactic H II regions observed with the 10.4 m Gran Telescopio Canarias telescope and located at Galactocentric distances (RG) from 4 to 8 kpc. The distances of the objects have been revised using Gaia DR2 parallaxes. We determine the electron temperature for all the nebulae, which allows a precise computation of their ionic abundances. We have included published data of an additional sample of Galactic H II regions, providing a final data set of 42 objects. The shape of the radial gradients of O/H and N/H is linear and constant, discarding any substantial change of the slope, at least for RG between 4 and 17 kpc. The small dispersion of the O/H and N/H values with respect to the computed gradients implies the absence of significant azimuthal variations of the chemical abundances, at least in the quadrant covered by our observations. We find an almost flat N/O versus O/H diagram relation. This result is not observed in other nearby spiral galaxy except M31. Finally, we compare our computed gradients with those obtained using far infrared (FIR) spectra. We confirm the significant offset in the N/O distribution between the optical and FIR observations. Possible explanations involve ionization correction factors and the strong dependence on density of the abundance determinations based on FIR lines.


title: Regularizing Action Policies for Smooth Control with Reinforcement Learning

abstract: A critical problem with the practical utility of controllers trained with deep Reinforcement Learning (RL) is the notable lack of smoothness in the actions learned by the RL policies. This trend often presents itself in the form of control signal oscillation and can result in poor control, high power consumption, and undue system wear. We introduce Conditioning for Action Policy Smoothness (CAPS), an effective yet intuitive regularization on action policies, which offers consistent improvement in the smoothness of the learned state-to-action mappings of neural network controllers, reflected in the elimination of high-frequency components in the control signal. Tested on a real system, improvements in controller smoothness on a quadrotor drone resulted in an almost 80% reduction in power consumption while consistently training flight-worthy controllers. Project website: http://ai.bu.edu/caps


title: An approximate closed formula for European Mortgage Options

abstract: The aim of this paper is to investigate the use of close formula approximation for pricing European mortgage options. Under the assumption of logistic duration and normal mortgage rates the underlying price at the option expiry is approximated by shifted lognormal or regular lognormal distribution by matching moments. Once the price function is approximated by lognormal distributions, the option price can be computed directly as an integration of the distribution function over the payoff at the option expiry by using Black-Scholes-Merton close formula. We will see that lower curvature levels correspond to positively skewness price distributions and in this case lognormal approximation leads to close parametric formula representation in terms of all model parameters. The proposed methodologies are tested against Monte Carlo approach under different market and contract parameters and the tests confirmed that the close form approximation have a very good accuracy.


title: A fine-grained parallelization of the immersed boundary method

abstract: We present new algorithms for the parallelization of Eulerian-Lagrangian interaction operations in the immersed boundary method. Our algorithms rely on two well-studied parallel primitives: key-value sort and segmented reduce. The use of these parallel primitives allows us to implement our algorithms on both graphics processing units (GPUs) and on other shared memory architectures. We present strong and weak scaling tests on problems involving scattered points and elastic structures. Our tests show that our algorithms exhibit near-ideal scaling on both multicore CPUs and GPUs.


title: Analysis of multi-registry kidney exchange program with individual rationality constraints

abstract: Kidney exchange programs have been developed to overcome the compatibility challenges for the patient with incompatible donors in kidney transplantation. A registry of such incompatible donor-recipient pairs is created, and compatibility is achieved through the exchange of donors. Single-center kidney exchange registries do not have a significant pool size for achieving the full benefits of the kidney exchange program. Thus multi-registry exchange transplants seem to be a natural way forward. There are various challenges in multi-registry exchange programs, such as different constraints for each registry, varying bound on cycle lengths, and data sharing among registries. Researchers have proposed different merging mechanisms with a sequential merger, full merger with fair allocations, and full merger without any restrictions. To form a stable multi-registry sharing mechanism, registries should be individually rational both in the long run as well as in the short run. In this paper, an Integer Programming model for a multi-registry exchange program has been proposed with individual rationality constraints for each registry. A simulation study is conducted to compare the benefits of multi-registry programs with individual registry allocations on Indian data, and multiple analyses have been done. It was concluded that registries with lower arrival rates could benefit 7-9\% in terms of the number of transplants and 4-6\% in terms of quality of matches in a typical kidney exchange registry. Registry with more hard to match patients will benefit significantly larger than registries with more easy to match patients.


title: Large Deviations of Cancer Recurrence Timing

abstract: We study large deviation events in the timing of disease recurrence. In particular, we are interested in modeling cancer treatment failure due to mutation-induced drug resistance. We first present a two-type branching process model of this phenomenon, where an initial population of cells that are sensitive to therapy can produce mutants that are resistant to the therapy. In this model, we investigate two random times, the recurrence time and the crossover time. Recurrence time is defined as the first time that the population size of mutant cells exceeds a given proportion of the initial population size of drug-sensitive cells. Crossover time is defined as the first time that the resistant cell population dominates the total population. We establish convergence in probability results for both recurrence and crossover time. We then develop expressions for the large deviations rate of early recurrence and early crossover events. We characterize how the large deviation rates and rate functions depend on the initial size of the mutant cell population. We finally look at the large deviations rate of early recurrence conditioned on the number of mutant clones present at recurrence in the special case of a deterministically decaying sensitive population. We find that if recurrence occurs before the predicted law of large numbers limit then there will likely be an increase in the number of clones present at recurrence time.


title: Forecasts of the trend in global-mean temperature to 2100 arising from the scenarios of first-difference CO2 and peak fossil fuel

abstract: Two future scenarios that are not explicitly in the range of scenarios (the Representative Concentration Pathway scenarios) utilised by the IPCC. These two scenarios are the emissions trend under peak fossil fuel (for example, Mohr et al., 2015); and the climate sensitivity determinable from the relationship between first-difference CO2 and temperature recently shown by Leggett and Ball (2015). This paper provides forecasts of a global surface temperature trajectory to 2100 resulting from the effect of these two scenarios. The time-series models developed both displayed high statistical significance and converged in their forecasts, so adding to the potential robustness of the findings. Under the effect of the combination of the peak fossil fuel and first-difference CO2 scenarios, we found that temperature is forecast to continue to rise, but only gently, until around 2023 where it reaches a level slightly higher than at present; and from then to display an also gentle steady decrease. It is shown that this trajectory is markedly lower than that generated by the IPCC business-as-usual level-of-CO2 (RCP8.5) model (Riahi et al. 2011). These lower results are evidence that the climate problem may require less future preventative action than is presently being considered necessary. If so, it is stressed that the same evidence is support for the case that the peak fossil fuel problem would require ameliorative action. This globally required action is the same as it would have been for climate (as embodied in the Paris Agreement on Climate Change of the 21st Conference of the Parties of the UNFCCC in Paris adopted in 2015): the rapid transition to a predominantly renewable global energy system.


title: D$^2$IM-Net: Learning Detail Disentangled Implicit Fields from Single Images

abstract: We present the first single-view 3D reconstruction network aimed at recovering geometric details from an input image which encompass both topological shape structures and surface features. Our key idea is to train the network to learn a detail disentangled reconstruction consisting of two functions, one implicit field representing the coarse 3D shape and the other capturing the details. Given an input image, our network, coined D$^2$IM-Net, encodes it into global and local features which are respectively fed into two decoders. The base decoder uses the global features to reconstruct a coarse implicit field, while the detail decoder reconstructs, from the local features, two displacement maps, defined over the front and back sides of the captured object. The final 3D reconstruction is a fusion between the base shape and the displacement maps, with three losses enforcing the recovery of coarse shape, overall structure, and surface details via a novel Laplacian term.


title: Epistemic Logic of Know-Who

abstract: The paper suggests a definition of "know who" as a modality using Grove-Halpern semantics of names. It also introduces a logical system that describes the interplay between modalities "knows who", "knows", and "for all agents". The main technical result is a completeness theorem for the proposed system.


title: Inferring urban social networks from publicly available data

abstract: The emergence of social networks and the definition of suitable generative models for synthetic yet realistic social graphs are widely studied problems in the literature. By not being tied to any real data, random graph models cannot capture all the subtleties of real networks and are inadequate for many practical contexts -- including areas of research, such as computational epidemiology, which are recently high on the agenda. At the same time, the so-called contact networks describe interactions, rather than relationships, and are strongly dependent on the application and on the size and quality of the sample data used to infer them. To fill the gap between these two approaches, we present a data-driven model for urban social networks, implemented and released as open source software. Given a territory of interest, and only based on widely available aggregated demographic and social-mixing data, we construct an age-stratified and geo-referenced synthetic population whose individuals are connected by "strong ties" of two types: intra-household (e.g., kinship) or friendship. While household links are entirely data-driven, we propose a parametric probabilistic model for friendship, based on the assumption that distances and age differences play a role, and that not all individuals are equally sociable. The demographic and geographic factors governing the structure of the obtained network, under different configurations, are thoroughly studied through extensive simulations focused on three Italian cities of different size.


title: Statistical tools for a better optical model

abstract: Background: Modern statistical tools provide the ability to compare the information content of observables and provide a path to explore which experiments would be most useful to give insight into and constrain theoretical models. Purpose: In this work we study three such tools in the context of nuclear reactions with the goal of constraining the optical potential. Method: The three statistical tools examined are: i) the principal component analysis; ii) the sensitivity analysis based on derivatives; and iii) the Bayesian evidence. We first apply these tools to a toy-model case, comparing the form of the imaginary part of the optical potential. Then we consider two different reaction observables, elastic angular distributions and polarization data for reactions on 48Ca at two different beam energies. Results: For the toy-model case, we find significant discrimination power in the sensitivities and the Bayesian evidence, showing clearly that the volume imaginary term is more useful to describe scattering at higher energies. When comparing between elastic cross sections and polarization data using realistic optical models, sensitivity studies indicate that both observables are roughly equally sensitive but the variability of the optical model parameters is strongly angle dependent. The Bayesian evidence shows some variability between the two observables, but the Bayes factor obtained is not sufficient to discriminate between angular distributions and polarization. Conclusions: From the cases considered, we conclude that in general elastic scattering angular distributions have similar impact in constraining the optical potential parameters compared to the polarization data. The angular ranges for the optimum experimental constraints can vary significantly with the observable considered.


title: The extension of Elementary Numerical Operator Theory: Simply Generic, Right-Stochastically Non-Gaussian, Sub-Canonically H-Positive Definite Monoids

abstract: The concept of polytopes was a milestone in elementary integral group theory. In this article, we will extend the concept of contra-linearly stochastic lines to H-Positive Definite Monoids. We will show that by adding the Non-Gaussian assumption, the problem of ordinary Elementary Numerical Operator will have Sub-Canonically properties. Moreover, we will discuss the possibility to classify associative, bijective, and conditionally covariant subalgebras. The prove is given in each section.


title: A multi-period and bi-objective approach for locating ambulances: a case study in Belo Horizonte, Brazil

abstract: This work aims to apply the Facility Location Problem in the Emergency Medical Service (EMS) of Belo Horizonte, Brazil. The objective is to improve two previous optimization models from literature to handle base locations and ambulances allocation/relocation problems. The proposed multi-period models introduce the concept of relocation to the local EMS, which allows ambulances to move among bases in different periods to raise the system coverage. This paper also proposes a bi-objective approach aiming to minimize the number of bases and maximize the coverage of demands, which is solved using the epsilon-constraint method. Results show that coverage levels increase by up to 31% for the deterministic approach and up to 24% for the probabilistic approach. Ambulances' relocation optimization might improve coverage levels by up to 21% for the deterministic approach and up to 15% for the probabilistic approach when change scenarios from static to multi-period. Also, since the multi-period model solutions result in installing a larger number of bases, the bi-objective approach is a powerful tool for the decision-maker. Bi-objective results suggest modest increments in the objective function when the stations' number exceeds 28 for the deterministic approach. The probabilistic approach increments in the objective function start to narrow above 30 installed stations.


title: How to Train your Quadrotor: A Framework for Consistently Smooth and Responsive Flight Control via Reinforcement Learning

abstract: We focus on the problem of reliably training Reinforcement Learning (RL) models (agents) for stable low-level control in embedded systems and test our methods on a high-performance, custom-built quadrotor platform. A common but often under-studied problem in developing RL agents for continuous control is that the control policies developed are not always smooth. This lack of smoothness can be a major problem when learning controllers %intended for deployment on real hardware as it can result in control instability and hardware failure. Issues of noisy control are further accentuated when training RL agents in simulation due to simulators ultimately being imperfect representations of reality - what is known as the reality gap. To combat issues of instability in RL agents, we propose a systematic framework, `REinforcement-based transferable Agents through Learning' (RE+AL), for designing simulated training environments which preserve the quality of trained agents when transferred to real platforms. RE+AL is an evolution of the Neuroflight infrastructure detailed in technical reports prepared by members of our research group. Neuroflight is a state-of-the-art framework for training RL agents for low-level attitude control. RE+AL improves and completes Neuroflight by solving a number of important limitations that hindered the deployment of Neuroflight to real hardware. We benchmark RE+AL on the NF1 racing quadrotor developed as part of Neuroflight. We demonstrate that RE+AL significantly mitigates the previously observed issues of smoothness in RL agents. Additionally, RE+AL is shown to consistently train agents that are flight-capable and with minimal degradation in controller quality upon transfer. RE+AL agents also learn to perform better than a tuned PID controller, with better tracking errors, smoother control and reduced power consumption.


title: A Simulation Study to Evaluate the Performance of the Cauchy Proximal Operator in Despeckling SAR Images of the Sea Surface

abstract: The analysis of ocean surface is widely performed using synthetic aperture radar (SAR) imagery as it yields information for wide areas under challenging weather conditions, during day or night, etc. Speckle noise constitutes however the main reason for reduced performance in applications such as classification, ship detection, target tracking and so on. This paper presents an investigation into the despeckling of SAR images of the ocean that include ship wake structures, via sparse regularisation using the Cauchy proximal operator. We propose a closed-form expression for calculating the proximal operator for the Cauchy prior, which makes it applicable in generic proximal splitting algorithms. In our experiments, we simulate SAR images of moving vessels and their wakes. The performance of the proposed method is evaluated in comparison to the L1 and TV norm regularisation functions. The results show a superior performance of the proposed method for all the utilised images generated.


title: Hiding in the Particles: When Return-Oriented Programming Meets Program Obfuscation

abstract: Largely known for attack scenarios, code reuse techniques at a closer look reveal properties that are appealing also for program obfuscation. We explore the popular return-oriented programming paradigm under this light, transforming program functions into ROP chains that coexist seamlessly with the surrounding software stack. We show how to build chains that can withstand popular static and dynamic deobfuscation approaches, evaluating the robustness and overheads of the design over common programs. The results suggest a significant amount of computational resources would be required to carry a deobfuscation attack for secret finding and code coverage goals.


title: DeCoAR 2.0: Deep Contextualized Acoustic Representations with Vector Quantization

abstract: Recent success in speech representation learning enables a new way to leverage unlabeled data to train speech recognition model. In speech representation learning, a large amount of unlabeled data is used in a self-supervised manner to learn a feature representation. Then a smaller amount of labeled data is used to train a downstream ASR system using the new feature representations. Based on our previous work DeCoAR and inspirations from other speech representation learning, we propose DeCoAR 2.0, a Deep Contextualized Acoustic Representation with vector quantization. We introduce several modifications over the DeCoAR: first, we use Transformers in encoding module instead of LSTMs; second, we introduce a vector quantization layer between encoder and reconstruction modules; third, we propose an objective that combines the reconstructive loss with vector quantization diversity loss to train speech representations. Our experiments show consistent improvements over other speech representations in different data-sparse scenarios. Without fine-tuning, a light-weight ASR model trained on 10 hours of LibriSpeech labeled data with DeCoAR 2.0 features outperforms the model trained on the full 960-hour dataset with filterbank features.


title: Understanding Spectral Graph Neural Network

abstract: Graph neural networks have developed by leaps and bounds in recent years due to the restriction of traditional convolutional filters on non-Euclidean structured data. Spectral graph theory mainly studies fundamental graph properties using algebraic methods to analyze the spectrum of the adjacency matrix or Laplacian matrix of a graph, which lays the foundation of graph convolutional neural networks. This report is more than notes and self-contained which comes from my Ph.D. first-year report literature review part, it illustrates how the graph convolutional neural network model is motivated by spectral graph theory, and discusses the major spectral-based models associated with their fundamentals. The practical applications of the graph convolutional neural networks defined in the spectral domain are also reviewed.


title: Lie automorphisms of incidence algebras

abstract: Let $X$ be a finite connected poset and $K$ a field. We give a full description of the Lie automorphisms of the incidence algebra $I(X,K)$. In particular, we show that they are in general not proper.


title: Protective Policy Transfer

abstract: Being able to transfer existing skills to new situations is a key capability when training robots to operate in unpredictable real-world environments. A successful transfer algorithm should not only minimize the number of samples that the robot needs to collect in the new environment, but also prevent the robot from damaging itself or the surrounding environment during the transfer process. In this work, we introduce a policy transfer algorithm for adapting robot motor skills to novel scenarios while minimizing serious failures. Our algorithm trains two control policies in the training environment: a task policy that is optimized to complete the task of interest, and a protective policy that is dedicated to keep the robot from unsafe events (e.g. falling to the ground). To decide which policy to use during execution, we learn a safety estimator model in the training environment that estimates a continuous safety level of the robot. When used with a set of thresholds, the safety estimator becomes a classifier for switching between the protective policy and the task policy. We evaluate our approach on four simulated robot locomotion problems and a 2D navigation problem and show that our method can achieve successful transfer to notably different environments while taking the robot's safety into consideration.


title: Exploiting the Dual-Tree Complex Wavelet Transform for Ship Wake Detection in SAR Imagery

abstract: In this paper, we analyse synthetic aperture radar (SAR) images of the sea surface using an inverse problem formulation whereby Radon domain information is enhanced in order to accurately detect ship wakes. This is achieved by promoting linear features in the images. For the inverse problem-solving stage, we propose a penalty function, which combines the dual-tree complex wavelet transform (DT-CWT) with the non-convex Cauchy penalty function. The solution to this inverse problem is based on the forward-backward (FB) splitting algorithm to obtain enhanced images in the Radon domain. The proposed method achieves the best results and leads to significant improvement in terms of various performance metrics, compared to state-of-the-art ship wake detection methods. The accuracy of detecting ship wakes in SAR images with different frequency bands and spatial resolution reaches more than 90%, which clearly demonstrates an accuracy gain of 7% compared to the second-best approach.


title: Data-driven subgrid-scale modeling of forced Burgers turbulence using deep learning with generalization to higher Reynolds numbers via transfer learning

abstract: Developing data-driven subgrid-scale (SGS) models for large eddy simulations (LES) has received substantial attention recently. Despite some success, particularly in a priori (offline) tests, challenges have been identified that include numerical instabilities in a posteriori (online) tests and generalization (i.e., extrapolation) of trained data-driven SGS models, for example to higher Reynolds numbers. Here, using the stochastically forced Burgers turbulence as the test-bed, we show that deep neural networks trained using properly pre-conditioned (augmented) data yield stable and accurate a posteriori LES models. Furthermore, we show that transfer learning enables accurate/stable generalization to a flow with 10x higher Reynolds number.


title: Spin contribution to the perihelion precession in binary systems like OJ287: higher order corrections

abstract: Higher order corrections are obtained for the perihelion precession in binary systems like OJ287, Sagittarius A*-S2 and H1821+643 using both the Schwarzschild metric and the Kerr metric to take into account the spin effect. The corrections are performed considering the third root of the motion equation and developing the expansion in terms of $\epsilon \equiv r_s/\left(a(1-e^2)\right)$ and $\epsilon^{*} \equiv \left(1- \frac{2 \alpha E'}{cJ}\right) \epsilon $.The results are compared with those obtained in a previous paper.


title: Cooperative Location Privacy in Vehicular Networks: Why Simple Mix-zones are not Enough

abstract: Vehicular communications disclose rich information about the vehicles and their whereabouts. Pseudonymous authentication secures communication while enhancing user privacy. To enhance location privacy, cryptographic mix-zones were proposed to facilitate vehicles covertly transition to new ephemeral credentials. The resilience to (syntactic and semantic) pseudonym linking (attacks) highly depends on the geometry of the mix-zones, mobility patterns, vehicle density, and arrival rates. We introduce a tracking algorithm for linking pseudonyms before and after a cryptographically protected mix-zone. Our experimental results show that an eavesdropper, leveraging standardized vehicular communication messages and road layout, could successfully link 73% of pseudonyms during non-rush hours and 62% of pseudonyms during rush hours after vehicles change their pseudonyms in a mix-zone. To mitigate such inference attacks, we present a novel cooperative mix-zone scheme that enhances user privacy regardless of the vehicle mobility patterns, vehicle density, and arrival rate to the mix-zone. A subset of vehicles, termed relaying vehicles, are selected to be responsible for emulating non-existing vehicles. Such vehicles cooperatively disseminate decoy traffic without affecting safety-critical operations: with 50% of vehicles as relaying vehicles, the probability of linking pseudonyms (for the entire interval) drops from 68% to 18%. On average, this imposes 28 ms extra computation overhead, per second, on the Roadside Units (RSUs) and 4.67 ms extra computation overhead, per second, on the (relaying) vehicle side; it also introduces 1.46 KB/sec extra communication overhead by (relaying) vehicles and 45 KB/sec by RSUs for the dissemination of decoy traffic. Thus, user privacy is enhanced at the cost of low computation and communication overhead.


title: Avoiding The Double Descent Phenomenon of Random Feature Models Using Hybrid Regularization

abstract: We demonstrate the ability of hybrid regularization methods to automatically avoid the double descent phenomenon arising in the training of random feature models (RFM). The hallmark feature of the double descent phenomenon is a spike in the regularization gap at the interpolation threshold, i.e. when the number of features in the RFM equals the number of training samples. To close this gap, the hybrid method considered in our paper combines the respective strengths of the two most common forms of regularization: early stopping and weight decay. The scheme does not require hyperparameter tuning as it automatically selects the stopping iteration and weight decay hyperparameter by using generalized cross-validation (GCV). This also avoids the necessity of a dedicated validation set. While the benefits of hybrid methods have been well-documented for ill-posed inverse problems, our work presents the first use case in machine learning. To expose the need for regularization and motivate hybrid methods, we perform detailed numerical experiments inspired by image classification. In those examples, the hybrid scheme successfully avoids the double descent phenomenon and yields RFMs whose generalization is comparable with classical regularization approaches whose hyperparameters are tuned optimally using the test data. We provide our MATLAB codes for implementing the numerical experiments in this paper at https://github.com/EmoryMLIP/HybridRFM.


title: Selection Games on Hyperspaces

abstract: In this paper we connect selection principles on a topological space to corresponding selection principles on one of its hyperspaces. We unify techniques and generalize theorems from the known results about selection principles for common hyperspace constructions. This includes results of Lj.D.R. Ko\v{c}inac, Z. Li, and others. We use selection games to generalize selection principles and we work with strategies of various strengths for these games. The selection games we work with are primarily abstract versions of the selection principles of Rothberger, Menger, and Hurewicz type, as well as games of countable fan tightness and selective separability. The hyperspace constructions that we work with are the Vietoris and Fell topologies, both upper and full, generated by ideals of closed sets. Using a new technique we are able to extend straightforward connections between topological constructs to connections between selection games related to those constructs. This extension process works regardless of the length of the game, the kind of selection being performed, or the strength of the strategy being considered.


title: Kinetic Modeling of Magnetospheres

abstract: This paper presents the state of the art of kinetic modeling techniques for simulating plasma kinetic dynamics in magnetospheres. We describe the critical numerical techniques for enabling large-scale kinetic simulations of magnetospheres: parameter scaling, implicit Particle-in-Cell schemes, and fluid-kinetic coupling. We show an application of these techniques to study particle acceleration and heating in asymmetric magnetic reconnection in the Ganymede magnetosphere.


title: Adaptive Histogram-Based Gradient Boosted Trees for Federated Learning

abstract: Federated Learning (FL) is an approach to collaboratively train a model across multiple parties without sharing data between parties or an aggregator. It is used both in the consumer domain to protect personal data as well as in enterprise settings, where dealing with data domicile regulation and the pragmatics of data silos are the main drivers. While gradient boosted tree implementations such as XGBoost have been very successful for many use cases, its federated learning adaptations tend to be very slow due to using cryptographic and privacy methods and have not experienced widespread use. We propose the Party-Adaptive XGBoost (PAX) for federated learning, a novel implementation of gradient boosting which utilizes a party adaptive histogram aggregation method, without the need for data encryption. It constructs a surrogate representation of the data distribution for finding splits of the decision tree. Our experimental results demonstrate strong model performance, especially on non-IID distributions, and significantly faster training run-time across different data sets than existing federated implementations. This approach makes the use of gradient boosted trees practical in enterprise federated learning.


title: Detection of chaotic behavior in time series

abstract: Deterministic chaos is phenomenon from nonlinear dynamics and it belongs to greatest advances of twentieth-century science. Chaotic behavior appears apart of mathematical equations also in wide range in observable nature, so as in there originating time series. Chaos in time series resembles stochastic behavior, but apart of randomness it is totally deterministic and therefore chaotic data can provide us useful information. Therefore it is essential to have methods, which are able to detect chaos in time series, moreover to distinguish chaotic data from stochastic one. Here we present and discuss the performance of standard and machine learning methods for chaos detection and its implementation on two well known simple chaotic discrete dynamical systems - Logistic map and Tent map, which fit to the most of the definitions of chaos.


title: Exoplanets through extreme optics: from PLATO to SHARK-NIR

abstract: In the last years, the Observatory of Padova (Istituto Nazionale di Astrofisica - Osservatorio Astronomico di Padova) and the University of Padova have been involved massively in projects dedicated to the exoplanets search, both ground, and space-based. The activities concerning my Ph.D. have been exploited both in the framework of the space projects and in the field of ground instrumentation. PLATO, the acronym of "PLAnetary Transits and Oscillations of stars", is a ESA mission o with the target to detect and characterize exoplanets utilizing their transit on a bright star. The overall instrumental layout consists of a multi-telescope concept instrument, composed by several tens of telescope units, for which it has developed an all-refractive optical solution. These devices are characterized by a very large Field of View (>20 degrees on one side) with an optical quality that fits most of the energy into a single sensor pixel. I participated in the Assembly, Integration, and Verification (AIV) of the Telescope Optical Unit prototype, to validate the AIV procedure and the telescope optical performance in-flight conditions. SHARK-NIR, the acronym of "System for coronagraphy with High order Adaptive optics from R to K band - Near-Infrared", is an instrument designed to search and characterize the young exo-planetary system and star-forming regions in the NIR domain, in coronagraphic direct imaging and spectroscopic mode. It has been selected for the 2nd generation Large Binocular Telescope (LBT) instruments, and it will take advantage of the excellent performance of the LBT eXtreme Adaptive Optics (XAO) correction, necessary for SHARK-NIR to achieve the best possible coronagraphic performance, which is mandatory to detect faint planets orbiting around bright stars. Concerning SHARK, my activity has been performing optical alignment and qualification of the instrument.


title: Ruin Probabilities for a Sparre Andersen Model with Investments

abstract: We study a Sparre Andersen model in which the business activity of the company is described by a compound renewal process with drift assuming that the capital reserves are invested in a risky asset. The price of the latter is assumed to evolve according to a geometric L\'evy process. We prove that the asymptotic behavior of the ruin probability depends to a large extent only on the properties of the price process.


title: High throughput spatially sensitive single-shot quantitative phase microscopy

abstract: High space-bandwidth product with high spatial phase sensitivity is indispensable for a single-shot quantitative phase microscopy (QPM) system. It opens avenue for widespread applications of QPM in the field of biomedical imaging. Temporally low coherence length light sources are generally implemented to achieve high spatial phase sensitivity in QPM at the cost of either reduced temporal resolution or smaller field of view (FOV). On the contrary, high temporal coherence light sources like lasers are capable of exploiting the full FOV of the QPM systems at the expense of less spatial phase sensitivity. In the present work, we employed pseudo-thermal light source (PTLS) in QPM which overcomes the limitations of conventional light sources. The capabilities of PTLS over conventional light sources are systematically studied and demonstrated on various test objects like USAF resolution chart and thin optical waveguide (height ~ 8 nm). The spatial phase sensitivity of QPM in case of PTLS is measured to be equivalent to that for white light source. The high-speed and large FOV capabilities of PTLS based QPM is demonstrated by high-speed imaging of live sperm cells that is limited by the camera speed and by imaging extra-ordinary large FOV phase imaging on histopathology placenta tissue samples.


title: Clustered Sparse Channel Estimation for Massive MIMO Systems by Expectation Maximization-Propagation (EM-EP)

abstract: We study the problem of downlink channel estimation in multi-user massive multiple input multiple output (MIMO) systems. To this end, we consider a Bayesian compressive sensing approach in which the clustered sparse structure of the channel in the angular domain is employed to reduce the pilot overhead. To capture the clustered structure, we employ a conditionally independent identically distributed Bernoulli-Gaussian prior on the sparse vector representing the channel, and a Markov prior on its support vector. An expectation propagation (EP) algorithm is developed to approximate the intractable joint distribution on the sparse vector and its support with a distribution from an exponential family. The approximate distribution is then used for direct estimation of the channel. The EP algorithm assumes that the model parameters are known a priori. Since these parameters are unknown, we estimate these parameters using the expectation maximization (EM) algorithm. The combination of EM and EP referred to as EM-EP algorithm is reminiscent of the variational EM approach. Simulation results show that the proposed EM-EP algorithm outperforms several recently-proposed algorithms in the literature.


title: A new approach to the Dyson rank conjectures

abstract: In 1944 Dyson defined the rank of a partition as the largest part minus the number of parts, and conjectured that the residue of the rank mod 5 divides the partitions of 5n+4 into five equal classes. This gave a combinatorial explanation of Ramanujan's famous partition congruence mod 5. He made an analogous conjecture for the rank mod 7 and the partitions of 7n+5. In 1954 Atkin and Swinnerton-Dyer proved Dyson's rank conjectures by constructing several Lambert-series identities basically using the theory of elliptic functions. In 2016 the author gave another proof using the theory of weak harmonic Maass forms. In this paper we describe a new and more elementary approach using Hecke-Rogers series.


title: Focus optimization in a Computational Confocal Microscope

abstract: We consider the numerical optimization of performance for a computational extension of a confocal microscope. Using a system where the pinhole detector is replaced with a detector array, we seek to exploit this additional information for each point in the scan. We derive an optimal estimate of the light at focus which minimizes the contribution of out-of-focus light. We estimate the amount of improvement that would be theoretically possible in point-scanning and line-scanning systems and demonstrate with simulation. We find that even with a large degree of regularization, a significant improvement is possible, especially for line-scanning systems.


title: TabTransformer: Tabular Data Modeling Using Contextual Embeddings

abstract: We propose TabTransformer, a novel deep tabular data modeling architecture for supervised and semi-supervised learning. The TabTransformer is built upon self-attention based Transformers. The Transformer layers transform the embeddings of categorical features into robust contextual embeddings to achieve higher prediction accuracy. Through extensive experiments on fifteen publicly available datasets, we show that the TabTransformer outperforms the state-of-the-art deep learning methods for tabular data by at least 1.0% on mean AUC, and matches the performance of tree-based ensemble models. Furthermore, we demonstrate that the contextual embeddings learned from TabTransformer are highly robust against both missing and noisy data features, and provide better interpretability. Lastly, for the semi-supervised setting we develop an unsupervised pre-training procedure to learn data-driven contextual embeddings, resulting in an average 2.1% AUC lift over the state-of-the-art methods.


title: Convolutional LSTM Neural Networks for Modeling Wildland Fire Dynamics

abstract: As the climate changes, the severity of wildland fires is expected to worsen. Models that accurately capture fire propagation dynamics greatly help efforts for understanding, responding to and mitigating the damages caused by these fires. Machine learning techniques provide a potential approach for developing such models. The objective of this study is to evaluate the feasibility of using a Convolutional Long Short-Term Memory (ConvLSTM) recurrent neural network to model the dynamics of wildland fire propagation. The machine learning model is trained on simulated wildfire data generated by a mathematical analogue model. Three simulated datasets are analyzed, each with increasing degrees of complexity. The simplest dataset includes a constant wind direction as a single confounding factor, whereas the most complex dataset includes dynamic wind, complex terrain, spatially varying moisture content and heterogenous vegetation density distributions. We examine how effective the ConvLSTM can learn the fire-spread dynamics over consecutive time steps. It is shown that ConvLSTMs can capture local fire transmission events, as well as the overall fire dynamics, such as the rate at which the fire spreads. Finally, we demonstrate that ConvLSTMs outperform other network architectures that have previously been used to model similar wildland fire dynamics.


title: The physical and the geometrical properties of simulated cold HI structures

abstract: The objective of this paper is to help shedding some light on the nature and the properties of the cold structures formed via thermal instability in the magnetized atomic interstellar medium. To this end, we searched for clumps formed in forced (magneto)hydrodynamic simulations with an initial magnetic field ranging from 0 to 8.3$\mu$G. We statistically analyzed, through the use of Kernel Density Estimations, the physical and the morphological properties of a sample containing $\sim 1500$ clumps, as well as the relative alignments between the main direction of clumps and the internal velocity and magnetic field. The density ($n\sim 50-200$cm$^{-3}$), the thermal pressure ($P_{th}/k\sim 4.9\times 10^3-10^4$K cm$^{-3}$), the mean magnetic field ($\sim 3-11$$\mu$G ), and the sonic Mach number of the selected clumps have values comparable to those reported in observations. We find, however, that the cloud sample can not be described by a single regime concerning their pressure balance and their Alf\'enic Mach number. We measured the morphological properties of clumps mainly through the asphericity and the prolatness, which appear to be more sensitive than the aspect ratios. From this analysis we find that the presence of magnetic field, even if it is weak, does qualitatively affect the morphology of the clumps by increasing the probability of having highly aspherical and highly plolate clumps by a factor of two, that is by producing more filamentary clumps. Finally, we find that the angle between the main direction of the clumps and the local magnetic field lies between $\sim\pi/4-\pi/2$ and shifts to more perpendicular alignments as the intensity of this field increases, while the relative direction between the local density structure and the local magnetic field transits from parallel to perpendicular.


title: Linear embedding of nonlinear dynamical systems and prospects for efficient quantum algorithms

abstract: The simulation of large nonlinear dynamical systems, including systems generated by discretization of hyperbolic partial differential equations, can be computationally demanding. Such systems are important in both fluid and kinetic computational plasma physics. This motivates exploring whether a future error-corrected quantum computer could perform these simulations more efficiently than any classical computer. We describe a method for mapping any finite nonlinear dynamical system to an infinite linear dynamical system (embedding) and detail three specific cases of this method that correspond to previously-studied mappings. Then we explore an approach for approximating the resulting infinite linear system with finite linear systems (truncation). Using a number of qubits only logarithmic in the number of variables of the nonlinear system, a quantum computer could simulate truncated systems to approximate output quantities if the nonlinearity is sufficiently weak. Other aspects of the computational efficiency of the three detailed embedding strategies are also discussed.


title: Mind the Gap: Cake Cutting With Separation

abstract: We study the problem of fairly allocating a divisible resource, also known as cake cutting, with an additional requirement that the shares that different agents receive should be sufficiently separated from one another. This captures, for example, constraints arising from social distancing guidelines. While it is sometimes impossible to allocate a proportional share to every agent under the separation requirement, we show that the well-known criterion of maximin share fairness can always be attained. We then provide algorithmic analysis of maximin share fairness in this setting -- for instance, the maximin share of an agent cannot be computed exactly by any finite algorithm, but can be approximated with an arbitrarily small error. In addition, we consider the division of a pie (i.e., a circular cake) and show that an ordinal relaxation of maximin share fairness can be achieved. We also prove that an envy-free or equitable allocation that allocates the maximum amount of resource exists under separation.


title: Cortex: Harnessing Correlations to Boost Query Performance

abstract: Databases employ indexes to filter out irrelevant records, which reduces scan overhead and speeds up query execution. However, this optimization is only available to queries that filter on the indexed attribute. To extend these speedups to queries on other attributes, database systems have turned to secondary and multi-dimensional indexes. Unfortunately, these approaches are restrictive: secondary indexes have a large memory footprint and can only speed up queries that access a small number of records, and multi-dimensional indexes cannot scale to more than a handful of columns. We present Cortex, an approach that takes advantage of correlations to extend the reach of primary indexes to more attributes. Unlike prior work, Cortex can adapt itself to any existing primary index, whether single or multi-dimensional, to harness a broad variety of correlations, such as those that exist between more than two attributes or have a large number of outliers. We demonstrate that on real datasets exhibiting these diverse types of correlations, Cortex matches or outperforms traditional secondary indexes with $5\times$ less space, and it is $2-8\times$ faster than existing approaches to indexing correlations.


title: Faster Policy Learning with Continuous-Time Gradients

abstract: We study the estimation of policy gradients for continuous-time systems with known dynamics. By reframing policy learning in continuous-time, we show that it is possible construct a more efficient and accurate gradient estimator. The standard back-propagation through time estimator (BPTT) computes exact gradients for a crude discretization of the continuous-time system. In contrast, we approximate continuous-time gradients in the original system. With the explicit goal of estimating continuous-time gradients, we are able to discretize adaptively and construct a more efficient policy gradient estimator which we call the Continuous-Time Policy Gradient (CTPG). We show that replacing BPTT policy gradients with more efficient CTPG estimates results in faster and more robust learning in a variety of control tasks and simulators.


title: Coordinated Frequency and Voltage Regulation of Grid-Following and Grid-Forming Inverters

abstract: In a purely inverter-based microgrid, both grid-forming (GFM) and grid-following (GFL) inverters will have a crucial role to play in frequency/voltage regulation and maintaining power sharing through their grid support capabilities. Consequently, the coordination between these two technologies becomes increasingly important for optimal system performance. However, the existing work does not consider GFL's potential to participate in a secondary control in coordination with GFM, thus not able to utilize the full capability of inverter resources. In this paper, we show that it is possible to fully coordinate the GFL and GFM inverters to achieve accurate power sharing, frequency/voltage regulation, and circulating var mitigation in networked microgrids even without the support of any synchronous generators or the bulk power system. We use the leader-follower consensus framework to develop a GFM-GFL coordination control. The effectiveness of the proposed coordination is verified under different disturbances and communication degradation. In addition, we find that the proposed fully-coordinated secondary control outperforms other approaches such as un-coordinated and partially-coordinated secondary controls, in aspects of load sharing and frequency and voltage regulation. Overall, this study emphasizes the need and benefits of GFL-GFM coordination in microgrids.


title: Computing Machinery and Knowledge

abstract: The purpose of this paper is to discuss the possibilities for computing machinery, or AI agents, to know and to possess knowledge. This is done mainly from a virtue epistemology perspective and definition of knowledge. However, this inquiry also shed light on the human condition, what it means for a human to know, and to possess knowledge. The paper argues that it is possible for an AI agent to know and examines this from both current state-of-the-art in artificial intelligence as well as from the perspective of what the future AI development might bring in terms of superintelligent AI agents.


title: Proposal for room-temperature quantum repeaters with nitrogen-vacancy centers and optomechanics

abstract: We propose a quantum repeater architecture that can operate under ambient conditions. Our proposal builds on recent progress towards non-cryogenic spin-photon interfaces based on nitrogen-vacancy centers, which have excellent spin coherence times even at room temperature, and optomechanics, which allows to avoid phonon-related decoherence and also allows the emitted photons to be in the telecom band. We apply the photon number decomposition method to quantify the fidelity and the efficiency of entanglement established between two remote electron spins. We describe how the entanglement can be stored in nuclear spins and extended to long distances via quasi-deterministic entanglement swapping operations involving the electron and nuclear spins. We furthermore propose schemes to achieve high-fidelity readout of the spin states at room temperature using the spin-optomechanics interface. Our work shows that long-distance quantum networks made of solid-state components that operate at room temperature are within reach of current technological capabilities.


title: Compact Two-Dimensional Magneto-Optical Trap for Ultracold Atom Setups

abstract: We report on the design, implementation, and performance of a compact two-dimensional magneto-optical trap (2D MOT) for cesium. In a small-volume vacuum chamber, the setup uses cesium dispensers in close proximity to the trapping region of the 2D MOT and operates at low vapor pressures in the $10^{-9}$ torr range. We achieve a cold atom flux of $4 \times 10^8$ atoms/s that is comparable to the performance of more complex atomic sources. The setup is simple to construct and can be adapted to a broad range of atomic species.


title: Flow-plane decorrelations in heavy-ion collisions with multiple-plane cumulants

abstract: The azimuthal correlations between local flow planes at different (pseudo)rapidities ($\eta$) may reveal important details of the initial nuclear matter density distributions in heavy-ion collisions. Extensive experimental measurements of a factorization ratio ($r_2$) and its derivative ($F_2$) have shown evidence of the longitudinal flow-plane decorrelation. However, nonflow effects also affect this observable and prevent a quantitative understanding of the phenomenon. In this paper, to distinguish decorrelation and nonflow effects, we propose a new cumulant observable, $T_2$, which largely suppresses nonflow. The technique sensitivity to different initial-state scenarios and nonflow effects are tested with a simple Monte Carlo model, and in the end, the method is applied to events simulated by a multiphase transport model (AMPT) for Au+Au collisions at $\sqrt{s_{\rm NN}} =200$ GeV. We also emphasize that a distinct decorrelation signal requires not only the right sign of an observable, but also its proper dependence on the $\eta$-window of the reference flow plane, to be consistent with the pertinent decorrelation picture.


title: Yelp Review Rating Prediction: Machine Learning and Deep Learning Models

abstract: We predict restaurant ratings from Yelp reviews based on Yelp Open Dataset. Data distribution is presented, and one balanced training dataset is built. Two vectorizers are experimented for feature engineering. Four machine learning models including Naive Bayes, Logistic Regression, Random Forest, and Linear Support Vector Machine are implemented. Four transformer-based models containing BERT, DistilBERT, RoBERTa, and XLNet are also applied. Accuracy, weighted F1 score, and confusion matrix are used for model evaluation. XLNet achieves 70% accuracy for 5-star classification compared with Logistic Regression with 64% accuracy.


title: Parameter Estimation with Dense and Convolutional Neural Networks Applied to the FitzHugh-Nagumo ODE

abstract: Machine learning algorithms have been successfully used to approximate nonlinear maps under weak assumptions on the structure and properties of the maps. We present deep neural networks using dense and convolutional layers to solve an inverse problem, where we seek to estimate parameters of a FitzHugh-Nagumo model, which consists of a nonlinear system of ordinary differential equations (ODEs). We employ the neural networks to approximate reconstruction maps for model parameter estimation from observational data, where the data comes from the solution of the ODE and takes the form of a time series representing dynamically spiking membrane potential of a biological neuron. We target this dynamical model because of the computational challenges it poses in an inference setting, namely, having a highly nonlinear and nonconvex data misfit term and permitting only weakly informative priors on parameters. These challenges cause traditional optimization to fail and alternative algorithms to exhibit large computational costs. We quantify the prediction errors of model parameters obtained from the neural networks and investigate the effects of network architectures with and without the presence of noise in observational data. We generalize our framework for neural network-based reconstruction maps to simultaneously estimate ODE parameters and parameters of autocorrelated observational noise. Our results demonstrate that deep neural networks have the potential to estimate parameters in dynamical models and stochastic processes, and they are capable of predicting parameters accurately for the FitzHugh-Nagumo model.


title: Mathematical modeling the wildfire propagation in a Randers space

abstract: The devastating effects of wildfires on the wildlife and their impact on human lives and properties are undeniable. This shows the importance of studying the spread of wildfire, predicting its behavior and presenting more reliable models for its propagation. Here, by using the validity of the Huygens` envelope principle for wavefronts in Randers spaces, we present some models for the propagation of wildfire in an n-dimensional smooth manifold under the presence of wind. In the models, trajectories of fire particles are tracked and the equations that give the wavefront at each time are provided. Furthermore, we determine the paths and points of great importance in the process of wildfire management, called strategic paths and points. Finally, we consider two examples of spreading the wildfire in some agricultural land or woodland, for the sake of illustration.