title: Why are Megaprojects, Including Nuclear Power Plants, Delivered Overbudget and Late? Reasons and Remedies

abstract: In the first section, this report analyses Nuclear Power Plants (NPPs) in the context of megaprojects, explaining why they are often delivered over budget and late. In the second section, the report discusses how Small Modular Reactors (SMRs) might address these issues. Megaprojects are extremely risky and often implemented after a sub-optimal phase of project planning leading to underestimations of the costs and overestimation of short-term benefits. When considering adherence to schedule and budget, often megaprojects might be considered a failure, and optimism bias, strategic mis-rapresentation, complexity, poor planning, poor risk allocation, poor scope management are all reasons to explain their over budget and delay. For megaprojects, especially in the nuclear field, a key strategy to achieve good performances appears to be the standardization. This standardization needs to be twofold: (i) technical standardisation, i.e. the construction of very similar design over and over, and (ii) the project delivery chain standardisation, i.e. the same stakeholders involved in the delivery of a project that is replicable multiple times. Under this perspective, given their size and standardisation potential, SMRs, might be a suitable class of NPP for several countries. Yet, if the economy of scale is the only driver considered, SMRs are hardly competitive with large NPPs (or even with gas or coal power plants). However, a fleet of standard SMRs might balance the lack of economy of scale with the economy of multiples, and the delivery of several standardised SMR projects might be the key to achieve good project management performances in the nuclear sector. However, the deployment of SMRs faces a number of challenges from several perspectives, such as the licencing, supply chain and financing ones. These challenges might be enormous, but so are the potential rewards too.


title: A Novel Hybrid Islanding Detection Method for Inverter-based DG

abstract: A novel method for achieving a better performance using the combination of the available passive and active methods has been proposed. The algorithm detects the islanding in proper time by using harmonic detection, the average rate of change of voltage and shifting power generation. Harmonic detection in this method decreases process time and also differentiates between islanding and other power systems events. For harmonic detection, extended Kalman filter has been used. Besides, the reliability of the method increases using the average rate of change of the voltage. The proposed method uses a strategy for decreasing the non-detection zone. In this strategy, minimum and maximum average rates of change of voltage limits are defined to improve the security of the system. Therefore, three main specifications of a proper method, reliability, security and time of process are achievable by the combination of these passive and active methods. By applying different power system events under different power conditions, the proposed method has been verified in Simulink software.


title: On the chemical abundances of Miras in clusters: V1 in the metal-rich globular NGC 5927

abstract: We present the first spectroscopic abundance determination of iron, alpha-elements (Si, Ca and Ti) and sodium for the Mira variable V1 in the metal-rich globular cluster NGC 5927. We use high-resolution (R~ 28,000), high signal-to-noise ratio (~200) spectra collected with WINERED, a near-infrared (NIR) spectrograph covering simultaneously the wavelength range 0.91--1.35 micron. The effective temperature and the surface gravity at the pulsation phase of the spectroscopic observation were estimated using both optical (V) and NIR time-series photometric data. We found that the Mira is metal-rich ([Fe/H]=-0.55 \pm 0.15) and moderately alpha-enhanced ([alpha/Fe]=0.15 \pm 0.01, sigma=0.2). These values agree quite well with the mean cluster abundances based on high-resolution optical spectra of several cluster red giants available in the literature ([Fe/H]=-0.47 \pm 0.06, [alpha/Fe]=+0.24 \pm 0.05). We also found a Na abundance of +0.35 \pm 0.20 that is higher than the mean cluster abundance based on optical spectra (+0.18 \pm 0.13). However, the lack of similar spectra for cluster red giants and that of corrections for departures from local-thermodynamical equilibrium prevents us from establishing whether the difference is intrinsic or connected with multiple populations. These findings indicate a strong similarity between optical and NIR metallicity scales in spite of the difference in the experimental equipment, data analysis and in the adopted spectroscopic diagnostics.


title: A note concerning the modular valued von Neumann interaction operator

abstract: The modular valued operator $\widehat{V}_m$ of the von Neumann interaction operator for a projector is defined. The properties of $\widehat{V}_m$ are discussed and contrasted with those of the standard modular value of a projector. The associated notion of a faux qubit is introduced and its possible utility in quantum computation is noted. An experimental implementation of $\widehat{V}_m$ is also highlighted.


title: The search for radio emission from exoplanets using LOFAR beam-formed observations: Jupiter as an exoplanet

