title: Fast magnetic reconnection and the ideal evolution of a magnetic field

abstract: Regardless of how small non-ideal effects may be, phenomena associated with changes in magnetic field line connections are frequently observed to occur on an Alfv\'enic time scale. Since it is mathematically impossible for magnetic field line connections to change when non-ideal effects are identically zero, an ideal evolution must naturally lead to states of unbounded sensitivity to non-ideal effects. That such an evolution is natural is demonstrated using Lagrangian coordinates based on the flow velocity of the magnetic field lines. The Lagrangian representation of an evolving magnetic field is highly constrained when neither the magnetic field strength nor the forces exerted by the magnetic field increase exponentially with time. The development of a state of fast reconnection consistent with these constraints (1) requires a three-dimensional evolution, (2) has an exponentially increasing sensitivity to non-ideal effects, and (3) has a parallel current density, which lies in exponentially thinning but exponentially widening ribbons, with a magnitude that is limited to a slow growth. The implication is that exponential growth in sensitivity is the cause of fast magnetic reconnection when non-ideal effects are sufficiently small. The growth of the non-ideal effect of the resistivity multiplied by the parallel current density is far too slow to be competitive.


title: Comparison between super-hydrophobic, liquid infused and rough surfaces: a DNS study

abstract: Direct Numerical Simulations of two superposed fluids in a channel with a textured surface on the lower wall have been carried out. A parametric study varying the viscosity ratio between the two fluids has been performed to mimic both {\bf idealised} super hydrophobic and liquid infused surfaces and assess its effect on the frictional, form and total drag for three different textured geometries: longitudinal square bars, transversal square bars and staggered cubes. The interface between the two fluids is assumed to be slippery in the streamwise and spanwise directions and not deformable in the vertical direction, corresponding to the ideal case of infinite surface tension. To identify the role of the fluid-fluid interface, an extra set of simulations with a single fluid has been carried out and compared to the results obtained with two fluids of same viscosity separated by the interface. The drag and the maximum wall-normal velocity fluctuations were found to be highly correlated for all the surface configurations, whether they reduce or increase the drag. This implies that the structure of the near-wall turbulence is dominated by the total shear and not by the local boundary condition of super-hydrophobic, liquid--infused or rough surfaces.


title: KIC 10001893: A pulsating sdB star with multiple trapped modes

abstract: KIC 10001893 is a V1093 Her type pulsating subdwarf-B star, which was observed extensively by the Kepler spacecraft. It was a part of the survey searching for compact pulsators in the Kepler field of view. An amplitude spectrum of the star demonstrates a rich content of g-modes between 102 and 496 {\mu}Hz as well as a few p-modes above 2000 {\mu}Hz. In total, we found 110 frequencies. The g-mode region contains 104 frequencies, while the p-mode region contains just six, altogether indicating the hybrid nature of KIC 10001893. The main goal of our analysis was to identify the detected modes and to find some features, which will significantly help modeling of the star. We found no multiplets, which points at a pole-on orientation, however, we defined modal degrees and relative radial orders using asymptotic period spacing. Eventually, we assigned 32 dipole l = 1 and 18 quadrupole l = 2 modes. The most remarkable feature we found are trapped modes, which are clearly seen in a reduce period diagram. It is the first time that three trapped modes are detected in one pulsating sdB star. Since the more trapped modes we find, the better sounding of the stellar interior we can do, this feature provides important constraints on the physical structure of the star. Mode trapping is likely caused by the He-H transition region and therefore it provides crucial constraints for making realistic theoretical models of hot subdwarfs.


title: A Probe Towards Understanding GAN and VAE Models

abstract: This project report compares some known GAN and VAE models proposed prior to 2017. There has been significant progress after we finished this report. We upload this report as an introduction to generative models and provide some personal interpretations supported by empirical evidence. Both generative adversarial network models and variational autoencoders have been widely used to approximate probability distributions of data sets. Although they both use parametrized distributions to approximate the underlying data distribution, whose exact inference is intractable, their behaviors are very different. We summarize our experiment results that compare these two categories of models in terms of fidelity and mode collapse. We provide a hypothesis to explain their different behaviors and propose a new model based on this hypothesis. We further tested our proposed model on MNIST dataset and CelebA dataset.


title: Hund's metal regimes and orbital selective Mott transitions in three band systems

abstract: We analyze the electronic properties of interacting crystal field split three band systems. Using a rotationally invariant slave boson approach we analyze the behavior of the electronic mass renormalization as a function of the intralevel repulsion $U$, the Hund's coupling $J$, the crystal field splitting, and the number of electrons per site $n$. We first focus on the case in which two of the bands are identical and the levels of the third one are shifted by $\Delta>0$ with respect to the former. We find an increasing quasiparticle mass differentiation between the bands, for system away from half-filling ($n=3$), as the Hubbard interaction $U$ is increased. This leads to orbital selective Mott transitions where either the higher energy band (for $4>n>3$) or the lower energy degenerate bands ($2<n<3$) become insulating for $U$ larger than a critical interaction $U_{c}(n)$. Away from the half-filled case $|n-3|\gtrsim 0.3$ there is a wide range of parameters for $U<U_c(n)$ where the system presents a Hund's metal phase with the physics dominated by the local high spin multiplets. Finally, we study the fate of the $n=2$ Hund's metal as the energy splitting between orbitals is increased for different possible crystal distortions. We find a strong sensitivity of the Hund's metal regime to crystal fields due to the opposing effects of $J$ and the crystal field splittings on the charge distribution between the bands.


title: Split Regression Modeling

abstract: Sparse methods are the standard approach to obtain interpretable models with high prediction accuracy. Alternatively, algorithmic ensemble methods can achieve higher prediction accuracy at the cost of loss of interpretability. However, the use of blackbox methods has been heavily criticized for high-stakes decisions and it has been argued that there does not have to be a trade-off between accuracy and interpretability. To combine high accuracy with interpretability, we generalize best subset selection to best split selection. Best split selection constructs a small number of sparse models learned jointly from the data which are then combined in an ensemble. Best split selection determines the models by splitting the available predictor variables among the different models when fitting the data. The proposed methodology results in an ensemble of sparse and diverse models that each provide a possible explanation for the relationship between the predictors and the response. The high computational cost of best split selection motivates the need for computational tractable approximations. We evaluate a method developed by Christidis et al. (2020) which can be seen as a multi-convex relaxation of best split selection.


title: Rejecting cosmic background for exclusive neutrino interaction studies with Liquid Argon TPCs; a case study with the MicroBooNE detector

abstract: Cosmic ray (CR) interactions can be a challenging source of background for neutrino oscillation and cross-section measurements in surface detectors. We present methods for CR rejection in measurements of charged-current quasielastic-like (CCQE-like) neutrino interactions, with a muon and a proton in the final state, measured using liquid argon time projection chambers (LArTPCs). Using a sample of cosmic data collected with the MicroBooNE detector, mixed with simulated neutrino scattering events, a set of event selection criteria is developed that produces an event sample with minimal contribution from CR background. Depending on the selection criteria used a purity between 50% and 80% can be achieved with a signal selection efficiency between 50% and 25%, with higher purity coming at the expense of lower efficiency. While using a specific dataset from the MicroBooNE detector and selection criteria values optimized for CCQE-like events, the concepts presented here are generic and can be adapted for various studies of exclusive {\nu}{\mu} interactions in LArTPCs.

categories: physics.ins-det hep-ex nucl-ex physics.data-an


title: Periodic codings of Bratteli-Vershik systems

abstract: We develop conditions for the coding of a Bratteli-Vershik system according to initial path segments to be periodic, equivalently for a constructive symbolic recursive scheme corresponding to a cutting and stacking process to produce a periodic sequence. This is a step toward understanding when a Bratteli-Vershik system can be essentially faithfully represented by means of a natural coding as a subshift on a finite alphabet.


title: Relaxing LHC constraints on the $W_R$ mass

abstract: We study mass bounds of the $W_R$ gauge boson in generic left-right symmetric models. Assuming that the gauge bosons couple universally to quarks and leptons, we allow different gauge couplings $g_R \ne g_L$ and mass mixing, $V_{CKM}^L \ne V_{CKM}^R$ in the left and right sectors. Imposing constraints from collider experiments and $K^0$, $B_d$, $B_s$ physics, we investigate scenarios where $W_R$ is lighter, or heavier than the right handed neutrino $\nu_R$. In these scenarios, $W_R$ mass bounds can be considerably relaxed, while $Z_R$ mass bounds are much more stringent. In the case where $M_{W_R} \le M_{\nu_R}$, the experimental constraints come from $W_R \to tb $ and $W_R \to jj$ channels, while if $M_{W_R} \ge M_{\nu_R}$, the dominant constraints come from $W_R \to \ell \ell jj $. The observed (expected) limits in the two-dimensional ($M_{W_R}$, $M_{\nu_R}$) mass plane excluded at 95\% confidence level extend to approximately $M_{W_R}$= 3.1 (3.3) TeV in the $ee$ channel and 3.3 (3.4) TeV in the ($\mu\mu$) channel, for a large range of right-handed neutrino masses up to $M_{\nu_R}$= 2.1 (2.1) TeV in the $ee$ channel and 2.6 (2.5) in the ($\mu\mu$) channel, representing a significant relaxation of the mass bounds.


title: All-Sky Measurement of the Anisotropy of Cosmic Rays at 10 TeV and Mapping of the Local Interstellar Magnetic Field

abstract: We present the first full-sky analysis of the cosmic ray arrival direction distribution with data collected by the HAWC and IceCube observatories in the Northern and Southern hemispheres at the same median primary particle energy of 10 TeV. The combined sky map and angular power spectrum largely eliminate biases that result from partial sky coverage and holds a key to probe into the propagation properties of TeV cosmic rays through our local interstellar medium and the interaction between the interstellar and heliospheric magnetic fields. From the map we determine the horizontal dipole components of the anisotropy $\delta_{0h} = 9.16 \times 10^{-4}$ and $\delta_{6h} = 7.25 \times 10^{-4}~(\pm0.04 \times 10^{-4})$. In addition, we infer the direction ($229.2\pm 3.5^\circ$ RA , $11.4\pm 3.0^\circ$ Dec.) of the interstellar magnetic field from the boundary between large scale excess and deficit regions from which we estimate the missing corresponding vertical dipole component of the large scale anisotropy to be $\delta_N \sim -3.97 ^{+1.0}_{-2.0} \times 10^{-4}$.


