title: Switching EEG Headsets Made Easy: Reducing Offline Calibration Effort Using Active Weighted Adaptation Regularization

abstract: Electroencephalography (EEG) headsets are the most commonly used sensing devices for Brain-Computer Interface. In real-world applications, there are advantages to extrapolating data from one user session to another. However, these advantages are limited if the data arise from different hardware systems, which often vary between application spaces. Currently, this creates a need to recalibrate classifiers, which negatively affects people's interest in using such systems. In this paper, we employ active weighted adaptation regularization (AwAR), which integrates weighted adaptation regularization (wAR) and active learning, to expedite the calibration process. wAR makes use of labeled data from the previous headset and handles class-imbalance, and active learning selects the most informative samples from the new headset to label. Experiments on single-trial event-related potential classification show that AwAR can significantly increase the classification accuracy, given the same number of labeled samples from the new headset. In other words, AwAR can effectively reduce the number of labeled samples required from the new headset, given a desired classification accuracy, suggesting value in collating data for use in wide scale transfer-learning applications.


title: PowerAlert: An Integrity Checker using Power Measurement

abstract: We propose PowerAlert, an efficient external integrity checker for untrusted hosts. Current attestation systems suffer from shortcomings in requiring complete checksum of the code segment, being static, use of timing information sourced from the untrusted machine, or use of timing information with high error (network round trip time). We address those shortcomings by (1) using power measurements from the host to ensure that the checking code is executed and (2) checking a subset of the kernel space over a long period of time. We compare the power measurement against a learned power model of the execution of the machine and validate that the execution was not tampered. Finally, power diversifies the integrity checking program to prevent the attacker from adapting. We implement a prototype of PowerAlert using Raspberry pi and evaluate the performance of the integrity checking program generation. We model the interaction between PowerAlert and an attacker as a game. We study the effectiveness of the random initiation strategy in deterring the attacker. The study shows that \power forces the attacker to trade-off stealthiness for the risk of detection, while still maintaining an acceptable probability of detection given the long lifespan of stealthy attacks.


title: Subconvex equidistribution of cusp forms: reduction to Eisenstein observables

abstract: Let $\pi$ traverse a sequence of cuspidal automorphic representations of GL(2) with large prime level, unramified central character and bounded infinity type. For G either of the groups GL(1) or PGL(2), let H(G) denote the assertion that subconvexity holds for G-twists of the adjoint $L$-function of $\pi$, with polynomial dependence upon the conductor of the twist. We show that H(GL(1)) implies H(PGL(2)). In geometric terms, H(PGL(2)) corresponds roughly to an instance of arithmetic quantum unique ergodicity with a power savings in the error term, H(GL(1)) to the special case in which the relevant sequence of measures is tested against an Eisenstein series.


title: Active subspaces of airfoil shape parameterizations

abstract: Design and optimization benefit from understanding the dependence of a quantity of interest (e.g., a design objective or constraint function) on the design variables. A low-dimensional active subspace, when present, identifies important directions in the space of design variables; perturbing a design along the active subspace associated with a particular quantity of interest changes that quantity more, on average, than perturbing the design orthogonally to the active subspace. This low-dimensional structure provides insights that characterize the dependence of quantities of interest on design variables. Airfoil design in a transonic flow field with a parameterized geometry is a popular test problem for design methodologies. We examine two particular airfoil shape parameterizations, PARSEC and CST, and study the active subspaces present in two common design quantities of interest, transonic lift and drag coefficients, under each shape parameterization. We mathematically relate the two parameterizations with a common polynomial series. The active subspaces enable low-dimensional approximations of lift and drag that relate to physical airfoil properties. In particular, we obtain and interpret a two-dimensional approximation of both transonic lift and drag, and we show how these approximation inform a multi-objective design problem.


title: Roche tomography of CVs

abstract: Roche tomography is a technique designed for mapping the intensity distribution over the surface of binary star components that fill (or nearly fill) their Roche lobes and so are both rotationally and tidally distorted. It builds on and develops the successful techniques used for mapping rapidly rotating single stars, which have revealed the presence of starspots. It has so far been applied to only a small number of cataclysmic variables, and a few other related systems. This paper reviews the technique and presents maps of spots on most of these systems. The review has been written with the help of contributions of material from Chris Watson and Colin Hill of Queen's University, Belfast.


title: $G$-Strands on symmetric spaces

abstract: We study the $G$-strand equations that are extensions of the classical chiral model of particle physics in the particular setting of broken symmetries described by symmetric spaces. These equations are simple field theory models whose configuration space is a Lie group, or in this case a symmetric space. In this class of systems, we derive several models that are completely integrable on finite dimensional Lie group $G$ and we treat in more details examples with symmetric space $SU(2)/S^1$ and $SO(4)/SO(3)$. The later model simplifies to an apparently new integrable $9$ dimensional system. We also study the $G$-strands on the infinite dimensional group of diffeomorphisms, which gives, together with the Sobolev norm, systems of $1+2$ Camassa-Holm equations. The solutions of these equations on the complementary space related to the Witt algebra decomposition are the odd function solutions.


title: Randomized Dynamic Mode Decomposition