abstract: $\textit{Context.}$ The magnetized Solar System planets are strong radio emitters and theoretical studies suggest that the radio emission from nearby exoplanets in close-in orbits could reach intensity levels $10^{3}-10^{7}$ times higher than Jupiter's decametric emission. Detection of exoplanets in the radio domain would open up a brand new field of research, however, currently there are no confirmed detections at radio frequencies. $\textit{Aims.}$ We investigate the radio emission from Jupiter, scaled such that it mimics emission coming from an exoplanet, with low-frequency beam-formed observations using LOFAR. The goals are to define a set of observables that can be used as a guideline in the search for exoplanetary radio emission and to measure effectively the sensitivity limit for LOFAR beam-formed observations. $\textit{Methods.}$ We observe "Jupiter as an exoplanet" by dividing a LOFAR observation of Jupiter by a down-scaling factor and adding this observation to beam-formed data of the "sky background". Then we run this artificial dataset through our total intensity (Stokes-I) and circular polarization (Stokes-V) processing and post-processing pipelines and determine up to which down-scaling factor Jupiter is still detected in the dataset. $\textit{Results.}$ We find that exoplanetary radio bursts can be detected at 5 pc if the circularly polarized flux is $10^5$ times stronger than the typical level of Jupiter's radio bursts during active emission events ($\sim4\times10^{5}$ Jy). Equivalently, circularly polarized radio bursts can be detected up to a distance of 20 pc (encompassing the known exoplanets 55 Cnc, Tau Bo\"{o}tis, and Upsilon Andromedae) assuming the level of emission is $10^{5}$ times stronger than the peak flux of Jupiter's decametric burst emission ($\sim6\times10^{6}$ Jy).


title: Transport through a network of topological states in twisted bilayer graphene

abstract: We explore a network of electronic quantum valley Hall (QVH) states in the moir\'e crystal of minimally twisted bilayer graphene. In our transport measurements we observe Fabry-P\'erot and Aharanov-Bohm oscillations which are robust in magnetic fields ranging from 0 to 8T, in strong contrast to more conventional 2D systems where trajectories in the bulk are bent by the Lorentz force. This persistence in magnetic field and the linear spacing in density indicate that charge carriers in the bulk flow in topologically protected, one dimensional channels. With this work we demonstrate coherent electronic transport in a lattice of topologically protected states.


title: High-temperature structural phase transition and infrared dielectric features of La2CoMnO6

abstract: Temperature-dependent FAR-infrared reflectivity spectra of partially ordered magnetodielectric La2CoMnO6 is presented, from room temperature up to 675 K. A clear first-ordered structural phase transition (SPT) from a monoclinic structure with P2_1/n symmetry to a rhombohedral phase with R-3 symmetry was evidenced from the behaviour of polar phonon modes at TC~590 K. The temperature dependences of the transversal and longitudinal phonon branches, dielectric strengths, and damping of the strongest dielectric modes confirm the significant contribution of the phonon modes on the SPT, and revealed an important lattice anharmonicity, particularly for the low frequency modes. In addition, these investigations showed that structural ordering does not inhibit the SPT, and provided valuable information towards the polar phonons, their implications on intrinsic dielectric constant in double perovskites and in related compounds.


title: The Updated BaSTI Stellar Evolution Models and Isochrones: I. Solar Scaled Calculations

abstract: We present an updated release of the BaSTI (a Bag of Stellar Tracks and Isochrones) stellar model and isochrone library for a solar scaled heavy element distribution. The main input physics changed from the previous BaSTI release include the solar metal mixture, electron conduction opacities, a few nuclear reaction rates, bolometric corrections, and the treatment of the overshooting efficiency for shrinking convective cores. The new model calculations cover a mass range between 0.1 and 15 Msun, 22 initial chemical compositions between [Fe/H]=-3.20 and +0.45, with helium to metal enrichment ratio dY /dZ=1.31. The isochrones cover an age range between 20 Myr and 14.5 Gyr, take consistently into account the pre-main sequence phase, and have been translated to a large number of popular photometric systems. Asteroseismic properties of the theoretical models have also been calculated. We compare our isochrones with results from independent databases and with several sets of observations, to test the accuracy of the calculations. All stellar evolution tracks, asteroseismic properties and isochrones are made available through a dedicated Web site.


title: Degrees of extensionality in the theory of B\"ohm trees and Sall\'e's conjecture

abstract: The main observational equivalences of the untyped lambda-calculus have been characterized in terms of extensional equalities between B\"ohm trees. It is well known that the lambda-theory H*, arising by taking as observables the head normal forms, equates two lambda-terms whenever their B\"ohm trees are equal up to countably many possibly infinite eta-expansions. Similarly, two lambda-terms are equal in Morris's original observational theory H+, generated by considering as observable the beta-normal forms, whenever their B\"ohm trees are equal up to countably many finite eta-expansions. The lambda-calculus also possesses a strong notion of extensionality called "the omega-rule", which has been the subject of many investigations. It is a longstanding open problem whether the equivalence B-omega obtained by closing the theory of B\"ohm trees under the omega-rule is strictly included in H+, as conjectured by Sall\'e in the seventies. In this paper we demonstrate that the two aforementioned theories actually coincide, thus disproving Sall\'e's conjecture. The proof technique we develop for proving the latter inclusion is general enough to provide as a byproduct a new characterization, based on bounded eta-expansions, of the least extensional equality between B\"ohm trees. Together, these results provide a taxonomy of the different degrees of extensionality in the theory of B\"ohm trees.


title: Robustness in Consensus Networks

abstract: We consider the problem of robustness in large consensus networks that occur in many areas such as distributed optimization. Robustness, in this context, is the scaling of performance measures, e.g. H2-norm, as a function of network dimension. We provide a formal framework to quantify the relation between such performance scaling and the convergence speed of the network. Specifically, we provide upper and lower bounds for the convergence speed in terms of robustness and discuss how these bounds scale with the network topology. The main contribution of this work is that we obtain tight bounds, that hold regardless of network topology. The work here also encompasses some results in convergence time analysis in previous literature.


title: Broyden's method for nonlinear eigenproblems

abstract: Broyden's method is a general method commonly used for nonlinear systems of equations, when very little information is available about the problem. We develop an approach based on Broyden's method for nonlinear eigenvalue problems. Our approach is designed for problems where the evaluation of a matrix vector product is computationally expensive, essentially as expensive as solving the corresponding linear system of equations. We show how the structure of the Jacobian matrix can be incorporated into the algorithm to improve convergence. The algorithm exhibits local superlinear convergence for simple eigenvalues, and we characterize the convergence. We show how deflation can be integrated and combined such that the method can be used to compute several eigenvalues. A specific problem in machine tool milling, coupled with a PDE is used to illustrate the approach. The simulations are done in the julia programming language, and are provided as publicly available module for reproducability.


title: Parametric amplification and squeezing with an ac- and dc-voltage biased superconducting junction

abstract: We theoretically investigate a near-quantum-limited parametric amplifier based on the nonlinear dynamics of quasiparticles flowing through a superconducting-insulator-superconducting junction. Photon-assisted tunneling, resulting from the combination of dc- and ac-voltage bias, gives rise to a strong parametric interaction for the electromagnetic modes reflected by the junction coupled to a transmission line. We show phase-sensitive and phase-preserving amplification, together with single- and two-mode squeezing. For an aluminum junction pumped at twice the center frequency, $\omega_0/2\pi=6$~GHz, we predict narrow-band phase-sensitive amplification of microwaves signals to more than 20 dB, and broadband phase-preserving amplification of 20 dB over a 1.2 GHz 3-dB bandwidth. We also predict single- and two-mode squeezing reaching more than -12 dB over 5.3 GHz 3-dB bandwidth. Moreover, with a simple impedance matching circuit, we demonstrate 3 dB bandwidth reaching 4.3 GHz for 20 dB of gain. A key feature of the device is that its performance can be controlled in-situ with the applied dc- and ac-voltage biases.


title: Using Semi-Supervised Learning for Predicting Metamorphic Relations

abstract: Software testing is difficult to automate, especially in programs which have no oracle, or method of determining which output is correct. Metamorphic testing is a solution this problem. Metamorphic testing uses metamorphic relations to define test cases and expected outputs. A large amount of time is needed for a domain expert to determine which metamorphic relations can be used to test a given program. Metamorphic relation prediction removes this need for such an expert. We propose a method using semi-supervised machine learning to detect which metamorphic relations are applicable to a given code base. We compare this semi-supervised model with a supervised model, and show that the addition of unlabeled data improves the classification accuracy of the MR prediction model.


title: The chemical connection between damped Lyman-\alpha\ systems and Local Group dwarf galaxies

abstract: Abundances of the volatile elements S and Zn have now been measured in around 80 individual stars in the Sculptor dwarf spheroidal galaxy, covering the metallicity range $-2.4\leq\text{[Fe/H]}\leq-0.9$. These two elements are of particular interest as they are not depleted onto dust in gas, and their ratio, [S/Zn], has thus commonly been used as a proxy for [$\alpha$/Fe] in Damped Lyman-$\alpha$ systems. The S abundances in Sculptor are similar to other $\alpha$-elements in this galaxy, consistent with S being mainly created in core-collapse supernovae, but also having some contribution from supernovae Type Ia. However, our results show that Zn and Fe do not trace all the same nucleosynthetic production channels. In particular, (contrary to Fe) Zn is not significantly produced by supernovae Type Ia. Thus, [S/Zn] cannot be reliably used as a proxy for [$\alpha$/Fe]. We propose [O/S] as a function of [S/H] as a possible alternative. At higher metallicities, the values of [S/Zn] measured in Damped Lyman-$\alpha$ systems are inconsistent with those in local dwarf galaxies, and are more compatible with the Milky Way disk. Low-metallicity Damped Lyman-$\alpha$ systems are, however, consistent with the most metal-poor stars in Local Group dwarf spheroidal galaxies. Assuming that the dust depletions of S and Zn are negligible, our comparison indicates that the star formation histories of Damped Lyman-$\alpha$ systems are on average different from both the Milky Way and the Sculptor dwarf spheroidal galaxy.


title: Short-range antiferromagnetic interaction and spin-phonon coupling in La2CoMnO6 double perovskite

abstract: Weak antiferromagnetic (AF) interaction in the ferromagnetic (FM) partially ordered La2CoMnO6 (LCMO) was detected by Raman spectroscopy by monitoring spin-phonon coupling. Because of the sensibility to probe short-range disorder and lattice modifications, the Raman spectroscopy showed to be an useful tool to indicate less remarkable magnetic transitions in LCMO compound. Apart from the expected spin-phonon coupling due to the long-range FM superexchange (Tc~230 K), phonons parameters pointed out an additional spin-phonon coupling related to the short-range AF interaction at Tc~135 K, which was not detected from the magnetic bulk response. These results reinforce the Raman spectroscopy as a powerful technique to detect antisite disorder into A2B'B"O6 magnetic double perovskites, whose magnetic properties are driven by superexchange interactions.


title: Asymmetric skyrmion Hall effect in systems with a hybrid Dzyaloshinskii-Moriya interaction

abstract: We examine the current-induced dynamics of a skyrmion that is subject to both structural and bulk inversion asymmetry. There arises a hybrid type of Dzyaloshinskii-Moriya interaction (DMI) which is in the form of a mixture of interfacial and bulk DMIs. Examples include crystals with symmetry classes C$_n$ as well as magnetic multilayers composed of a ferromagnet with a noncentrosymmetric crystal and a nonmagnet with strong spin-orbit coupling. As a striking result, we find that, in systems with a hybrid DMI, the spin-orbit-torque-induced skyrmion Hall angle is asymmetric for the two different skyrmion polarities ($\pm 1$ given by out-of-plane core magnetization), even allowing one of them to be tuned to zero. We propose several experimental ways to achieve the necessary straight skyrmion motion (with zero Hall angle) for racetrack memories, even without antiferromagnetic interactions or any interaction with another magnet. Our results can be understood within a simple picture by using a global spin rotation which maps the hybrid DMI model to an effective model containing purely interfacial DMI. The formalism directly reveals the effective spin torque and effective current that result in qualitatively different dynamics. Our work provides a way to utilize symmetry breaking to eliminate detrimental phenomena as hybrid DMI eliminates the skyrmion Hall angle.


title: Gamma-ray burst models in light of the GRB 170817A - GW170817 connection

abstract: For the first time, a short gamma-ray burst (GRB) was unambiguously associated with a gravitational wave (GW) observation from a binary neutron star (NS) merger. This allows us to link the details of the central engine properties to GRB emission models. We find that photospheric models (both dissipative and non-dissipative variants) have difficulties accounting for the observations. Internal shocks give the most natural account of the observed peak energy, viewing angle and total energy. We also show that a simple external shock model can reproduce the observed GRB pulse with parameters consistent with those derived from the afterglow modeling. We find a simple cocoon shock breakout model is in mild tension with the observed spectral evolution, however it cannot be excluded based on gamma-ray data alone. Future joint observations of brighter GRBs will pose even tighter constraints on prompt emission models.


title: Bayesian Incremental Learning for Deep Neural Networks

abstract: In industrial machine learning pipelines, data often arrive in parts. Particularly in the case of deep neural networks, it may be too expensive to train the model from scratch each time, so one would rather use a previously learned model and the new data to improve performance. However, deep neural networks are prone to getting stuck in a suboptimal solution when trained on only new data as compared to the full dataset. Our work focuses on a continuous learning setup where the task is always the same and new parts of data arrive sequentially. We apply a Bayesian approach to update the posterior approximation with each new piece of data and find this method to outperform the traditional approach in our experiments.


title: A folded model for compositional data analysis

abstract: A folded type model is developed for analyzing compositional data. The proposed model involves an extension of the $\alpha$-transformation for compositional data and provides a new and flexible class of distributions for modeling data defined on the simplex sample space. Despite its rather seemingly complex structure, employment of the EM algorithm guarantees efficient parameter estimation. The model is validated through simulation studies and examples which illustrate that the proposed model performs better in terms of capturing the data structure, when compared to the popular logistic normal distribution, and can be advantageous over a similar model without folding.


title: Strongly enhanced tunneling at total charge neutrality in double bilayer graphene-WSe$_2$ heterostructures

abstract: We report the experimental observation of strongly enhanced tunneling between graphene bilayers through a WSe$_2$ barrier when the graphene bilayers are populated with carriers of opposite polarity and equal density. The enhanced tunneling increases sharply in strength with decreasing temperature, and the tunneling current exhibits a vertical onset as a function of interlayer voltage at a temperature of 1.5 K. The strongly enhanced tunneling at overall neutrality departs markedly from single-particle model calculations that otherwise match the measured tunneling current-voltage characteristics well, and suggests the emergence of a many-body state with condensed interbilayer excitons when electrons and holes of equal densities populate the two layers.


title: Achievability Bounds for T-Fold Irregular Repetition Slotted ALOHA Scheme in the Gaussian MAC

abstract: We address the problem of uncoordinated massive random-access in the Gaussian multiple access channel (MAC). The performance of low-complexity T-fold irregular repetition slotted ALOHA (IRSA) scheme is investigated and achievability bounds are derived. The main difference of this scheme in comparison to IRSA is as follows: any collisions of order up to T can be resolved with some probability of error introduced by noise. In order to optimize the parameters of the scheme we combine the density evolution method (DE) proposed by G. Liva and a finite length random coding bound for the Gaussian MAC proposed by Y. Polyanskiy. As energy efficiency is of critical importance for massive machine-type communication (mMTC), then our main goal is to minimize the energy-per-bit required to achieve the target packet loss ratio (PLR). We consider two scenarios: (a) the number of active users is fixed; (b) the number of active users is a Poisson random variable.


title: Modeling of Supersonic Radiative Marshak waves using Simple Models and Advanced Simulations

abstract: We study the problem of radiative heat (Marshak) waves using advanced approximate approaches. Supersonic radiative Marshak waves that are propagating into a material are radiation dominated (i.e. hydrodynamic motion is negligible), and can be described by the Boltzmann equation. However, the exact thermal radiative transfer problem is a nontrivial one, and there still exists a need for approximations that are simple to solve. The discontinuous asymptotic $P_1$ approximation, which is a combination of the asymptotic $P_1$ and the discontinuous asymptotic diffusion approximations, was tested in previous work via theoretical benchmarks. Here we analyze a fundamental and typical experiment of a supersonic Marshak wave propagation in a low-density $\mathrm{SiO_2}$ foam cylinder, embedded in gold walls. First, we offer a simple analytic model, that grasps the main effects dominating the physical system. We find the physics governing the system to be dominated by a simple, one-dimensional effect, based on the careful observation of the different radiation temperatures that are involved in the problem. The model is completed with the main two-dimensional effect which is caused by the loss of energy to the gold walls. Second, we examine the validity of the discontinuous asymptotic $P_1$ approximation, comparing to exact simulations with good accuracy. Specifically, the heat front position as a function of the time is reproduced perfectly in compare to exact Boltzmann solutions.


title: Multigap superconductivity in RbCa2Fe4As4F2 investigated using muSR measurements

abstract: The superconducting properties of the recently discovered double Fe$_2$As$_2$ layered high-$T_c$ superconductor RbCa$_2$Fe$_4$As$_4$F$_2$ with $T_c\approx$ 30~K have been investigated using magnetization, heat capacity, transverse-field (TF) and zero-field (ZF) muon-spin rotation/relaxation ($\mu$SR) measurements. Our low field magnetization measurements and heat capacity (C$_p$) reveal an onset of bulk superconductivity with $T_{\bf c}\sim$ 30.0(4) K. Furthermore, the heat capacity exhibits a jump at $T_{\bf c}$ of $\Delta$C$_p$/$T_{\bf c}$=94.6 (mJ/mole-K$^2$) and no clear effect of applied magnetic fields was observed on C$_p$(T) up to 9 T between 2 K and 5 K. Our analysis of the TF-$\mu$SR results shows that the temperature dependence of the magnetic penetration depth is better described by a two-gap model, either isotropic $s$+$s$-wave or $s$+$d$-wave than a single gap isotropic $s$-wave or $d$-wave model for the superconducting gap. The presence of two superconducting gaps in RbCa$_2$Fe$_4$As$_4$F$_2$ suggests a multiband nature of the superconductivity, which is consistent with the multigap superconductivity observed in other Fe-based superconductors, including ACa$_2$Fe$_4$As$_4$F$_2$ (A=K and Cs). Furthermore, from our TF-$\mu$SR study we have estimated an in-plane penetration depth $\lambda_{\mathrm{ab}}$$(0)$ =231.5(3) nm, superconducting carrier density $n_s = 7.45 \times 10^{26}~ $m$^{-3}$, and carrier's effective-mass $m^*$ = 2.45\textit{m}$_{e}$. Our ZF $\mu$SR measurements do not reveal a clear sign of time reversal symmetry breaking at $T_{\bf c}$, but the temperature dependent relaxation between 150 K and 1.2 K might indicate the presence of spin-fluctuations. The results of our present study have been compared with those reported for other Fe pnictide superconductors.


title: Performance Comparison of two novel Relay-Assisted Hybrid FSO / RF Communication Systems

abstract: In this manuscript, two novel multi-hop relay-assisted hybrid Free Space Optical / Radio Frequency (FSO / RF) communication systems are presented and compared. In these structures, RF and FSO links, at each hop, are parallel and send data simultaneously. This is the first time that in a multihop hybrid FSO / RF structure, Detect and Forward protocol is used. In the first structure, at each hop, received signals with higher Signal to Noise Ratio (SNR) is selected. But in the second structure, at each hop, received FSO and RF signals are separately detected and forwarded and selection is done only at the last hop. Considering FSO link in Negative Exponential atmospheric turbulence and RF link in Rayleigh fading, for the first time, closed-form expressions are derived for Outage Probability (P_out) and Bit Error Rate (BER) of the proposed structures. MATLAB simulations are provided to verify derived expressions. The main motivation of this work is to answer this question that how much is the difference of selection at each hop and selection at the last hop. Results indicate that the structure with selection at each hop has better performance than the structure with selection at the last hop. At different target Outage Probability, selection at the last hop consumes about two times (~3dB) more power than selection at each hop. Both structures are particularly suitable for long-range communications. However, selection at each hop, in the cost of more complexity, is recommended for applications which have problem with supplying the required power for communication.


title: Minimal group determinants and the Lind-Lehmer problem for dihedral groups

abstract: We find the minimal non-trivial integer variable group determinant for any dihedral group of order less than $3.79\times 10^{47}$. We think of this as the Lind-Lehmer problem for the dihedral group. We give a complete description of the determinants for some dihedral groups including $D_{2p}$ and $D_{4p}$.


title: NGC 4051: Black hole mass and photon index-mass accretion rate correlation

abstract: We present a discovery of the correlation between the X-ray spectral (photon) index and mass accretion rate observed in AGN NGC 4051. We analyzed spectral transition episodes observed in NGC 4051 using XMM/Newton, Suzaku and RXTE. We applied a scaling technique for a black hole (BH) mass evaluation which uses a correlation between the photon index and normalization of the seed (disk) component, which is proportional to a mass accretion rate. We developed an analytical model that shows the spectral (photon) index of the BH emergent spectrum undergoes an evolution from lower to higher values depending on a mass accretion rate in the accretion disk. We considered Cygnus X-1 and GRO~J1550-564 as reference sources for which distances, inclination angles and the BH masses are evaluated by dynamical measurements. Application of the scaling technique for the photon index-mass accretion rate correlation provides an estimate of the black hole mass in NGC 4051 to be more than 6x10^5 solar masses.


title: Fertilizer usage and cadmium in soils, crops and food

abstract: Phosphate fertilizers were first implicated by Schroeder and Balassa in 1963 for increasing the Cd concentration in cultivated soils and crops. This suggestion has become a part of the accepted paradigm on soil toxicity. Consequently, stringent fertilizer control programs to monitor Cd have been launched. Attempts to link Cd toxicity and fertilizers to chronic diseases are common. A re-assessment of this "accepted" paradigm is timely, given the larger body of data available today. The data show that both the input and output of Cd per hectare from fertilizers are negligibly small compared to the total amount of Cd/hectare usually present in the soil itself. Calculations based on current agricultural practices are used to show that it will take about 18 centuries to double the ambient soil-cadmium level, and about 8 centuries to double the soil-fluoride level, even after neglecting leaching and other removal effects. Hence the concern of long-term agriculture should be the depletion of available phosphate fertilizers, rather than the contamination of the soil by trace metals or fluoride. This conclusion is confirmed by showing that the claimed correlations between fertilizer input and cadmium accumulation in crops are not robust. Alternative scenarios that explain the data are examined. Thus soil acidulation on fertilizer loading, and the effect of magnesium, zinc, and fluoride ions contained in fertilizers are considered using recent Cd$^{2+}$, Mg$^{2+}$ and F$^-$ ion-association theories. The protective role of ions like Zn, Se, Fe, etc., is emphasized, and the question of cadmium toxicity in the presence of other ions is considered. These help to clarify and rectify difficulties found in the standard point of view. This analysis does not modify the accepted views on Cd contamination by airborne delivery, smoking, and industrial activity, or P-contamination causing algal blooms.


title: Exclusive and diffractive $\mu^+ \mu^-$ production in $pp$ collisions at the LHC

abstract: In this paper we estimate the production of dimuons ($\mu^+ \mu^-$) in exclusive photon -- photon ($\gamma \gamma$) and diffractive Pomeron - Pomeron $(I\!P I\!P)$, Pomeron - Reggeon $(I\!P I\!R)$ and Reggeon - Reggeon $(I\!R I\!R)$ interactions in $pp$ collisions at the LHC energy. The invariant mass, rapidity and tranverse momentum distributions are calculated using the Forward Physics Monte Carlo (FPMC), which allows to obtain realistic predictions for the dimuon production with two leading intact hadrons. In particular, predictions taking into account the CMS and LHCb acceptances are presented. Moreover, the contribution of the single diffraction for the dimuon production also is estimated. Our results demonstrate that the experimental separation of these different mechanisms is feasible. In particular, the events characterized by pairs with large squared transverse momentum are dominated by diffractive interactions, which allows to investigate the underlying assumptions present in the description of these processes.


title: Influence of non-linearity of medium on the laser induced filamentation instability in magnetized plasma

abstract: The effects of the non-linearity of the medium on the growth rate of filamentation instability in a magnetized plasma interacting with an intense laser pulse, is investigated. The non-linearity of the medium, modeled by Kerr non-linearity, is an important factor, which determines the rate of instability growth. Sensitivity of the rate of filamentation growth, to the Kerr non-linear coefficient could be adjusted by the external magnetic field and laser intensity.


title: An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations. Part I: Theory and Numerical Verification

abstract: The first paper of this series presents a discretely entropy stable discontinuous Galerkin (DG) method for the resistive magnetohydrodynamics (MHD) equations on three-dimensional curvilinear unstructured hexahedral meshes. Compared to other fluid dynamics systems such as the shallow water equations or the compressible Navier-Stokes equations, the resistive MHD equations need special considerations because of the divergence-free constraint on the magnetic field. For instance, it is well known that for the symmetrization of the ideal MHD system as well as the continuous entropy analysis a non-conservative term proportional to the divergence of the magnetic field, typically referred to as the Powell term, must be included. As a consequence, the mimicry of the continuous entropy analysis in the discrete sense demands a suitable DG approximation of the non-conservative terms in addition to the ideal MHD terms. This paper focuses on the resistive MHD equations: Our first contribution is a proof that the resistive terms are symmetric and positive-definite when formulated in entropy space as gradients of the entropy variables. This enables us to show that the entropy inequality holds for the resistive MHD equations. This continuous analysis is the key for our DG discretization and guides the path for the construction of an approximation that discretely mimics the entropy inequality, typically termed entropy stability. Our second contribution is a detailed derivation and analysis of the discretization on three-dimensional curvilinear meshes. The discrete analysis relies on the summation-by-parts property, which is satisfied by the DG spectral element method (DGSEM) with Legendre-Gauss-Lobatto (LGL) nodes. Although the divergence-free constraint is included in the non-conservative terms, the resulting method has no particular treatment of the magnetic field divergence errors...


title: A joint decomposition method for global optimization of multiscenario nonconvex mixed-integer nonlinear programs

abstract: This paper proposes a joint decomposition method that combines La- grangian decomposition and generalized Benders decomposition, to efficiently solve multiscenario nonconvex mixed-integer nonlinear programming (MINLP) problems to global optimality, without the need for explicit branch and bound search. In this approach, we view the variables coupling the scenario dependent variables and those causing nonconvexity as complicating variables. We systemat- ically solve the Lagrangian decomposition subproblems and the generalized Ben- ders decomposition subproblems in a unified framework. The method requires the solution of a difficult relaxed master problem, but the problem is only solved when necessary. Enhancements to the method are made to reduce the number of the relaxed master problems to be solved and ease the solution of each relaxed master problem. We consider two scenario-based, two-stage stochastic nonconvex MINLP problems that arise from integrated design and operation of process net- works in the case study, and we show that the proposed method can solve the two problems significantly faster than state-of-the-art global optimization solvers.


title: Mass functions from the excursion set model

abstract: Aims. We aim to study the stochastic evolution of the smoothed overdensity $\delta$ at scale $S$ of the form $\delta(S) = \int_{0}^S K(S,u)\mathrm{d}W(u)$, where $K$ is a kernel and $\mathrm{d}W$ is the usual Wiener process. Methods. For a Gaussian density field, smoothed by the top-hat filter, in real space, we used a simple kernel that gives the correct correlation between scales. A Monte Carlo procedure was used to construct random walks and to calculate first crossing distributions and consequently mass functions for a constant barrier. Results. We show that the evolution considered here improves the agreement with the results of N-body simulations relative to analytical approximations which have been proposed from the same problem by other authors. In fact, we show that an evolution which is fully consistent with the ideas of the excursion set model, describes accurately the mass function of dark matter haloes for values of $\nu \leq 1$ and underestimates the number of larger haloes. Finally, we show that a constant threshold of collapse, lower than it is usually used, it is able to produce a mass function which approximates the results of N-body simulations for a variety of redshifts and for a wide range of masses. Conclusions. A mass function in good agreement with N-body simulations can be obtained analytically using a lower than usual constant collapse threshold.


title: Coconut: Threshold Issuance Selective Disclosure Credentials with Applications to Distributed Ledgers

abstract: Coconut is a novel selective disclosure credential scheme supporting distributed threshold issuance, public and private attributes, re-randomization, and multiple unlinkable selective attribute revelations. Coconut integrates with blockchains to ensure confidentiality, authenticity and availability even when a subset of credential issuing authorities are malicious or offline. We implement and evaluate a generic Coconut smart contract library for Chainspace and Ethereum; and present three applications related to anonymous payments, electronic petitions, and distribution of proxies for censorship resistance. Coconut uses short and computationally efficient credentials, and our evaluation shows that most Coconut cryptographic primitives take just a few milliseconds on average, with verification taking the longest time (10 milliseconds).


title: On a class of solutions to the generalized KdV type equation

abstract: We consider the IVP associated to the generalized KdV equation with low degree of non-linearity \begin{equation*} \partial_t u + \partial_x^3 u \pm |u|^{\alpha}\partial_x u = 0,\; x,t \in \mathbb{R},\;\alpha \in (0,1). \end{equation*} By using an argument similar to that introduced by Cazenave and Naumkin [2] we establish the local well-posedness for a class of data in an appropriate weighted Sobolev space. Also, we show that the solutions obtained satisfy the propagation of regularity principle proven in [3] in solutions of the $k$-generalized KdV equation.


title: Cooperative Robot Localization Using Event-triggered Estimation

abstract: This paper describes a novel communication-spare cooperative localization algorithm for a team of mobile unmanned robotic vehicles. Exploiting an event-based estimation paradigm, robots only send measurements to neighbors when the expected innovation for state estimation is high. Since agents know the event-triggering condition for measurements to be sent, the lack of a measurement is thus also informative and fused into state estimates. The robots use a Covariance Intersection (CI) mechanism to occasionally synchronize their local estimates of the full network state. In addition, heuristic balancing dynamics on the robots' CI-triggering thresholds ensure that, in large diameter networks, the local error covariances remains below desired bounds across the network. Simulations on both linear and nonlinear dynamics/measurement models show that the event-triggering approach achieves nearly optimal state estimation performance in a wide range of operating conditions, even when using only a fraction of the communication cost required by conventional full data sharing. The robustness of the proposed approach to lossy communications, as well as the relationship between network topology and CI-based synchronization requirements, are also examined.


title: Electron-Phonon Systems on a Universal Quantum Computer

abstract: We present an algorithm that extends existing quantum algorithms for simulating fermion systems in quantum chemistry and condensed matter physics to include bosons in general and phonons in particular. We introduce a qubit representation for the low-energy subspace of phonons which allows an efficient simulation of the evolution operator of the electron-phonon systems. As a consequence of the Nyquist-Shannon sampling theorem, the phonons are represented with exponential accuracy on a discretized Hilbert space with a size that increases linearly with the cutoff of the maximum phonon number. The additional number of qubits required by the presence of phonons scales linearly with the size of the system. The additional circuit depth is constant for systems with finite-range electron-phonon and phonon-phonon interactions and linear for long-range electron-phonon interactions. Our algorithm for a Holstein polaron problem was implemented on an Atos Quantum Learning Machine (QLM) quantum simulator employing the Quantum Phase Estimation method. The energy and the phonon number distribution of the polaron state agree with exact diagonalization results for weak, intermediate and strong electron-phonon coupling regimes.


title: A novel hybrid FSO / RF communication system with receive diversity

abstract: In mobile communication system, due to the limitations of mobile device such as low power supply as well as small size, most of the processing should be done at the Base Station. Using multi-receive structure at the Base Station really helps better recovery of the original signal by combining different received signals. In this paper, for the first time, receive diversity is used in single-hop hybrid Free Space Optical / Radio Frequency (FSO / RF) communication system. Also it is the first time that a single-hop hybrid FSO / RF system is investigated at saturate atmospheric turbulence regime. For the first time, closed-form expression is derived for Outage Probability of the proposed system and verified through MATLAB simulation. Results indicate a significant improvement in the performance of the proposed structure compared with common FSO and RF systems with receive diversity. Therefore it can be concluded that although the proposed structure requires a complex receiver, but addition of this complexity could significantly reduce processing or power consumption required for performance maintenance of the system.


title: New regime of plasma wake field acceleration in the extreme blowout regime

abstract: Three dimensional particle in cell simulations are used for studying proton driven plasma wake-field acceleration that uses a high-energy proton bunch to drive a plasma wake-field for electron beam acceleration. A new parameter regime was found which generates essentially constant electric field that is three orders magnitudes larger than that of AWAKE design, i.e. of the order of $2 \times 10^{3}$ GV/m. This is achieved in the the extreme blowout regime, when number density of the driving proton bunch exceeds plasma electron number density 100 times.


title: Reimagine Procrastination: Music Preference and Health Habits as Factors on Self-Perceived Procrastination of Young People

abstract: As the buzzword phenomenon, procrastination holds a continued need for a comprehensive examination of its nature and the associated factors. The presented study explores the potential relationship between music taste, life style and the youngsters' procrastination through quantitative modelling. To handle the big set of survey statistics and the uncertainty caused by the data missingness, the combined methods of factor analysis, multiple imputation (MI) and Ordered logit regression are employed. The result reveals that the music preference for Hip-hop, AlternativeR and Opera have a significant effect on procrastination. Concerning the living habits, the eating habit and local authority (city/rural) also yield strong connection to the self-perceived procrastination. Implications for this procrastination research is discussed.


title: Devon: Deformable Volume Network for Learning Optical Flow

abstract: State-of-the-art neural network models estimate large displacement optical flow in multi-resolution and use warping to propagate the estimation between two resolutions. Despite their impressive results, it is known that there are two problems with the approach. First, the multi-resolution estimation of optical flow fails in situations where small objects move fast. Second, warping creates artifacts when occlusion or dis-occlusion happens. In this paper, we propose a new neural network module, Deformable Cost Volume, which alleviates the two problems. Based on this module, we designed the Deformable Volume Network (Devon) which can estimate multi-scale optical flow in a single high resolution. Experiments show Devon is more suitable in handling small objects moving fast and achieves comparable results to the state-of-the-art methods in public benchmarks.


title: Crystal graphs for shifted tableaux

abstract: We define crystal operators on semistandard shifted tableaux, giving a new proof that Schur $P$-functions are Schur positive. We define a queer crystal operator to construct a connected queer crystal on semistandard shifted tableaux of a given shape, providing a new proof that products of Schur $P$-functions are Schur $P$-positive. We also give a rectification map from shifted tableaux to Young tableaux that commutes with the crystal operators and provides a dual algorithm to shifted insertion.


title: Non-relativistic Arbitrary l-states of Quarkonium through Asymptotic Iteration Method

abstract: The energy eigenvalues with any l-states and mass of heavy quark- antiquark system (quarkonium) are obtained by using Asymptotic Iteration Method in the view of non-relativistic quantum chromodynamics, in which the quarks are considered as spinless for easiness, and are bounded by Cornell potential. A semi-analytical formula for energy eigenvalues and mass is achieved via the method in scope of the perturbation theory. The accuracy of this formula is checked by comparing the eigenvalues with the ones numerically obtained in this study, and with exact ones in literature. Furthermore, semi-analytical formula is applied to some meson systems for comparing the masses with the experimental data.


title: Quality Assurance of Bioinformatics Software: A Case Study of Testing a Biomedical Text Processing Tool Using Metamorphic Testing

abstract: Bioinformatics software plays a very important role in making critical decisions within many areas including medicine and health care. However, most of the research is directed towards developing tools, and little time and effort is spent on testing the software to assure its quality. In testing, a test oracle is used to determine whether a test is passed or failed during testing, and unfortunately, for much of bioinformatics software, the exact expected outcomes are not well defined. Thus, the main challenge associated with conducting systematic testing on bioinformatics software is the oracle problem. Metamorphic testing (MT) is a technique used to test programs that face the oracle problem. MT uses metamorphic relations (MRs) to determine whether a test has passed or failed and specifies how the output should change according to a specific change made to the input. In this work, we use MT to test LingPipe, a tool for processing text using computational linguistics, often used in bioinformatics for bio-entity recognition from biomedical literature. First, we identify a set of MRs for testing any bio-entity recognition program. Then we develop a set of test cases that can be used to test LingPipe's bio-entity recognition functionality using these MRs. To evaluate the effectiveness of this testing process, we automatically generate a set of faulty versions of LingPipe. According to our analysis of the experimental results, we observe that our MRs can detect the majority of these faulty versions, which shows the utility of this testing technique for quality assurance of bioinformatics software.


title: Controlling magnetism in 2D CrI3 by electrostatic doping

abstract: The atomic thickness of two-dimensional (2D) materials provides a unique opportunity to control material properties and engineer new functionalities by electrostatic doping. Electrostatic doping has been demonstrated to tune the electrical and optical properties of 2D materials in a wide range, as well as to drive the electronic phase transitions. The recent discovery of atomically thin magnetic insulators has opened up the prospect of electrical control of magnetism and new devices with unprecedented performance. Here we demonstrate control of the magnetic properties of monolayer and bilayer CrI3 by electrostatic doping using a dual-gate field-effect device structure. In monolayer CrI3, doping significantly modifies the saturation magnetization, coercive force and Curie temperature, showing strengthened (weakened) magnetic order with hole (electron) doping. Remarkably, in bilayer CrI3 doping drastically changes the interlayer magnetic order, causing a transition from an antiferromagnetic ground state in the pristine form to a ferromagnetic ground state above a critical electron density. The result reveals a strongly doping-dependent interlayer exchange coupling, which enables robust switching of magnetization in bilayer CrI3 by small gate voltages.


title: Hydrodynamic fluctuations in quasi-two dimensional diffusion

abstract: We study diffusion of colloids on a fluid-fluid interface using particle simulations and fluctuating hydrodynamics. Diffusion on a two-dimensional interface with three-dimensional hydrodynamics is known to be anomalous, with the collective diffusion coefficient diverging like the inverse of the wavenumber. This unusual collective effect arises because of the compressibility of the fluid flow in the plane of the interface, and leads to a nonlinear nonlocal convolution term in the diffusion equation for the ensemble-averaged concentration. We extend the previous hydrodynamic theory to account for a species/color labeling of the particles, as necessary to model experiments based on fluorescent techniques. We study the magnitude and dynamics of density and color density fluctuations using a novel Brownian dynamics algorithm, as well as fluctuating hydrodynamics theory and simulation. We find that hydrodynamic coupling between a single tagged particle and collective density fluctuations leads to a reduction of the long-time self-diffusion coefficient, even for an ideal gas of non-interacting particles. Using linearized fluctuating hydrodynamics theory, we show that for diffusion on a fluid-fluid interface, nonequilibrium fluctuations of the total density are small compared to the equilibrium fluctuations, but fluctuations of color density are giant and exhibit a spectrum that decays as the inverse cubed power of the wavenumber. We confirm these predictions through Brownian dynamics simulations of diffusive mixing with two indistinguishable species. We also examine nonequilibrium fluctuations in systems with two-dimensional hydrodynamics, such as thin smectic films in vacuum. We find that nonequilibrium fluctuations are colossal and comparable in magnitude to the mean, and can be accurately modeled using numerical solvers for the nonlinear equations of fluctuating hydrodynamics.


title: Tien Shan data on inelastic proton-air cross section at 10 PeV

abstract: New data on the absolute value on the inelastic proton-air cross section at 10 PeV are presented. Conclusion is made about the growth with energy of the inelastic proton-air cross section according to comparisons with experimental data were obtained at the Tien Shan complex array on various components of extensive air showers: hadrons, Cherenkov light and electrons with many different calculated models of cosmic rays interactions at the atmosphere. The analysis showed that the rise conforms to 7-9 % per one order of the energy from 0.2 TeV (accelerators with fixed targets) to 10 PeV (EAS cosmic rays). That corresponds to around 350 mb at 1 PeV and 380 mb at 10 PeV of primary cosmic rays. These data corresponds better to the QGSJET-II-04 model version.


title: Magnetize Topological Surface States of Bi2Se3 with a CrI3 Monolayer

abstract: To magnetize surfaces of topological insulators without damaging their topological feature is a crucial step for the realization of the quantum anomalous Hall effect (QAHE), and still remains as a challenging task. Through density functional calculations, we found that adsorption of a two-dimensional van der Waals (2D-vdW) ferromagnetic insulator CrI3 monolayer can create a sizable spin splitting at the Dirac point of the topological surface states of Bi2Se3 films. Furthermore, general rules that connect different quantum and topological parameters are established through model analyses. This work provides a useful guideline for the realization of QAHE at high temperature in heterostructures of 2D-vdW magnetic monolayers and topological insulators.


title: Positive temperature dynamics on Gelfand-Tsetlin patterns restricted by wall

abstract: The thesis focuses on processes on symplectic Gelfand-Tsetlin patterns. In chapter 4, a process with dynamics inspired by the Berele correspondence [Ber86] is presented. It is proved that the shape of the pattern is a Doob $h$-transform of independent random walks with $h$ given by the symplectic Schur function. This is followed by an extension to a $q$-weighted version. This randomised version has itself a branching structure and is related to a $q$-deformation of the $so_{2n+1}$-Whittaker functions. In chapter 5, we present a fully randomised process. This process $q$-deforms a process proposed in [WW09]. In chapter 7 we prove the convergence of the $q$-deformation of the $so_{2n+1}$-Whittaker functions to the classical $so_{2n+1}$-Whittaker functions when $q \to 1$. Finally, in chapter 8 we turn our interest to the continuous setting and construct a process on patterns which contains a positive temperature analogue of the Dyson's Brownian motion of type $B/C$. The processes obtained are $h$-transforms of Brownian motions killed at a continuous rate that depends on their distance from the boundary of the Weyl chamber of type $B/C$, with $h$ related with the $so_{2n+1}$-Whittaker functions.


title: Comparison of the Magnetic properties of Mn3Fe2Si3O12 as a crystalline garnet and as a glass

abstract: The crystalline garnet Mn3Fe2Si3O12 and an amorphous phase of the same nominal composition are synthesized at high pressure. The magnetic properties of the two forms are reported. Both phases order antiferromagnetically. The crystalline phase exhibits a Curie-Weiss theta of -47.2 K, with a sharp ordering transition at 12 K. The glassy phase exhibits a larger antiferromagnetic Curie-Weiss theta, of -83.0 K, with a broad ordering transition observed at 2.5 K. Both phases can be classified as magnetically frustrated, although the amorphous phase shows a much higher degree of frustration. The amorphous phase exhibits spin-glass behavior and is determined to have an actual composition of Mn3Fe2Si3O13.


title: Fault Detection Effectiveness of Source Test Case Generation Strategies for Metamorphic Testing

abstract: Metamorphic testing is a well known approach to tackle the oracle problem in software testing. This technique requires the use of source test cases that serve as seeds for the generation of follow-up test cases. Systematic design of test cases is crucial for the test quality. Thus, source test case generation strategy can make a big impact on the fault detection effectiveness of metamorphic testing. Most of the previous studies on metamorphic testing have used either random test data or existing test cases as source test cases. There has been limited research done on systematic source test case generation for metamorphic testing. This paper provides a comprehensive evaluation on the impact of source test case generation techniques on the fault finding effectiveness of metamorphic testing. We evaluated the effectiveness of line coverage, branch coverage, weak mutation and random test generation strategies for source test case generation. The experiments are conducted with 77 methods from 4 open source code repositories. Our results show that by systematically creating source test cases, we can significantly increase the fault finding effectiveness of metamorphic testing. Further, in this paper we introduce a simple metamorphic testing tool called "METtester" that we use to conduct metamorphic testing on these methods.


title: Statistical Software for Psychology: Comparing Development Practices Between CRAN and Other Communities

abstract: Different communities rely heavily on software, but use quite different software development practices. {\bf Objective}: We wanted to measure the state of the practice in the area of statistical software for psychology to understand how it compares to best practices. {\bf Method}: We compared and ranked 30 software tools with respect to adherence to best software engineering practices on items that could be measured by end-users. {\bf Results} We found that R packages use quite good practices, that while commercial packages were quite usable, many aspects of their development is too opaque to be measures, and that research projects vary a lot in their practices. {\bf Conclusion} We recommend that more organizations adopt practices similar to those used by CRAN to facilitate success, even for small teams. We also recommend close coupling of source code and documentation, to improve verifiability.


title: Mean Field Approximations to a Queueing System with Threshold-Based Workload Control Scheme

abstract: In this paper, motivated by considerations of server utilization and energy consumptions in cloud computing, we investigate a homogeneous queueing system with a threshold-based workload control scheme. In this system, a virtual machine will be turned off when there are no tasks in its buffer upon the completion of a service by the machine, and turned on when the number of tasks in its buffer reaches a pre-set threshold value. Due to complexity of this system, we propose approximations to system performance measures by mean field limits. An iterative algorithm is suggested for the solution to the mean field limit equations. In addition, numerical and simulation results are presented to justify the proposed approximation method and to provide a numerical analysis on the impact of the system performances by system parameters.


title: Controlling stability and transport of magnetic microswimmers by an external field

abstract: We investigate the hydrodynamic stability and transport of magnetic microswimmers in an external field using a kinetic theory framework. Combining linear stability analysis and nonlinear 3D continuum simulations, we show that for sufficiently large activity and magnetic field strengths, a homogeneous polar steady state is unstable for both puller and pusher swimmers. This instability is caused by the amplification of anisotropic hydrodynamic interactions due to the external alignment and leads to a partial depolarization and a reduction of the average transport speed of the swimmers in the field direction. Notably, at higher field strengths a reentrant hydrodynamic stability emerges where the homogeneous polar state becomes stable and a transport efficiency identical to that of active particles without hydrodynamic interactions is restored.


title: Lorentz-violating contributions to the nuclear Schiff moment and nuclear EDM

abstract: In the context of an atom endowed with nuclear electric dipole moment (EDM), we consider the effects on the Schiff moment of $CPT$-even Lorentz-violating (LV) terms that modify the Coulomb potential. First, we study the modifications on the Schiff moment when the nucleus interacts with the electronic cloud by means of a Coulomb potential altered only by the $P$-even LV components. Next, by supposing the existence of an additional intrinsic LV EDM generated by other LV sources, we assess the corrections to the Schiff moment when the interaction nucleus-electrons runs mediated by a Coulomb potential modified by both the $P$-odd and $P$-even LV components. We then use known estimates and EDM measurements to discuss upper bounds on the new Schiff moment components and the possibility of an intrisic nuclear EDM component ascribed to LV effects.


title: Free complete Wasserstein algebras

abstract: We present an algebraic account of the Wasserstein distances $W_p$ on complete metric spaces, for $p \geq 1$. This is part of a program of a quantitative algebraic theory of effects in programming languages. In particular, we give axioms, parametric in $p$, for algebras over metric spaces equipped with probabilistic choice operations. The axioms say that the operations form a barycentric algebra and that the metric satisfies a property typical of the Wasserstein distance $W_p$. We show that the free complete such algebra over a complete metric space is that of the Radon probability measures with finite moments of order $p$, equipped with the Wasserstein distance as metric and with the usual binary convex sums as operations.


title: Strain Fields in Repulsive Colloidal Crystals

abstract: The concept of a local linear elastic strain field is commonly used in the metallurgical research community to approximate the collective effect of atomic displacements around crystalline defects. Here we show that the elastic strain field approximation is a useful tool in colloidal systems. For colloidal crystals with repulsive particle interaction potentials, given similar mechanical properties, sharper potentials lead to: 1) free energies of deformation dominated by entropy, 2) lower variance in strain field fluctuations, 3) increased tension-compression asymmetry near dislocation core regions, and 4) smaller windows of applicability of the linear elastic approximation. We show that the window of linear behavior for entropic colloidal crystals is broadened for pressures at which the inter-particle separation sufficiently exceeds the range of steep repulsive interactions.