title: Crosstalk Impacts on Homogeneous Weakly-Coupled Multicore Fiber Based IM/DD System

abstract: We numerically discussed crosstalk impacts on homogeneous weakly-coupled multicore fiber based intensity modulation/direct-detection (IM/DD) systems taking into account mean crosstalk power fluctuation, walk-off between cores, laser frequency offset, and laser linewidth.


title: Solutions to Diophantine Equation of Erdos-Straus Conjecture

abstract: In number theory, the Erdos-Straus conjecture states that for all n >=2, the rational number 4/n can be expressed as the sum of three unit fractions. Paul Erdos and Ernst G. Straus formulated the conjecture in 1948. The restriction that the three unit fractions be positive is essential to the difficulty of the problem, for if negative values were allowed the problem could always be solved. This paper presents an explicit solutions to this conjecture for all n >=2 excepting some n such that n=1(mod8).


title: Derived category of projectivization and generalized linear duality

abstract: In this note, we generalize the linear duality between vector subbundles (or equivalently quotient bundles) of dual vector bundles to coherent quotients $V \twoheadrightarrow \mathscr{L}$ considered in arXiv:1811.12525, in the framework of Kuznetsov's homological projective duality (HPD). As an application, we obtain a generalized version of the fundamental theorem of HPD for the $\mathbb{P}(\mathscr{L})$-sections and the respective dual sections of a given HPD pair.


title: \texttt{GooStats}: A GPU-based framework for multi-variate analysis in particle physics

abstract: \texttt{GooStats} is a software framework that provides a flexible environment and common tools to implement multi-variate statistical analysis. The framework is built upon the \texttt{CERN ROOT}, \texttt{MINUIT} and \texttt{GooFit} packages. Running a multi-variate analysis in parallel on graphics processing units yields a huge boost in performance and opens new possibilities. The design and benchmark of \texttt{GooStats} are presented in this article along with illustration of its application to statistical problems.


title: Ablation of a Robot's Brain: Neural Networks Under a Knife

abstract: It is still not fully understood exactly how neural networks are able to solve the complex tasks that have recently pushed AI research forward. We present a novel method for determining how information is structured inside a neural network. Using ablation (a neuroscience technique for cutting away parts of a brain to determine their function), we approach several neural network architectures from a biological perspective. Through an analysis of this method's results, we examine important similarities between biological and artificial neural networks to search for the implicit knowledge locked away in the network's weights.


title: Measurements of UHECR Mass Composition by Telescope Array

abstract: Telescope Array (TA) has recently published results of nearly nine years of $X_{\mathrm{max}}$ observations providing its highest statistics measurement of ultra high energy cosmic ray (UHECR) mass composition to date for energies exceeding $10^{18.2}$ eV. This analysis measured agreement of observed data with results expected for four different single elements. Instead of relying only on the first and second moments of $X_{\mathrm{max}}$ distributions, we employ a morphological test of agreement between data and Monte Carlo to allow for systematic uncertainties in data and in current UHECR hadronic models. Results of this latest analysis and implications of UHECR composition observed by TA are presented. TA can utilize different analysis methods to understand composition as both a crosscheck on results and as a tool to understand systematics affecting $X_{\mathrm{max}}$ measurements. The different analysis efforts underway at TA to understand composition are also discussed.


title: Scaling of intrinsic domain wall magneto-resistance with confinement in electromigrated nanocontacts

abstract: In this work we study the evolution of intrinsic domain wall magnetoresistance (DWMR) with domain wall confinement. Clean permalloy notched half-ring nanocontacts are fabricated using a special ultra-high vacuum electromigration procedure to tailor the size of the wire in-situ and through the resulting domain wall confinement we tailor the domain wall width from a few tens of nm down to a few nm. Through measurements of the dependence of the resistance with respect to the applied field direction we extract the contribution of a single domain wall to the MR of the device, as a function of the domain wall width in the confining potential at the notch. In this size range, an intrinsic positive MR is found, which dominates over anisotropic MR, as confirmed by comparison to micromagnetic simulations. Moreover, the MR is found to scale monotonically with the size of the domain wall, $\delta_{DW}$, as 1/$\delta_{DW}^b$, with $b=2.31\pm 0.39 $. The experimental result is supported by quantum-mechanical transport simulations based on ab-initio density functional theory calculations.


title: Vortex lattices and broken time reversal symmetry in the topological superconductor UPt3

abstract: The topological superconductor UPt3, has three distinct vortex phases, a strong indication of its unconventional character. Using small-angle neutron scattering we have probed the vortex lattice in the UPt3 B phase with the magnetic field along the crystal c-axis. We find a difference in the vortex lattice configuration depending on the sign of the magnetic field relative to the field direction established upon entering the B phase at low temperature in a field sweep, showing that the vortices in this material posses an internal degree of freedom. This observation is facilitated by the discovery of a field driven non-monotonic vortex lattice rotation, driven by competing effects of the superconducting gap distortion and the vortex-core structure. From our bulk measurements we infer that the superconducting order parameter in the UPt3 B phase breaks time reversal symmetry and exhibits chiral symmetry with respect to the c-axis.


title: Dose-response modeling in high-throughput cancer drug screenings: An end-to-end approach

abstract: Personalized cancer treatments based on the molecular profile of a patient's tumor are an emerging and exciting class of treatments in oncology. As genomic tumor profiling is becoming more common, targeted treatments to specific molecular alterations are gaining traction. To discover new potential therapeutics that may apply to broad classes of tumors matching some molecular pattern, experimentalists and pharmacologists rely on high-throughput, in-vitro screens of many compounds against many different cell lines. We propose a hierarchical Bayesian model of how cancer cell lines respond to drugs in these experiments and develop a method for fitting the model to real-world high-throughput screening data. Through a case study, the model is shown to capture nontrivial associations between molecular features and drug response, such as requiring both wild type TP53 and overexpression of MDM2 to be sensitive to Nutlin-3(a). In quantitative benchmarks, the model outperforms a standard approach in biology, with ~20% lower predictive error on held out data. When combined with a conditional randomization testing procedure, the model discovers biomarkers of therapeutic response that recapitulate known biology and suggest new avenues for investigation. All code for the paper is publicly available at https://github.com/tansey/deep-dose-response.


title: Bayesian Sparsification of Gated Recurrent Neural Networks

abstract: Bayesian methods have been successfully applied to sparsify weights of neural networks and to remove structure units from the networks, e. g. neurons. We apply and further develop this approach for gated recurrent architectures. Specifically, in addition to sparsification of individual weights and neurons, we propose to sparsify preactivations of gates and information flow in LSTM. It makes some gates and information flow components constant, speeds up forward pass and improves compression. Moreover, the resulting structure of gate sparsity is interpretable and depends on the task. Code is available on github: https://github.com/tipt0p/SparseBayesianRNN


title: The influence of lattice termination on the edge states of the quantum spin Hall insulator monolayer $1T'$-WTe$_2$

abstract: We study the influence of sample termination on the electronic properties of the novel quantum spin Hall insulator monolayer $1T'$-WTe$_2$. For this purpose, we construct an accurate, minimal 4-orbital tight-binding model with spin-orbit coupling by employing a combination of density-functional theory calculations, symmetry considerations, and fitting to experimental data. Based on this model, we compute energy bands and 2-terminal conductance spectra for various ribbon geometries with different terminations, with and without magnetic field. Because of the strong electron-hole asymmetry we find that the edge Dirac point is buried in the bulk bands for most edge terminations. In the presence of a magnetic field, an in-gap edge Dirac point leads to exponential suppression of conductance as an edge Zeeman gap opens, whereas the conductance stays at the quantized value when the Dirac point is buried in the bulk bands. Finally, we find that disorder in the edge termination drastically changes this picture: the conductance of a sufficiently rough edge is uniformly suppressed for all energies in the bulk gap regardless of the orientation of the edge.


title: Background Discrimination for Neutrinoless Double Beta Decay in Liquid Xenon Using Cherenkov Light

abstract: Neutrinoless double beta decays in liquid xenon produce a significant amount of Cherenkov light, with a photon number and angular distribution that distinguishes these events from common backgrounds. A GEANT4 simulation was used to simulate Cherenkov photon production and measurement in a liquid xenon detector and a multilayer perceptron was used to analyze the resulting distributions to classify events based on their Cherenkov photons. Our results show that a modest improvement in the sensitivity of neutrinoless double beta decay searches is possible using this technique, but the kinematics of the neutrinoless double beta decay and electron scattering in liquid xenon substantially limit this approach.


title: A Finite Strain Constitutive Model Considering Transformation Induced Plasticity for Shape Memory Alloys under Cyclic Loading

abstract: Many engineering applications of Shape Memory Alloys (SMAs) involve passing back and forth through phase transformation many times. Repeated phase transformation develops permanent deformations originating from the significant distortion that phase transformation induces at the austenite-martensite interfaces and grain boundaries. This distortion drives dislocation activity resulting in an observable macroscopic Transformation-Induced Plastic (TRIP) deformation, which occurs at effective stress levels much lower than the plastic yield limit of the material. TRIP strains may accumulate up to 20% during the lifetime of an SMA component and thus a finite strain constitutive model is required for simulating their response. In this work, a 3-D finite strain model is developed based on logarithmic strain which is the only strain measure whose rate is equal to the stretching, that describes phase transformation and TRIP deformation in a thermodynamically consistent setting. The model is implemented in Abaqus finite element software through a user-defined material subroutine (UMAT). Boundary value problems such as strip and a torque tube under both pseudoelastic and actuation cyclic loadings are performed to test the capabilities of the newly proposed model.


title: A Finite Strain Constitutive Model for Martensitic Transformation in Shape Memory Alloys Based on Logarithmic Strain