abstract: This paper presents a randomized algorithm for computing the near-optimal low-rank dynamic mode decomposition (DMD). Randomized algorithms are emerging techniques to compute low-rank matrix approximations at a fraction of the cost of deterministic algorithms, easing the computational challenges arising in the area of `big data'. The idea is to derive a small matrix from the high-dimensional data, which is then used to efficiently compute the dynamic modes and eigenvalues. The algorithm is presented in a modular probabilistic framework, and the approximation quality can be controlled via oversampling and power iterations. The effectiveness of the resulting randomized DMD algorithm is demonstrated on several benchmark examples of increasing complexity, providing an accurate and efficient approach to extract spatiotemporal coherent structures from big data in a framework that scales with the intrinsic rank of the data, rather than the ambient measurement dimension. For this work we assume that the dynamics of the problem under consideration is evolving on a low-dimensional subspace that is well characterized by a fast decaying singular value spectrum.


title: Convergence radius of perturbative Lindblad driven non-equilibrium steady states

abstract: We address the problem of analyzing the radius of convergence of perturbative expansion of non-equilibrium steady states of Lindblad driven spin chains. A simple formal approach is developed for systematically computing the perturbative expansion of small driven systems. We consider the paradigmatic model of an open $XXZ$ spin 1/2 chain with boundary supported ultralocal Lindblad dissipators and treat two different perturbative cases: (i) expansion in system-bath coupling parameter and (ii) expansion in driving (bias) parameter. In the first case (i) we find that the radius of convergence quickly shrinks with increasing the system size, while in the second case (ii) we find that the convergence radius is always larger than $1$, and in particular it approaches $1$ from above as we change the anisotropy from easy plane ($XY$) to easy axis (Ising) regime.


title: Spatial Filtering for EEG-Based Regression Problems in Brain-Computer Interface (BCI)

abstract: Electroencephalogram (EEG) signals are frequently used in brain-computer interfaces (BCIs), but they are easily contaminated by artifacts and noises, so preprocessing must be done before they are fed into a machine learning algorithm for classification or regression. Spatial filters have been widely used to increase the signal-to-noise ratio of EEG for BCI classification problems, but their applications in BCI regression problems have been very limited. This paper proposes two common spatial pattern (CSP) filters for EEG-based regression problems in BCI, which are extended from the CSP filter for classification, by making use of fuzzy sets. Experimental results on EEG-based response speed estimation from a large-scale study, which collected 143 sessions of sustained-attention psychomotor vigilance task data from 17 subjects during a 5-month period, demonstrate that the two proposed spatial filters can significantly increase the EEG signal quality. When used in LASSO and k-nearest neighbors regression for user response speed estimation, the spatial filters can reduce the root mean square estimation error by 10.02-19.77%, and at the same time increase the correlation to the true response speed by 19.39-86.47%.


title: Recent Advances in Local Energy Trading in the Smart Grid Based on Game--Theoretic Approaches

abstract: The global move towards efficient energy consumption and production has led to remarkable advancements in the design of the smart grid infrastructure. Local energy trading is one way forward. It typically refers to the transfer of energy from an entity of the smart grid surplus energy to one with a deficit. In this paper, we present a detailed review of the recent advances in the application of game--theoretic methods to local energy trading scenarios. An extensive description of a complete game theory-based energy trading framework is presented. It includes a taxonomy of the methods and an introduction to the smart grid architecture with a focus on renewable energy generation and energy storage. Finally, we present a critical evaluation of the current shortcomings and identify areas for future research.


title: Point-contact spectroscopy of tantalum, molybdenum, and tungsten

abstract: Electron-phonon interaction in superconducting transition metals (Ta, Mo, W) is investigated by the method of point-contact spectroscopy. The spectra obtained are compared with model densities of the phonon states $F(\omega)$ computed from experiments on neutron scattering, and with the tunnel electron-phonon interaction function $g(\omega)$ for Ta. Close agreement between the appropriate spectra is obtained for Ta and Mo, while a significant discrepancy is observed for W, which is the appearance of a blurred maximum instead of the two expected separated peaks on the point-contact spectra. The absolute intensities of a number of point-contact spectra $G(eV)$ and electron-phonon coupling constants $\lambda$ are calculated. For all three metals the values obtained turned out to be considerably below those expected. The possible reasons for this disagreement are discussed.


title: Topological first-order vortices in a gauged CP(2) model

abstract: We study time-independent radially symmetric first-order solitons in a CP(2) model interacting with an Abelian gauge field whose dynamics is controlled by the usual Maxwell term. In this sense, we develop a consistent first-order framework verifying the existence of a well-defined lower bound for the corresponding energy. We saturate such a lower bound by focusing on those solutions satisfying a particular set of coupled first-order differential equations. We solve these equations numerically using appropriate boundary conditions giving rise to regular structures possessing finite-energy. We also comment the main features these configurations exhibit. Moreover, we highlight that, despite the different solutions we consider for an auxiliary function $\beta \left( r\right) $ labeling the model (therefore splitting our investigation in two a priori distinct branches), all resulting scenarios engender the very same phenomenology, being physically equivalent.


title: The Geometry of Strong Koszul Algebras

abstract: Koszul algebras with quadratic Groebner bases, called strong Koszul algebras, are studied. We introduce affine algebraic varieties whose points are in one-to-one correspondence with certain strong Koszul algebras and we investigate the connection between the varieties and the algebras.


title: Connection probabilities for conformal loop ensembles

abstract: The goal of the present paper is to explain, based on properties of the conformal loop ensembles CLE$_\kappa$ (both with simple and non-simple loops, i.e., for the whole range $\kappa \in (8/3, 8)$) how to derive the connection probabilities in conformal rectangles for a conditioned version of CLE$_\kappa$ which can be interpreted as a CLE$_{\kappa}$ with wired/free/wired/free boundary conditions on four boundary arcs (the wired parts being viewed as portions of to-be-completed loops). In particular, in the case of a conformal square, we prove that the probability that the two wired sides hook up so that they create one single loop is equal to $1/(1 - 2 \cos (4 \pi / \kappa ))$. Comparing this with the corresponding connection probabilities for discrete O(N) models for instance indicates that if a dilute O(N) model (respectively a critical FK(q)-percolation model on the square lattice) has a non-trivial conformally invariant scaling limit, then necessarily this scaling limit is CLE$_\kappa$ where $\kappa$ is the value in $(8/3, 4]$ such that $-2 \cos (4 \pi / \kappa )$ is equal to $N$ (resp. the value in $[4,8)$ such that $-2 \cos (4\pi / \kappa)$ is equal to $\sqrt {q}$). Our arguments and computations build on the one hand on Dub\'edat's SLE commutation relations (as developed and used by Dub\'edat, Zhan or Bauer-Bernard-Kyt\"ol\"a) and on the other hand, on the construction and properties of the conformal loop ensembles and their relation to Brownian loop-soups, restriction measures, and the Gaussian free field, as recently derived in works with Sheffield and with Qian.


title: Self-regulation in continuum population models

abstract: We study the Markov dynamics of an infinite birth-and-death system of point entities placed in $\mathbb{R}^d$, in which the constituents disperse and die, also due to competition. Assuming that the dispersal and competition kernels are just continuous and integrable we prove that the evolution of states of this model preserves their sub-Poissonicity, and hence the local self-regulation (suppression of clustering) takes place. Upper bounds for the correlation functions of all orders are also obtained for both long and short dispersals, and for all values of the intrinsic mortality rate $m\geq 0$.


title: Gravitational Coupling from Active Gravity

abstract: We attempt to construct a gravitational coupling by pre-selecting an energy-momentum tensor as the source for gravitational field. The energy-momentum tensor we take is a recently derived new expression motivated by joint localization of energy and momentum in quantum measurement. This energy-momentum tensor differs from the traditional canonical and symmetric ones, and the theory we obtain is of an Einstein-Cartan type, but derived from a minimal coupling of a Lagrangian with second-derivative, and leads to additional interaction between torsion and matter, including the scalar field. For the scalar field, the theory can also be derived in the Riemann space-time by a non-minimal coupling. Our study gives hint on more general tests of general relativistic effects.


title: Constraints on Pulsed Emission Model for Repeating FRB 121102

abstract: Recent localization of the repeating Fast Radio Burst (FRB) 121102 revealed the distance of its host galaxy and luminosities of the bursts. We investigated constraints on the young neutron star (NS) model, that (a) the FRB intrinsic luminosity is supported by the spin-down energy, and (b) the FRB duration is shorter than the NS rotation period. In the case of a circular cone emission geometry, conditions (a) and (b) determine the NS parameters within very small ranges, compared with that from only condition (a) discussed in previous works. Anisotropy of the pulsed emission does not affect the area of the allowed parameter region by virtue of condition (b). The determined parameters are consistent with those independently limited by the properties of the possible persistent radio counterpart and the circumburst environments such as surrounding materials. Since the NS in the allowed parameter region is older than the spin-down timescale, the hypothetical GRP-like model expects a rapid radio flux decay of $\lesssim1$ Jy within a few years as the spin-down luminosity decreases. The continuous monitoring will give a hint of discrimination of the models. If no flux evolution will be seen, we need to consider an alternative model, e.g., the magnetically powered flare.


title: Smoothing expansion rate data to reconstruct cosmological matter perturbations

abstract: The existing degeneracy between different dark energy and modified gravity cosmologies at the background level may be broken by analysing quantities at the perturbative level. In this work, we apply a non-parametric smoothing (NPS) method to reconstruct the expansion history of the Universe ($H(z)$) from model-independent cosmic chronometers and high-$z$ quasar data. Assuming a homogeneous and isotropic flat universe and general relativity (GR) as the gravity theory, we calculate the non-relativistic matter perturbations in the linear regime using the $H(z)$ reconstruction and realistic values of $\Omega_{m0}$ and $\sigma_8$ from Planck and WMAP-9 collaborations. We find a good agreement between the measurements of the growth rate and $f\sigma_8(z)$ from current large-scale structure observations and the estimates obtained from the reconstruction of the cosmic expansion history. Considering a recently proposed null test for GR using matter perturbations, we also apply the NPS method to reconstruct $f\sigma_8(z)$. For this case, we find a $\sim 2\sigma$ tension (good agreement) with the standard relativistic cosmology when the Planck (WMAP-9) priors are used.


title: Two-Higgs-Doublet-Portal Dark-Matter Models in Light of Direct Search and LHC Data

abstract: We explore simple Higgs-portal models of dark matter (DM) with spin 1/2, 3/2, and 1, respectively, applying to them constraints from the LUX and PandaX-II direct detection experiments and from LHC measurements on the 125-GeV Higgs boson. With only one Higgs doublet, we find that the spin-1/2 DM having a purely scalar effective coupling to the doublet is viable only in a narrow range of mass near the Higgs pole, whereas the vector DM is still allowed if its mass is also close to the Higgs pole or exceeds 1.4 TeV, both in line with earlier analyses. Moreover, the spin-3/2 DM is in a roughly similar situation to the spin-1/2 DM, but has surviving parameter space which is even more restricted. We also consider the two-Higgs-doublet extension of each of the preceding models, assuming that the expanded Yukawa sector is that of the two-Higgs-doublet model of type II. We show that in these two-Higgs-doublet-portal models significant portions of the DM mass regions excluded in the simplest scenarios by direct search bounds can be reclaimed due to suppression of the effective DM interactions with nucleons at some ratios of the $CP$-even Higgs bosons' couplings to the up and down quarks. The regained parameter space contains areas which can yield a DM-nucleon scattering cross-section that is far less than its current experimental limit or even goes below the neutrino-background floor.


title: EAC-Net: A Region-based Deep Enhancing and Cropping Approach for Facial Action Unit Detection

abstract: In this paper, we propose a deep learning based approach for facial action unit detection by enhancing and cropping the regions of interest. The approach is implemented by adding two novel nets (layers): the enhancing layers and the cropping layers, to a pretrained CNN model. For the enhancing layers, we designed an attention map based on facial landmark features and applied it to a pretrained neural network to conduct enhanced learning (The E-Net). For the cropping layers, we crop facial regions around the detected landmarks and design convolutional layers to learn deeper features for each facial region (C-Net). We then fuse the E-Net and the C-Net to obtain our Enhancing and Cropping (EAC) Net, which can learn both feature enhancing and region cropping functions. Our approach shows significant improvement in performance compared to the state-of-the-art methods applied to BP4D and DISFA AU datasets.


title: The Word Problem of $\mathbb{Z}^n$ Is a Multiple Context-Free Language

abstract: The \emph{word problem} of a group $G = \langle \Sigma \rangle$ can be defined as the set of formal words in $\Sigma^*$ that represent the identity in $G$. When viewed as formal languages, this gives a strong connection between classes of groups and classes of formal languages. For example, Anisimov showed that a group is finite if and only if its word problem is a regular language, and Muller and Schupp showed that a group is virtually-free if and only if its word problem is a context-free language. Above this, not much was known, until Salvati showed recently that the word problem of $\mathbb{Z}^2$ is a multiple context-free language, giving first such example. We generalize Salvati's result to show that the word problem of $\mathbb{Z}^n$ is a multiple context-free language for any $n$.


title: Measurement of $W^{\pm}$ single spin asymmetries and $W$ cross section ratio in polarized $p+p$ collisions at $\sqrt{s}=510$ GeV at STAR

abstract: We present the preliminary results of measurements of single spin asymmetries, $A_L$ for $W^\pm$ boson production in longitudinally polarized p+p collisions at $\sqrt{s}=510$ GeV and measurements of cross section ratios, $\sigma _{W^+} / \sigma _{W^-}$, for $W^+,W^-$ boson production in $p+p$ collisions at $\sqrt{s}$ = 500 and 510 GeV. The asymmetry measurements are based on 246.2 pb$^{-1}$ of data taken in the RHIC 2013 run and the cross section ratio measurements are based on 102 pb$^{-1}$ of data taken during RHIC 2011 and 2012 runs by the STAR experiment. While the asymmetry results are shown as a function of the decay lepton pseudorapidity, $\eta_e$, the cross section results are shown as a function of both $\eta_e$ and W boson rapidity, $y_W$ in the mid-rapidity region ($|\eta_e|<1$). At this kinematics, $W^\pm$ single spin asymmetries provide a theoretically clean probe of the proton's polarized quark and antiquark distributions and the W cross section ratio provides sensitivity to unpolarized sea quark distributions at the scale of the W mass. The asymmetry results are consistent with recently published STAR $A^{W^\pm}_L$ results based on the data collected during RHIC 2011 and 2012 runs which showed a preference for a sizable, positive up antiquark polarization in the range $0.05 < x < 0.2$. The new preliminary results can be considered as the most precise results of $A^{W^\pm}_L$ in the world to date, with uncertainties reduced by 40\% in comparison to the published results.


title: Equivariant Weiss Calculus and Loop Spaces of Stiefel Manifolds

abstract: In the mid 1980s, Steve Mitchell and Bill Richter produced a filtration of the Stiefel manifolds $O(V; W)$ and $U(V; W)$ of orthogonal and unitary, respectively, maps $V \to V \oplus W$ stably split as a wedge sum of Thom spaces defined over Grassmanians. Additionally, they produced a similar filtration for loops on $SU(V)$, with a similar splitting. A few years later, Michael Crabb made explicit the equivariance of the Stiefel manifold splittings and conjectured that the splitting of the loop space was equivariant as well. However, it has long been unknown whether the loop space of the real Stiefel manifold (or even the special case of $\Omega SO_n$) has a similar splitting. Here, inspired by the work of Greg Arone that made use of Weiss' orthogonal calculus to generalize the results of Mitchell and Richter, we obtain an $\mathbb{Z}/2 \mathbb{Z}$-equivariant splitting theorem using an equivariant version of Weiss calculus. In particular, we show that $\Omega U(V; W)$ has an equivariant stable splitting when $dim(W) > 0$. By considering the (geometric) fixed points of this loop space, we also obtain, as a corollary, a stable splitting of the space $\Omega(U(V; W),O(V_{\mathbb{R}}; W_{\mathbb{R}}))$ of paths in $U(V; W)$ from $I$ to a point of $O(V_{\mathbb{R}}; W_{\mathbb{R}})$ as well. In particular, by setting $W = \mathbb{C}$, this gives us a stable splitting of $\Omega(SU_n / SO_n)$.


title: Double-$\beta$ Decay Matrix Elements from Lattice Quantum Chromodynamics

abstract: A lattice quantum chromodynamics (LQCD) calculation of the nuclear matrix element relevant to the $nn\to ppee\overline{\nu}_e\overline{\nu}_e$ transition is described in detail, expanding on the results presented in Ref. [1]. This matrix element, which involves two insertions of the weak axial current, is an important input for phenomenological determinations of double-$\beta$ decay rates of nuclei. From this exploratory study, performed using unphysical values of the quark masses, the long-distance deuteron-pole contribution to the matrix element is separated from shorter-distance hadronic contributions. This polarizability, which is only accessible in double-weak processes, cannot be constrained from single-$\beta$ decay of nuclei, and is found to be smaller than the long-distance contributions in this calculation, but non-negligible. In this work, technical aspects of the LQCD calculations, and of the relevant formalism in the pionless effective field theory, are described. Further calculations of the isotensor axial polarizability, in particular near and at the physical values of the light-quark masses, are required for precise determinations of both two-neutrino and neutrinoless double-$\beta$ decay rates in heavy nuclei.


title: Microscopic origin of ideal conductivity in integrable quantum models

abstract: Non-ergodic dynamical systems display anomalous transport properties. A prominent example are integrable quantum systems, whose exceptional property are diverging DC conductivities. In this Letter, we explain the microscopic origin of ideal conductivity by resorting to the thermodynamic particle content of a system. Using group-theoretic arguments we rigorously resolve the long-standing controversy regarding the nature of spin and charge Drude weights in the absence of chemical potentials. In addition, by employing a hydrodynamic description, we devise an efficient computational method to calculate exact Drude weights from the stationary currents generated in an inhomogeneous quench from bi-partitioned initial states. We exemplify the method on the anisotropic Heisenberg model at finite temperatures for the entire range of anisotropies, accessing regimes which are out of reach with other approaches. Quite remarkably, spin Drude weight and asymptotic spin current rates reveal a completely discontinuous (fractal) dependence on the anisotropy parameter.


title: Decoupling mixed finite elements on hierarchical triangular grids for parabolic problems

abstract: In this paper, we propose a numerical method for the solution of time-dependent flow problems in mixed form. Such problems can be efficiently approximated on hierarchical grids, obtained from an unstructured coarse triangulation by using a regular refinement process inside each of the initial coarse elements. If these elements are considered as subdomains, we can formulate a non-overlapping domain decomposition method based on the lowest-order Raviart-Thomas elements, properly enhanced with Lagrange multipliers on the boundaries of each subdomain (excluding the Dirichlet edges). A suitable choice of mixed finite element spaces and quadrature rules yields a cell-centered scheme for the pressures with a local 10-point stencil. The resulting system of differential-algebraic equations is integrated in time by the Crank-Nicolson method, which is known to be a stiffly accurate scheme. As a result, we obtain independent subdomain linear systems that can be solved in parallel. The behaviour of the algorithm is illustrated on a variety of numerical experiments.


title: Evolution of low mass stars: lithium problem and alpha-enhanced tracks and isochrones (Ph.D. thesis)

abstract: Precise studies on the Galactic bulge, halo, thick disk, and globular clusters require stellar models with alpha enhancement and various He contents. For this purpose we complement existing PARSEC (PAdova-TRieste Stellar Evolution Code) models, which are based on the solar partition, with alpha-enhanced partitions. We collect detailed measurements on the metal mixture and He abundance for the two populations of 47Tuc from literature, and calculate stellar tracks and isochrones with these chemical compositions ([alpha/Fe]~0.4 and 0.2). By fitting precise CMD with HST data from low main sequence till horizontal branch, we calibrate free parameters that are important for the evolution of low mass stars. The new calibration and alpha-enhanced mixtures are applied to other metallicities. The new models reproduce the RGB bump observations much better than previous models. PARSEC is able to predict any chemical pattern in stellar evolution. Li is one of the most intriguing ones. Li abundance derived in POP II dwarfs is about three times lower than the standard Big Bang Nucleosynthesis predicts. This disagreement is referred as the cosmological lithium problem. We reconsider the stellar Li evolution from pre-main sequence to the end of main sequence by introducing the effects of convective overshooting and residual accretion. $^7$Li could be significantly depleted by convective overshooting in the pre-main sequence and then partially restored in the stellar atmosphere by a tail of matter accretion which could be regulated by EUV photo-evaporation. By considering the conventional nuclear burning and microscopic diffusion along the main sequence, we reproduce the Spite plateau for $m_0=0.62-0.80M_{\odot}$ stars, and the Li declining branch for lower mass dwarfs, with a wide range of metallicities (Z=0.00001 to Z=0.0005), starting from the primordial abundance A(Li)=2.72.


title: Trapped imbalanced fermionic superfluids in one dimension: A variational approach

abstract: We propose and analyze a variational wave function for a population-imbalanced one-dimensional Fermi gas that allows for Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) type pairing correlations among the two fermion species, while also accounting for the harmonic confining potential. In the strongly interacting regime, we find large spatial oscillations of the order parameter, indicative of an FFLO state. The obtained density profiles versus imbalance are consistent with recent experimental results as well as with theoretical calculations based on combining Bethe ansatz with the local density approximation. Although we find no signature of the FFLO state in the densities of the two fermion species, we show that the oscillations of the order parameter appear in density-density correlations, both in-situ and after free expansion. Furthermore, above a critical polarization, the value of which depends on the interaction, we find the unpaired Fermi-gas state to be energetically more favorable.


title: Comprehensive Survey of Evolutionary Morphological Soft Robotic Systems

abstract: Evolutionary robotics aims to automatically design autonomous adaptive morphological robots that can evolve to accomplish a specific task while adapting to environmental changes. Soft robotics have demonstrated the feasibility of evolutionary robotics for the synthesis of robots control and morphology. The motivation of developing evolutionary soft computing techniques to that can generate task oriented structures for morphological robots makes the domain of soft robotics worthy of serious investigation and research, and hence this article summarizes an important volume of research for a computational and software architecture perspective. This paper reviews the literature and discusses various aspects of evolutionary robotics including the application on morphological soft robots to allow self assembly, self reconfiguration, self repair, and self reproduction. Then, major milestones are outlined along with important morphological soft robotic prototypes due to their importance in the field. Finally, the current state of the art in the field is assessed.


title: Quasiparticle energy in a strongly interacting homogeneous Bose-Einstein condensate

abstract: Using two-photon Bragg spectroscopy, we study the energy of particle-like excitations in a strongly interacting homogeneous Bose-Einstein condensate, and observe dramatic deviations from Bogoliubov theory. In particular, at large scattering length $a$ the shift of the excitation resonance from the free-particle energy changes sign from positive to negative. For an excitation with wavenumber $q$, this sign change occurs at $a \approx 4/(\pi q)$, in agreement with the Feynman energy relation and the static structure factor expressed in terms of the two-body contact. For $a \gtrsim 3/q$ we also see a breakdown of this theory, and better agreement with calculations based on the Wilson operator product expansion. Neither theory explains our observations across all interaction regimes, inviting further theoretical efforts.

categories: cond-mat.quant-gas cond-mat.stat-mech physics.atom-ph quant-ph


title: A bijective proof of Macdonald's reduced word formula

abstract: We give a bijective proof of Macdonald's reduced word identity using pipe dreams and Little's bumping algorithm. This proof extends to a principal specialization due to Fomin and Stanley. Such a proof has been sought for over 20 years. Our bijective tools also allow us to solve a problem posed by Fomin and Kirillov from 1997 using work of Wachs, Lenart, Serrano and Stump. These results extend earlier work by the third author on a Markov process for reduced words of the longest permutation.


title: A Generalization of Permanent Inequalities and Applications in Counting and Optimization

abstract: A polynomial $p\in\mathbb{R}[z_1,\dots,z_n]$ is real stable if it has no roots in the upper-half complex plane. Gurvits's permanent inequality gives a lower bound on the coefficient of the $z_1z_2\dots z_n$ monomial of a real stable polynomial $p$ with nonnegative coefficients. This fundamental inequality has been used to attack several counting and optimization problems. Here, we study a more general question: Given a stable multilinear polynomial $p$ with nonnegative coefficients and a set of monomials $S$, we show that if the polynomial obtained by summing up all monomials in $S$ is real stable, then we can lowerbound the sum of coefficients of monomials of $p$ that are in $S$. We also prove generalizations of this theorem to (real stable) polynomials that are not multilinear. We use our theorem to give a new proof of Schrijver's inequality on the number of perfect matchings of a regular bipartite graph, generalize a recent result of Nikolov and Singh, and give deterministic polynomial time approximation algorithms for several counting problems.

categories: cs.DS cs.DM cs.IT math.CO math.IT math.PR


title: Inflection-point inflation in hyper-charge oriented U(1)$_X$ model

abstract: Inflection-point inflation is an interesting possibility to realize a successful slow-roll inflation when inflation is driven by a single scalar field with its value during inflation below the Planck mass ($\phi_I \lesssim M_{Pl}$). In order for a renormalization group (RG) improved effective $\lambda \phi^4$ potential to develop an inflection-point, the running quartic coupling $\lambda(\phi)$ must exhibit a minimum with an almost vanishing value in its RG evolution, namely $\lambda(\phi_I) \simeq 0$ and $\beta_{\lambda}(\phi_I) \simeq 0$, where $\beta_{\lambda}$ is the beta-function of the quartic coupling. In this paper, we consider the inflection-point inflation in the context of the minimal gauged U(1)$_X$ extended Standard Model (SM), which is a generalization of the minimal U(1)$_{B-L}$ model, and is constructed as a linear combination of the SM U(1)$_Y$ and U(1)$_{B-L}$ gauge symmetries. We identify the U(1)$_X$ Higgs field with the inflaton field. For a successful inflection-point inflation to be consistent with the current cosmological observations, the mass ratios among the U(1)$_X$ gauge boson, the right-handed neutrinos and the U(1)$_X$ Higgs boson are fixed. Focusing on the case that the extra U(1)$_X$ gauge symmetry is mostly aligned along the SM U(1)$_Y$ direction, we investigate a consistency between the inflationary predictions and the latest LHC Run-2 results on the search for a narrow resonance with the di-lepton final state. %In addition, the inflection-point inflation provides a unique prediction for the running of the spectral index $\alpha \simeq - 2.7 \times %10^{-3}\left(\frac{60}{N}\right)^2$ ($N$ is the e-folding number), which can be tested in the near future.


title: cellGPU: massively parallel simulations of dynamic vertex models

abstract: Vertex models represent confluent tissue by polygonal or polyhedral tilings of space, with the individual cell interacting via force laws that depend on both the geometry of the cells and the topology of the tessellation. This dependence on the connectivity of the cellular network introduces several complications to performing molecular-dynamics-like simulations of vertex models, and in particular makes parallelizing the simulations difficult. cellGPU addresses this difficulty and lays the foundation for massively parallelized, GPU-based simulations of these models. This article discusses its implementation for a pair of two-dimensional models, and compares the typical performance that can be expected between running cellGPU entirely on the CPU versus its performance when running on a range of commercial and server-grade graphics cards. By implementing the calculation of topological changes and forces on cells in a highly parallelizable fashion, cellGPU enables researchers to simulate time- and length-scales previously inaccessible via existing single-threaded CPU implementations.


title: Stark profiles modeling of radiation lines originating from atomic autoionizing states in dense plasmas, solid state matter under short intense XUV/X-ray free electron laser irradiation

abstract: In the present paper we propose to update the standard calculation of Stark line profiles for autoionizing atomic states in dense plasmas, i.e. warm dense matter and strongly coupled plasmas. This is motivated by the importance of taking into account the effect of the electric microfield generated by the plasma charge constituents in the calculation of the atomic populations of autoionizing states which leads to the modification of the calculation of Stark line profiles in the frame of the standard theory. This is due to the properties of the autoionizing atomic states that are far from Boltzmann distribution even in dense plasmas bringing serious doubts on the statistical approach used in the standard theory of Stark profile calculations. We discuss the importance of the present analysis in view of the radiation emission originating from dense plasmas created by the interaction of the X-ray Free Electron Lasers (XFEL's) with solid density matter.


title: Investigating the temperature dependence of the specific shear viscosity of QCD matter with dilepton radiation

abstract: This work reports on investigations of the effects on the evolution of viscous hydrodynamics and on the flow coefficients of thermal dileptons, originating from a temperature-dependent specific shear viscosity $\eta/s (T)$ at temperatures beyond 180 MeV formed at the Relativistic Heavy-Ion Collider (RHIC). We show that the elliptic flow of thermal dileptons can resolve the magnitude of $\eta/s$ at the high temperatures, where partonic degrees of freedom become relevant, whereas discriminating between different specific functional forms will likely not be possible at RHIC using this observable.


title: Nonequilibrium steady states and transient dynamics of conventional superconductors under phonon driving

abstract: We perform a systematic analysis of the influence of phonon driving on the superconducting Holstein model coupled to heat baths by studying both the transient dynamics and the nonequilibrium steady state (NESS) in the weak and strong electron-phonon coupling regimes. Our study is based on the nonequilibrium dynamical mean-field theory, and for the NESS we present a Floquet formulation adapted to electron-phonon systems. The analysis of the phonon propagator suggests that the effective attractive interaction can be strongly enhanced in a parametric resonant regime because of the Floquet side bands of phonons. While this may be expected to enhance the superconductivity (SC), our fully self-consistent calculations, which include the effects of heating and nonthermal distributions, show that the parametric phonon driving generically results in a suppression or complete melting of the SC order. In the strong coupling regime, the NESS always shows a suppression of the SC gap, the SC order parameter and the superfluid density as a result of the driving, and this tendency is most prominent at the parametric resonance. Using the real-time nonequilibrium DMFT formalism, we also study the dynamics towards the NESS, which shows that the heating effect dominates the transient dynamics, and SC is weakened by the external modulations, in particular at the parametric resonance. In the weak coupling regime, we find that the SC fluctuations above the transition temperature are generally weakened under the driving. The strongest suppression occurs again around the parametric resonances because of the efficient energy absorption.


title: Time-Series Analysis of Broadband Photometry of Neptune from K2

abstract: We report here on our search for excess power in photometry of Neptune collected by the K2 mission that may be due to intrinsic global oscillations of the planet Neptune. To conduct this search, we developed new methods to correct for instrumental effects such as intrapixel variability and gain variations. We then extracted and analyzed the time-series photometry of Neptune from 49 days of nearly continuous broadband photometry of the planet. We find no evidence of global oscillations and place an upper limit of $\sim$5 ppm at 1000 \uhz\ for the detection of a coherent signal. With an observed cadence of 1-minute and point-to-point scatter less than 0.01\%, the photometric signal is dominated by reflected light from the Sun, which is in turn modulated by atmospheric variability of Neptune at the 2\% level. A change in flux is also observed due to the increasing distance between Neptune and the K2 spacecraft, and solar variability with convection-driven solar p modes present.


title: Dark Matter "Transporting" Mechanism Explaining Positron Excesses

abstract: We propose a novel mechanism to explain the positron excesses, which are observed by satellite-based telescopes including PAMELA and AMS-02, in dark matter (DM) scenarios. The novelty behind the proposal is that it makes direct use of DM around the Galactic Center where DM populates most densely, allowing us to avoid tensions from cosmological and astrophysical measurements. The key ingredients of this mechanism include DM annihilation into unstable states with a very long laboratory-frame life time and their "retarded" decay near the Earth to electron-positron pair(s) possibly with other (in)visible particles. We argue that this sort of explanation is not in conflict with relevant constraints from big bang nucleosynthesis and cosmic microwave background. Regarding the resultant positron spectrum, we provide a generalized source term in the associated diffusion equation, which can be readily applicable to any type of two-"stage" DM scenarios wherein production of Standard Model particles occurs at completely different places from those of DM annihilation. We then conduct a data analysis with the recent AMS-02 data to validate our proposal.


title: Molecular gas at supernova local environments unveiled by EDGE

abstract: CO observations allow estimations of the gas content of molecular clouds, which trace the reservoir of cold gas fuelling star formation, as well as to determine extinction via H$_2$ column density, N(H$_2$). Here, we studied millimetric and optical properties at 26 supernovae (SNe) locations of different types in a sample of 23 nearby galaxies by combining molecular $^{12}$C$^{16}$O (J = 1 $\rightarrow$ 0) resolved maps from the EDGE survey and optical Integral Field Spectroscopy from the CALIFA survey. We found an even clearer separation between type II and type Ibc SNe in terms of molecular gas than what we found in the optical using H$\alpha$ emission as a proxy for current SF rate, which reinforces the fact that SNe Ibc are more associated with SF-environments. While A$_V$ at SN locations is similar for SNe II and SNe Ibc, and higher compared to SNe Ia, N(H$_2$) is significantly higher for SNe Ibc than for SNe II and SNe Ia. When compared to alternative extinction estimations directly from SN photometry and spectroscopy, we find that our SNe Ibc have also redder color excess but showed standard Na I D absorption pseudo-equivalent widths ($\sim$1 \AA). In some cases we find no extinction when estimated from the environment, but high amounts of extinction when measured from SN observations, which suggests that circumstellar material or dust sublimation may be playing a role. This work serves as a benchmark for future studies combining last generation millimeter and optical IFS instruments to reveal the local environmental properties of extragalactic SNe.


title: The global and small-scale magnetic fields of fully convective, rapidly spinning M dwarf pair GJ65A and B

abstract: The nearby M dwarf binary GJ65AB, also known as BL Cet and UV Cet, is a unique benchmark for investigation of dynamo-driven activity of low-mass stars. Magnetic activity of GJ65 was repeatedly assessed by indirect means, such as studies of flares, photometric variability, X-ray and radio emission. Here we present a direct analysis of large-scale and local surface magnetic fields in both components. Interpreting high-resolution circular polarisation spectra (sensitive to a large-scale field geometry) we uncovered a remarkable difference of the global stellar field topologies. Despite nearly identical masses and rotation rates, the secondary exhibits an axisymmetric, dipolar-like global field with an average strength of 1.3 kG while the primary has a much weaker, more complex and non-axisymmetric 0.3 kG field. On the other hand, an analysis of the differential Zeeman intensification (sensitive to the total magnetic flux) shows the two stars having similar magnetic fluxes of 5.2 and 6.7 kG for GJ65A and B, respectively, although there is an evidence that the field strength distribution in GJ65B is shifted towards a higher field strength compared to GJ65A. Based on these complementary magnetic field diagnostic results we suggest that the dissimilar radio and X-ray variability of GJ65A and B is linked to their different global magnetic field topologies. However, this difference appears to be restricted to the upper atmospheric layers but does not encompass the bulk of the stars and has no influence on the fundamental stellar properties.


title: Disentangling Heavy Flavor at Colliders

abstract: We propose two new analysis strategies for studying charm and beauty quarks at colliders. The first strategy is aimed at testing the kinematics of heavy-flavor quarks within an identified jet. Here, we use the SoftDrop jet-declustering algorithm to identify two subjets within a large-radius jet, using subjet flavor tagging to test the heavy-quark splitting functions of QCD. For subjets containing a $J / \psi$ or $\Upsilon$, this declustering technique can also help probe the mechanism for quarkonium production. The second strategy is aimed at isolating heavy-flavor production from gluon splitting. Here, we introduce a new FlavorCone algorithm, which smoothly interpolates from well-separated heavy-quark jets to the gluon-splitting regime where jets overlap. Because of its excellent ability to identify charm and beauty hadrons, the LHCb detector is ideally suited to pursue these strategies, though similar measurements should also be possible at ATLAS and CMS. Together, these SoftDrop and FlavorCone studies should clarify a number of aspects of heavy-flavor physics at colliders, and provide crucial information needed to improve heavy-flavor modeling in parton-shower generators.


title: Elliptic Genera of 2d (0,2) Gauge Theories from Brane Brick Models

abstract: We compute the elliptic genus of abelian 2d (0,2) gauge theories corresponding to brane brick models. These theories are worldvolume theories on a single D1-brane probing a toric Calabi-Yau 4-fold singularity. We identify a match with the elliptic genus of the non-linear sigma model on the same Calabi-Yau background, which is computed using a new localization formula. The matching implies that the quantum effects do not drastically alter the correspondence between the geometry and the 2d (0,2) gauge theory. In theories whose matter sector suffers from abelian gauge anomaly, we propose an ansatz for an anomaly cancelling term in the integral formula for the elliptic genus. We provide an example in which two brane brick models related to each other by Gadde-Gukov-Putrov triality give the same elliptic genus.


title: Electroweak Kaluza-Klein Dark Matter

abstract: In models with universal extra dimensions (UED), the lightest Kaluza-Klein excitation of neutral electroweak gauge bosons is a stable, weakly interacting massive particle and thus is a candidate for dark matter thanks to Kaluza-Klein parity. We examine concrete model realizations of such dark matter in the context of non-minimal UED extensions. The boundary localized kinetic terms for the electroweak gauge bosons lead to a non-trivial mixing among the first Kaluza-Klein excitations of the ${\rm SU}(2)_W$ and ${\rm U}(1)_Y$ gauge bosons and the resultant low energy phenomenology is rich. We investigate implications of various experiments including low energy electroweak precision measurements, direct and indirect detection of dark matter particles and direct collider searches at the LHC. Notably, we show that the electroweak Kaluza-Klein dark matter can be as heavy as 2.4 TeV, which is significantly higher than $1.3$ TeV as is indicated as an upper bound in the minimal UED model.


title: Second-order chiral kinetic theory: Chiral magnetic and pseudomagnetic waves

abstract: The consistent chiral kinetic theory accurate to the second order in electromagnetic and pseudoelectromagnetic fields is derived for a relativistic matter with two Weyl fermions. By making use of such a framework, the properties of longitudinal collective excitations, which include both chiral magnetic and chiral pseudomagnetic waves, are studied. It is shown that the proper treatment of dynamical electromagnetism transforms these gapless waves into chiral (pseudo)magnetic plasmons, whose Langmuir (plasma) gap receives corrections quadratic in both magnetic and pseudomagnetic fields.


title: What can the programming language Rust do for astrophysics?

abstract: The astrophysics community uses different tools for computational tasks such as complex systems simulations, radiative transfer calculations or big data. Programming languages like Fortran, C or C++ are commonly present in these tools and, generally, the language choice was made based on the need for performance. However, this comes at a cost: safety. For instance, a common source of error is the access to invalid memory regions, which produces random execution behaviors and affects the scientific interpretation of the results. In 2015, Mozilla Research released the first stable version of a new programming language named Rust. Many features make this new language attractive for the scientific community, it is open source and it guarantees memory safety while offering zero-cost abstraction. We explore the advantages and drawbacks of Rust for astrophysics by re-implementing the fundamental parts of Mercury-T, a Fortran code that simulates the dynamical and tidal evolution of multi-planet systems.


title: Glassy quantum dynamics in translation invariant fracton models

abstract: We investigate relaxation in the recently discovered "fracton" models and discover that these models naturally host glassy quantum dynamics in the absence of quenched disorder. We begin with a discussion of "type I" fracton models, in the taxonomy of Vijay, Haah, and Fu. We demonstrate that in these systems, the mobility of charges is suppressed exponentially in the inverse temperature. We further demonstrate that when a zero temperature type I fracton model is placed in contact with a finite temperature heat bath, the approach to equilibrium is a logarithmic function of time over an exponentially wide window of time scales. Generalizing to the more complex "type II" fracton models, we find that the charges exhibit subdiffusion upto a relaxation time that diverges at low temperatures as a super-exponential function of inverse temperature. This behaviour is reminiscent of "nearly localized" disordered systems, but occurs with a translation invariant three-dimensional Hamiltonian. We also conjecture that fracton models with conserved charge may support a phase which is a thermal metal but a charge insulator.

categories: cond-mat.stat-mech cond-mat.dis-nn cond-mat.str-el quant-ph


title: The Potassium abundance in the globular clusters NGC104, NGC6752 and NGC6809

abstract: We derived Potassium abundances in red giant branch stars in the Galactic globular clusters NGC104 (144 stars), NGC6752 (134 stars) and NGC6809 (151 stars) using high-resolution spectra collected with FLAMES at the ESO - Very Large Telescope. In the considered samples we do not find significant intrinsic spreads in [K/Fe] (confirming the previous findings by Carretta et al.), at variance with the cases of the massive clusters NGC2419 and NGC2808. Additionally, marginally significant [K/Fe]-[O/Fe] anti-correlations are found in NGC104 and NGC6809, and [K/Fe]-[Na/Fe] correlations are found in NGC104 and NGC6752. No evidence of [K/Fe]-[Mg/Fe] anti-correlation are found. The results of our analysis are consistent with a scenario in which the process leading to the multi-populations in globular clusters implies also enrichment in the K abundance, the amplitude of the associated [K/Fe] enhancement becoming measurable only in stars showing the most extreme effects of O and Mg depletion. Stars enhanced in [K/Fe] have been found so far only in clusters harbouring some Mg-poor stars, while the other globulars, without a Mg-poor sub-population, show small or null [K/Fe] spreads.


abstract: The non-perturbative structure of the photon and gluon propagators plays an important role in governing the dynamics of quantum electrodynamics (QED) and quantum chromodynamics (QCD) respectively. Although it is often assumed that these interacting field propagators can be decomposed into longitudinal and transverse components, as for the free case, it turns out that in general this is not possible. Moreover, the non-abelian gauge symmetry of QCD permits the momentum space gluon propagator to contain additional singular terms involving derivatives of $\delta(p)$, the appearance of which is related to confinement. Despite the possibility of the failure of the transverse-longitudinal decomposition for the photon and gluon propagators, and the appearance of singular terms in the gluon propagator, the Slavnov-Taylor identity nevertheless remains preserved.


title: Detecting transient rate-tipping using Steklov averages and Lyapunov vectors

abstract: A wide variety of physical systems ranging from the firing of neurons to eutrophication of lakes to the presence of Arctic summer sea ice exhibit a phenomenon known as tipping. In mathematical models, tipping can be caused by bifurcations, noise, and the rate at which parameters are changing in time [2]. Because traditional methods in dynamical systems are usually concerned with the long-term behavior of the system, these methods are not always able to detect the transient dynamics characteristic of rate-tipping. In this paper, we consider one- and two-dimensional dynamical systems with nonautonomous parameters that exhibit rate-tipping, as defined as not tracking the evolution of stable equilibria (QSEs) in the corresponding autonomous systems. We find that nonautonomous stability spectra in the form of Steklov averages and their derivatives appear to be correlated with transient rate-tipping in systems with unique QSEs or with parameters that change at a constant rate. Furthermore, for systems in two dimensions and higher, comparison of the angle between leading Lyapunov vectors of different trajectories admits a possible criterion for detecting rate-tipping. Our heuristic results add to the body of work dedicated to studying and understanding the phenomenon of rate-tipping.


title: Disorder induced transitions in resonantly driven Floquet Topological Insulators

abstract: We investigate the effects of disorder in Floquet topological insulators (FTIs) occurring in semiconductor quantum wells. Such FTIs are induced by resonantly driving a transition between the valence and conduction band. We show that when disorder is added, the topological nature of such FTIs persists as long as there is a mobility gap at the resonant quasi-energy. For strong enough disorder, this gap closes and all the states become localized as the system undergoes a transition to a trivial insulator. Interestingly, the effects of disorder are not necessarily adverse: we show that in the same quantum well, disorder can also induce a transition from a trivial to a topological system, thereby establishing a Floquet Topological Anderson Insulator (FTAI). We identify the conditions on the driving field necessary for observing such a transition.


title: A versatile, high-power 460nm laser system for Rydberg excitation of ultracold potassium

abstract: We present a versatile laser system which provides more than 1.5W of narrowband light, tunable in the range from 455-463 nm. It consists of a commercial Titanium-Sapphire laser which is frequency doubled using resonant cavity second harmonic generation and stabilized to an external reference cavity. We demonstrate a wide wavelength tuning range combined with a narrow linewidth and low intensity noise. This laser system is ideally suited for atomic physics experiments such as two-photon excitation of Rydberg states of potassium atoms with principal quantum numbers n > 18. To demonstrate this we perform two-photon spectroscopy on ultracold potassium gases in which we observe an electromagnetically induced transparency resonance corresponding to the 35s1/2 state and verify the long-term stability of the laser system. Additionally, by performing spectroscopy in a magneto-optical trap we observe strong loss features corresponding to the excitation of s, p, d and higher-l states accessible due to a small electric Field.


title: Lipschitz regularity of solutions to two-phase free boundary problems

abstract: We prove Lipschitz continuity of viscosity solutions to a class of two-phase free boundary problems governed by fully nonlinear operators.


title: Electron paramagnetic resonance and photochromism of $\mathrm{N}_{3}\mathrm{V}^{0}$ in diamond

abstract: The defect in diamond formed by a vacancy surrounded by three nearest-neighbor nitrogen atoms and one carbon atom, $\mathrm{N}_{3}\mathrm{V}$, is found in $\approx98\%$ of natural diamonds. Despite $\mathrm{N}_{3}\mathrm{V}^{0}$ being the earliest electron paramagnetic resonance spectrum observed in diamond, to date no satisfactory simulation of the spectrum for an arbitrary magnetic field direction has been produced due to its complexity. In this work, $\mathrm{N}_{3}\mathrm{V}^{0}$ is identified in $^{15}\mathrm{N}$-doped synthetic diamond following irradiation and annealing. The $\mathrm{^{15}N}_{3}\mathrm{V}^{0}$ spin Hamiltonian parameters are revised and used to refine the parameters for $\mathrm{^{14}N}_{3}\mathrm{V}^{0}$, enabling the latter to be accurately simulated and fitted for an arbitrary magnetic field direction. Study of $\mathrm{^{15}N}_{3}\mathrm{V}^{0}$ under excitation with green light indicates charge transfer between $\mathrm{N}_{3}\mathrm{V}$ and $\mathrm{N_s}$. It is argued that this charge transfer is facilitated by direct ionization of $\mathrm{N}_{3}\mathrm{V}^{-}$, an as-yet unobserved charge state of $\mathrm{N}_{3}\mathrm{V}$.


title: The dusty tori of nearby QSOs as constrained by high-resolution mid-IR observations

abstract: We present mid-infrared (MIR, 7.5-13.5 $\mu$m) imaging and spectroscopy observations obtained with the CanariCam (CC) instrument on the 10.4m Gran Telescopio CANARIAS for a sample of 20 nearby, MIR bright and X-ray luminous QSOs. We find that for the majority of QSOs the MIR emission is unresolved at angular scales nearly 0.3 arcsec, corresponding to physical scales $<600$ pc. We find that the higher-spatial resolution CC spectra have similar shapes to those obtained with Spitzer/IRS, and hence we can assume that the spectra are not heavily contaminated by extended emission in the host galaxy. We thus take advantage of the higher signal to noise Spitzer/IRS spectra, as a fair representation of the nuclear emission, to decompose it into a combination of active galactic nuclei (AGN), polycyclic aromatic hydrocarbon (PAH) and stellar components. In most cases the AGN is the dominant component, with a median contribution of 85 per cent of the continuum light at MIR (5-15 $\mu$m) within the IRS slit. This IR AGN emission is well reproduced by clumpy torus models. We find evidence for significant differences in the parameters that describe the dusty tori of QSOs when compared with the same parameters of Seyfert 1 and 2 nuclei. In particular, we find a lower number of clouds ($N_{0}<12$), steeper radial distribution of clouds ($q=1.5-3.0$), and clouds that are less optically thick ($\tau_{V}<100$) than in Seyfert 1, which could be attributed to dusty structures that have been partially evaporated and piled up by the higher radiation field in QSOs. We find that the combination of the angular width $\sigma_{torus}$, viewing angle $i$, and number of clouds along the equatorial line $N_{0}$, produces large escape probabilities ($P_{esc} > 2$ per cent) and low geometrical covering factors ($f_{2}<0.6$), as expected for AGN with broad lines in their optical spectra.


title: The maximum isotropic energy of gamma-ray bursts

abstract: The most energetic gamma-ray bursts (GRBs) are remarkable sources releasing huge amounts of energy on short timescales. Their prompt emission, which usually lasts few seconds, is so bright that it is visible across the whole observable universe. Studying these extreme events may provide clues on the nature of GRB progenitors and on the physical processes at work in relativistic jets. In this paper, we study the bright end of the isotropic energy distribution of long GRBs. We use two samples of long GRBs with redshift detected by Fermi/GBM or Konus-Wind, two instruments which measure the spectral shape and the energetics of the prompt emission accurately. We focus on GRBs within a range of redshifts z = 1 -- 5, a volume that contains a large number of energetic GRBs, and we propose a simple method to reconstruct the bright end of the GRB energy distribution from the observed one. We find that the GRB energy distribution cannot be described by a simple power law but requires a strong cutoff above $1-3 \times 10^{54}$ erg. We attribute this feature to an intrinsic limit on the energy per unit of solid angle radiated by gamma-ray bursts.


title: Precise deviations for Hawkes processes

abstract: Hawkes process is a class of simple point processes with self-exciting and clustering properties. Hawkes process has been widely applied in finance, neuroscience, social networks, criminology, seismology, and many other fields. In this paper, we study precise deviations for Hawkes processes for large time asymptotics, that strictly extends and improves the existing results in the literature. Numerical illustrations will also be provided.


title: Vortex motion around a circular cylinder above a plane

abstract: The study of vortex flows around solid obstacles is of considerable interest from both a theoretical and practical perspective. One geometry that has attracted renewed attention recently is that of vortex flows past a circular cylinder placed above a plane wall, where a stationary recirculating eddy can form in front of the cylinder, in contradistinction to the usual case (without the plane boundary) for which a vortex pair appears behind the cylinder. Here we analyze the problem of vortex flows past a cylinder near a wall through the lenses of the point-vortex model. By conformally mapping the fluid domain onto an annular region in an auxiliary complex plane, we compute the vortex Hamiltonian analytically in terms of certain special functions related to elliptic theta functions. A detailed analysis of the equilibria of the model is then presented. The location of the equilibrium in front of the cylinder is shown to be in qualitative agreement with the experimental findings. We also show that a topological transition occurs in phase space as the parameters of the systems are varied


title: Decay photons from the ALP burst of type-II supernovae

abstract: We determine limits from SN 1987A on massive axion-like particles (ALPs) with masses in the 10 keV - 100 MeV range and purely coupled to two photons. ALPs produced in the core collapse escape from the star and decay into photons that can be observed as a delayed and diffuse burst. We discuss the time and angular distribution of such a signal. Looking into the future we also estimate the possible improvements if the red supergiant Betelgeuse explodes in a supernova event.


title: The parabolic Anderson model on Riemann surfaces

abstract: We show well-posedness for the parabolic Anderson model on $2$-dimensional closed Riemannian manifolds. To this end we extend the notion of regularity structures to curved space, and explicitly construct the minimal structure required for this equation. A central ingredient is the appropriate re-interpretation of the polynomial model, which we build up to any order.


title: A monitoring and diagnostic approach for stochastic textured surfaces

abstract: We develop a supervised-learning-based approach for monitoring and diagnosing texture-related defects in manufactured products characterized by stochastic textured surfaces that satisfy the locality and stationarity properties of Markov random fields. Examples of stochastic textured surface data include images of woven textiles; image or surface metrology data for machined, cast, or formed metal parts; microscopy images of material microstructure samples; etc. To characterize the complex spatial statistical dependencies of in-control samples of the stochastic textured surface, we use rather generic supervised learning methods, which provide an implicit characterization of the joint distribution of the surface texture. We propose two spatial moving statistics, which are computed from residual errors of the fitted supervised learning model, for monitoring and diagnosing local aberrations in the general spatial statistical behavior of newly manufactured stochastic textured surface samples in a statistical process control context. We illustrate the approach using images of textile fabric samples and simulated 2-D stochastic processes, for which the algorithm successfully detects local defects of various natures. Supplemental discussions, results, data and computer codes are available online.


title: AC optimal power flow in the presence of renewable sources and uncertain loads

abstract: The increasing penetration of renewable energy resources, paired with the fact that load can vary significantly, introduce a high degree of uncertainty in the behavior of modern power grids. Given that classical dispatch solutions are "rigid," their performance in such an uncertain environment is in general far from optimal. For this reason, in this paper, we consider AC optimal power flow (AC-OPF) problems in the presence of uncertain loads and (uncertain) renewable energy generators. The goal of AC-OPF design is to guarantee that controllable generation is dispatched at minimum cost, while satisfying constraints on generation and transmission for "almost all" realizations of the uncertainty. We propose an approach based on a randomized technique recently developed, named "scenario with certificates", which allows to tackle the problem without assuming any a-priori dependence of the voltages in the network on the uncertain generators/loads. The proposed solution can exploit the usually available probabilistic description of the uncertainty and variability, and provides solutions with a-priori probabilistic guarantees on the risk of violating the constraints on generation and transmission.


title: Comparative benchmarking of cloud computing vendors with High Performance Linpack

abstract: We present a comparative analysis of the maximum performance achieved by the Linpack benchmark on compute intensive hardware publicly available from multiple cloud providers. We study both performance within a single compute node, and speedup for distributed memory calculations with up to 32 nodes or at least 512 computing cores. We distinguish between hyper-threaded and non-hyper-threaded scenarios and estimate the performance per single computing core. We also compare results with a traditional supercomputing system for reference. Our findings provide a way to rank the cloud providers and demonstrate the viability of the cloud for high performance computing applications.


title: Network-Cognizant Voltage Droop Control for Distribution Grids

abstract: This paper examines distribution systems with a high integration of distributed energy resources (DERs) and addresses the design of local control methods for real-time voltage regulation. Particularly, the paper focuses on proportional control strategies where the active and reactive output-powers of DERs are adjusted in response to (and proportionally to) local changes in voltage levels. The design of the voltage-active power and voltage-reactive power characteristics leverages suitable linear approximations of the AC power-flow equations and is network-cognizant; that is, the coefficients of the controllers embed information on the location of the DERs and forecasted non-controllable loads/injections and, consequently, on the effect of DER power adjustments on the overall voltage profile. A robust approach is pursued to cope with uncertainty in the forecasted non-controllable loads/power injections. Stability of the proposed local controllers is analytically assessed and numerically corroborated.


title: The Price of Selection in Differential Privacy

abstract: In the differentially private top-$k$ selection problem, we are given a dataset $X \in \{\pm 1\}^{n \times d}$, in which each row belongs to an individual and each column corresponds to some binary attribute, and our goal is to find a set of $k \ll d$ columns whose means are approximately as large as possible. Differential privacy requires that our choice of these $k$ columns does not depend too much on any on individual's dataset. This problem can be solved using the well known exponential mechanism and composition properties of differential privacy. In the high-accuracy regime, where we require the error of the selection procedure to be to be smaller than the so-called sampling error $\alpha \approx \sqrt{\ln(d)/n}$, this procedure succeeds given a dataset of size $n \gtrsim k \ln(d)$. We prove a matching lower bound, showing that a dataset of size $n \gtrsim k \ln(d)$ is necessary for private top-$k$ selection in this high-accuracy regime. Our lower bound is the first to show that selecting the $k$ largest columns requires more data than simply estimating the value of those $k$ columns, which can be done using a dataset of size just $n \gtrsim k$.


title: Chemical evolution of the Galactic bulge as traced by microlensed dwarf and subgiant stars. VI. Age and abundance structure of the stellar populations in the central sub-kpc of the Milky Way

abstract: We present a detailed elemental abundance study of 90 F and G dwarf, turn-off and subgiant stars in the Galactic bulge. Based on high-resolution spectra acquired during gravitational microlensing events, stellar ages and abundances for 11 elements (Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Zn, Y and Ba) have been determined. We find that the Galactic bulge has a wide metallicity distribution with significant peaks at [Fe/H]=-1.09, -0.63, -0.20, +0.12, +0.41. We also find a high fraction of intermediate-age to young stars: at [Fe/H]>0 more than 35 % are younger than 8 Gyr. For [Fe/H]<-0.5 most stars are 10 Gyr or older. We have also identified several episodes when significant star formation in the bulge happened: 3, 6, 8, and 12 Gyr ago. We further find that the "knee" in the alpha-element abundance trends of the sub-solar metallicity bulge is located at about 0.1 dex higher [Fe/H] than in the local thick disk. The Galactic bulge has complex age and abundance properties that appear to be tightly connected to the main Galactic stellar populations. In particular, the peaks in the metallicity distribution, the star formation episodes, and the abundance trends, show similarities with the properties of the Galactic thin and thick disks. At the same time there are additional components not seen outside the bulge region, and that most likely can be associated with the Galactic bar. For instance, the star formation rate appears to have been slightly faster in the bulge than in the local thick disk, which most likely is an indication of the denser stellar environment closer to the Galactic centre. Our results strengthen the observational evidence that support the idea of a secular origin for the Galactic bulge, formed out of the other main Galactic stellar populations present in the central regions of our Galaxy.


title: Trace Properties from Separation Logic Specifications

abstract: We propose a formal approach for relating abstract separation logic library specifications with the trace properties they enforce on interactions between a client and a library. Separation logic with abstract predicates enforces a resource discipline that constrains when and how calls may be made between a client and a library. Intuitively, this can enforce a protocol on the interaction trace. This intuition is broadly used in the separation logic community but has not previously been formalised. We provide just such a formalisation. Our approach is based on using wrappers which instrument library code to induce execution traces for the properties under examination. By considering a separation logic extended with trace resources, we prove that when a library satisfies its separation logic specification then its wrapped version satisfies the same specification and, moreover, maintains the trace properties as an invariant. Consequently, any client and library implementation that are correct with respect to the separation logic specification will satisfy the trace properties.


title: Bayesian Multiscale Finite Element Methods. Modeling missing subgrid information probabilistically

abstract: In this paper, we develop a Bayesian multiscale approach based on a multiscale finite element method. Because of scale disparity in many multiscale applications, computational models can not resolve all scales. Various subgrid models are proposed to represent un-resolved scales. Here, we consider a probabilistic approach for modeling un-resolved scales using the Multiscale Finite Element Method (cf., [1, 2]). By representing dominant modes using the Generalized Multiscale Finite Element, we propose a Bayesian framework, which provides multiple inexpensive (computable) solutions for a deterministic problem. These approximate probabilistic solutions may not be very close to the exact solutions and, thus, many realizations are needed. In this way, we obtain a rigorous probabilistic description of approximate solutions. In the paper, we consider parabolic and wave equations in heterogeneous media. In each time interval, the domain is divided into subregions. Using residual information, we design appropriate prior and posterior distributions. The likelihood consists of the residual minimization. To sample from the resulting posterior distribution, we consider several sampling strategies. The sampling involves identifying important regions and important degrees of freedom beyond permanent basis functions, which are used in residual computation. Numerical results are presented. We consider two sampling algorithms. The first algorithm uses sequential sampling and is inexpensive. In the second algorithm, we perform full sampling using the Gibbs sampling algorithm, which is more accurate compared to the sequential sampling. The main novel ingredients of our approach consist of: defining appropriate permanent basis functions and the corresponding residual; setting up a proper posterior distribution; and sampling the posteriors.


title: Experimental Constraint on an Exotic Spin- and Velocity-Dependent Interaction in the Sub-meV Range of Axion Mass with a Spin-Exchange Relaxation-Free Magnetometer

abstract: We conducted a search for an exotic spin- and velocity-dependent interaction for polarized electrons with an experimental approach based on a high-sensitivity spin-exchange relaxation-free (SERF) magnetometer, which serves as both a source of polarized electrons and a magnetic-field sensor. The experiment aims to sensitively detect magnetic-fieldlike effects from the exotic interaction between the polarized electrons in a SERF vapor cell and unpolarized nucleons of a closely located solid-state mass. We report experimental results on the interaction with 82 h of data averaging, which sets an experimental limit on the coupling strength around $10^{-19}$ for the axion mass $m_a \lesssim 10^{-3}$ eV, within the important axion window.


title: Angular Momentum in Disk Wind Revealed in the Young Star MWC349A

abstract: Disk winds are thought to play a critical role in star birth. As winds extract excess angular momentum from accretion disks, matter in the disk can be transported inward to the star to fuel mass growth. However, the observational evidence of wind carrying angular momentum has been very limited. We present Submillimeter Array (SMA) observations of the young star MWC349A in the H26$\alpha$ and H30$\alpha$ recombination lines. The high signal-to-noise ratios made possible by the maser emission process allow us to constrain the relative astrometry of the maser spots to a milli-arcsecond precision. Previous observations of the H30$\alpha$ line with the SMA and the Plateau de Bure interferometer (PdBI) showed that masers are distributed in the disk and wind. Our new high resolution observations of the H26$\alpha$ line reveal differences in spatial distribution from that of the H30$\alpha$ line. H26$\alpha$ line masers in the disk are excited in a thin annulus with a radius of about 25 AU, while the H30$\alpha$ line masers are formed in a slightly larger annulus with a radius of 30 AU. This is consistent with expectations for maser excitation in the presence of an electron density variation of approximately $R^{-4}$. In addition, the H30$\alpha$ and H26$\alpha$ line masers arise from different parts in the wind. This difference is also expected from the maser theory. The wind component of both masers exhibits line-of-sight velocities that closely follow a Keplerian law. This result provides strong evidence that the disk wind extracts significant angular momentum and thereby facilitating mass accretion in the young star.


title: Asymptotic approximations of the solution to a boundary-value problem in a thin aneurysm-type domain

abstract: A nonuniform Neumann boundary-value problem is considered for the Poisson equation in a thin $3D$ aneurysm-type domain that consists of thin curvilinear cylinders that are joined through an aneurysm of diameter $\mathcal{O}(\varepsilon).$ A rigorous procedure is developed to construct the complete asymptotic expansion for the solution as the parameter $\varepsilon \to 0.$ The asymptotic expansion consists of a regular part that is located inside of each cylinder, a boundary-layer part near the base of each cylinder, and an inner part discovered in a neighborhood of the aneurysm. Terms of the inner part of the asymptotics are special solutions of boundary-value problems in an unbounded domain with different outlets at infinity. It turns out that they have polynomial growth at infinity. By matching these parts, we derive the limit problem $(\varepsilon =0)$ in the corresponding graph and a recurrence procedure to determine all terms of the asymptotic expansion. Energetic and uniform pointwise estimates are proved. These estimates allow us to observe the impact of the aneurysm.


title: Scalability Analysis of the RADAR Decision Support Tool

abstract: This report presents a theoretical complexity analysis and empirical scalability analysis of the Requirements and Architecture Decision Analyser (RADAR).


title: Elastic Resource Management with Adaptive State Space Partitioning of Markov Decision Processes

abstract: Modern large-scale computing deployments consist of complex applications running over machine clusters. An important issue in these is the offering of elasticity, i.e., the dynamic allocation of resources to applications to meet fluctuating workload demands. Threshold based approaches are typically employed, yet they are difficult to configure and optimize. Approaches based on reinforcement learning have been proposed, but they require a large number of states in order to model complex application behavior. Methods that adaptively partition the state space have been proposed, but their partitioning criteria and strategies are sub-optimal. In this work we present MDP_DT, a novel full-model based reinforcement learning algorithm for elastic resource management that employs adaptive state space partitioning. We propose two novel statistical criteria and three strategies and we experimentally prove that they correctly decide both where and when to partition, outperforming existing approaches. We experimentally evaluate MDP_DT in a real large scale cluster over variable not-encountered workloads and we show that it takes more informed decisions compared to static and model-free approaches, while requiring a minimal amount of training data.


title: Effective spin physics in two-dimensional cavity QED arrays

abstract: We investigate a strongly correlated system of light and matter in two-dimensional cavity arrays. We formulate a Jaynes-Cummings Hamiltonian for two-level atoms coupled to cavity modes and driven by an external laser field which reduces to an effective spin Hamiltonian in the dispersive regime. In one dimension we provide exact analytical solution. In two dimensions, we perform mean-field study and large scale quantum Monte Carlo simulations of both the Jaynes-Cummings and the effective spin models. We discuss the phase diagram and the parameter regime which gives rise to frustrated interactions between the spins. We provide quantitative description of the phase transitions and correlation properties featured by the system and we discuss graph-theoretical properties of the ground states in terms of graph colorings using P\'{o}lya's enumeration theorem.


title: Computability and Complexity of Unconventional Computing Devices

abstract: We discuss some claims that certain UCOMP devices can perform hypercomputation (compute Turing-uncomputable functions) or perform super-Turing computation (solve NP-complete problems in polynomial time). We discover that all these claims rely on the provision of one or more unphysical resources.


title: Trigonometric integrators for quasilinear wave equations

abstract: Trigonometric time integrators are introduced as a class of explicit numerical methods for quasilinear wave equations. Second-order convergence for the semi-discretization in time with these integrators is shown for a sufficiently regular exact solution. The time integrators are also combined with a Fourier spectral method into a fully discrete scheme, for which error bounds are provided without requiring any CFL-type coupling of the discretization parameters. The proofs of the error bounds are based on energy techniques and on the semiclassical G\aa rding inequality.


title: Fixing an error in Caponnetto and de Vito (2007)

abstract: The seminal paper of Caponnetto and de Vito (2007) provides minimax-optimal rates for kernel ridge regression in a very general setting. Its proof, however, contains an error in its bound on the effective dimensionality. In this note, we explain the mistake, provide a correct bound, and show that the main theorem remains true.


title: Neutron Stars in Screened Modified Gravity: Chameleon vs Dilaton

abstract: We consider the scalar field profile around relativistic compact objects such as neutron stars for a range of modified gravity models with screening mechanisms of the chameleon and Damour-Polyakov types. We focus primarily on inverse power law chameleons and the environmentally dependent dilaton as examples of both mechanisms. We discuss the modified Tolman-Oppenheimer-Volkoff equation and then implement a relaxation algorithm to solve for the scalar profiles numerically. We find that chameleons and dilatons behave in a similar manner and that there is a large degeneracy between the modified gravity parameters and the neutron star equation of state. This is exemplified by the modifications to the mass-radius relationship for a variety of model parameters.


title: Graded multiplications on iterated bar constructions

abstract: We define a bar construction endofunctor on the category of commutative augmented monoids $A$ of a symmetric monoidal category $\mathcal{V}$ endowed with a left adjoint monoidal functor $F:s\mathbf{Set}\to \mathcal{V}$. To do this, we need to carefully examine the monoidal properties of the well-known (reduced) simplicial bar construction $B_\bullet(1,A,1)$. We define a geometric realization $|-|$ with respect to the image under $F$ of the canonical cosimplicial simplicial set. This guarantees good monoidal properties of $|-|$: it is monoidal, and given a left adjoint monoidal functor $G:\mathcal{V}\to \mathcal{W}$, there is a monoidal transformation $|G-|\Rightarrow G|-|$. We can then consider $BA=|B_\bullet A|$ and the iterations $B^nA$. We establish the existence of a graded multiplication on these objects, provided the category $\mathcal{V}$ is cartesian and $A$ is a ring object. The examples studied include simplicial sets and modules, topological spaces, chain complexes and spectra.


title: Minimal complexes of cotorsion flat modules

abstract: Let R be a commutative noetherian ring. We give criteria for a complex of cotorsion flat R-modules to be minimal, in the sense that every self homotopy equivalence is an isomorphism. To do this, we exploit Enochs' description of the structure of cotorsion flat R-modules. More generally, we show that any complex built from covers in every degree (or envelopes in every degree) is minimal, as well as give a partial converse to this in the context of cotorsion pairs. As an application, we show that every R-module is isomorphic in the derived category over R to a minimal semi-flat complex of cotorsion flat R-modules.


title: Totally acyclic complexes and locally Gorenstein rings

abstract: A commutative noetherian ring with a dualizing complex is Gorenstein if and only if every acyclic complex of injective modules is totally acyclic. We extend this characterization, which is due to Iyengar and Krause, to arbitrary commutative noetherian rings, i.e. we remove the assumption about a dualizing complex. In this context Gorenstein, of course, means locally Gorenstein at every prime.


title: Modification of microstructure and micromagnetic properties in Gd-Fe thin films by rapid thermal processing

abstract: Impact of rapid thermal processing (RTP) on microstructure and magnetic properties of Gd-Fe thin films have been investigated with a special emphasis to magnetic microstructure. 100 nm thick amorphous Gd-Fe film shows elongated stripe domains with characteristic feature size of 122 nm, which signifies the development of perpendicular magnetic anisotropy (PMA) in this system. RTP at 550^oC for different time intervals viz. 5, 10, 15, 20 minutes induces the crystallization of Fe over the amorphous Gd-Fe matrix. Cross-sectional imaging with transmission electron microscope reveals the presence of Fe nanocrystal and the reacted film-substrate interface, which has also been verified by Rutherford backscattering spectrometry. With increasing RTP time, the characteristic size and the threshold height of the topographic feature increases. While the magnetic contrast drops down for the processed films, the increasing threshold height of the nanocrystalline grains results a topography dominated mixed phase in the magnetic force microscopy images. The magnetization measurements clearly indicate a re-orientation of the magnetization direction from perpendicular to the plane of the surface. To unravel the effect of anisotropy modification on domain structure, 3D micromagnetic simulations have been performed, which showed that with 50% reduction of original anisotropy, the magnetization tries to orient in the plane of the surface, whereas an in-plane (along xz plane) tilt angle of 30^o has been estimated to be the threshold for the decrease in perpendicular component of magnetization in the investigated Gd-Fe thin film system.


title: New Hermite-Hadamard type integral inequalities for convex functions and theirs applications

abstract: In this paper, we establish (presumably new type) integral inequalities for convex functions via the Hermite--Hadamard's inequalities. As applications, we apply these new inequalities to construct inequalities involving special means of real numbers, some error estimates for the formula midpoint are given. Finally, new inequalities for some special and $q-$special functions are also pointed out.


title: Incompressible fluid problems on embedded surfaces: Modeling and variational formulations

abstract: Governing equations of motion for a viscous incompressible material surface are derived from the balance laws of continuum mechanics. The surface is treated as a time-dependent smooth orientable manifold of codimension one in an ambient Euclidian space. We use elementary tangential calculus to derive the governing equations in terms of exterior differential operators in Cartesian coordinates. The resulting equations can be seen as the Navier-Stokes equations posed on an evolving manifold. We consider a splitting of the surface Navier-Stokes system into coupled equations for the tangential and normal motions of the material surface. We then restrict ourselves to the case of a geometrically stationary manifold of codimension one embedded in $\Bbb{R}^n$. For this case, we present new well-posedness results for the simplified surface fluid model consisting of the surface Stokes equations. Finally, we propose and analyze several alternative variational formulations for this surface Stokes problem, including constrained and penalized formulations, which are convenient for Galerkin discretization methods.


title: A line-smeared treatment of opacities for the spectra and light curves from macronovae

abstract: Gravitational wave observations need accompanying electromagnetic signals to accurately determine the sky positions of the sources. The ejecta of neutron star mergers are expected to produce such electromagnetic transients, called macronovae. Characteristics of the ejecta include large velocity gradients and the presence of heavy $r$-process elements, which pose significant challenges to the accurate calculation of radiative opacities and radiation transport. For example, these opacities include a dense forest of bound-bound features arising from near-neutral lanthanide and actinide elements. Here we investigate the use of fine-structure, line-smeared opacities that preserve the integral of the opacity over frequency, which is an alternative approach to the expansion-opacity formalism. The use of area-preserving line profiles produces frequency-dependent opacities that are one to two orders of magnitude greater than those obtained with the use of expansion opacities in the Sobolev approximation. Opacities are calculated for four $r$-process elements (cerium, neodymium, samarium and uranium) using fully and semi-relativistic methods, as well as different amounts of configuration interaction in the atomic structure, in order to test the sensitivity of the emission to the underlying atomic physics. We determine the effect of these new opacities on simulated spectra and broad-band light curves by applying a multi-dimensional ray-trace method to the ejecta predicted from 3-dimensional merger calculations. With our substantially larger opacities the simulations yield slightly lower luminosities that peak in the mid-IR. These results suggest that those radioactively powered transients that are related to the very heaviest $r$-process material are more difficult to observe than previously believed.


abstract: Solar flares are complicated physical phenomena that are observable in a broad range of the electromagnetic spectrum, from radiowaves to $\gamma$-rays. For a more comprehensive understanding of flares, it is necessary to perform a combined multi-wavelength analysis using observations from many satellites and ground-based observatories. For efficient data search, integration of different flare lists and representation of observational data, we have developed an Interactive Multi-Instrument Database of Solar Flares (https://solarflare.njit.edu/). The web accessible database is fully functional and allows the user to search for uniquely-identified flare events based on their physical descriptors and availability of observations by a particular set of instruments. Currently, the data from three primary flare lists (GOES, RHESSI and HEK) and a variety of other event catalogs (Hinode, Fermi GBM, Konus-Wind, OVSA flare catalogs, CACTus CME catalog, Filament eruption catalog) and observing logs (IRIS and Nobeyama coverage) are integrated, and an additional set of physical descriptors (temperature and emission measure) is provided along with an observing summary, data links, and multi-wavelength light curves for each flare event since January, 2002. We envision that this new tool will allow researchers to significantly speed up the search of events of interest for statistical and case studies.


title: The 2HWC HAWC Observatory Gamma Ray Catalog

abstract: We present the first catalog of TeV gamma-ray sources realized with the recently completed High Altitude Water Cherenkov Observatory (HAWC). It is the most sensitive wide field-of-view TeV telescope currently in operation, with a 1-year survey sensitivity of ~5-10% of the flux of the Crab Nebula. With an instantaneous field of view >1.5 sr and >90% duty cycle, it continuously surveys and monitors the sky for gamma ray energies between hundreds GeV and tens of TeV. HAWC is located in Mexico at a latitude of 19 degree North and was completed in March 2015. Here, we present the 2HWC catalog, which is the result of the first source search realized with the complete HAWC detector. Realized with 507 days of data and represents the most sensitive TeV survey to date for such a large fraction of the sky. A total of 39 sources were detected, with an expected contamination of 0.5 due to background fluctuation. Out of these sources, 16 are more than one degree away from any previously reported TeV source. The source list, including the position measurement, spectrum measurement, and uncertainties, is reported. Seven of the detected sources may be associated with pulsar wind nebulae, two with supernova remnants, two with blazars, and the remaining 23 have no firm identification yet.


title: Short term unpredictability of high Reynolds number turbulence --- rough dependence on initial data

abstract: Short term unpredictability is discovered numerically for high Reynolds number fluid flows under periodic boundary conditions. Furthermore, the abundance of the short term unpredictability is also discovered. These discoveries support our theory that fully developed turbulence is constantly driven by such short term unpredictability.

categories: physics.flu-dyn math.AP nlin.CD physics.ao-ph