title: Gaussian Random Number Generator: Implemented in FPGA for Quantum Key Distribution

abstract: Quantum Key Distribution is the process of using quantum communication to establish a shared key between two parties. It has been demonstrated the unconditional security and effective communication of quantum communication system can be guaranteed by an excellent Gaussian random number generator with high speed and an extended random period. In this paper, we propose to construct the Gaussian random number generator using Field-Programmable Gate Array (FPGA) which is able to process large data in high speed. We also compare three algorithms ofGRNgeneration: Box-Muller algorithm, polarization decision algorithm, and central limit algorithm. We demonstrate that the polarization decision algorithm implemented inFPGArequires less computing resources and also produces a high-quality Gaussian random number, through the null hypothesis test.


title: On the Statistical Challenges of Echo State Networks and Some Potential Remedies

abstract: Echo state networks are powerful recurrent neural networks. However, they are often unstable and shaky, making the process of finding an good ESN for a specific dataset quite hard. Obtaining a superb accuracy by using the Echo State Network is a challenging task. We create, develop and implement a family of predictably optimal robust and stable ensemble of Echo State Networks via regularizing the training and perturbing the input. Furthermore, several distributions of weights have been tried based on the shape to see if the shape of the distribution has the impact for reducing the error. We found ESN can track in short term for most dataset, but it collapses in the long run. Short-term tracking with large size reservoir enables ESN to perform strikingly with superior prediction. Based on this scenario, we go a further step to aggregate many of ESNs into an ensemble to lower the variance and stabilize the system by stochastic replications and bootstrapping of input data.