abstract: Shape Memory Alloys (SMAs) are materials with the ability to recover apparently permanent deformation under specific thermomechanical loading. The majority of constitutive models for SMAs are developed based on the infinitesimal strain theory. However, such assumption may not be proper in the presence of geometric discontinuities, such as cracks, and repeated cycling loading that has been reported to induce irrecoverable strains up to 20% due to transformation induced plasticity. In addition to finite strains, SMA-based devices may also undergo large rotations. Thus, it is indispensable to develop a constitutive model based on the finite strain to provide accurate predictions of these actuators response. A three-dimensional phenomenological constitutive model for SMAs considering finite strains and finite rotations is proposed in this work. This model utilizes the logarithmic strain as the strain measure that is the strain measure whose logarithmic rate in a corotating material frame is equal to the rate of deformation tensor. In the proposed model, the martensitic volume fraction and the second-order logarithmic transformation strain tensor are chosen as the internal state variables associated with the inelastic transformation process. Numerical simulations considering basic SMAs component geometries such as a bar, a beam, and a torque tube are performed to test the capabilities of the proposed model under both mechanically and thermally induced phase transformation. The presented model formulation will be extended in future work for the incorporation of transformation-induced plasticity.


title: Higher Moment Estimation for Elliptically-distributed Data: Is it Necessary to Use a Sledgehammer to Crack an Egg?

abstract: Multivariate elliptically-contoured distributions are widely used for modeling economic and financial data. We study the problem of estimating moment parameters of a semi-parametric elliptical model in a high-dimensional setting. Such estimators are useful for financial data analysis and quadratic discriminant analysis. For low-dimensional elliptical models, efficient moment estimators can be obtained by plugging in an estimate of the precision matrix. Natural generalizations of the plug-in estimator to high-dimensional settings perform unsatisfactorily, due to estimating a large precision matrix. Do we really need a sledgehammer to crack an egg? Fortunately, we discover that moment parameters can be efficiently estimated without estimating the precision matrix in high-dimension. We propose a marginal aggregation estimator (MAE) for moment parameters. The MAE only requires estimating the diagonal of covariance matrix and is convenient to implement. With mild sparsity on the covariance structure, we prove that the asymptotic variance of MAE is the same as the ideal plug-in estimator which knows the true precision matrix, so MAE is asymptotically efficient. We also extend MAE to a block-wise aggregation estimator (BAE) when estimates of diagonal blocks of covariance matrix are available. The performance of our methods is validated by extensive simulations and an application to financial returns.


title: Self-propelled particle in a nonconvex external potential: Persistent limit in one dimension

abstract: Equilibrium mapping techniques for nonaligning self-propelled particles have made it possible to predict the density profile of an active ideal gas in a wide variety of external potentials, however they fail when the self-propulsion is very persistent and the potential is nonconvex, which is precisely when the most uniquely active phenomena occur. Here we show how to predict the density profile of a 1D active Ornstein-Uhlenbeck particle in an arbitrary external potential in the persistent limit and discuss the consequences of the potential's nonconvexity on the structure of the solution, including the central role of the potential's inflection points and the nonlocal dependence of the density profile on the potential.


title: KeV Scale Frozen-in Self-Interacting Fermionic Dark Matter

abstract: We present a model in which the dark matter particle is frozen-in at MeV scale. In this model, the mediator between the standard model sector and the dark sector can automatically provide a self-interaction for dark matter. The interaction strength is naturally to be the in the region in favor of the cluster mass deficit anomaly. Due to the self-scattering, the Lyman-$\alpha$ constraint can be relaxed to $m_D \gtrsim 2 $ keV. In this region the self-interaction and the Fermi pressure both play roles on forming a dark matter core at the center of the dwarf galaxies.


title: A Three-dimensional Constitutive Model for Polycrystalline Shape Memory Alloys Under Large Strains Combined With Large Rotations

abstract: Shape Memory Alloys (SMAs) has been widely aware of working as actuators for active/smart morphing structures in the engineering industry. Because of the high actuation energy density of SMAs, compared to other active materials, structures integrated with SMA-based actuators has high advantage in terms of trade-offs between overall structure weight, integrity and functionality. The majority of available constitutive models for SMAs are developed within the infinitesimal strain regime. However, it was reported that particular SMAs can generate transformation strains nearly up to 8%-10%, for which the adopted infinitesimal strain assumption is no longer appropriate. Furthermore, industry applications may require SMA actuators undergo large rotation deformation at work. Combining the above two facts, a constitutive model for SMAs developed on a finite deformation framework is required. A three-dimensional constitutive model for SMAs considering large strains with large rotations is proposed in this work. This model utilizes the logarithmic strain as a finite strain measure for large deformation analysis so that its rate form hypo-elastic constitutive relation can be consistently integrated to deliver a free energy based hyper-elastic constitutive relation. The martensitic volume fraction and the second-order transformation strain tensor are chosen as the internal state variables to characterize the inelastic response exhibited by polycrystalline SMAs. The presented formulation and its numerical implementation scheme can be extended in future work for the incorporation of other inelastic phenomena such as transformation-induced plasticity, viscoplasticity and creep under large deformations.


title: The maximal current carried by a normal-superconducting interface in the absence of magnetic field

abstract: Modeling a normal/superconducting interface, we consider a semi-infinite wire whose edge is adjacent to a normal magnetic metal, assuming asymptotic convergence, away from the boundary, to the purely superconducting state. We obtain that the maximal current which can be carried by the interface diminishes in the small normal conductivity limit.


title: Search for power-efficient wide-range reversible resistance modulation of $VO_2$ single crystals

abstract: The abrupt metal insulator transition in $VO_2$ is attracting considerable interest from both fundamental and applicative angles. We report on DC I-V characteristics measured on $VO_2$ single crystals in the two-probe configuration at several ambient temperatures below the insulator-metal transition. The insulator-mixed-metal-insulator transition is induced by Joule heating above ambient temperature in the range of negative differential resistivity (NDR). In this range the stability of V(I) is governed by the load resistance $R_L$. Steady state I(V) is obtained for $R_L> |dV/dI|_{max}$ in the NDR regime. For $R_L< |dV/dI|_{max}$ there is switching between initial and final steady states associated with peaks in the Joule power, that are higher the lower $R_L$ is. The peaks caused by steep switching are superfluous and damaging the samples. On the other hand, the large $R_L$ needed for steady state is the main power consumer in the circuit at high currents. The present work is motivated by the need to avoid damaging switching in the NDR regime while reducing the power consumption in the circuit. It is shown here that large resistance modulation can be obtained under steady state conditions with reduced power consumption by increasing the ambient temperature of the device above room temperature.


title: Gyrokinetic continuum simulations of plasma turbulence in the Texas Helimak

abstract: The first gyrokinetic simulations of plasma turbulence in the Texas Helimak device, a simple magnetized torus, are presented. The device has features similar to the scrape-off layer region of tokamaks, such as bad-curvature-driven instabilities and sheath boundary conditions on the end plates, which are included in these simulations. Comparisons between simulations and measurements from the experiment show similarities, including equilibrium profiles and fluctuation amplitudes that approach experimental values, but also some important quantitative differences. Both experimental and simulation results exhibit turbulence statistics that are characteristic of blob transport.


title: On gravitational lensing by symmetric and asymmetric wormholes

abstract: We discuss the peculiarities of gravitational lensing by spherically symmetric wormholes if they are not symmetric with respect to their throats. It is noticed, in particular, that wormholes always contain the so-called photon spheres, near which the photon deflection angles can be arbitrarily large, but, in general, the throat is such a sphere only for symmetric wormholes. In some cases, photons from outside can cross the throat and return back from a neighborhood of a photon sphere if the latter is located beyond the throat. Two families of generally asymmetric wormhole configurations are considered as examples: (1) anti-Fisher wormholes with a massless phantom scalar field as a source of gravity, and (2) wormholes with a zero Ricci scalar that may be interpreted as vacuum configurations in a brane world. For these metrics, the photon effective potentials and deflection angles are found and discussed.


title: Exploring Embedding Methods in Binary Hyperdimensional Computing: A Case Study for Motor-Imagery based Brain-Computer Interfaces

abstract: Key properties of brain-inspired hyperdimensional (HD) computing make it a prime candidate for energy-efficient and fast learning in biosignal processing. The main challenge is however to formulate embedding methods that map biosignal measures to a binary HD space. In this paper, we explore variety of such embedding methods and examine them with a challenging application of motor imagery brain-computer interface (MI-BCI) from electroencephalography (EEG) recordings. We explore embedding methods including random projections, quantization based thermometer and Gray coding, and learning HD representations using end-to-end training. All these methods, differing in complexity, aim to represent EEG signals in binary HD space, e.g. with 10,000 bits. This leads to development of a set of HD learning and classification methods that can be selectively chosen (or configured) based on accuracy and/or computational complexity requirements of a given task. We compare them with state-of-the-art linear support vector machine (SVM) on an NVIDIA TX2 board using the 4-class BCI competition IV-2a dataset as well as a new 3-class dataset. Compared to SVM, results on 3-class dataset show that simple thermometer embedding achieves moderate average accuracy (79.56% vs. 82.67%) with 26.8$\times$ faster training time and 22.3$\times$ lower energy; on the other hand, switching to end-to-end training with learned HD representations wipes out these training benefits while boosting the accuracy to 84.22% (1.55% higher than SVM). Similar trend is observed on the 4-class dataset where SVM achieves on average 74.29%: the thermometer embedding achieves 89.9$\times$ faster training time and 58.7$\times$ lower energy, but a lower accuracy (67.09%) than the learned representation of 72.54%.


title: Bias-Variance Trade-off and Model Selection for Proton Radius Extractions

abstract: Intuitively, a scientist might assume that a more complex regression model will necessarily yield a better predictive model of experimental data. Herein, we disprove this notion in the context of extracting the proton charge radius from charge form factor data. Using a Monte Carlo study, we show that a simpler regression model can in certain cases be the better predictive model. This is especially true with noisy data where the complex model will fit the noise instead of the physical signal. Thus, in order to select the appropriate regression model to employ, a clear technique should be used such as the Akaike information criterion or Bayesian information criterion, and ideally selected previous to seeing the results. Also, to ensure a reasonable fit, the scientist should also make regression quality plots, such as residual plots, and not just rely on a single criterion such as reduced chi2. When we apply these techniques to low four-momentum transfer cross section data, we find a proton radius that is consistent with the muonic Lamb shift results. While presented for the case of proton radius extraction, these concepts are applicable in general and can be used to illustrate the necessity of balancing bias and variance when building a regression model and validating results, ideas that are at the heart of modern machine learning algorithms.