title: SufiSent - Universal Sentence Representations Using Suffix Encodings

abstract: Computing universal distributed representations of sentences is a fundamental task in natural language processing. We propose a method to learn such representations by encoding the suffixes of word sequences in a sentence and training on the Stanford Natural Language Inference (SNLI) dataset. We demonstrate the effectiveness of our approach by evaluating it on the SentEval benchmark, improving on existing approaches on several transfer tasks.


title: Gravitational Wave Opacity from Gauge Field Dark Energy

abstract: We show that astrophysical gravitational waves can undergo an anomalous modulation when propagating through cosmic gauge field dark energy. A sufficiently strong effect, dependent on the gauge field energy density, would appear as a redshift-dependent opacity, thereby impacting the use of gravitational wave standard sirens to constrain the expansion history of the Universe. We investigate a particular model of cosmic gauge field dark energy and show that at early times it behaves like dark radiation, whereas a novel interaction causes it to drive cosmic acceleration at late times. Joint constraints on the cosmological scenario due to type 1a supernovae, baryon acoustic oscillations, and cosmic microwave background data are presented. In view of these constraints, we show that standard siren luminosity distances in the redshift range 0.5 < z < 1.5 would systematically dim by up to 1%, which may be distinguishable by third-generation gravitational wave detectors.


title: Stochastic Variance-Reduced Cubic Regularization for Nonconvex Optimization

abstract: Cubic regularization (CR) is an optimization method with emerging popularity due to its capability to escape saddle points and converge to second-order stationary solutions for nonconvex optimization. However, CR encounters a high sample complexity issue for finite-sum problems with a large data size. %Various inexact variants of CR have been proposed to improve the sample complexity. In this paper, we propose a stochastic variance-reduced cubic-regularization (SVRC) method under random sampling, and study its convergence guarantee as well as sample complexity. We show that the iteration complexity of SVRC for achieving a second-order stationary solution within $\epsilon$ accuracy is $O(\epsilon^{-3/2})$, which matches the state-of-art result on CR types of methods. Moreover, our proposed variance reduction scheme significantly reduces the per-iteration sample complexity. The resulting total Hessian sample complexity of our SVRC is ${\Oc}(N^{2/3} \epsilon^{-3/2})$, which outperforms the state-of-art result by a factor of $O(N^{2/15})$. We also study our SVRC under random sampling without replacement scheme, which yields a lower per-iteration sample complexity, and hence justifies its practical applicability.


title: Polynomial convolutions in max-plus algebra

abstract: Recently, in a work that grew out of their exploration of interlacing polynomials, Marcus, Spielman and Srivastava and then Marcus studied certain combinatorial polynomial convolutions. These convolutions preserve real-rootedness and capture expectations of characteristic polynomials of unitarily invariant random matrices, thus providing a link to free probability. We explore analogues of these types of convolutions in the setting of max-plus algebra. In this setting the max-permanent replaces the determinant, the maximum is the analogue of the expected value and real-rootedness is replaced by full canonical form. Our results resemble those of Marcus et al., however, in contrast to the classical setting we obtain an exact and simple description of all roots.