title: The polylog quotient and the Goncharov quotient in computational Chabauty-Kim theory I

abstract: Polylogarithms are those multiple polylogarithms that factor through a certain quotient of the de Rham fundamental group of the thrice punctured line known as the polylogarithmic quotient. Building on work of Dan-Cohen, Wewers, and Brown, we push the computational boundary of our explicit motivic version of Kim's method in the case of the thrice punctured line over an open subscheme of $\mathrm{Spec}\,\mathbb{Z}$. To do so, we develop a greatly refined version of the algorithm of Dan-Cohen tailored specifically to this case, and we focus attention on the polylogarithmic quotient. This allows us to restrict our calculus with motivic iterated integrals to the so-called depth-one part of the mixed Tate Galois group studied extensively by Goncharov. We also discover an interesting consequence of the symmetry-breaking nature of the polylog quotient that forces us to symmetrize our polylogarithmic version of Kim's conjecture. In this first part of a two-part series, we focus on a specific example, which allows us to verify an interesting new case of Kim's conjecture.


title: A Fourier Frame for the Middle-Third Cantor Measure

abstract: In this paper we show that if $\mu$ is any locally and uniformly $\alpha$-dimensional measure supported on a $\alpha$-quasi-regular set $E$, then $L^2(\mu)$ admits a frame of exponentials. In particular, for the uniform middle third Cantor measure, $\mu_C,$ our result shows that there exists a countable set $\Lambda$ such that $\{e^{2\pi i t \lambda}\}_{\lambda \in \Lambda}$ is a frame for $L^2(\mu_C)$ (i.e. the measure $\mu_C$ admits a generalized spectrum), answering an old outstanding question about the existence of a frame of exponentials for the space $L^2(\mu_C)$.


title: Inequality in Death from Social Conflicts: A Gini & Kolkata indices-based Study