title: On the scaling of polynomial features for representation matching

abstract: In many neural models, new features as polynomial functions of existing ones are used to augment representations. Using the natural language inference task as an example, we investigate the use of scaled polynomials of degree 2 and above as matching features. We find that scaling degree 2 features has the highest impact on performance, reducing classification error by 5% in the best models.


title: Periodicity in Data Streams with Wildcards

abstract: We investigate the problem of detecting periodic trends within a string $S$ of length $n$, arriving in the streaming model, containing at most $k$ wildcard characters, where $k=o(n)$. A wildcard character is a special character that can be assigned any other character. We say $S$ has wildcard-period $p$ if there exists an assignment to each of the wildcard characters so that in the resulting stream the length $n-p$ prefix equals the length $n-p$ suffix. We present a two-pass streaming algorithm that computes wildcard-periods of $S$ using $\mathcal{O}(k^3\,\mathsf{polylog}\,n)$ bits of space, while we also show that this problem cannot be solved in sublinear space in one pass. We then give a one-pass randomized streaming algorithm that computes all wildcard-periods $p$ of $S$ with $p<\frac{n}{2}$ and no wildcard characters appearing in the last $p$ symbols of $S$, using $\mathcal{O}(k^3\mathsf{polylog}\, n)$ space.


title: Core Emergence in a Massive Infrared Dark Cloud: A Comparison Between Mid-IR Extinction and 1.3 mm Emission

abstract: Stars are born from dense cores in molecular clouds. Observationally, it is crucial to capture the formation of cores in order to understand the necessary conditions and rate of the star formation process. The {\it Atacama Large Mm/sub-mm Array} (ALMA) is extremely powerful for identifying dense gas structures, including cores, at mm wavelengths via their dust continuum emission. Here we use ALMA to carry out a survey of dense gas and cores in the central region of the massive ($\sim10^5\:M_\odot$) Infrared Dark Cloud (IRDC) G28.37+0.07. The observation consists of a mosaic of 86 pointings of the 12m-array and produces an unprecedented view of the densest structures of this IRDC. In this first paper about this data set, we focus on a comparison between the 1.3 mm continuum emission and a mid-infrared (MIR) extinction map of the IRDC. This allows estimation of the "dense gas" detection probability function (DPF), i.e., as a function of the local mass surface density, $\Sigma$, for various choices of thresholds of mm continuum emission to define "dense gas". We then estimate the dense gas mass fraction, $f_{\rm dg}$, in the central region of the IRDC and, via extrapolation with the DPF and the known $\Sigma$ probability distribution function, to the larger-scale surrounding regions, finding values of about 5\% to 15\% for the fiducial choice of threshold. We argue that this observed dense gas is a good tracer of the protostellar core population and, in this context, estimate a star formation efficiency per free-fall time in the central IRDC region of $\epsilon_{\rm ff}\sim$10\%, with approximately a factor of two systematic uncertainties.


title: C*-Algebras for partial product systems over N

abstract: We define partial product systems over N. They generalise product systems over N and Fell bundles over Z. We define Toeplitz C*-algebras and relative Cuntz-Pimsner algebras for them and show that the section C*-algebra of a Fell bundle over Z is a relative Cuntz-Pimsner algebra. We describe the gauge-invariant ideals in the Toeplitz C*-algebra.


title: Perturbation theories behind thermal mode spectroscopy for high-accuracy measurement of thermal diffusivity of solids

abstract: Thermal mode spectroscopy (TMS) has been recently proposed for accurately measuring thermal diffusivity of solids from a temperature decay rate of a specific thermal mode selected by three- dimensional (anti)nodal information [Phys. Rev. Lett., 117, 195901 (2016)]. In this paper, we find out the following advantages of TMS by use of perturbation analyses. First, TMS is applicable to the measurement of high thermal diffusivity with a small size specimen. Second, it is less affected by thermally resistive films on a specimen in the sense that the resistance at the interface does not affect the first-order correction of thermal diffusivity. Third, it can perform doubly accurate measurement of the thermal diffusivity specified at a thermal equilibrium state even if the diffusivity depends on temperature in the sense that the measurement can be performed within tiny temperature difference from the given state and that the decay rate of the slowest decaying mode is not affected by the dependence.


title: Scalable Label Propagation for Multi-relational Learning on the Tensor Product of Graphs

abstract: Multi-relational learning on knowledge graphs infers high-order relations among the entities across the graphs. This learning task can be solved by label propagation on the tensor product of the knowledge graphs to learn the high-order relations as a tensor. In this paper, we generalize a widely used label propagation model to the normalized tensor product graph, and propose an optimization formulation and a scalable Low-rank Tensor-based Label Propagation algorithm (LowrankTLP) to infer multi-relations for two learning tasks, hyperlink prediction and multiple graph alignment. The optimization formulation minimizes the upper bound of the noisy tensor estimation error for multiple graph alignment, by learning with a subset of the eigen-pairs in the spectrum of the normalized tensor product graph. We also provide a data-dependent transductive Rademacher bound for binary hyperlink prediction. We accelerate LowrankTLP with parallel tensor computation which enables label propagation on a tensor product of 100 graphs each of size 1000 in less than half hour in the simulation. LowrankTLP was also applied to predicting the author-paper-venue hyperlinks in publication records, alignment of segmented regions across up to 26 CT-scan images and alignment of protein-protein interaction networks across multiple species. The experiments demonstrate that LowrankTLP indeed well approximates the original label propagation with better scalability and accuracy.


title: Fast Nonconvex Deconvolution of Calcium Imaging Data