abstract: Human deaths caused by individual man-made conflicts (e.g., wars, armed-conflicts, terrorist-attacks etc.) occur unequally across the events (conflicts) and such inequality (in deaths) have been studied here using Lorenz curve and values of the inequality indices Gini ($g$) and Kolkata ($k$) have been estimated from it. The data are taken from various well-known databases maintained by some Universities and Peace Research Institutes. The inequality measures for man-made conflicts are found to have very high values ($g$ = $0.82 ~\pm~ $0.02, $k$ = $0.84~ \pm~ $0.02), which is rarely seen in economic (income or wealth) inequality measures across the world ($g \leq 0.4$, $k \leq 0.6$; presumably because of various welfare measures). We also investigated the inequalities in human deaths from natural disasters (like earthquakes, floods, etc.). Interestingly, we observe that the social inequality measures ($g$ and $k$ values) from man-made conflicts compare well with those of academic centers (inequality in citations; found in earlier studies) of different institutions of the world, while those for natural disasters can be even higher. We discuss about the `similarity classes' of social inequality (similar higher values of $g$ and $k$ indices) for man-made competitive societies like academic institutions and man-made social conflicts, and connect our observations with that of the growing recent trend of economic inequality across the world (with rapid disappearance of welfare strategies).


title: FPETS : Fully Parallel End-to-End Text-to-Speech System

abstract: End-to-end Text-to-speech (TTS) system can greatly improve the quality of synthesised speech. But it usually suffers form high time latency due to its auto-regressive structure. And the synthesised speech may also suffer from some error modes, e.g. repeated words, mispronunciations, and skipped words. In this paper, we propose a novel non-autoregressive, fully parallel end-to-end TTS system (FPETS). It utilizes a new alignment model and the recently proposed U-shape convolutional structure, UFANS. Different from RNN, UFANS can capture long term information in a fully parallel manner. Trainable position encoding and two-step training strategy are used for learning better alignments. Experimental results show FPETS utilizes the power of parallel computation and reaches a significant speed up of inference compared with state-of-the-art end-to-end TTS systems. More specifically, FPETS is 600X faster than Tacotron2, 50X faster than DCTTS and 10X faster than Deep Voice3. And FPETS can generates audios with equal or better quality and fewer errors comparing with other system. As far as we know, FPETS is the first end-to-end TTS system which is fully parallel.


title: Two Relaxation Time Lattice Boltzmann Method Coupled to Fast Fourier Transform Poisson Solver: Application to Electroconvective Flow

abstract: Electroconvective flow between two infinitely long parallel electrodes is investigated via a multiphysics computational model. The model solves for spatiotemporal flow properties using two-relaxation-time Lattice Boltzmann Method for fluid and charge transport coupled to Fast Fourier Transport Poisson solver for the electric potential. The segregated model agrees with the previous analytical and numerical results providing a robust approach for modeling electrohydrodynamic flows.


title: All-2D Material Inkjet-Printed Capacitors: Towards Fully-Printed Integrated Circuits

abstract: A well-defined insulating layer is of primary importance in the fabrication of passive (e.g. capacitors) and active (e.g. transistors) components in integrated circuits. One of the most widely known 2-Dimensional (2D) dielectric materials is hexagonal boron nitride (hBN). Solution-based techniques are cost-effective and allow simple methods to be used for device fabrication. In particular, inkjet printing is a low-cost, non-contact approach, which also allows for device design flexibility, produces no material wastage and offers compatibility with almost any surface of interest, including flexible substrates. In this work we use water-based and biocompatible graphene and hBN inks to fabricate all-2D material and inkjet-printed capacitors. We demonstrate an areal capacitance of 2.0 \pm 0.3 nF cm^(-2) for a dielectric thickness of \sim 3 \mu m and negligible leakage currents, averaged across more than 100 devices. This gives rise to a derived dielectric constant of 6.1 \pm 1.7. The inkjet printed hBN dielectric has a breakdown field of 1.9 \pm 0.3 MV cm^(-1). Fully printed capacitors with sub-/mu m hBN layer thicknesses have also been demonstrated. The capacitors are then exploited in two fully printed demonstrators: a resistor-capacitor (RC) low-pass filter and a graphene-based field effect transistor.


title: Application of the method of multiple scales to unravel energy exchange in nonlinear locally resonant metamaterials

abstract: In this paper, the effect of weak nonlinearities in 1D locally resonant metamaterials is investigated via the method of multiple scales. Commonly employed to the investigate the effect of weakly nonlinear interactions on the free wave propagation through a phononic structure or on the dynamic response of a Duffing oscillator, the method of multiple scales is here used to investigate the forced wave propagation through locally resonant metamaterials. The perturbation approach reveals that energy exchange may occur between propagative and evanescent waves induced by quadratic nonlinear local interaction.


title: Helium load on W-O coatings grown by pulsed laser deposition

abstract: W-O deposits with complex morphologies and significant He contents will growth on the surface of plasma facing components exposed to He discharges. To mimic the re/co-deposition process, W-O-He coatings were produced by implanting He+ ions on W films grown by pulsed laser deposition (PLD). The use of appropriate PLD experimental parameters such as Ar or He background atmospheres induces the deposition of porous or amorphous-like W-O structures, respectively. After multiple ion implantation stages using 150 keV, 100 keV and 50 keV incident He+ ion beams with a total fluence of 5 x 1017 ion/cm^2, significant amounts of He were identified in porous coatings by Rutherford backscattering (RBS). Time-of-flight elastic recoil detection (TOF-ERDA) measurements showed that most of the implanted He was already released from the porous coatings five month after implantation while for the case of amorphous ones the He content remains significant at deeper layers and smoothly decrease towards the surface, as result of a different morphology and nanostructure. The proposed method involving PLD and ion implantation seems adequate to produce W-O-He reference samples for fusion investigations.


title: Optimal error estimates for analytic continuation in the upper half-plane

abstract: Analytic functions in the Hardy class $H^2$ over the upper half-plane $\mathbb{H}_+$ are uniquely determined by their values on any curve $\Gamma$ lying in the interior or on the boundary of $\mathbb{H}_+$. The goal of this paper is to provide a quantitative version of this statement. Given that $f$ from a unit ball in $H^2$ is small on $\Gamma$ (say, its $L^2$ norm is of order $\epsilon$), how does this affect the magnitude of $f$ at a point $z$ away from the curve? When $\Gamma \subset \partial \mathbb{H}_+$, we give a sharp upper bound on $|f(z)|$ of the form $\epsilon^\gamma$, with an explicit exponent $\gamma=\gamma(z) \in (0,1)$ and describe the maximizer function attaining the upper bound. When $\Gamma \subset \mathbb{H}_+$ we give an implicit sharp upper bound in terms of a solution of an integral equation on $\Gamma$. We conjecture and give evidence that this bound also behaves like $\epsilon^\gamma$ for some $\gamma=\gamma(z) \in (0,1)$. These results can also be transplanted to other domains conformally equivalent to the upper half-plane.


title: On the 1/f spectrum in the solar wind and its connection with magnetic compressibility

abstract: We discuss properties of Alfv\'enic fluctuations with large amplitude in plasmas characterised by low magnetic field compression. We note that in such systems power laws can not develop with arbitrarily steep slopes at large scales, i.e. when $|\delta \bf{B}|$ becomes of the order of the background field $|\bf{B}|$. In such systems there is a scale $l_0$ at which the spectrum has to break due to the condition of weak compressibility. A very good example of this dynamics is offered by solar wind fluctuations in Alfv\'enic fast streams, characterised by the property of constant field magnitude. We show here that the distribution of $\delta B=|\delta \bf{B}|$ in the fast wind displays a strong cut-off at $\delta B/|{\bf B}|\lesssim2$, as expected for fluctuations bounded on a sphere of radius $B=|{\bf B}|$. This is also associated with a saturation of the rms of the fluctuations at large scales and introduces a specific length $l_0$ above which the amplitude of the fluctuations becomes independent on the scale $l$. Consistent with that, the power spectrum at $l>l_0$ is characterised by a -1 spectral slope, as expected for fluctuations that are scale-independent. Moreover, we show that the spectral break between the 1/f and inertial range in solar wind spectra indeed corresponds to the scale $l_0$ at which $\left<\delta B/B\right>\sim1$. Such a simple model provides a possible alternative explanation of magnetic spectra observed in interplanetary space, also pointing out the inconsistency for a plasma to simultaneously maintain $|\bf{B}|\sim$const. at arbitrarily large scales and satisfy a Kolmogorov scaling.


title: Network Coding with Random Packet-Index Assignment for Data Collection Networks

abstract: This paper considers a data collection network consisting of uncoordinated, heterogeneous, and possibly mobile devices. These devices use medium and short-range radio technologies, which require multi-hop communication to deliver data to some sink. While numerous techniques from managed networks can be adapted, one of the most efficient (from the energy and spectrum use perspective) is \emph{network coding} (NC). NC is well suited to networks with mobility and unreliability, however, practical NC requires a precise identification of individual packets that have been mixed together. In a purely decentralized system, this requires either conveying identifiers in headers along with coded information as in COPE, or integrating a more complex protocol in order to efficiently identify the sources (participants) and their payloads. This paper presents a novel solution, Network Coding with Random Packet Index Assignment (NeCoRPIA), where packet indices in NC headers are selected in a decentralized way, by simply choosing them randomly. Traditional network decoding techniques apply directly when all original packets have different indices. When this is not the case, i.e., in case of collisions of indices, a specific decoding algorithm is proposed. A theoretical analysis of its performance in terms of complexity and decoding error probability is described. Simulation results match well the theoretical results. Comparisons of NeCoRPIA header lengths with those of a COPE-based NC protocol are also provided.


title: Multilayer Network Model of Movie Script

abstract: Network models have been increasingly used in the past years to support summarization and analysis of narratives, such as famous TV series, books and news. Inspired by social network analysis, most of these models focus on the characters at play. The network model well captures all characters interactions, giving a broad picture of the narration's content. A few works went beyond by introducing additional semantic elements, always captured in a single layer network. In contrast, we introduce in this work a multilayer network model to capture more elements of the narration of a movie from its script: people, locations, and other semantic elements. This model enables new measures and insights on movies. We demonstrate this model on two very popular movies.


title: Convergence of a Relaxed Variable Splitting Method for Learning Sparse Neural Networks via $\ell_1, \ell_0$, and transformed-$\ell_1$ Penalties

abstract: Sparsification of neural networks is one of the effective complexity reduction methods to improve efficiency and generalizability. We consider the problem of learning a one hidden layer convolutional neural network with ReLU activation function via gradient descent under sparsity promoting penalties. It is known that when the input data is Gaussian distributed, no-overlap networks (without penalties) in regression problems with ground truth can be learned in polynomial time at high probability. We propose a relaxed variable splitting method integrating thresholding and gradient descent to overcome the lack of non-smoothness in the loss function. The sparsity in network weight is realized during the optimization (training) process. We prove that under $\ell_1, \ell_0$; and transformed-$\ell_1$ penalties, no-overlap networks can be learned with high probability, and the iterative weights converge to a global limit which is a transformation of the true weight under a novel thresholding operation. Numerical experiments confirm theoretical findings, and compare the accuracy and sparsity trade-off among the penalties.


title: Why ReLU networks yield high-confidence predictions far away from the training data and how to mitigate the problem

abstract: Classifiers used in the wild, in particular for safety-critical systems, should not only have good generalization properties but also should know when they don't know, in particular make low confidence predictions far away from the training data. We show that ReLU type neural networks which yield a piecewise linear classifier function fail in this regard as they produce almost always high confidence predictions far away from the training data. For bounded domains like images we propose a new robust optimization technique similar to adversarial training which enforces low confidence predictions far away from the training data. We show that this technique is surprisingly effective in reducing the confidence of predictions far away from the training data while maintaining high confidence predictions and test error on the original classification task compared to standard training.


title: Stochastic Gradient Descent for Spectral Embedding with Implicit Orthogonality Constraint

abstract: In this paper, we propose a scalable algorithm for spectral embedding. The latter is a standard tool for graph clustering. However, its computational bottleneck is the eigendecomposition of the graph Laplacian matrix, which prevents its application to large-scale graphs. Our contribution consists of reformulating spectral embedding so that it can be solved via stochastic optimization. The idea is to replace the orthogonality constraint with an orthogonalization matrix injected directly into the criterion. As the gradient can be computed through a Cholesky factorization, our reformulation allows us to develop an efficient algorithm based on mini-batch gradient descent. Experimental results, both on synthetic and real data, confirm the efficiency of the proposed method in term of execution speed with respect to similar existing techniques.


title: Classes of operators related to m-isometric operators

abstract: Isometries played a pivotal role in the development of operator theory, in particular with the theory of contractions and polar decompositions and has been widely studied due to its fundamental importance in the theory of stochastic processes, the intrinsic problem of modeling the general contractive operator via its isometric dilation and many other areas in applied mathematics. In this paper we present some properties of n-quasi-(m;C)-isometric operators. We show that a power of a n-quasi-(m;C)-isometric operator is again a n-quasi-(m;C)-isometric operator and some products and tens


title: On the sign recovery by LASSO, thresholded LASSO and thresholded Basis Pursuit Denoising

abstract: Basis Pursuit (BP), Basis Pursuit DeNoising (BPDN), and LASSO are popular methods for identifying important predictors in the high-dimensional linear regression model, i.e. when the number of rows of the design matrix X is smaller than the number of columns. By definition, BP uniquely recovers the vector of regression coefficients b if there is no noise and the vector b has the smallest L1 norm among all vectors s such that Xb=Xs (identifiability condition). Furthermore, LASSO can recover the sign of b only under a much stronger irrepresentability condition. Meanwhile, it is known that the model selection properties of LASSO can be improved by hard-thresholding its estimates. This article supports these findings by proving that thresholded LASSO, thresholded BPDN and thresholded BP recover the sign of b in both the noisy and noiseless cases if and only if b is identifiable and large enough. In particular, if X has iid Gaussian entries and the number of predictors grows linearly with the sample size, then these thresholded estimators can recover the sign of b when the signal sparsity is asymptotically below the Donoho-Tanner transition curve. This is in contrast to the regular LASSO, which asymptotically recovers the sign of b only when the signal sparsity tends to 0. Numerical experiments show that the identifiability condition, unlike the irrepresentability condition, does not seem to be affected by the structure of the correlations in the $X$ matrix.


title: Surface Melting and Breakup of Metal Nanowires: Theory and Molecular Dynamics Simulations

abstract: We consider the surface melting of metal nanowires by solving a phenomenological two-parabola Landau model and by conducting molecular dynamics simulations of nickel and aluminium nanowires. The model suggests that surface melting will precede bulk melting when the spreading parameter $\Delta \gamma$ for the melt in contact with the solid surface is positive (i.e. if the melt wets or partially wets the surface) and the wire is sufficiently thick, as is the case for planar surfaces and sufficiently large nanoparticles. Surface melting does not occur if $\Delta \gamma$ is negative. We test this model, which assumes the surface energies of the wire are isotropic, using molecular dynamics simulations. For nickel, we observe the onset of anisotropic surface melting associated with each of the two surface facets present, but this gives way to uniform surface melting and the solid melts radially until the solid core eventually breaks up. For aluminium, while we observe complete surface melting of one facet, the lowest energy surface remains partially dry even up to the point where the melt completely penetrates the solid core.


title: Training Set Camouflage

abstract: We introduce a form of steganography in the domain of machine learning which we call training set camouflage. Imagine Alice has a training set on an illicit machine learning classification task. Alice wants Bob (a machine learning system) to learn the task. However, sending either the training set or the trained model to Bob can raise suspicion if the communication is monitored. Training set camouflage allows Alice to compute a second training set on a completely different -- and seemingly benign -- classification task. By construction, sending the second training set will not raise suspicion. When Bob applies his standard (public) learning algorithm to the second training set, he approximately recovers the classifier on the original task. Training set camouflage is a novel form of steganography in machine learning. We formulate training set camouflage as a combinatorial bilevel optimization problem and propose solvers based on nonlinear programming and local search. Experiments on real classification tasks demonstrate the feasibility of such camouflage.


title: Theoretical Evaluation of Electronic Density-of-states and Transport Effects on Field Emission from $n$-type Ultrananocrystalline Diamond Films

abstract: In the nitrogen-incorporated ultrananocrystalline diamond ((N)UNCD) films, representing an $n$-type highly conductive two-phase material comprised of $sp^3$ diamond grains and $sp^2$-rich graphitic grain boundaries, the current is carried by a high concentration of mobile electrons within the large-volume grain boundary networks. Fabricated in a simple thin-film planar form, (N)UNCD was found to be an efficient field emitter capable of emitting a significant amount of charge starting at the applied electric field as low as a few V/$\mu$m which makes it a promising material for designing electron sources. Despite the semimetallic conduction, field emission (FE) characteristics of this material demonstrate a strong deviation from the Fowler-Nordheim law in a high-current-density regime when (N)UNCD field emitters switch from a diode-like to resistor-like behavior. Such phenomenon resembles the current-density saturation effect in conventional semiconductors. In the present paper, we adapt the formalism developed for conventional semiconductors to study current-density saturation in (N)UNCD field emitters. We provide a comprehensive theoretical investigation of ($i$) the influence of partial penetration of the electric field into the material, ($ii$) transport effects (such as electric-field-dependent mobility), and ($iii$) features of a complex density-of-states structure (position and shape of $\pi-\pi^*$ bands, controlling the concentration of charge carriers) on the FE characteristics of (N)UNCD. We show that the formation of the current-density saturation plateau can be explained by the limited supply of electrons within the impurity $\pi-\pi^*$ bands and decreasing electron mobility in high electric field. Theoretical calculations are consistent with experiment.


title: Processing Transactions in a Predefined Order

abstract: In this paper we provide a high performance solution to the problem of committing transactions while enforcing a predefined order. We provide the design and implementation of three algorithms, which deploy a specialized cooperative transaction execution model. This model permits the propagation of written values along the chain of ordered transactions. We show that, even in the presence of data conflicts, the proposed algorithms are able to outperform single-threaded execution, and other baseline and specialized state-of-the-art competitors (e.g., STMLite). The maximum speedup achieved in micro benchmarks, STAMP, PARSEC and SPEC200 applications is in the range of 4.3x -- 16.5x.


title: A Consideration of Cosmic Evolution from the Points of View of the Inflationary and Cyclic Theories

abstract: This study reviews the advances made in inflationary theory, especially regarding the seeming disparity between inflation energy and dark energy, and their significance to cosmic evolution as a whole. I attempt to connect the two sources of expansion and thereby enhance the predictive capacity of the consensus model. I also attempt to contrast its strengths and weaknesses with those of the cyclic theory developed by Turok and Steinhardt, particularly where cosmic evolution over larger time scales is concerned. Furthermore, I endeavor to provide a physical meaning to the existence and workings of dark energy as well as dark matter in the cyclic model (as compared to their simply being extra parameters in the inflation picture), so that each observed component of the universe is given significance.


title: A Novel PN junction between Mechanically Exfoliated \b{eta}-Ga2O3 and p-GaN

abstract: Several pn junctions were constructed from mechanically exfoliated ultrawide bandgap (UWBG) beta-phase gallium oxide (\b{eta}-Ga2O3) and p-type gallium nitride (GaN). The mechanical exfoliation process, which is described in detail, is similar to that of graphene and other 2D materials. Atomic force microscopy (AFM) scans of the exfoliated \b{eta}-Ga2O3 flakes show very smooth surfaces with average roughness of 0.647 nm and transmission electron microscopy (TEM) scans reveal flat, clean interfaces between the \b{eta}-Ga2O3 flakes and p-GaN. The device showed a rectification ratio around 541.3 (V+5/V-5). Diode performance improved over the temperature range of 25{\deg}C and 200{\deg}C, leading to an unintentional donor activation energy of 135 meV. As the thickness of exfoliated \b{eta}-Ga2O3 increases, ideality factors decrease as do the diode turn on voltages, tending toward an ideal threshold voltage of 3.2 V as determined by simulation. This investigation can help increase study of novel devices between mechanically exfoliated \b{eta}-Ga2O3 and other materials.


title: A 6.7 GHz methanol maser survey of the central molecular zone

abstract: The Central Molecular Zone (CMZ) spans the inner ~450 pc (3 degrees) of our Galaxy. This region is defined by its enhanced molecular emission and contains 5% of the entire Galaxy's molecular gas mass. However, the number of detected star forming sites towards the CMZ may be low for the amount of molecular gas that is present, and improved surveys of star formation indicators can help clarify this. With the Karl G Jansky Very Large Array (VLA), we conducted a blind survey of 6.7 GHz methanol masers spanning the inner 3deg x 40arcmin (450 pc x 100 pc) of the Galaxy. We detected 43 methanol masers towards 28 locations, 16 of which are new detections. The velocities of most of these masers are consistent with being located within the CMZ. A majority of the detected methanol masers are distributed towards positive Galactic longitudes, similar to 2/3 of the molecular gas mass distributed at positive Galactic longitudes. The 6.7 GHz methanol maser is an excellent indicator of high mass (>8 solar mass) star formation, with new detections indicating sites of massive star formation in the CMZ.


title: Revisiting Iterative Relevance Feedback for Document and Passage Retrieval

abstract: As more and more search traffic comes from mobile phones, intelligent assistants, and smart-home devices, new challenges (e.g., limited presentation space) and opportunities come up in information retrieval. Previously, an effective technique, relevance feedback (RF), has rarely been used in real search scenarios due to the overhead of collecting users' relevance judgments. However, since users tend to interact more with the search results shown on the new interfaces, it becomes feasible to obtain users' assessments on a few results during each interaction. This makes iterative relevance feedback (IRF) techniques look promising today. IRF has not been studied systematically in the new search scenarios and its effectiveness is mostly unknown. In this paper, we re-visit IRF and extend it with RF models proposed in recent years. We conduct extensive experiments to analyze and compare IRF with the standard top-k RF framework on document and passage retrieval. Experimental results show that IRF is at least as effective as the standard top-k RF framework for documents and much more effective for passages. This indicates that IRF for passage retrieval has huge potential.


title: Strength and Cohesive Behavior of Thermoset Polymers at the Microscale: A Size-Effect Study

abstract: This study investigated, experimentally and numerically, the fracturing behavior of thermoset polymer structures featuring cracks and sharp u-notches. It is shown that, even for cases in which the sharpness of the notch would suggest otherwise, the failure behavior of cracked and pre-notched specimens is substantially different, the failure loads of the former configuration being about three times lower than the latter one. To capture this interesting behavior a two-scale cohesive model is proposed. The model is in excellent agreement with the experimental data and its predictions allow to conclude that (a) residual plastic stresses cannot explain the very high failure loads of notched structures; (b) the strength of the polymer at the microscale can be from six to ten times larger than the values measured from conventional tests whereas the fracture energy at the microscale can be about forty times lower; (c) the pre-notched specimens investigated in this work failed when the stress at the tip reached the microscale strength whereas the cracked specimens failed when the energy release rate reached the total fracture energy of the material. The foregoing considerations are of utmost importance for the design of microelectronic devices or polymer matrix composites for which the main damage mechanisms are governed by the strength and cohesive behavior at the microscale.


title: Constraining Scatter in the Stellar Mass--Halo Mass Relation for Haloes Less Massive than the Milky Way

abstract: Most galaxies are hosted by massive, invisible dark matter haloes, yet little is known about the scatter in the stellar mass--halo mass relation for galaxies with host halo masses $M_{h}\le 10^{11}M_{\odot}$. Using mock catalogues based on dark matter simulations, we find that two observable signatures are sensitive to scatter in the stellar mass--halo mass relation even at these mass scales; i.e., conditional stellar mass functions and velocity distribution functions for neighbouring galaxies. We compute these observables for 179,373 galaxies in the Sloan Digital Sky Survey (SDSS) with stellar masses $M_{\ast} > 10^9 M_{\odot}$ and redshifts 0.01 $< z <$ 0.307. We then compare to mock observations generated from the $\textit{Bolshoi-Planck}$ dark matter simulation for stellar mass--halo mass scatters ranging from 0 to 0.6 dex. The observed results are consistent with simulated results for most values of scatter ($<$0.6 dex), and SDSS statistics are insufficient to provide firm constraints. However, this method could provide much tighter constraints on stellar mass--halo mass scatter in the future if applied to larger data sets, especially the anticipated Dark Energy Spectroscopic Instrument Bright Galaxy Survey. Constraining the value of scatter could have important implications for galaxy formation and evolution.


title: Graphs that are cospectral for the distance Laplacian

abstract: The distance matrix $\mathcal{D}(G)$ of a graph $G$ is the matrix containing the pairwise distances between vertices, and the distance Laplacian matrix is $\mathcal{D}^L(G)=T(G)-\mathcal{D}(G)$, where $T(G)$ is the diagonal matrix of row sums of $\mathcal{D}(G)$. We establish several general methods for producing $\mathcal{D}^L$-cospectral graphs that can be used to construct infinite families. We provide examples showing that various properties are not preserved by $\mathcal{D}^L$-cospectrality, including examples of $\mathcal{D}^L$-cospectral strongly regular and circulant graphs. We establish that the absolute values of coefficients of the distance Laplacian characteristic polynomial are decreasing, i.e., $|\delta^L_{1}|\geq \dots \geq |\delta^L_{n}|$ where $\delta^L_{k}$ is the coefficient of $x^k$.


title: EDRIXS: An open source toolkit for simulating spectra of resonant inelastic x-ray scattering

abstract: Resonant inelastic x-ray scattering (RIXS) has become a very powerful experimental technique to probe a broad range of intrinsic elementary excitations, for example, from low energy phonons and (bi-)magnons to high energy $d$-$d$, charge-transfer and plasmons excitations in strongly correlated electronic systems. Due to the complexity of the RIXS cross-section and strong core-hole effects, theoretical simulation of the experimental RIXS spectra is still a difficult task which hampers the understanding of RIXS spectra and the development of the RIXS technique. In this paper, we present an open source toolkit (dubbed EDRIXS) to facilitate the simulations of RIXS spectra of strongly correlated materials based on exact diagonalization (ED) of certain model Hamiltonians. The model Hamiltonian can be from a single atom, small cluster or Anderson impurity model, with model parameters from density functional theory plus Wannier90 or dynamical mean-field theory calculations. The spectra of x-ray absorption spectroscopy (XAS) and RIXS are then calculated using Krylov subspace techniques. This toolkit contains highly efficient ED, XAS and RIXS solvers written in modern Fortran 90 language and a convenient Python library used to prepare inputs and set up calculations. We first give a short introduction to RIXS spectroscopy, and then we discuss the implementation details of this toolkit. Finally, we show three examples to demonstrate its usage.


title: Detecting unseen visual relations using analogies

abstract: We seek to detect visual relations in images of the form of triplets t = (subject, predicate, object), such as "person riding dog", where training examples of the individual entities are available but their combinations are unseen at training. This is an important set-up due to the combinatorial nature of visual relations : collecting sufficient training data for all possible triplets would be very hard. The contributions of this work are three-fold. First, we learn a representation of visual relations that combines (i) individual embeddings for subject, object and predicate together with (ii) a visual phrase embedding that represents the relation triplet. Second, we learn how to transfer visual phrase embeddings from existing training triplets to unseen test triplets using analogies between relations that involve similar objects. Third, we demonstrate the benefits of our approach on three challenging datasets : on HICO-DET, our model achieves significant improvement over a strong baseline for both frequent and unseen triplets, and we observe similar improvement for the retrieval of unseen triplets with out-of-vocabulary predicates on the COCO-a dataset as well as the challenging unusual triplets in the UnRel dataset.


title: Effectiveness of Hierarchical Softmax in Large Scale Classification Tasks

abstract: Typically, Softmax is used in the final layer of a neural network to get a probability distribution for output classes. But the main problem with Softmax is that it is computationally expensive for large scale data sets with large number of possible outputs. To approximate class probability efficiently on such large scale data sets we can use Hierarchical Softmax. LSHTC datasets were used to study the performance of the Hierarchical Softmax. LSHTC datasets have large number of categories. In this paper we evaluate and report the performance of normal Softmax Vs Hierarchical Softmax on LSHTC datasets. This evaluation used macro f1 score as a performance measure. The observation was that the performance of Hierarchical Softmax degrades as the number of classes increase.


title: Gravitational Wave Emission from 3D Explosion Models of Core-Collapse Supernovae with Low and Normal Explosion Energies

abstract: Understanding gravitational wave emission from core-collapse supernovae will be essential for their detection with current and future gravitational wave detectors. This requires a sample of waveforms from modern 3D supernova simulations reaching well into the explosion phase, where gravitational wave emission is expected to peak. However, recent waveforms from 3D simulations with multi-group neutrino transport do not reach far into the explosion phase, and some are still obtained from non-exploding models. We therefore calculate waveforms up to 0.9\,s after bounce using the neutrino hydrodynamics code \textsc{CoCoNuT-FMT}. We consider two models with low and normal explosion energy, namely explosions of an ultra-stripped progenitor with an initial helium star mass of $3.5\,M_{\odot}$, and of an $18\,M_{\odot}$ single star. Both models show gravitational wave emission from the excitation of surface g-modes in the proto-neutron star with frequencies between $\mathord{\sim}800\,\mathrm{Hz}$ and 1000\,Hz at peak emission. The peak amplitudes are about $6\, \mathrm{cm}$ and $10\, \mathrm{cm}$, respectively, which is somewhat higher than in most recent 3D models of the pre-explosion or early explosion phase. Using a Bayesian analysis, we determine the maximum detection distances for our models in simulated Advanced LIGO, Advanced Virgo, and Einstein Telescope design sensitivity noise. The more energetic $18 M_\odot$ explosion will be detectable to about $17.5 \,\mathrm{kpc}$ by the LIGO/Virgo network and to about $180\, \mathrm{kpc}$ with the Einstein Telescope.


title: The Tersoff potential for extreme environment

abstract: A novel modification of the Tersoff potential for Si is presented. The modification improves the transferability of the Tersoff potential for liquid states without the change of original parameters and with no alteration of bulk properties. Also, the modification introduces a correction term for high-pressure states. The modification is meaningful considering that by high energy irradiations local liquid structures and unstable high-pressure manifolds may occur, therefore an interatomic potential must have an acceptable reliability on high thermal/pressure situations to simulate such phenomenon. Particularly, in the modification, a new screening function replaces the radial cutoff function and the bond order function is slightly changed. Also, a repulsive energy function is replaced by a correction function within a specific pair distance.


title: Pay Voice: Point of Sale Recognition for Visually Impaired People

abstract: Millions of visually impaired people depend on relatives and friends to perform their everyday tasks. One relevant step towards self-sufficiency is to provide them with means to verify the value and operation presented in payment machines. In this work, we developed and released a smartphone application, named Pay Voice, that uses image processing, optical character recognition (OCR) and voice synthesis to recognize the value and operation presented in POS and PIN pad machines, and thus informing the user with auditive and visual feedback. The proposed approach presented significant results for value and operation recognition, especially for POS, due to the higher display quality. Importantly, we achieved the key performance indicators, namely, more than 80% of accuracy in a real-world scenario, and less than $5$ seconds of processing time for recognition. Pay Voice is publicly available on Google Play and App Store for free.


title: Constrained inference through posterior projections

abstract: Bayesian approaches are appealing for constrained inference problems by allowing a probabilistic characterization of uncertainty, while providing a computational machinery for incorporating complex constraints in hierarchical models. However, the usual Bayesian strategy of placing a prior on the constrained space and conducting posterior computation with Markov chain Monte Carlo algorithms is often intractable. An alternative is to conduct inference for a less constrained posterior and project samples to the constrained space through a minimal distance mapping. We formalize and provide a unifying framework for such posterior projections. For theoretical tractability, we initially focus on constrained parameter spaces corresponding to closed and convex subsets of the original space. We then consider non-convex Stiefel manifolds. We provide a general formulation of projected posteriors in a Bayesian decision-theoretic framework. We show that asymptotic properties of the unconstrained posterior are transferred to the projected posterior, leading to asymptotically correct credible intervals. We demonstrate numerically that projected posteriors can have better performance that competitor approaches in real data examples.


title: The N2O-CS2 dimer is cross-shaped

abstract: The infrared spectrum of the cross-shaped van der Waals complex N2O-CS2 is observed in the region of the N2O nu1 fundamental band (~2220 cm-1) using a tuneable diode laser to probe a pulsed supersonic slit jet expansion. Both 14N- and 15N-substituted species are studied. Analysis of their spectra establishes that this dimer has a cross-shaped structure, similar to its isoelectronic cousin CO2-CS2. This is the first spectroscopic observation of N2O-CS2, and the molecular parameters determined here should be useful for detection of its pure rotational microwave spectrum.


title: An Incentive-Aware Job Offloading Control Framework for Mobile Edge Computing

abstract: This paper considers a scenario in which an access point (AP) is equipped with a mobile edge server of finite computing power, and serves multiple resource-hungry mobile users by charging users a price. Pricing provides users with incentives in offloading. However, existing works on pricing are based on abstract concave utility functions (e.g, the logarithm function), giving no dependence on physical layer parameters. To that end, we first introduce a novel utility function, which measures the cost reduction by offloading as compared with executing jobs locally. Based on this utility function we then formulate two offloading games, with one maximizing individual's interest and the other maximizing the overall system's interest. We analyze the structural property of the games and admit in closed form the Nash Equilibrium and the Social Equilibrium, respectively. The proposed expressions are functions of the user parameters such as the weights of computational time and energy, the distance from the AP, thus constituting an advancement over prior economic works that have considered only abstract functions. Finally, we propose an optimal pricing-based scheme, with which we prove that the interactive decision-making process with self-interested users converges to a Nash Equilibrium point equal to the Social Equilibrium point.


title: FREE -- Fine-grained Scheduling for Reliable and Energy Efficient Data Collection in LoRaWAN

abstract: LoRaWAN promises to provide wide-area network access to low-cost devices that can operate for up to 10 years on a single 1000 mAh battery. This makes LoRaWAN particularly suited to data collection applications (e.g. monitoring applications), where device lifetime is a key performance metric. However, when supporting a large number of devices, LoRaWAN suffers from a scalability issue due to the high collision probability of its Aloha-based MAC layer. The performance worsens further when using acknowledged transmissions due to the duty cycle restriction at the gateway. For this, we propose FREE, a fine-grained scheduling scheme for reliable and energy-efficient data collection in LoRaWAN. FREE takes advantage of applications that do not have hard delay requirements on data delivery by supporting synchronized bulk data transmission. This means data is buffered for transmission in scheduled time slots instead of transmitted straight away. FREE allocates spreading factors, transmission powers, frequency channels, time slots, and schedules slots in frames for LoRaWAN end-devices. As a result, FREE overcomes the scalability problem of LoRaWAN by eliminating collisions and grouping acknowledgments. We evaluate the performance of FREE versus different legacy LoRaWAN configurations. The numerical results show that FREE scales well and achieves almost 100% data delivery and the device lifetime is estimated to over 10 years independent of traffic type and network size. Comparing to poor scalability, low data delivery and device lifetime of fewer than 2 years for acknowledged data traffic in the standard LoRaWAN configurations.


title: Data Security and Privacy Protection Data Security and Privacy Protection in Public Cloud

abstract: This paper discusses about the challenges, advantages and shortcomings of existing solutions in data security and privacy in public cloud computing. As in cloud computing, oceans of data will be stored. Data stored in public cloud would face both outside attacks and inside attacks since public cloud provider themselves are untrusted. Conventional encryption could be used for storage, however most data in cloud needs further computation. Decryption before computation will cause large overheads for data operation and lots of inconvenience. Thus, efficient methods to protect data security as well as privacy for large amount of data in cloud are necessary. In the paper, different mechanisms to protect data security and privacy in public cloud are discussed. A data security and privacy enabled multi-cloud architecture is proposed.


title: On uniform K-stability of pairs

abstract: In this paper, we discuss stable pairs, which were first studied by S. Paul, and give a proof for a result I learned from him. As a consequence, we will show that the K-stability implies the CM-stability.


title: Field Angle Tuned Metamagnetism and Lifschitz Transitions in UPt3

abstract: Strongly correlated electronic systems can harbor a rich variety of quantum spin states. Understanding and controlling such spin states in quantum materials is of great current interest. Focusing on the simple binary system UPt3 with ultrasound (US) as a probe we identify clear signatures in field sweeps demarkating new high field spin states. Magnetostriction (MS) measurements performed up to 65 T also show signatures at the same fields confirming these state transitions. At the very lowest temperatures (<200 mK) we also observe magneto-acoustic quantum oscillations which for theta = 90{\deg} and vicinity abruptly become very strong in the 24.8-30 T range. High resolution magnetization measurements for this same angle reveal a continuous variation of the magnetization implying the subtle nature of the implied transitions. With B rotated away from the c-axis, the US signatures occur at nearly the same field. These state transitions merge with the separate sequence of the well known metamagnetic transition which commences at 20 T for theta = 0{\deg} but moves to higher fields as 1/cos(\theta). This merge, suggesting a tricritical behavior, occurs at \theta ~ 51{\deg} from the ab-plane. This is an unique off-symmetry angle where the length change is precisely zero due to the anisotropic nature of MS in UPt3 for all magnetic field values.


title: Dynamic Programming with Recursive Preferences: Optimality and Applications

abstract: This paper provides new conditions for dynamic optimality in discrete time and uses them to establish fundamental dynamic programming results for several commonly used recursive preference specifications. These include Epstein-Zin preferences, risk-sensitive preferences, narrow framing models and recursive preferences with sensitivity to ambiguity. The results obtained for these applications include (i) existence of optimal policies, (ii) uniqueness of solutions to the Bellman equation, (iii) a complete characterization of optimal policies via Bellman's principle of optimality, and (iv) a globally convergent method of computation via value function iteration.


title: Quantum Reflection and Transmission in Ring Systems with Double Y-Junctions: Occurrence of Perfect Reflection

abstract: We consider the scattering problems of a quantum particle in a system with a single Y-junction and in ring systems with double Y-junctions. We provide new formalism for such quantum mechanical problems. Based on a path integral approach, we find compact formulas for probability amplitudes in the ring systems. We also discuss quantum reflection and transmission in the ring systems under scale-invariant junction conditions. It is remarkable that perfect reflection can occur in an anti-symmetric ring system, in contrast with the one-dimensional quantum systems having singular nodes of degree 2.


title: Ordered and convex geometric trees with linear extremal function

abstract: The extremal functions $ex_{\rightarrow}(n,F)$ and $ex_{\cir}(n,F)$ for ordered and convex geometric acyclic graphs $F$ have been extensively investigated by a number of researchers. Basic questions are to determine when $ex_{\rightarrow}(n,F)$ and $ex_{\cir}(n,F)$ are linear in $n$, the latter posed by Bra\ss-K\'arolyi-Valtr in 2003. In this paper, we answer both these questions for every tree $F$. We give a forbidden subgraph characterization for a family $\cal T$ of ordered trees with $k$ edges, and show that $ex_{\rightarrow}(n,T) = (k - 1)n - {k \choose 2}$ for all $n \geq k + 1$ when $T \in {\cal T}$ and $ex_{\rightarrow}(n,T) = \Omega(n\log n)$ for $T \not\in {\cal T}$. We also describe the family of the convex geometric trees with linear Tur\' an number and show that for every convex geometric tree $F$ not in this family, $ex_{\cir}(n,F)= \Omega(n\log \log n)$.


title: Modeling of the sunspot-associated microwave emission using a new method of DEM inversion

abstract: We developed a method to compute the temperature and density structure along the line of sight by inversion of the differential emission measure (DEM), under the assumptions of stratification and hydrostatic equilibrium. We applied this method to the DEM obtained from AIA observations and used the results, together with potential extrapolations of the photosheric magnetic field, to compute the microwave emission of three sunspots, which we compared with observations from the RATAN-600 radio telescope and the Nobeyama Radioheliograph (NoRH). Our DEM based models reproduced very well the observations of the moderate-size spot on October 2011 and within 25% the data of a similar sized spot on March 2016, but predicted too low values for the big spot of April 14, 2016. The latter was better fitted by a constant conductive flux atmospheric model which, however, could not reproduce the peak brightness temperature of $4.7\times10^6$ K and the shape of the source at the NoRH frequency. We propose that these deviations could be due to low intensity non-thermal emission associated to a moving pore and to an opposite polarity light bridge. We also found that the double structure of the big spot at high RATAN-600 frequencies could be interpreted in terms of the variation of the angle between the magnetic field and the line of sight along the sunspot.


title: A Commute in Data: The comma2k19 Dataset

abstract: comma.ai presents comma2k19, a dataset of over 33 hours of commute in California's 280 highway. This means 2019 segments, 1 minute long each, on a 20km section of highway driving between California's San Jose and San Francisco. The dataset was collected using comma EONs that have sensors similar to those of any modern smartphone including a road-facing camera, phone GPS, thermometers and a 9-axis IMU. Additionally, the EON captures raw GNSS measurements and all CAN data sent by the car with a comma grey panda. Laika, an open-source GNSS processing library, is also introduced here. Laika produces 40% more accurate positions than the GNSS module used to collect the raw data. This dataset includes pose (position + orientation) estimates in a global reference frame of the recording camera. These poses were computed with a tightly coupled INS/GNSS/Vision optimizer that relies on data processed by Laika. comma2k19 is ideal for development and validation of tightly coupled GNSS algorithms and mapping algorithms that work with commodity sensors.


title: Type $\theta$ Stokes' Theorem for Chains

abstract: We present type $\theta$ Stokes' theorem for type $\theta$ $k$-chains which extends the fundamental theorem of calculus in higher dimensions.


title: Numerical calculation of overlap distribution of $(2+1)$-dimensional directed polymer in random media

abstract: We investigate $(2+1)$-dimensional discretized directed polymers in Gaussian random media. By numerically calculating the probability distribution function of overlap between two independent and identical systems on a common random potential, we show that there is no replica symmetry breaking. We also show that while the mean-squared displacement of one polymer end exhibits superdiffusive behavior, the mean-squared relative distance of two polymer ends exhibits subdiffusive behavior.


title: Transport through periodically driven systems: Green's function approach formulated within frequency domain

abstract: The steady-state electronic transport across periodically driven systems can be efficiently addressed using Landauer-B\"{u}ttiker formalism. The time-dependent nonequilibrium Green's function theory then may be adapted for developing direct and universal calculation schemes. Here we propose an alternative scheme to carry out the calculations. The idea is based on treating the transport problem in frequency domain and designing Green's functions for the corresponding Hamiltonians rewritten in Floquet-Sambe formalism. We show that within our approach the expressions for time-periodic currents and densities essentially replicate the well known formulas from time-independent theory. The results are then simplified for easier implementation in numeric computations.


title: Lost Waterways: Clues from Digitized Historical Maps of Manila and Other Philippine Cities

abstract: We search for lost bodies of water in the cities of Manila, Tacloban, Iloilo, Cebu, Davao, and Naga by aligning their digitized Spanish-era and American-era maps to Google maps. These vanished ancient waterways can either become flooding hazards in case of extreme weather events, or liquefaction hazards, in case of earthquakes. Digitized historical maps of the cities were georectified, overlaid on current Google maps, and checked for potential missing bodies of water. Inspection through field visits and interviews with locals were conducted to verify the actual status of suspected sites. The validation identified lost, found, and even new bodies of water. There was also evidence of affected buildings, rainless flooding, and a new normal for the meaning of flooding among frequently inundated residents.


title: Frequency-resolved photon-electronic spectroscopy for excited state population detection

abstract: Atomic excitation to excited states in strong laser field is the key to high-order harmonic generation below ionization threshold, yet remains unclear mainly due to the lack of proper detection methods. We propose a frequency-resolved photon-electron spectroscopy technique to reconstruct population of excited states with the second delayed laser pulse. The technique utilizes Fourier transformation to separate ionization from different excited states to different positions on the spectrum. With the advantage of separation, we provide a scheme to reconstruct populations on different excited states after the first pulse. The scheme is validated by high-precision population reconstruction of helium and hydrogen atoms.


title: On Stacked Denoising Autoencoder based Pre-training of ANN for Isolated Handwritten Bengali Numerals Dataset Recognition

abstract: This work attempts to find the most optimal parameter setting of a deep artificial neural network (ANN) for Bengali digit dataset by pre-training it using stacked denoising autoencoder (SDA). Although SDA based recognition is hugely popular in image, speech and language processing related tasks among the researchers, it was never tried in Bengali dataset recognition. For this work, a dataset of 70000 handwritten samples were used from (Chowdhury and Rahman, 2016) and was recognized using several settings of network architecture. Among all these settings, the most optimal setting being found to be five or more deeper hidden layers with sigmoid activation and one output layer with softmax activation. We proposed the optimal number of neurons that can be used in the hidden layer is 1500 or more. The minimum validation error found from this work is 2.34% which is the lowest error rate on handwritten Bengali dataset proposed till date.


title: Simulation and inference algorithms for stochastic biochemical reaction networks: from basic concepts to state-of-the-art

abstract: Stochasticity is a key characteristic of intracellular processes such as gene regulation and chemical signalling. Therefore, characterising stochastic effects in biochemical systems is essential to understand the complex dynamics of living things. Mathematical idealisations of biochemically reacting systems must be able to capture stochastic phenomena. While robust theory exists to describe such stochastic models, the computational challenges in exploring these models can be a significant burden in practice since realistic models are analytically intractable. Determining the expected behaviour and variability of a stochastic biochemical reaction network requires many probabilistic simulations of its evolution. Using a biochemical reaction network model to assist in the interpretation of time course data from a biological experiment is an even greater challenge due to the intractability of the likelihood function for determining observation probabilities. These computational challenges have been subjects of active research for over four decades. In this review, we present an accessible discussion of the major historical developments and state-of-the-art computational techniques relevant to simulation and inference problems for stochastic biochemical reaction network models. Detailed algorithms for particularly important methods are described and complemented with MATLAB implementations. As a result, this review provides a practical and accessible introduction to computational methods for stochastic models within the life sciences community.


title: A Halo Occupation Interpretation Of Quasars At $z\sim1.5$ Using Very Small Scale Clustering Information

abstract: We combine the most precise small scale ($< 100\, \rm h^{-1}kpc$) quasar clustering constraintsto date with recent measurements at large scales ($> 1\, \rm h^{-1}Mpc$) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) to better constrain the satellite fraction of quasars at $z\sim 1.5$ in the halo occupation formalism. We build our Halo Occupation Distribution (HOD) framework based on commonly used analytic forms for the one and two-halo terms with two free parameters: the minimum halo mass that hosts a central quasar and the fraction of satellite quasars that are within one halo. Inspired by recent studies that propose a steeper density profile for the dark matter haloes that host quasars, we explore HOD models at kiloparsec scales and best-fit parameters for models with $10\times$ higher concentration parameter. We find that an HOD model with a satellite fraction of $f_{\rm sat} = 0.071_{-0.004}^{+0.009}$ and minimum mass of $\rm M_{m} = 2.31_{-0.38}^{+0.41} \times 10^{12}\, \, \rm h^{-1} M_{\odot}$ for the host dark matter haloes best describes quasar clustering (on all scales) at $z \sim 1.5$. Our results are marginally inconsistent with earlier work that studied brighter quasars, hinting at a luminosity-dependence to the one-halo term.


title: Groundstates for Choquard equations with the upper critical exponent

abstract: In this paper, an autonomous Choquard equation with the upper critical exponent is considered. By using the Poho\v{z}aev constraint method, the subcritical approximation method and the compactness lemma of Strauss, a groundstate solution in $H^1(\mathbb{R}^N)$ which is positive and radially symmetric is obtained. The result here extends and complements the earlier theorems.