abstract: Calcium imaging data promises to transform the field of neuroscience by making it possible to record from large populations of neurons simultaneously. However, determining the exact moment in time at which a neuron spikes, from a calcium imaging data set, amounts to a non-trivial deconvolution problem which is of critical importance for downstream analyses. While a number of formulations have been proposed for this task in the recent literature, in this paper we focus on a formulation recently proposed in Jewell and Witten (2017) which has shown initial promising results. However, this proposal is slow to run on fluorescence traces of hundreds of thousands of timesteps. Here we develop a much faster online algorithm for solving the optimization problem of Jewell and Witten (2017) that can be used to deconvolve a fluorescence trace of 100,000 timesteps in less than a second. Furthermore, this algorithm overcomes a technical challenge of Jewell and Witten (2017) by avoiding the occurrence of so-called "negative" spikes. We demonstrate that this algorithm has superior performance relative to existing methods for spike deconvolution on calcium imaging datasets that were recently released as part of the spikefinder challenge (http://spikefinder.codeneuro.org/). Our C++ implementation, along with R and python wrappers, is publicly available on Github at https://github.com/jewellsean/FastLZeroSpikeInference.


title: How to Subvert Backdoored Encryption: Security Against Adversaries that Decrypt All Ciphertexts

abstract: We study secure and undetectable communication in a world where governments can read all encrypted communications of citizens. We consider a world where the only permitted communication method is via a government-mandated encryption scheme, using government-mandated keys. Citizens caught trying to communicate otherwise (e.g., by encrypting strings which do not appear to be natural language plaintexts) will be arrested. The one guarantee we suppose is that the government-mandated encryption scheme is semantically secure against outsiders: a perhaps advantageous feature to secure communication against foreign entities. But what good is semantic security against an adversary that has the power to decrypt? Even in this pessimistic scenario, we show citizens can communicate securely and undetectably. Informally, there is a protocol between Alice and Bob where they exchange ciphertexts that look innocuous even to someone who knows the secret keys and thus sees the corresponding plaintexts. And yet, in the end, Alice will have transmitted her secret message to Bob. Our security definition requires indistinguishability between unmodified use of the mandated encryption scheme, and conversations using the mandated encryption scheme in a modified way for subliminal communication. Our topics may be thought to fall broadly within the realm of steganography: the science of hiding secret communication in innocent-looking messages, or cover objects. However, we deal with the non-standard setting of adversarial cover object distributions (i.e., a stronger-than-usual adversary). We leverage that our cover objects are ciphertexts of a secure encryption scheme to bypass impossibility results which we show for broader classes of steganographic schemes. We give several constructions of subliminal communication schemes based on any key exchange protocol with random messages (e.g., Diffie-Hellman).


title: Generic Coreset for Scalable Learning of Monotonic Kernels: Logistic Regression, Sigmoid and more

abstract: Coreset (or core-set) is a small weighted \emph{subset} $Q$ of an input set $P$ with respect to a given \emph{monotonic} function $f:\mathbb{R}\to\mathbb{R}$ that \emph{provably} approximates its fitting loss $\sum_{p\in P}f(p\cdot x)$ to \emph{any} given $x\in\mathbb{R}^d$. Using $Q$ we can obtain approximation of $x^*$ that minimizes this loss, by running \emph{existing} optimization algorithms on $Q$. In this work we provide: (i) A lower bound which proves that there are sets with no coresets smaller than $n=|P|$ for general monotonic loss functions. (ii) A proof that, under a natural assumption that holds e.g. for logistic regression and the sigmoid activation functions, a small coreset exists for \emph{any} input $P$. (iii) A generic coreset construction algorithm that computes such a small coreset $Q$ in $O(nd+n\log n)$ time, and (iv) Experimental results which demonstrate that our coresets are effective and are much smaller in practice than predicted in theory.


title: Artinian algebras and Jordan type

abstract: The Jordan type of an element $\ell$ of the maximal ideal of an Artinian k-algebra A acting on an A-module M of k-dimension n, is the partition of n given by the Jordan block decomposition of the multiplication map $m_\ell$ on M. In general the Jordan type has more information than whether the pair $(\ell,M)$ is strong or weak Lefschetz. We develop basic properties of the Jordan type and their loci for modules over graded or local Artinian algebras. We as well study the relation of generic Jordan type of $A$ to the Hilbert function of $A$. We introduce and study a finer invariant, the Jordan degree type. In our last sections we give an overview of topics such as the Jordan types for Nagata idealizations, for modular tensor products, and for free extensions, including examples and some new results. We as well propose open problems.


title: Interpreting Neural Network Judgments via Minimal, Stable, and Symbolic Corrections

abstract: We present a new algorithm to generate minimal, stable, and symbolic corrections to an input that will cause a neural network with ReLU activations to change its output. We argue that such a correction is a useful way to provide feedback to a user when the network's output is different from a desired output. Our algorithm generates such a correction by solving a series of linear constraint satisfaction problems. The technique is evaluated on three neural network models: one predicting whether an applicant will pay a mortgage, one predicting whether a first-order theorem can be proved efficiently by a solver using certain heuristics, and the final one judging whether a drawing is an accurate rendition of a canonical drawing of a cat.


title: Universal Growth in Production Economies

abstract: We consider a simple variant of the von Neumann model of an expanding economy, in which multiple producers make goods according to their production function. The players trade their goods at the market and then use the bundles acquired for the production in the next round. We study a simple decentralized dynamic---known as proportional response---in which players update their bids proportionally to how useful the investments were in the past round. We show this dynamic leads to growth of the economy in the long term (whenever growth is possible) but also creates unbounded inequality, i.e. very rich and very poor players emerge. We analyze several other phenomena, such as how the relation of a player with others influences its development and the Gini index of the system. One of the key technical findings is that the players learn a global feature of the network (the optimal cycle) in a decentralized way, while interacting locally with their direct neighbors. We obtain this by studying the volume in the resulting dynamical system and showing that the volume of each cycle expands or contracts by a constant factor in each round.


title: Subspace Methods for 3-Parameter Eigenvalue Problems

abstract: We propose subspace methods for 3-parameter eigenvalue problems. Such problems arise when separation of variables is applied to separable boundary value problems; a particular example is the Helmholtz equation in ellipsoidal and paraboloidal coordinates. While several subspace methods for 2-parameter eigenvalue problems exist, their extensions to three parameter setting seem to be challenging. An inherent difficulty is that, while for 2-parameter eigenvalue problems we can exploit a relation to Sylvester equations to obtain a fast Arnoldi type method, such a relation does not seem to exist when there are three or more parameters. Instead, we introduce a subspace iteration method with projections onto generalized Krylov subspaces that are constructed from scratch at every iteration using certain Ritz vectors as the initial vectors. Another possibility is a Jacobi--Davidson type method for three or more parameters, which we generalize from its 2-parameter counterpart. For both approaches, we introduce a selection criterion for deflation that is based on the angles between left and right eigenvectors. The Jacobi--Davidson approach is devised to locate eigenvalues close to a prescribed target, yet it often also performs well when eigenvalues are sought based on the proximity of one of the components to a prescribed target. The subspace iteration method is devised specifically for the latter task. The proposed approaches are suitable especially for problems where the computation of several eigenvalues is required with high accuracy. Matlab implementations of both methods have been made available in the package MultiParEig.


title: Dynamics of Viscous Entrapped Saturated Zones in Partially Wetted Porous Media

abstract: As a typical multiphase fluid flow process, drainage in porous media is of fundamental interest in nature and industrial applications. During drainage processes in unsaturated soils and porous media in general, saturated clusters, in which a liquid phase fully occupies the pore space between solid grains, affect the relative permeability and effective stress of the system. In this study, we experimentally studied drainage processes in unsaturated granular media as a model porous system. The distribution of saturated clusters is analysed by an optical imaging method under different drainage conditions, in which pore-scale information from Voronoi and Delaunay tessellation was used to characterise the topology of saturated cluster distributions. By employing statistical analyses, the observed spatial and temporal information of multiphase flow and fluid entrapment in porous media are described. The results indicate that the distributions of both the crystallised cell size and pore size are positively correlated to the spatial and temporal distribution of saturated cluster sizes. The saturated cluster size is found to follow a lognormal distribution, in which the generalised Bond number correlates negatively to the scale parameter and positively to the shape parameter. These findings can be used to connect pore-scale behaviour with overall hydro-mechanical characteristics in partially saturated porous media, using both the degree of saturation and generalised Bond number.


title: Canonical heights on hyper-K\"ahler varieties and the Kawaguchi-Silverman conjecture

abstract: The Kawaguchi--Silverman conjecture predicts that if $f\colon X \dashrightarrow X$ is a dominant rational-self map of a projective variety over $\overline{\mathbb{Q}}$, and $P$ is a $\overline{\mathbb{Q}}$-point of $X$ with Zariski-dense orbit, then the dynamical and arithmetic degrees of $f$ coincide: $\lambda_1(f) = \alpha_f(P)$. We prove this conjecture in several higher-dimensional settings, including all endomorphisms of non-uniruled smooth projective threefolds with degree larger than $1$, and all endomorphisms of hyper-K\"ahler varieties in any dimension. In the latter case, we construct a canonical height function associated to any automorphism $f\colon X \to X$ of a hyper-K\"ahler variety defined over $\overline{\mathbb{Q}}$.


title: 3LC: Lightweight and Effective Traffic Compression for Distributed Machine Learning

abstract: The performance and efficiency of distributed machine learning (ML) depends significantly on how long it takes for nodes to exchange state changes. Overly-aggressive attempts to reduce communication often sacrifice final model accuracy and necessitate additional ML techniques to compensate for this loss, limiting their generality. Some attempts to reduce communication incur high computation overhead, which makes their performance benefits visible only over slow networks. We present 3LC, a lossy compression scheme for state change traffic that strikes balance between multiple goals: traffic reduction, accuracy, computation overhead, and generality. It combines three new techniques---3-value quantization with sparsity multiplication, quartic encoding, and zero-run encoding---to leverage strengths of quantization and sparsification techniques and avoid their drawbacks. It achieves a data compression ratio of up to 39--107X, almost the same test accuracy of trained models, and high compression speed. Distributed ML frameworks can employ 3LC without modifications to existing ML algorithms. Our experiments show that 3LC reduces wall-clock training time of ResNet-110--based image classifiers for CIFAR-10 on a 10-GPU cluster by up to 16--23X compared to TensorFlow's baseline design.


title: Quantum contextuality implies a logic that does not obey the principle of bivalence

abstract: In the paper, a value assignment for projection operators relating to a quantum system is equated with assignment of truth-values to the propositions associated with these operators. In consequence, the Kochen-Specker theorem (its localized variant, to be exact) can be treated as the statement that a logic of those projection operators does not obey the principle of bivalence. This implies that such a logic has a gappy (partial) semantics or many-valued semantics.


title: Glassy, Gardner-like Phenomenology in Minimally Polydisperse Crystalline Systems

abstract: We report on a non-equilibrium phase of matter, the minimally disordered crystal phase, which we find exists between the maximally amorphous glasses and the ideal crystal. Even though these near crystals appear highly ordered, they display glassy and jamming features akin to those observed in amorphous solids. Structurally, they exhibit a power-law scaling in their probability distribution of weak forces and small interparticle gaps as well as a flat density of vibrational states. Dynamically, they display anomalous aging above a characteristic pressure. Quantitatively this disordered crystal phase has much in common with the Gardner-like phase seen in maximally disordered solids. Near crystals should be amenable to experimental realizations in commercially-available particulate systems and are to be indispensable in verifying the theory of amorphous materials.


title: Chromospheric heating due to cancellation of quiet Sun internetwork fields

abstract: The heating of the solar chromosphere remains one of the most important questions in solar physics. Our current understanding is that small-scale internetwork (IN) magnetic fields play an important role as a heating agent. Indeed, cancellations of IN magnetic elements in the photosphere can produce transient brightenings in the chromosphere and transition region. These bright structures might be the signature of energy release and plasma heating, probably driven by magnetic reconnection of IN field lines. Although single events are not expected to release large amounts of energy, their global contribution to the chromosphere may be significant due to their ubiquitous presence in quiet Sun regions. In this paper we study cancellations of IN elements and analyze their impact on the energetics and dynamics of the quiet Sun atmosphere. We use high resolution, multiwavelength, coordinated observations obtained with the Interface Region Imaging Spectrograph (IRIS) and the Swedish 1-m Solar Telescope (SST) to identify cancellations of IN magnetic flux patches and follow their evolution. We find that, on average, these events live for ~3 minutes in the photosphere and ~12 minutes in the chromosphere and/or transition region. Employing multi-line inversions of the Mg II h & k lines we show that cancellations produce clear signatures of heating in the upper atmospheric layers. However, at the resolution and sensitivity accessible to the SST, their number density still seems to be one order of magnitude too low to explain the global chromospheric heating.


title: A lateral-type spin-photodiode based on Fe/x-AlOx/p-InGaAs junctions with a refracting-facet side window

abstract: A lateral-type spin-photodiode having a refracting facet on a side edge of the device is proposed and demonstrated at room temperature. The light shed horizontally on the side of the device is refracted and introduced directly into a thin InGaAs active layer under the spin-detecting Fe contact in which spin-polarized carriers are generated and injected into the Fe contact through a crystalline AlOx tunnel barrier. Experiments have been carried out with a circular polarization spectrometry set up, through which helicity-dependent photocurrent component, dI, is obtained with the conversion efficiency F ~ 0.4 %, where F is the ratio between dI and total photocurrent Iph. This value is the highest reported so far for pure lateral-type spin-photodiodes. It is discussed through analysis with a model consisting of drift-diffusion and quantum tunneling equations that a factor that limits the F value is unoccupied spin-polarized density-of-states of Fe in energy region into which spin-polarized electrons in a semiconductor are injected.


title: The "quantum" Turan problem for operator systems

abstract: Let V be a linear subspace of M_n(C) which contains the identity matrix and is stable under Hermitian transpose. A "quantum k-clique" for V is a rank k orthogonal projection P in M_n(C) for which dim(PVP) = k^2, and a "quantum k-anticlique" is a rank k orthogonal projection for which dim(PVP) = 1. We give upper and lower bounds both for the largest dimension of V which would ensure the existence of a quantum k-anticlique, and for the smallest dimension of V which would ensure the existence of a quantum k-clique.

categories: math.OA math.CO math.FA math.RA quant-ph


title: A Mixed Mimetic Spectral Element Model of the Rotating Shallow Water Equations on the Cubed Sphere

abstract: In a previous article [J. Comp. Phys. $\mathbf{357}$ (2018) 282-304], the mixed mimetic spectral element method was used to solve the rotating shallow water equations in an idealized geometry. Here the method is extended to a smoothly varying, non-affine, cubed sphere geometry. The differential operators are encoded topologically via incidence matrices due to the use of spectral element edge functions to construct tensor product solution spaces in $H(\mathrm{rot})$, $H(\mathrm{div})$ and $L_2$. These incidence matrices commute with respect to the metric terms in order to ensure that the mimetic properties are preserved independent of the geometry. This ensures conservation of mass, vorticity and energy for the rotating shallow water equations using inexact quadrature on the cubed sphere. The spectral convergence of errors are similarly preserved on the cubed sphere, with the generalized Piola transformation used to construct the metric terms for the physical field quantities.


title: Galois representations and ordinary reduction

abstract: We provide conditions on the p-adic Galois representation of a smooth proper variety over a complete nonarchimedean extension of Q_p to have (potentially) good ordinary reduction.


title: PTL-separability and closures for WQOs on words

abstract: We introduce a flexible class of well-quasi-orderings (WQOs) on words that generalizes the ordering of (not necessarily contiguous) subwords. Each such WQO induces a class of piecewise testable languages (PTLs) as Boolean combinations of upward closed sets. In this way, a range of regular language classes arises as PTLs. Moreover, each of the WQOs guarantees regularity of all downward closed sets. We consider two problems. First, we study which (perhaps non-regular) language classes permit a decision procedure to decide whether two given languages are separable by a PTL with respect to a given WQO. Second, we want to effectively compute downward closures with respect to these WQOs. Our first main result that for each of the WQOs, under mild assumptions, both problems reduce to the simultaneous unboundedness problem (SUP) and are thus solvable for many powerful system classes. In the second main result, we apply the framework to show decidability of separability of regular languages by $\mathcal{B}\Sigma_1[<, \mathsf{mod}]$, a fragment of first-order logic with modular predicates.


title: Investigating Rumor News Using Agreement-Aware Search

abstract: Recent years have witnessed a widespread increase of rumor news generated by humans and machines. Therefore, tools for investigating rumor news have become an urgent necessity. One useful function of such tools is to see ways a specific topic or event is represented by presenting different points of view from multiple sources. In this paper, we propose Maester, a novel agreement-aware search framework for investigating rumor news. Given an investigative question, Maester will retrieve related articles to that question, assign and display top articles from agree, disagree, and discuss categories to users. Splitting the results into these three categories provides the user a holistic view towards the investigative question. We build Maester based on the following two key observations: (1) relatedness can commonly be determined by keywords and entities occurring in both questions and articles, and (2) the level of agreement between the investigative question and the related news article can often be decided by a few key sentences. Accordingly, we use gradient boosting tree models with keyword/entity matching features for relatedness detection, and leverage recurrent neural network to infer the level of agreement. Our experiments on the Fake News Challenge (FNC) dataset demonstrate up to an order of magnitude improvement of Maester over the original FNC winning solution, for agreement-aware search.