title: The Fermat-type equations x^5 + y^5 = 2z^p or 3z^p solved through Q-curves

abstract: We solve the Diophantine equations $x^5 + y^5 = dz^p$ with $d=2, 3$ for a set of prime numbers of density 1/4, 1/2, respectively. The method consists in relating a possible solution to another Diophantine equation and solving the later by using Q-curves and a generalized modular technique as in work of Ellenberg and Dieulefait-Jimenez along with some new techniques for eliminating newforms.


title: Thermodynamics of deformed AdS$_5$ model with a positive/negative quadratic correction in graviton-dilaton system

abstract: By solving the Einstein equations of the graviton coupling with a real scalar dilaton field, we establish a general framework to self-consistently solve the geometric background with black-hole for any given phenomenological holographic models. In this framwork, we solve the black-hole background, the corresponding dilaon field and the dilaton potential for the deformed AdS$_5$ model with a positive/negative quadratic correction. We systematically investigate the thermodynamical properties of the deformed AdS$_5$ model with a positive and negative quadratic correction, respectively, and compare with lattice QCD on the results of the equation of state, the heavy quark potential, the Polyakov loop and the spatial Wilson loop. We find that the bulk thermodynamical properties are not sensitive to the sign of the quadratic correction, and the results of both deformed holographic QCD models agree well with lattice QCD result for pure SU(3) gauge theory. However, the results from loop operators favor a positive quadratic correction, which agree well with lattice QCD result. Especially, the result from the Polyakov loop excludes the model with a negative quadratic correction in the warp factor of ${\rm AdS}_5$.


title: Exploring the top and bottom of the quantum control landscape

abstract: A controlled quantum system possesses a search landscape defined by the target physical objective as a function of the controls. This paper focuses on the landscape for the transition probability Pif between the states of a finite level quantum system. Traditionally, the controls are applied fields; here we extend the notion of control to also include the Hamiltonian structure, in the form of time independent matrix elements. Level sets of controls that produce the same transition probability value are shown to exist at the bottom Pif = 0.0 and top Pif = 1.0 of the landscape with the field and/or Hamiltonian structure as controls. We present an algorithm to continuously explore these level sets starting from an initial point residing at either extreme value of Pif . The technique can also identify control solutions that exhibit the desirable properties of (a) robustness at the top and (b) the ability to rapidly rise towards an optimal control from the bottom. Numerical simulations are presented to illustrate the varied control behavior at the top and bottom of the landscape for several simple model systems.


title: Time Evolution with the DMRG Algorithm: A Generic Implementation for Strongly Correlated Electronic Systems

abstract: A detailed description of the time-step-targetting time evolution method within the DMRG algorithm is presented. The focus of this publication is on the implementation of the algorithm, and on its generic application. The case of one-site excitations within a Hubbard model is analyzed as a test for the algorithm, using open chains and two-leg ladder geometries. The accuracy of the procedure in the case of the recently discussed holon-doublon photo excitations of Mott insulators is also analyzed. Performance and parallelization issues are discussed. In addition, the full open-source code is provided as supplementary material.


title: Uniformisation in dimension four: towards a conjecture of Iitaka

abstract: Let X be a compact K\"ahler manifold whose universal covering is $\mathbb C^n$. A conjecture of Iitaka claims that X is a torus, up to finite \'etale cover. We prove this conjecture in various cases in dimension four. We also show that in the projective case Iitaka's conjecture is a consequence of the non-vanishing conjecture.


title: The 14N(p,gamma)15O reaction studied with a composite germanium detector

abstract: The rate of the carbon-nitrogen-oxygen (CNO) cycle of hydrogen burning is controlled by the 14N(p,gamma)15O reaction. The reaction proceeds by capture to the ground states and several excited states in O-15. In order to obtain a reliable extrapolation of the excitation curve to astrophysical energy, fits in the R-matrix framework are needed. In an energy range that sensitively tests such fits, new cross section data are reported here for the four major transitions in the 14N(p,gamma)15O reaction. The experiment has been performed at the Laboratory for Underground Nuclear Astrophysics (LUNA) 400 kV accelerator placed deep underground in the Gran Sasso facility in Italy. Using a composite germanium detector, summing corrections have been considerably reduced with respect to previous studies. The cross sections for capture to the ground state and to the 5181, 6172, and 6792 keV excited states in O-15 have been determined at 359, 380, and 399 keV beam energy. In addition, the branching ratios for the decay of the 278 keV resonance have been remeasured.


title: Theoretical approaches to the physics of spectral line polarization

abstract: Due to the continuous developments in polarimetric instrumentation, which will become even more dramatic in the near future with the availability of new generation solar telescopes, we are now severely confronted with a variety of new detailed observations of high diagnostic potential, whose interpretation requires a firmly established theoretical framework. In this contribution, I review the fundamental physical processes that underlie the generation and transfer of polarized radiation in stellar atmospheres, and I discuss the present status of the theoretical schemes now available, pointing out their main successes and limitations. I also present some ideas about the theoretical improvements that I consider necessary to achieve a correct interpretation of the complex phenomenology shown by polarimetric observations, focusing particularly on the second solar spectrum, which can be considered as one of the most important test benches of the theory.


title: Implications of the Babinet Principle for Casimir Interactions

abstract: We formulate the Babinet Principle (BP) as a relation between the scattering amplitudes for electromagnetic waves, and combine it with multiple scattering techniques to derive new properties of Casimir forces. We show that the Casimir force exerted by a planar conductor or dielectric on a self- complementary perforated planar mirror is approximately half that on a uniform mirror independent of the distance between them. The BP suggests that Casimir edge effects are anomalously small, supporting results obtained earlier in special cases. Finally, we illustrate how the BP can be used to estimate Casimir forces between perforated planar mirrors.


title: The effect of primordial non-Gaussianity on the skeleton of cosmic shear maps

abstract: (abridged) We explore the imprints of deviations from Gaussian primordial density fluctuations on the skeleton of the large-scale matter distribution as mapped through cosmological weak lensing. We computed the skeleton length of simulated effective convergence maps covering $\sim 35$ sq. deg each, extracted from a suite of cosmological $n-$body runs with different levels of local primordial non-Gaussianity. The latter is expected to alter the structure formation process with respect to the fiducial Gaussian scenario, and thus to leave a signature on the cosmic web. We found that alterations of the initial conditions consistently modify both the cumulative and the differential skeleton length, although the effect is generically smaller than the cosmic variance and depends on the smoothing of the map prior to the skeleton computation. Nevertheless, the qualitative shape of these deviations is rather similar to their primordial counterparts, implying that skeleton statistics retain good memory of the initial conditions. We performed a statistical analysis in order to find out at what Confidence Level primordial non-Gaussianity could be constrained by the skeleton test on cosmic shear maps of the size we adopted. At 68.3% Confidence Level we found an error on the measured level of primordial non-Gaussianity of $\Delta f_\mathrm{NL}\sim 300$, while at 90% Confidence Level it is of $\Delta f_\mathrm{NL}\sim 500$. While these values by themselves are not competitive with the current constraints, weak lensing maps larger than those used here would have a smaller field-to-field variance, and thus would likely lead to tighter constraints. A rough estimate indicates $\Delta f_\mathrm{NL} \sim$ a few tens at 68.3% Confidence Level for an all-sky weak lensing survey.


title: Analysis of a fully packed loop model arising in a magnetic Coulomb phase

abstract: The Coulomb phase of spin ice, and indeed the Ic phase of water ice, naturally realise a fully-packed two-colour loop model in three dimensions. We present a detailed analysis of the statistics of these loops, which avoid themselves and other loops of the same colour, and contrast their behaviour to an analogous two-dimensional model. The properties of another extended degree of freedom are also addressed, flux lines of the emergent gauge field of the Coulomb phase, which appear as "Dirac strings" in spin ice. We mention implications of these results for related models, and experiments.


title: Benford's Law Detects Quantum Phase Transitions similarly as Earthquakes

abstract: A century ago, it was predicted that the first significant digit appearing in a data would be nonuniformly distributed, with the number one appearing with the highest frequency. This law goes by the name of Benford's law. It holds for data ranging from infectious disease cases to national greenhouse gas emissions. Quantum phase transitions are cooperative phenomena where qualitative changes occur in many-body systems at zero temperature. We show that the century-old Benford's law can detect quantum phase transitions, much like it detects earthquakes. Therefore, being certainly of very different physical origins, seismic activity and quantum cooperative phenomena may be detected by similar methods. The result has immediate implications in precise measurements in experiments in general, and for realizable quantum computers in particular. It shows that estimation of the first significant digit of measured physical observables is enough to detect the presence of quantum phase transitions in macroscopic systems.


title: Parameter Estimation for Hidden Markov Models with Intractable Likelihoods

abstract: Approximate Bayesian computation (ABC) is a popular technique for approximating likelihoods and is often used in parameter estimation when the likelihood functions are analytically intractable. Although the use of ABC is widespread in many fields, there has been little investigation of the theoretical properties of the resulting estimators. In this paper we give a theoretical analysis of the asymptotic properties of ABC based maximum likelihood parameter estimation for hidden Markov models. In particular, we derive results analogous to those of consistency and asymptotic normality for standard maximum likelihood estimation. We also discuss how Sequential Monte Carlo methods provide a natural method for implementing likelihood based ABC procedures.


title: Optical theorem for Aharonov-Bohm scattering

abstract: Quantum-mechanical scattering off a magnetic vortex is considered, and the optical theorem is derived. The vortex core is assumed to be impermeable to scattered particles, and its transverse size is taken into account. We show that the scattering Aharonov-Bohm effect is independent of the choice of boundary conditions from the variety of the Robin ones. The behaviour of the scattering amplitude in the forward direction is analyzed, and the persistence of the Fraunhofer diffraction in the short-wavelength limit is shown to be crucial for maintaining the optical theorem in the quasiclassical limit.


title: Precision cosmology and 7Li data

abstract: At variance from 2H, 3He and 4He abundances, 7Li abundance data yield an extra constraint to cosmological parameters, in top of those deriving from CMB, BAO, SNIa or H0 data. This constraint, often disregarded, would favor smaller $\Omega_b h^2$ values, also indicating a preference for Dark Energy state equations well in the phantom regime, simultaneously softening the upper limit on the sum of neutrino masses, up to ~1.6eV.


title: Infrared Spectral Energy Distribution of Galaxies in the AKARI All Sky Survey: Correlations with Galaxy Properties, and Their Physical Origin

abstract: We have studied the properties of more than 1600 low-redshift galaxies by utilizing high-quality infrared flux measurements of the AKARI All-Sky Survey and physical quantities based on optical and 21-cm observations. Our goal is to understand the physics determining the infrared spectral energy distribution (SED). The ratio of the total infrared luminosity L_TIR, to the star-formation rate (SFR) is tightly correlated by a power-law to specific SFR (SSFR), and L_TIR is a good SFR indicator only for galaxies with the largest SSFR. We discovered a tight linear correlation for normal galaxies between the radiation field strength of dust heating, estimated by infrared SED fits (U_h), and that of galactic-scale infrared emission (U_TIR ~ L_TIR/R^2), where R is the optical size of a galaxy. The dispersion of U_h along this relation is 0.3 dex, corresponding to 13% dispersion in the dust temperature. This scaling and the U_h/U_TIR ratio can be explained physically by a thin layer of heating sources embedded in a thicker, optically-thick dust screen. The data also indicate that the heated fraction of the total dust mass is anti-correlated to the dust column density, supporting this interpretation. In the large U_TIR limit, the data of circumnuclear starbursts indicate the existence of an upper limit on U_h, corresponding to the maximum SFR per gas mass of ~ 10 Gyr^{-1}. We find that the number of galaxies sharply drops when they become optically thin against dust-heating radiation, suggesting that a feedback process to galaxy formation (likely by the photoelectric heating) is working when dust-heating radiation is not self-shielded on a galactic scale. Implications are discussed for the M_HI-size relation, the Kennicutt-Schmidt relation, and galaxy formation in the cosmological context.


title: On the stability of self-similar solutions of 1D cubic Schrodinger equations

abstract: In this paper we will study the stability properties of self-similar solutions of 1-d cubic NLS equations with time-dependent coefficients of the form iu_t+u_{xx}+\frac{u}{2} (|u|^2-\frac{A}{t})=0, A\in \R (cubic). The study of the stability of these self-similar solutions is related, through the Hasimoto transformation, to the stability of some singular vortex dynamics in the setting of the Localized Induction Equation (LIE), an equation modeling the self-induced motion of vortex filaments in ideal fluids and superfluids. We follow the approach used by Banica and Vega that is based on the so-called pseudo-conformal transformation, which reduces the problem to the construction of modified wave operators for solutions of the equation iv_t+ v_{xx} +\frac{v}{2t}(|v|^2-A)=0. As a by-product of our results we prove that equation (cubic) is well-posed in appropriate function spaces when the initial datum is given by u(0,x)= z_0 \pv \frac{1}{x} for some values of z_0\in \C\setminus\{0\}, and A is adequately chosen. This is in deep contrast with the case when the initial datum is the Dirac-delta distribution.


title: Model-independent analysis of Higgs spin and CP properties in the process $e^+ e^- \to t \bar t \Phi$

abstract: In this paper we investigate methods to study the $t\bar{t}$ Higgs coupling. The spin and CP properties of a Higgs boson are analysed in a model-independent way in its associated production with a $t\bar{t}$ pair in high-energy $e^+e^-$ collisions. We study the prospects of establishing the CP quantum numbers of the Higgs boson in the CP-conserving case as well as those of determining the CP-mixing if CP is violated. We explore in this analysis the combined use of the total cross section and its energy dependence, the polarisation asymmetry of the top quark and the up-down asymmetry of the antitop with respect to the top-electron plane. We find that combining all three observables remarkably reduces the error on the determination of the CP properties of the Higgs Yukawa coupling. Furthermore, the top polarisation asymmetry and the ratio of cross sections at different collider energies are shown to be sensitive to the spin of the particle produced in association with the top quark pair.


title: Network Estimation and Packet Delivery Prediction for Control over Wireless Mesh Networks

abstract: Much of the current theory of networked control systems uses simple point-to-point communication models as an abstraction of the underlying network. As a result, the controller has very limited information on the network conditions and performs suboptimally. This work models the underlying wireless multihop mesh network as a graph of links with transmission success probabilities, and uses a recursive Bayesian estimator to provide packet delivery predictions to the controller. The predictions are a joint probability distribution on future packet delivery sequences, and thus capture correlations between successive packet deliveries. We look at finite horizon LQG control over a lossy actuation channel and a perfect sensing channel, both without delay, to study how the controller can compensate for predicted network outages.


title: Fractional variational calculus for nondifferentiable functions

abstract: We prove necessary optimality conditions, in the class of continuous functions, for variational problems defined with Jumarie's modified Riemann-Liouville derivative. The fractional basic problem of the calculus of variations with free boundary conditions is considered, as well as problems with isoperimetric and holonomic constraints.


title: Data augmentation for non-Gaussian regression models using variance-mean mixtures

abstract: We use the theory of normal variance-mean mixtures to derive a data-augmentation scheme for a class of common regularization problems. This generalizes existing theory on normal variance mixtures for priors in regression and classification. It also allows variants of the expectation-maximization algorithm to be brought to bear on a wider range of models than previously appreciated. We demonstrate the method on several examples, including sparse quantile regression and binary logistic regression. We also show that quasi-Newton acceleration can substantially improve the speed of the algorithm without compromising its robustness.


title: Spectral Risk Measures with an Application to Futures Clearinghouse Variation Margin Requirements

abstract: This paper applies an AR(1)-GARCH (1, 1) process to detail the conditional distributions of the return distributions for the S&P500, FT100, DAX, Hang Seng, and Nikkei225 futures contracts. It then uses the conditional distribution for these contracts to estimate spectral risk measures, which are coherent risk measures that reflect a user's risk-aversion function. It compares these to more familiar VaR and Expected Shortfall (ES) measures of risk, and also compares the precision and discusses the relative usefulness of each of these risk measures in setting variation margins that incorporate time-varying market conditions. The goodness of fit of the model is confirmed by a variety of backtests.


title: Exponential Spectral Risk Measures

abstract: Spectral risk measures are attractive risk measures as they allow the user to obtain risk measures that reflect their subjective risk-aversion. This paper examines spectral risk measures based on an exponential utility function, and finds that these risk measures have nice intuitive properties. It also discusses how they can be estimated using numerical quadrature methods, and how confidence intervals for them can be estimated using a parametric bootstrap. Illustrative results suggest that estimated exponential spectral risk measures obtained using such methods are quite precise in the presence of normally distributed losses.


title: Political protest Italian-style: The dissonance between the blogosphere and mainstream media in the promotion and coverage of Beppe Grillo's V-day

abstract: We analyze the organization, promotion and public perception of V-day, a political rally that took place on September 8, 2007, to protest against corruption in the Italian Parliament. Launched by blogger Beppe Grillo, and promoted via a word of mouth mobilization on the Italian blogosphere, V-day brought close to one million Italians in the streets on a single day, but was mostly ignored by mainstream media. This article is divided into two parts. In the first part, we analyze the volume and content of online articles published by both bloggers and mainstream news sources from June 14 (the day V-day was announced) until September 15, 2007 (one week after it took place) . We find that the success of V-day can be attributed to the coverage of bloggers and small-scale local news outlets only, suggesting a strong grassroots component in the organization of the rally. We also find a dissonant thematic relationship between content published by blogs and mainstream media: while the majority of blogs analyzed promote V-day, major mainstream media sources critique the methods of information production and dissemination employed by Grillo. Based on this finding, in the second part of the study, we explore the role of Grillo in the organization of the rally from a network analysis perspective. We study the interlinking structure of the V-day blogosphere network, to determine its structure, its levels of heterogeneity, and resilience. Our analysis contradicts the hypothesis that Grillo served as a top-down, broadcast-like source of information. Rather, we find that information about V-day was transferred across heterogeneous nodes in a moderately robust and resilient core network of blogs. We speculate that the organization of V-day represents the very first case, in Italian history, of a political demonstration developed and promoted primarily via the use of social media on the web.


title: Hedging Effectiveness under Conditions of Asymmetry

abstract: We examine whether hedging effectiveness is affected by asymmetry in the return distribution by applying tail specific metrics to compare the hedging effectiveness of short and long hedgers using crude oil futures contracts. The metrics used include Lower Partial Moments (LPM), Value at Risk (VaR) and Conditional Value at Risk (CVAR). Comparisons are applied to a number of hedging strategies including OLS and both Symmetric and Asymmetric GARCH models. Our findings show that asymmetry reduces in-sample hedging performance and that there are significant differences in hedging performance between short and long hedgers. Thus, tail specific performance metrics should be applied in evaluating hedging effectiveness. We also find that the Ordinary Least Squares (OLS) model provides consistently good performance across different measures of hedging effectiveness and estimation methods irrespective of the characteristics of the underlying distribution.


title: Margin setting with high-frequency data1

abstract: Both in practice and in the academic literature, models for setting margin requirements in futures markets classically use daily closing price changes. However, as well documented by research on high-frequency data, financial markets have recently shown high intraday volatility, which could bring more risk than expected. This paper tries to answer two questions relevant for margin committees in practice: is it right to compute margin levels based on closing prices and ignoring intraday dynamics? Is it justified to implement intraday margin calls? The paper focuses on the impact of intraday dynamics of market prices on daily margin levels. Daily margin levels are obtained in two ways: first, by using daily price changes defined with different time-intervals (say from 3 pm to 3 pm on the following trading day instead of traditional closing times); second, by using 5-minute and 1-hour price changes and scaling the results to one day. Our empirical analysis uses the FTSE 100 futures contract traded on LIFFE.


title: Mixed axion/neutralino cold dark matter in supersymmetric models

abstract: We consider supersymmetric (SUSY) models wherein the strong CP problem is solved by the Peccei-Quinn (PQ) mechanism with a concommitant axion/axino supermultiplet. We examine R-parity conserving models where the neutralino is the lightest SUSY particle, so that a mixture of neutralinos and axions serve as cold dark matter. The mixed axion/neutralino CDM scenario can match the measured dark matter abundance for SUSY models which typically give too low a value of the usual thermal neutralino abundance, such as models with wino-like or higgsino-like dark matter. The usual thermal neutralino abundance can be greatly enhanced by the decay of thermally-produced axinos to neutralinos, followed by neutralino re-annihilation at temperatures much lower than freeze-out. In this case, the relic density is usually neutralino dominated, and goes as \sim (f_a/N)/m_{axino}^{3/2}. If axino decay occurs before neutralino freeze-out, then instead the neutralino abundance can be augmented by relic axions to match the measured abundance. Entropy production from late-time axino decays can diminish the axion abundance, but ultimately not the neutralino abundance. In mixed axion/neutralino CDM models, it may be possible to detect both a WIMP and an axion as dark matter relics. We also discuss possible modifications of our results due to production and decay of saxions. In the appendices, we present expressions for the Hubble expansion rate and the axion and neutralino relic densities in radiation, matter and decaying-particle dominated universes.


title: Modeling Long Memory in REITs

abstract: One stylized feature of financial volatility impacting the modeling process is long memory. This paper examines long memory for alternative risk measures, observed absolute and squared returns for Daily REITs and compares the findings for a non- REIT equity index. The paper utilizes a variety of tests for long memory finding evidence that REIT volatility does display persistence, in contrast to the actual return series. Trading volume is found to be strongly associated with long memory. The results do however suggest differences in the findings with regard to REITs in comparison to the broader equity sector which may be due to relatively thin trading during the sample period.


title: Generic decompositions and semi-invariants for string algebras

abstract: We investigate the rings of semi-invariants for tame string algebras A(n) of non-polynomial growth. We are interested in dimension vectors of band modules. We use geometric technique related to the description of coordinate rings on varieties of complexes. The fascinating combinatorics emerges, showing that our rings of invariants are the rings of some toric varieties. We show that for $n\le 6$ the rings of semi-invariants are complete intersections but we show an example for $n=7$ that this is not the case in general.


title: Minimum Capital Requirement Calculations for UK Futures

abstract: Key to the imposition of appropriate minimum capital requirements on a daily basis requires accurate volatility estimation. Here, measures are presented based on discrete estimation of aggregated high frequency UK futures realisations underpinned by a continuous time framework. Squared and absolute returns are incorporated into the measurement process so as to rely on the quadratic variation of a diffusion process and be robust in the presence of fat tails. The realized volatility estimates incorporate the long memory property. The dynamics of the volatility variable are adequately captured. Resulting rescaled returns are applied to minimum capital requirement calculations.


title: Uncovering Volatility Dynamics in Daily REIT Returns

abstract: Using a time-varying approach, this paper examines the dynamics of volatility in the REIT sector. The results highlight the attractiveness and suitability of using GARCH based approaches in the modeling of daily REIT volatility. The paper examines the influencing factors on REIT volatility, documenting the return and volatility linkages between REIT sub-sectors and also examines the influence of other US equity series. The results contrast with previous studies of monthly REIT volatility. Linkages within the REIT sector and with related sectors such as value stocks are diminished, while the general influence of market sentiment, coming through the large cap indices is enhanced. This would indicate that on a daily basis general market sentiment plays a more fundamental role than more intuitive relationships within the capital markets.


title: Tail Behaviour of the Euro

abstract: This paper empirically analyses risk in the Euro relative to other currencies. Comparisons are made between a sub period encompassing the final transitional stage to full monetary union with a sub period prior to this. Stability in the face of speculative attack is examined using Extreme Value Theory to obtain estimates of tail exchange rate changes. The findings are encouraging. The Euro's common risk measures do not deviate substantially from other currencies. Also, the Euro is stable in the face of speculative pressure. For example, the findings consistently show the Euro being less risky than the Yen, and having similar inherent risk to the Deutsche Mark, the currency that it is essentially replacing.


title: $D$ to $K$ semileptonic decay form factors from HISQ light and charm quarks

abstract: We present a study of $D \rightarrow K, l \nu$ semileptonic decays on the lattice which employs the HISQ action for both the charm and the light quarks. We work with MILC unquenched $N_f = 2 + 1$ lattices and determine the scalar form factor $f_0(q^2)$. This form factor is obtained from a scalar current matrix element that does not require any operator matching. We find $f^{D \rightarrow K}_0(0) \equiv f^{D \rightarrow K}_+(0) = 0.747(19)$ in the chiral plus continuum limit and hereby improve the theory error on this quantity by a factor of $\sim$4 compared to previous lattice determinations. Combining the new theory result with recent experimental measurements of the product $f^{D \rightarrow K}_+(0) * |V_{cs}| $ from BaBar and CLEO-c leads to a very precise direct determination of the CKM matrix element $|V_{cs}| $, $|V_{cs}| = 0.961(11)(24)$, where the first error comes from experiment and the second is the lattice QCD theory error.


title: Nonlinear transport effects in mass separation by effusion

abstract: Generalizations of Onsager reciprocity relations are established for the nonlinear response coefficients of ballistic transport in the effusion of gaseous mixtures. These generalizations, which have been established on the basis of the fluctuation theorem for the currents, are here considered for mass separation by effusion. In this kinetic process, the mean values of the currents depend nonlinearly on the affinities or thermodynamic forces controlling the nonequilibrium constraints. These nonlinear transport effects are shown to play an important role in the process of mass separation. In particular, the entropy efficiency turns out to be significantly larger than it would be the case if the currents were supposed to depend linearly on the affinities.


title: Context-free ordinals

abstract: We consider context-free languages equipped with the lexicographic ordering. We show that when the lexicographic ordering of a context-free language is scattered, then its Hausdorff rank is less than $\omega^\omega$. As a corollary of this result we obtain that an ordinal is the order type of a well-ordered context-free language iff it is less than $\omega^{\omega^\omega}$.


title: High-Dimensional Topological Insulators with Quaternionic Analytic Landau Levels

abstract: We study the 3D topological insulators in the continuum by coupling spin-1/2 fermions to the Aharonov-Casher SU(2) gauge field. They exhibit flat Landau levels in which orbital angular momentum and spin are coupled with a fixed helicity. The 3D lowest Landau level wavefunctions exhibit the quaternionic analyticity as a generalization of the complex analyticity of the 2D case. Each Landau level contributes one branch of gapless helical Dirac modes to the surface spectra, whose topological properties belong to the Z2-class. The flat Landau levels can be generalized to an arbitrary dimension. Interaction effects and experimental realizations are also studied.


title: Linearly repetitive Delone sets are rectifiable

abstract: In this paper we prove that, for any integer $d>0$, every linearly repetitive Delone set in the Euclidean $d$-space $\RR^d$ is equivalent, up to a bi-Lipschitz homeomorphism, to the integer lattice $\ZZ^d$. In the particular case when the Delone set $X$ in $\RR^d$ comes from a primitive substitution tiling of $\RR^d$, we give a condition on the eigenvalues of the substitution matrix which implies the existence of a homeomorphism with bounded displacement from $X$ to the lattice lattice $\lambda\ZZ^d$ for some positive $\lambda$. This condition includes primitive Pisot substitution tilings but also concerns a much broader set of substitution tilings.


title: InAs/GaAs (211)B quantum dots with negligible FSS for the generation of entangled photons

abstract: Polarization-resolved single dot spectroscopy performed on (211)B InAs/GaAs quantum dots reveals that the fine structure splitting of the excitonic levels in these dots is much lower compared to the usual (100)-grown InAs dots. Time-resolved measurements confirm the high oscillator strength of these dots, and thus their good quantum efficiency at 4 K, comparable with that of (100) InAs/GaAs dots. Last, photon correlation measurements demonstrate single photon emission out of the excitonic optical transition of these dots. All these features make this novel dot system very promising for implementing solid-state entangled photon sources.


title: Simple model squirmers with tunable velocity

abstract: We present a scheme of self-propelling liquid droplets which closely mimics the locomotion of some protozoal organisms, so-called squirmers. In contrast to other schemes proposed earlier, locomotion paths are not self-avoiding, since the effect of the squirmer on the surrounding medium is weak. Our results suggest that not only the velocity, but also the mode of operation (i.e., the spherical harmonics of the flow field) can be controlled by appropriate variation of parameters.


title: Interference Channels with Rate-Limited Feedback

abstract: We consider the two-user interference channel with rate-limited feedback. Related prior works focus on the case where feedback links have infinite capacity, while no research has been done for the rate-limited feedback problem. Several new challenges arise due to the capacity limitations of the feedback links, both in deriving inner-bounds and outer-bounds. We study this problem under three different interference models: the El Gamal-Costa deterministic model, the linear deterministic model, and the Gaussian model. For the first two models, we develop an achievable scheme that employs three techniques: Han-Kobayashi message splitting, quantize-and-binning, and decode-and-forward. We also derive new outer-bounds for all three models and we show the optimality of our scheme under the linear deterministic model. In the Gaussian case, we propose a transmission strategy that incorporates lattice codes, inspired by the ideas developed in the first two models. For symmetric channel gains, we prove that the gap between the achievable sum-rate of the proposed scheme and our new outer-bounds is bounded by a constant number of bits, independent of the channel gains.


title: The mean value of Frobenius numbers with three arguments

abstract: We prove an asymptotic formula for the mean value of Frobenius numbers with three arguments. To prove this we use a new method invented by A. Ustinov, Rodseth's algorithm an bounds for exponential sums.


title: Two-Dimensional Arrays of RF Ion Traps with Addressable Interactions

abstract: We describe the advantages of 2-dimensional, addressable arrays of spherical Paul traps. They would provide for the ability to address and tailor the interaction strengths of trapped objects in 2D and could establish a valuable new tool for quantum information processing. Simulations of trapping ions are compared to first tests using printed circuit board trap arrays loaded with dust particles. Pair-wise interactions in the array are addressed by means of an adjustable radio-frequency (RF) electrode shared between trapping sites. By attenuating this RF electrode potential, neighboring pairs of trapped objects have their interaction strength increase and are moved closer to one another. In the limit of the adjustable electrode being held at RF ground, the two formerly spherical traps are merged into one linear Paul trap.


title: A geometric characterization of a sharp Hardy inequality

abstract: In this paper, we prove that the distance function of an open connected set in $\mathbb R^{n+1}$ with a $C^{2}$ boundary is superharmonic in the distribution sense if and only if the boundary is {\em weakly mean convex}. We then prove that Hardy inequalities with a sharp constant hold on {weakly mean convex} $C^{2}$ domains. Moreover, we show that the {weakly mean convexity} condition cannot be weakened. We also prove various improved Hardy inequalities on mean convex domains along the line of Brezis-Marcus \cite{BM}.


title: Some Generalized Harmonic Number Identities

abstract: Summation by parts is used to find the sum of a finite series of generalized harmonic numbers involving a specific polynomial or rational function. The Euler-Maclaurin formula for sums of powers is used to find the sums of some finite series of generalized harmonic numbers involving nonnegative integer powers, which can be used to evaluate the sums of the finite series of generalized harmonic numbers involving polynomials. Many examples and a computer program are provided.


title: Identification of Nonlinear Systems with Stable Limit Cycles via Convex Optimization

abstract: We propose a convex optimization procedure for black-box identification of nonlinear state-space models for systems that exhibit stable limit cycles (unforced periodic solutions). It extends the "robust identification error" framework in which a convex upper bound on simulation error is optimized to fit rational polynomial models with a strong stability guarantee. In this work, we relax the stability constraint using the concepts of transverse dynamics and orbital stability, thus allowing systems with autonomous oscillations to be identified. The resulting optimization problem is convex, and can be formulated as a semidefinite program. A simulation-error bound is proved without assuming that the true system is in the model class, or that the number of measurements goes to infinity. Conditions which guarantee existence of a unique limit cycle of the model are proved and related to the model class that we search over. The method is illustrated by identifying a high-fidelity model from experimental recordings of a live rat hippocampal neuron in culture.


title: How to display science since images have no mass

abstract: Education, science, in fact the whole society, extensively use images. Between us and the world are the visual displays. Screens, small and large, individual or not, are everywhere. Images are increasingly the 2D substrate of our virtual interaction with reality. However images will never support a complete representation of the reality. Three-dimensional representations will not change that. Images are primarily a spatial representation of our world dedicated to our sight. Key aspects such as energy and the associated forces are not spatially materialized. In classical physics, interaction description is based on Newton equations with trajectory and force as the dual central concepts. Images can in real time show all aspects of trajectories but not the associated dynamical aspects described by forces and energies. Contrary to the real world, the world of images opposes no constrain, nor resistance to our actions. Only the physical quantities, that do not contain mass in their dimension can be satisfactory represented by images. Often symbols such as arrows are introduced to visualize the force vectors.


title: Multifaceted Faculty Network Design and Management: Practice and Experience Report

abstract: We report on our experience on multidimensional aspects of our faculty's network design and management, including some unique aspects such as campus-wide VLANs and ghosting, security and monitoring, switching and routing, and others. We outline a historical perspective on certain research, design, and development decisions and discuss the network topology, its scalability, and management in detail; the services our network provides, and its evolution. We overview the security aspects of the management as well as data management and automation and the use of the data by other members of the IT group in the faculty.


title: Unification of Pulses in Long and Short Gamma-Ray Bursts: Evidence from Pulse Properties and their Correlations

abstract: We demonstrate that distinguishable gamma-ray burst pulses exhibit similar behaviors as evidenced by correlations among the observable pulse properties of duration, peak luminosity, fluence, spectral hardness, energy-dependent lag, and asymmetry. Long and Short burst pulses exhibit these behaviors, suggesting that a similar process is responsible for producing all GRB pulses. That these properties correlate in the observer's frame indicates that intrinsic correlations are strong enough to not be diluted into insignificance by the dispersion in distances and redshift. We show how all correlated pulse characteristics can be explained by hard-to-soft pulse evolution, and we demonstrate that "intensity tracking" pulses not having these properties are not single pulses; they instead appear to be composed of two or more overlapping hard-to-soft pulses. In order to better understand pulse characteristics, we recognize that hard-to-soft evolution provides a more accurate definition of a pulse than its intensity variation. This realization, coupled with the observation that pulses begin near-simultaneously across a wide range of energies, leads us to conclude that the observed pulse emission represents the energy decay resulting from an initial injection, and that one simple and as yet unspecified physical mechanism is likely to be responsible for all gamma-ray burst pulses regardless of the environment in which they form and, if GRBs originate from different progenitors, then of the progenitors that supply them with energy.


title: Efficient Algorithms for Optimal Control of Quantum Dynamics: The "Krotov'' Method unencumbered

abstract: Efficient algorithms for the discovery of optimal control designs for coherent control of quantum processes are of fundamental importance. One important class of algorithms are sequential update algorithms generally attributed to Krotov. Although widely and often successfully used, the associated theory is often involved and leaves many crucial questions unanswered, from the monotonicity and convergence of the algorithm to discretization effects, leading to the introduction of ad-hoc penalty terms and suboptimal update schemes detrimental to the performance of the algorithm. We present a general framework for sequential update algorithms including specific prescriptions for efficient update rules with inexpensive dynamic search length control, taking into account discretization effects and eliminating the need for ad-hoc penalty terms. The latter, while necessary to regularize the problem in the limit of infinite time resolution, i.e., the continuum limit, are shown to be undesirable and unnecessary in the practically relevant case of finite time resolution. Numerical examples show that the ideas underlying many of these results extend even beyond what can be rigorously proved.


title: Relaxation of spherical systems with long-range interactions: a numerical investigation

abstract: The process of relaxation of a system of particles interacting with long-range forces is relevant to many areas of Physics. For obvious reasons, in Stellar Dynamics much attention has been paid to the case of 1/r^2 force law. However, recently the interest in alternative gravities emerged, and significant differences with respect to Newtonian gravity have been found in relaxation phenomena. Here we begin to explore this matter further, by using a numerical model of spherical shells interacting with an 1/r^alpha force law obeying the superposition principle. We find that the virialization and phase-mixing times depend on the exponent alpha, with small values of alpha corresponding to longer relaxation times, similarly to what happens when comparing for N-body simulations in classical gravity and in Modified Newtonian Dynamics.


title: General Relationship Between the Entanglement Spectrum and the Edge State Spectrum of Topological Quantum States

abstract: We consider (2+1)-dimensional topological quantum states which possess edge states described by a chiral (1+1)-dimensional Conformal Field Theory (CFT), such as e.g. a general quantum Hall state. We demonstrate that for such states the reduced density matrix of a finite spatial region of the gapped topological state is a thermal density matrix of the chiral edge state CFT which would appear at the spatial boundary of that region. We obtain this result by applying a physical instantaneous cut to the gapped system, and by viewing the cutting process as a sudden "quantum quench" into a CFT, using the tools of boundary conformal field theory. We thus provide a demonstration of the observation made by Li and Haldane about the relationship between the entanglement spectrum and the spectrum of a physical edge state.


title: The VLBA Calibrator Search for the BeSSeL Survey

abstract: We present the results of a survey of radio continuum sources near the Galactic plane using the Very Long Baseline Array (VLBA). Our observations are designed to identify compact extragalactic sources of milliarcsecond size that can be used for parallax measurements in the Bar and Spiral Structure Legacy Survey. We selected point sources from the NVSS and CORNISH catalogs with flux densities above 30 mJy and within $1.5\degr$ of known maser targets. Of the 1529 sources observed, 199 were detected. For sources detected on 3 or more baselines, we determined accurate positions and evaluated their quality as potential calibrators. Most of the 1330 sources that were not detected with the VLBA are probably of extragalactic origin.


title: Composition-tuned smeared phase transitions

abstract: Phase transitions in random systems are smeared if individual spatial regions can order independently of the bulk system. In this paper, we study such smeared phase transitions (both classical and quantum) in substitutional alloys A$_{1-x}$B$_x$ that can be tuned from an ordered phase at composition $x=0$ to a disordered phase at $x=1$. We show that the ordered phase develops a pronounced tail that extends over all compositions $x<1$. Using optimal fluctuation theory, we derive the composition dependence of the order parameter and other quantities in the tail of the smeared phase transition. We also compare our results to computer simulations of a toy model, and we discuss experiments.


title: Motion of position-dependent mass as a damping-antidamping process: Application to the Fermi gas and to the Morse potential

abstract: The object of this paper is to investigate, classically and quantum mechanically, the relation existing between the position-dependent effective mass and damping-antidamping dynamics. The quantization of the equations of motion is carried out using the geometric interpretation of the motion, and we compare it with the one based on the ordering ambiguity scheme. Furthermore, we apply the obtained results to a Fermi gas of damped-antidamped particles, and we solve the Schr\"odinger equation for an exponentially increasing (decreasing) mass in the presence of the Morse potential.


title: Nobody but You: Sensor Selection for Voltage Regulation in Smart Grid

abstract: The increasing availability of distributed energy resources (DERs) and sensors in smart grid, as well as overlaying communication network, provides substantial potential benefits for improving the power system's reliability. In this paper, the problem of sensor selection is studied for the MAC layer design of wireless sensor networks for regulating the voltages in smart grid. The framework of hybrid dynamical system is proposed, using Kalman filter for voltage state estimation and LQR feedback control for voltage adjustment. The approach to obtain the optimal sensor selection sequence is studied. A sub- optimal sequence is obtained by applying the sliding window algorithm. Simulation results show that the proposed sensor selection strategy achieves a 40% performance gain over the baseline algorithm of the round-robin sensor polling.


title: Dimming of the 17th Century Sun

abstract: Reconstructions of total solar irradiance (TSI) rely mainly on linear relations between TSI variation and indices of facular area. When these are extrapolated to the prolonged 15th - 17th century Sp\"orer and Maunder solar activity minima, the estimated solar dimming is insufficient to explain the mid- millennial climate cooling of the Little Ice Age. We draw attention here to evidence that the relation departs from linearity at the lowest activity levels. Imaging photometry and radiometry indicate an increased TSI contribution per unit area from small network faculae by a factor of 2-4 compared to larger faculae in and around active regions. Even partial removal of this more TSI - effective network at prolonged minima could enable climatically significant solar dimming, yet be consistent with the weakened but persistent 11- yr cycle observed in Be 10 during the Maunder Minimum. The mechanism we suggest would not alter previous findings that increased solar radiative forcing is insufficient to account for 20th century global warming.


title: Anomaly-induced charges in nucleons

abstract: We show a novel charge structure of baryons in electromagnetic field due to the chiral anomaly. A key connection is to treat baryons as solitons of mesons. We use Skyrmions to calculate the charge distributions in a single nucleon and find an additional charge. We also perform calculations of charge distribution for classical multi-baryons with B=2, 3,...,8 and 17; they show amusing charge distributions.


title: SDRE Based Attitude Control Using Modified Rodriguez Parameters

abstract: The purpose of this paper is to present an application of the State Dependent Riccati Equation (SDRE) method to satellite attitude control where the satellite kinematics is modeled by Modified Rodriguez Parameters (MRP). The SDRE methodology is applicable on special forms of nonlinear systems where satellite model is one of the candidates. It is not easy to find an analytical solution from the SDRE. Thus point wise solutions are interpolated with respect to the operating conditions. The point wise solutions are obtained from the MATLAB algorithms which are derived from the positive definite solutions of the SDRE. The global stability analysis is difficult due to the nature of the methodology. The resultant attitude controllers outside the breakpoints (the selected operating conditions for interpolation) are suboptimal. The performance of the designs is examined by simulations on MATLAB - Simulink environment. The simulation results show that, the designed attitude controllers are working satisfactorily even in the presence of inertial uncertainties.


title: Left-Handed W Bosons at the LHC

abstract: The production of W bosons in association with jets is an important background to new physics at the LHC. Events in which the W carries large transverse momentum and decays leptonically lead to large missing energy and are of particular importance. We show that the left-handed nature of the W coupling, combined with valence quark domination at a pp machine, leads to a large left-handed polarization for both W^+ and W^- bosons at large transverse momenta. The polarization fractions are very stable with respect to QCD corrections. The leptonic decay of the W bosons translates the common left-handed polarization into a strong asymmetry in transverse momentum distributions between positrons and electrons, and between neutrinos and anti-neutrinos (missing transverse energy). Such asymmetries may provide an effective experimental handle on separating W + jets from top quark production, which exhibits very little asymmetry due to C invariance, and from various types of new physics.


title: Stellar Dynamics and Black Holes

abstract: Chandrasekhar's most important contribution to stellar dynamics was the concept of dynamical friction. I briefly review that work, then discuss some implications of Chandrasekhar's theory of gravitational encounters for motion in galactic nuclei.


title: On matrix variance inequalities

abstract: Olkin and Shepp (2005, J. Statist. Plann. Inference, vol. 130, pp. 351--358) presented a matrix form of Chernoff's inequality for Normal and Gamma (univariate) distributions. We extend and generalize this result, proving Poincare-type and Bessel-type inequalities, for matrices of arbitrary order and for a large class of distributions.


title: Numerical treatment of interfaces in Quantum Mechanics

abstract: In this article we develop a numerical scheme to deal with interfaces between touching numerical grids when solving Schr\"o{}dinger equation. In order to pass the information among grids we use the values of the fields only at the contact point between them. Surprisingly we obtain a convergent methods which is third order accurate with respect to the spatial resolution. In test cases, at the minimal resolution needed to describe correctly the waves, the error of this approximation is similar to that of a homogeneous (centered differences everywhere) scheme with three points stencil, that is a sixth order finite difference operator. The semi-discrete approximation preserves the norm and uses standard finite difference operators satisfying summation by parts. For the time integrator we use a semi-implicit IMEX Runge Kutta method.


title: Dissipation-induced pure Gaussian state

abstract: This paper provides some necessary and sufficient conditions for a generalMarkovian Gaussian master equation to have a unique pure steady state. The conditions are described by simple matrix equations; thus the so-called environment engineering problem for pure-Gaussian-state preparation can be straightforwardly dealt with in the linear algebraic framework. In fact, based on one of those conditions, for an arbitrary given pure Gaussian state,we obtain a complete parametrization of the Gaussian master equation having that state as a unique steady state; this leads to a systematic procedure for engineering a desired dissipative system.We demonstrate some examples including Gaussian cluster states.


title: Bianchi $VII_A$ solutions of quadratic gravity

abstract: It is believed that soon after the Planck time, Einstein's general relativity theory should be corrected to an effective quadratic theory. Numerical solutions for the anisotropic generalization of the Friedmann "open" model $H^ 3$ for this effective gravity are given. It must be emphasized that although numeric, these solutions are exact in the sense that they depend only on the precision of the machine. The solutions are identified asymptotically in a certain way. It is found solutions which asymptote de Sitter space, Riemann flat space and a singularity. The question of isotropisation of an initially anisotropic Universe is of great importance in the context of cosmology. Although isotropisation is not directly discussed in this present work, we show that sufficiently small anisotropies, do not increase indefinitely according to particular quadratic gravity theories. It can be understood as weak isotropisation, and we stress that this result is strongly dependent on initial conditions.


title: Complexity in human transportation networks: A comparative analysis of worldwide air transportation and global cargo ship movements

abstract: We present a comparative network theoretic analysis of the two largest global transportation networks: The worldwide air-transportation network (WAN) and the global cargoship network (GCSN). We show that both networks exhibit striking statistical similarities despite significant differences in topology and connectivity. Both networks exhibit a discontinuity in node and link betweenness distributions which implies that these networks naturally segragate in two different classes of nodes and links. We introduce a technique based on effective distances, shortest paths and shortest-path trees for strongly weighted symmetric networks and show that in a shortest-path-tree representation the most significant features of both networks can be readily seen. We show that effective shortest-path distance, unlike conventional geographic distance measures, strongly correlates with node centrality measures. Using the new technique we show that network resilience can be investigated more precisely than with contemporary techniques that are based on percolation theory. We extract a functional relationship between node characteristics and resilience to network disruption. Finally we discuss the results, their implications and conclude that dynamic processes that evolve on both networks are expected to share universal dynamic characteristics.


title: Logarithmic terms in entanglement entropies of 2D quantum critical points and Shannon entropies of spin chains

abstract: Universal logarithmic terms in the entanglement entropy appear at quantum critical points (QCPs) in one dimension (1D) and have been predicted in 2D at QCPs described by 2D conformal field theories. The entanglement entropy in a strip geometry at such QCPs can be obtained via the "Shannon entropy" of a 1D spin chain with open boundary conditions. The Shannon entropy of the XXZ chain is found to have a logarithmic term that implies, for the QCP of the square-lattice quantum dimer model, a logarithm with universal coefficient $\pm 0.25$. However, the logarithm in the Shannon entropy of the transverse-field Ising model, which corresponds to entanglement in the 2D Ising conformal QCP, is found to have a singular dependence on replica or R\'enyi index resulting from flows to different boundary conditions at the entanglement cut.


title: Using a Non-Commutative Bernstein Bound to Approximate Some Matrix Algorithms in the Spectral Norm

abstract: We focus on \emph{row sampling} based approximations for matrix algorithms, in particular matrix multipication, sparse matrix reconstruction, and \math{\ell_2} regression. For \math{\matA\in\R^{m\times d}} (\math{m} points in \math{d\ll m} dimensions), and appropriate row-sampling probabilities, which typically depend on the norms of the rows of the \math{m\times d} left singular matrix of \math{\matA} (the \emph{leverage scores}), we give row-sampling algorithms with linear (up to polylog factors) dependence on the stable rank of \math{\matA}. This result is achieved through the application of non-commutative Bernstein bounds. Keywords: row-sampling; matrix multiplication; matrix reconstruction; estimating spectral norm; linear regression; randomized


title: Ricci flat Kahler metrics with edge singularities

abstract: We construct Ricci flat Kahler metrics with cone singularities along a complex hypersurface. This construction is inspired in part by R. Mazzeo's program in the case of negative Einstein constant, and uses the linear theory developed recently by S. Donaldson.


title: The diffusion-induced nova scenario. CK Vul and PB 8 as possible observational counterparts

abstract: We propose a scenario for the formation of DA white dwarfs with very thin helium buffers. For these stars we explore the possible occurrence of diffusion-induced CNO- flashes, during their early cooling stage. In order to obtain very thin helium buffers, we simulate the formation of low mass remnants through an AGB final/late thermal pulse (AFTP/LTP scenario). Then we calculate the consequent white dwarf cooling evolution by means of a consistent treatment of element diffusion and nuclear burning. Based on physically sounding white dwarf models, we find that the range of helium buffer masses for these diffusion-induced novas to occur is significantly smaller than that predicted by the only previous study of this scenario. As a matter of fact, we find that these flashes do occur only in some low-mass (M < 0.6M) and low metallicity (Z_ZAMS <0.001) remnants about 10^6 - 10^7 yr after departing from the AGB. For these objects, we expect the luminosity to increase by about 4 orders of magnitude in less than a decade. We also show that diffusion-induced novas should display a very typical eruption lightcurve, with an increase of about a few magnitudes per year before reaching a maximum of M_V ~ -5 to -6. Our simulations show that surface abundances after the outburst are characterized by logNH/NHe ~ -0.15...0.6 and N>C>O by mass fractions. Contrary to previous speculations we show that these events are not recurrent and do not change substantially the final H-content of the cool (DA) white dwarf. (Abridged)


title: Hubble Space Telescope Observations of Main Belt Comet (596) Scheila

abstract: We present Hubble Space Telescope Observations of (596) Scheila during its recent dust outburst. The nucleus remained point-like with absolute magnitude V(1,1,0) = 8.85+/-0.02 in our data, equal to the pre-outburst value, with no secondary fragments of diameter =>100 m (for assumed albedos 0.04). We find a coma having a peak scattering cross-section ~2.2x10^4 km^2 corresponding to a mass in micron-sized particles of ~4x10^7 kg. The particles are deflected by solar radiation pressure on projected spatial scales ~2x10^4 km, in the sunward direction, and swept from the vicinity of the nucleus on timescales of weeks. The coma fades by ~30% between observations on UT 2010 December 27 and 2011 January 04. The observed mass loss is inconsistent with an origin either by rotational instability of the nucleus or by electrostatic ejection of regolith charged by sunlight. Dust ejection could be caused by the sudden but unexplained exposure of buried ice. However, the data are most simply explained by the impact, at ~5 km/s of a previously unknown asteroid ~35 m in diameter.


title: Combinatorics of embeddings

abstract: We offer the following explanation of the statement of the Kuratowski graph planarity criterion and of 6/7 of the statement of the Robertson-Seymour-Thomas intrinsic linking criterion. Let us call a cell complex 'dichotomial' if to every cell there corresponds a unique cell with the complementary set of vertices. Then every dichotomial cell complex is PL homeomorphic to a sphere; there exist precisely two 3-dimensional dichotomial cell complexes, and their 1-skeleta are K_5 and K_{3,3}; and precisely six 4-dimensional ones, and their 1-skeleta all but one graphs of the Petersen family. In higher dimensions n>2, we observe that in order to characterize those compact n-polyhedra that embed in S^{2n} in terms of finitely many "prohibited minors", it suffices to establish finiteness of the list of all (n-1)-connected n-dimensional finite cell complexes that do not embed in S^{2n} yet all their proper subcomplexes and proper cell-like combinatorial quotients embed there. Our main result is that this list contains the n-skeleta of (2n+1)-dimensional dichotomial cell complexes. The 2-skeleta of 5-dimensional dichotomial cell complexes include (apart from the three joins of the i-skeleta of (2i+2)-simplices) at least ten non-simplicial complexes.


title: Understanding the errors of SHAPE-directed RNA structure modeling

abstract: Single-nucleotide-resolution chemical mapping for structured RNA is being rapidly advanced by new chemistries, faster readouts, and coupling to computational algorithms. Recent tests have shown that selective 2'-hydroxyl acylation by primer extension (SHAPE) can give near-zero error rates (0-2%) in modeling the helices of RNA secondary structure. Here, we benchmark the method using six molecules for which crystallographic data are available: tRNA(phe) and 5S rRNA from Escherichia coli, the P4-P6 domain of the Tetrahymena group I ribozyme, and ligand-bound domains from riboswitches for adenine, cyclic di-GMP, and glycine. SHAPE-directed modeling of these highly structured RNAs gave an overall false negative rate (FNR) of 17% and a false discovery rate (FDR) of 21%, with at least one helix prediction error in five of the six cases. Extensive variations of data processing, normalization, and modeling parameters did not significantly mitigate modeling errors. Only one varation, filtering out data collected with deoxyinosine triphosphate during primer extension, gave a modest improvement (FNR = 12%, and FDR = 14%). The residual structure modeling errors are explained by the insufficient information content of these RNAs' SHAPE data, as evaluated by a nonparametric bootstrapping analysis. Beyond these benchmark cases, bootstrapping suggests a low level of confidence (<50%) in the majority of helices in a previously proposed SHAPE-directed model for the HIV-1 RNA genome. Thus, SHAPE-directed RNA modeling is not always unambiguous, and helix-by-helix confidence estimates, as described herein, may be critical for interpreting results from this powerful methodology.


title: Popping star clusters as building blocks of the Milky Way Thick Disc

abstract: It is widely believed that star clusters form with low star formation efficiencies. With the onset of stellar winds by massive stars or finally when the first super nova blows off, the residual gas is driven out of the embedded star cluster. Due to this fact a large amount, if not all, of the stars become unbound and disperse in the gravitational potential of the galaxy. In this context, Kroupa (2002) suggested a new mechanism for the emergence of thickened Galactic discs. Massive star clusters add kinematically hot components to the galactic field populations, building up in this way, the Galactic thick disc as well. In this work we perform, for the first time, numerical simulations to investigate this scenario for the formation of the galactic discs of the Milky Way. We find that a significant kinematically hot population of stars may be injected into the disk of a galaxy such that a thick disk emerges. For the MW the star clusters that formed the thick disk must have had masses of about 10^6 Msol.


title: Hubbard model description of silicon spin qubits: charge stability diagram and tunnel coupling in Si double quantum dots

abstract: We apply the recently introduced Hubbard model approach to quantitatively describe the experimental charge stability diagram and tunnel coupling of silicon double quantum dot systems. The results calculated from both the generalized Hubbard model and the microscopic theory are compared with existing experimental data, and excellent agreement between theory and experiment is found. The central approximation of our theory is a reduction of the full multi-electron multi-band system to an effective two-electron model, which is numerically tractable. In the microscopic theory we utilize the Hund-Mulliken approximation to the electron wave functions and compare the results calculated with two different forms of confinement potentials (biquadratic and Gaussian). We discuss the implications of our work for future studies.


title: Minimal flavour violation extensions of the seesaw

abstract: We analyze the most natural formulations of the minimal lepton flavour violation hypothesis compatible with a type-I seesaw structure with three heavy singlet neutrinos N, and satisfying the requirement of being predictive, in the sense that all LFV effects can be expressed in terms of low energy observables. We find a new interesting realization based on the flavour group $SU(3)_e\times SU(3)_{\ell+N}$ (being $e$ and $\ell$ respectively the SU(2) singlet and doublet leptons). An intriguing feature of this realization is that, in the normal hierarchy scenario for neutrino masses, it allows for sizeable enhancements of $\mu \to e$ transitions with respect to LFV processes involving the $\tau$ lepton. We also discuss how the symmetries of the type-I seesaw allow for a strong suppression of the N mass scale with respect to the scale of lepton number breaking, without implying a similar suppression for possible mechanisms of N production


title: Coordinates and 2MASS and OGLE identifications for all stars in Arp's 1965 finding chart for Baade's Window

abstract: Aims: We seek to provide 2MASS and OGLE identifications and coordinates for all stars in the finding chart published by Arp\,(1965). This chart covers the low extinction area around NGC 6522, also known as Baade's window, at coordinates (l,b)=(1.02,-3.92). Methods: A cross correlation, using numerical techniques, was performed between a scan of the original finding chart from Arp (1965) and 2MASS and OGLE-II images and stellar coordinates. Results: We provide coordinates for all stars in Arp's finding chart and 2MASS and OGLE identifications wherever possible. Two identifications in quadrant II do not appear in the original finding chart.


title: A large population of ultra-compact dwarf galaxies in the Hydra I cluster

abstract: We performed a large spectroscopic survey of compact, unresolved objects in the core of the Hydra I galaxy cluster (Abell 1060), with the aim of identifying ultra-compact dwarf galaxies (UCDs), and investigating the properties of the globular cluster (GC) system around the central cD galaxy NGC 3311. We obtained VIMOS medium resolution spectra of about 1200 candidate objects with apparent magnitudes 18.5 < V < 24.0 mag, covering both the bright end of the GC luminosity function and the luminosity range of all known UCDs. By means of spectroscopic redshift measurements, we identified 118 cluster members, from which 52 are brighter than M_V = -11.0 mag, and can therefore be termed UCDs. The brightest UCD in our sample has an absolute magnitude of M_V = -13.4 mag (corresponding to a mass of > 5 x 10^7 M_sun) and a half-light radius of 25 pc. This places it among the brightest and most massive UCDs ever discovered. Most of the GCs/UCDs are both spatially and dynamically associated to the central cD galaxy. The overall velocity dispersion of the GCs/UCDs is comparable to what is found for the cluster galaxies. However, when splitting the sample into a bright and a faint part, we observe a lower velocity dispersion for the bright UCDs/GCs than for the fainter objects. At a dividing magnitude of M_V = -10.75 mag, the dispersions differ by more than 200 km/s, and up to 300 km/s for objects within 5 arcmin around NGC 3311. We interpret these results in the context of different UCD formation channels, and conclude that interaction driven formation seems to play an important role in the centre of Hydra I.


title: Constraining the role of star cluster mergers in nuclear cluster formation: Simulations confront integral-field data

abstract: We present observations and dynamical models of the stellar nuclear clusters (NCs) at the centres of NGC 4244 and M33. We then compare these to an extensive set of simulations testing the importance of purely stellar dynamical mergers on the formation and growth of NCs. Mergers of star clusters are able to produce a wide variety of observed properties, including densities, structural scaling relations, shapes (including the presence of young discs) and even rapid rotation. Nonetheless, difficulties remain, most notably that the second order kinematic moment V_rms = (V^2 + sigma^2)^(1/2) of the models is too centrally peaked to match observations. This can be remedied by the merger of star clusters onto a pre-existing nuclear disc, but the line-of-sight velocity V is still more slowly rising than in NGC 4244. Our results therefore suggest that purely stellar dynamical mergers cannot form NCs, and that gas dissipation is a necessary ingredient for at least ~50% of a NC's mass. The negative vertical anisotropy found in NGC 4244 however requires at least 10% of the mass to be accreted as stars, since gas dissipation and in situ star formation leads to positive vertical anisotropy.


title: Measuring the dark matter equation of state

abstract: The nature of the dominant component of galaxies and clusters remains unknown. While the astrophysics community supports the cold dark matter (CDM) paradigm as a clue factor in the current cosmological model, no direct CDM detections have been performed. Faber and Visser 2006 have suggested a simple method for measuring the dark matter equation of state that combines kinematic and gravitational lensing data to test the widely adopted assumption of pressureless dark matter. Following this formalism, we have measured the dark matter equation of state for first time using improved techniques. We have found that the value of the equation of state parameter is consistent with pressureless dark matter within the errors. Nevertheless, the measured value is lower than expected because typically the masses determined with lensing are larger than those obtained through kinematic methods. We have tested our techniques using simulations and we have also analyzed possible sources of error that could invalidate or mimic our results. In the light of this result, we can now suggest that the understanding of the nature of dark matter requires a complete general relativistic analysis.


title: CN Morphology Studies of Comet 103P/Hartley 2

abstract: We report on narrowband CN imaging of Comet 103P/Hartley 2 obtained at Lowell Observatory on 39 nights from 2010 July until 2011 January. We observed two features, one generally to the north and the other generally to the south. The CN morphology varied during the apparition: no morphology was seen in July; in August and September the northern feature dominated and appeared as a mostly face-on spiral; in October, November, and December the northern and southern features were roughly equal in brightness and looked like more side-on corkscrews; in January the southern feature was dominant but the morphology was indistinct due to very low signal. The morphology changed smoothly during each night and similar morphology was seen from night to night. However, the morphology did not exactly repeat each rotation cycle, suggesting that there is a small non-principal axis rotation. Based on the repetition of the morphology, we find evidence that the fundamental rotation period was increasing: 16.7 hr from August 13-17, 17.2 hr from September 10-13, 18.2 hr from October 12-19, and 18.7 hr from October 31-November 7. We conducted Monte Carlo jet modeling to constrain the pole orientation and locations of the active regions based on the observed morphology. Our preliminary, self-consistent pole solution has an obliquity of 10 deg relative to the comet's orbital plane (i.e., it is centered near RA = 257 deg and Dec=+67 deg with an uncertainty around this position of about 15 deg) and has two mid-latitude sources, one in each hemisphere.


title: The correlation structure of dark matter halo properties

abstract: We investigate the correlation between nine different dark matter halo properties using a rank correlation analysis and a Principal Component Analysis for a sample of haloes spanning five orders of magnitude in mass. We consider mass and dimensionless measures of concentration, age, relaxedness, sphericity, triaxiality, substructure, spin, and environment, where the latter is defined in a way that makes it insensitive to mass. We find that concentration is the most fundamental property. Except for environment, all parameters are strongly correlated with concentration. Concentration, age, substructure, mass, sphericity and relaxedness can be considered a single family of parameters, albeit with substantial scatter. In contrast, spin, environment, and triaxiality are more independent, although spin does correlate strongly with substructure and both spin and triaxiality correlate substantially with concentration. Although mass sets the scale of a halo, all other properties are more sensitive to concentration.


title: Generalised massive gravity one-loop partition function and AdS/(L)CFT

abstract: The graviton 1-loop partition function is calculated for Euclidean generalised massive gravity (GMG) using AdS heat kernel techniques. We find that the results fit perfectly into the AdS/(L)CFT picture. Conformal Chern-Simons gravity, a singular limit of GMG, leads to an additional contribution in the 1-loop determinant from the conformal ghost. We show that this contribution has a nice interpretation on the conformal field theory side in terms of a semi-classical null vector at level two descending from a primary with conformal weights (3/2,-1/2).


title: On the Effective Description of Large Volume Compactifications

abstract: We study the reliability of the Two-Step moduli stabilization in the type-IIB Large Volume Scenarios with matter and gauge interactions. The general analysis is based on a family of N=1 Supergravity models with a factorizable Kaehler invariant function, where the decoupling between two sets of fields without a mass hierarchy is easily understood. For the Large Volume Scenario particular analyses are performed for explicit models, one of such developed for the first time here, finding that the simplified version, where the Dilaton and Complex structure moduli are regarded as frozen by a previous stabilization, is a reliable supersymmetric description whenever the neglected fields stand at their leading F-flatness conditions and be neutral. The terms missed by the simplified approach are either suppressed by powers of the Calabi-Yau volume, or are higher order operators in the matter fields, and then irrelevant for the moduli stabilization rocedure. Although the power of the volume suppressing such corrections depends on the particular model, up to the mass level it is independent of the modular weight for the matter fields. This at least for the models studied here but we give arguments to expect the same in general. These claims are checked through numerical examples. We discuss how the factorizable models present a context where despite the lack of a hierarchy with the supersymmetry breaking scale, the effective theory still has a supersymmetric description. This can be understood from the fact that it is possible to find vanishing solution for the auxiliary components of the fields being integrated out, independently of the remaining dynamics. Our results settle down the question on the reliability of the way the Dilaton and Complex structure are treated in type-IIB compactifications with large compact manifold volumina.


title: Resolvents and Seiberg-Witten representation for Gaussian beta-ensemble

abstract: The exact free energy of matrix model always obeys the Seiberg-Witten (SW) equations on a complex curve defined by singularities of the quasiclassical resolvent. The role of SW differential is played by the exact one-point resolvent. We show that these properties are preserved in generalization of matrix models to beta-ensembles. However, since the integrability and Harer-Zagier topological recursion are still unavailable for beta-ensembles, we need to rely upon the ordinary AMM/EO recursion to evaluate the first terms of the genus expansion. Consideration in this paper is restricted to the Gaussian model.


title: Evaluation of Polarization Mode Dispersion and Chromatic Dispersion using Quantum Interferometry

abstract: The use of quantum correlations between photons to measure polarization mode dispersion (PMD) and chromatic dispersion is investigated. Two types of apparatus are discussed which use coincidence counting of entangled photon pairs to allow sub-femtosecond resolution for measurement of both PMD and chromatic dispersion, as well as separation of even-order and odd-order chromatic effects in the PMD. Group delays can be measured with a resolution of order 0.1 fs, whereas attosecond resolution can be achieved for phase delays.


title: Constraints on Scalar Asymmetric Dark Matter from Black Hole Formation in Neutron Stars

abstract: We consider possibly observable effects of asymmetric dark matter (ADM) in neutron stars. Since dark matter does not self-annihilate in the ADM scenario, dark matter accumulates in neutron stars, eventually reaching the Chandrasekhar limit and forming a black hole. We focus on the case of scalar ADM, where the constraints from Bose-Einstein condensation and subsequent black hole formation are most severe due to the absence of Fermi degeneracy pressure. We also note that in some portions of this constrained parameter space, non-trivial effects from Hawking radiation can modify our limits. We find that for scalar ADM with mass between 100 keV and 10^5 GeV, the constraint from pulsars in globular clusters on the scattering cross-section with neutrons ranges from \sigma_n < 10^{-45} cm^2 to 10^{-52} cm}^2. In particular, for scalar ADM with mass between 1 GeV and 1 TeV (in the case where black hole evaporation due to Hawking radiation is unimportant), the constraint on the scattering cross-section is below what is reachable with ton scale direct detection experiments.


title: Denominator-preserving maps

abstract: Let F be a continuous injective map from an open subset of R^n to R^n. Assume that, for infinitely many k>1, F induces a bijection between the rational points of denominator k in the domain and those in the image (the denominator of (a_1/b_1,...,a_n/b_n) being the l.c.m. of b_1,...,b_n). Then F preserves the Lebesgue measure.


title: Telescope Bibliometrics 101

abstract: During recent years, bibliometric studies have become increasingly important in evaluating individual scientists, institutes, and entire observatories. In astronomy, often librarians are involved in maintaining publication databases and compiling statistics for their institutions. In this paper, we present a look behind the scenes to understand who is interested in bibliometric statistics, which methodologies astronomy librarians apply, and what kind of features next-generation bibliographies may include.


title: From convective to stellar dynamos

abstract: Convectively driven dynamos with rotation generating magnetic fields on scales large compared with the scale of the turbulent eddies are being reviewed. It is argued that such fields can be understood as the result of an alpha effect. Simulations in Cartesian domains show that such large-scale magnetic fields saturate on a time scale compatible with the resistive one, suggesting that the magnitude of the alpha effect is here still constrained by approximate magnetic helicity conservation. It is argued that, in the absence of shear and/or any other known large-scale dynamo effects, these simulations prove the existence of turbulent alpha^2-type dynamos. Finally, recent results are discussed in the context of solar and stellar dynamos.


title: Removable Matter-Power-Spectrum Covariance from Bias Fluctuations

abstract: We find a simple, accurate model for the covariance matrix of the real-space cosmological matter power spectrum on slightly nonlinear scales (k~0.1-0.8 h/Mpc at z=0), where off-diagonal matrix elements become substantial. The model includes a multiplicative, scale-independent modulation of the power spectrum. It has only one parameter, the variance (among realizations) of the variance of the nonlinear density field in cells, with little dependence on the cell size between 2-8 Mpc/h. Furthermore, we find that this extra covariance can be modeled out by instead measuring the power spectrum of (delta/sigma_cell), i.e. the ratio of the overdensity to its dispersion in cells a few Mpc in size. Dividing delta by sigma_cell essentially removes the non-Gaussian part of the covariance matrix, nearly diagonalizing it.


title: High Energy Neutrino Emission from the Earliest Gamma-Ray Bursts

abstract: We discuss the high energy neutrino emission from gamma-ray bursts resulting from the earliest generation (`population III') stars forming in the Universe, whose core collapses into a black hole. These gamma-ray bursts are expected to produce a highly relativistic, magnetically dominated jet, where protons can be accelerated to ultra-high energies. These interact with the photons produced by the jet, leading to ultra-high energy photo-meson neutrinos as well as secondary leptons and photons. The photon luminosity and the shock properties, and thus the neutrino spectrum, depend on the mass of the black holes as well as on the density of the surrounding external gas. We calculate the individual source neutrino spectral fluxes and the expected diffuse neutrino flux for various source parameters and evolution scenarios. Both the individual and diffuse signals appear detectable in the 1-300 PeV range with current and planned neutrino detectors such as IceCube and ARIANNA, provided the black hole mass is in excess of 30-100 solar masses. This provides a possible test for the debated mass of the progenitor stellar objects, as well as a probe for the early cosmological environment and the formation rate of the earliest structures.


title: Proof of the outage probability conjecture for MISO channels

abstract: In Telatar 1999, it is conjectured that the covariance matrices minimizing the outage probability for MIMO channels with Gaussian fading are diagonal with either zeros or constant values on the diagonal. In the MISO setting, this is equivalent to conjecture that the Gaussian quadratic forms having largest tale probability correspond to such diagonal matrices. We prove here the conjecture in the MISO setting.


title: Unicity conditions for low-rank matrix recovery

abstract: Low-rank matrix recovery addresses the problem of recovering an unknown low-rank matrix from few linear measurements. Nuclear-norm minimization is a tractible approach with a recent surge of strong theoretical backing. Analagous to the theory of compressed sensing, these results have required random measurements. For example, m >= Cnr Gaussian measurements are sufficient to recover any rank-r n x n matrix with high probability. In this paper we address the theoretical question of how many measurements are needed via any method whatsoever --- tractible or not. We show that for a family of random measurement ensembles, m >= 4nr - 4r^2 measurements are sufficient to guarantee that no rank-2r matrix lies in the null space of the measurement operator with probability one. This is a necessary and sufficient condition to ensure uniform recovery of all rank-r matrices by rank minimization. Furthermore, this value of $m$ precisely matches the dimension of the manifold of all rank-2r matrices. We also prove that for a fixed rank-r matrix, m >= 2nr - r^2 + 1 random measurements are enough to guarantee recovery using rank minimization. These results give a benchmark to which we may compare the efficacy of nuclear-norm minimization.

categories: math.NA cs.IT cs.SY math.IT math.OC math.PR


title: Universality of Phases in QCD and QCD-like Theories

abstract: We argue that the whole or the part of the phase diagrams of QCD and QCD-like theories should be universal in the large-N_c limit through the orbifold equivalence. The whole phase diagrams, including the chiral phase transitions and the BEC-BCS crossover regions, are identical between SU(N_c) QCD at finite isospin chemical potential and SO(2N_c) and Sp(2N_c) gauge theories at finite baryon chemical potential. Outside the BEC-BCS crossover region in these theories, the phase diagrams are also identical to that of SU(N_c) QCD at finite baryon chemical potential. We give examples of the universality in some solvable cases: (i) QCD and QCD-like theories at asymptotically high density where the controlled weak-coupling calculations are possible, (ii) chiral random matrix theories of different universality classes, which are solvable large-N (large volume) matrix models of QCD. Our results strongly suggest that the chiral phase transition and the QCD critical point at finite baryon chemical potential can be studied using sign-free theories, such as QCD at finite isospin chemical potential, in lattice simulations.


title: Spindown of massive rotating stars

abstract: Models of rapidly rotating massive stars at low metallicities show significantly different evolution and higher metal yields compared to non-rotating stars. We estimate the spin-down time-scale of rapid rotating non-convective stars supporting an alpha-Omega dynamo. The magnetic dynamo gives rise to mass loss in a magnetically controlled stellar wind and hence stellar spin down owing to loss of angular momentum. The dynamo is maintained by strong horizontal rotation-driven turbulence which dominates over the Parker instability. We calculate the spin-down time-scale and find that it could be relatively short, a small fraction of the main-sequence lifetime. The spin-down time-scale decreases dramatically for higher surface rotations suggesting that rapid rotators may only exhibit such high surface velocities for a short time, only a small fraction of their main-sequence lifetime.


title: Deformation of quotients on a product

abstract: We consider the general problem of deforming a surjective map of modules $f : E \to F$ over a coproduct sheaf of rings $B=B_1 \otimes_A B_2$ when the domain module $E = B_1 \otimes_A E_2$ is obtained via extension of scalars from a $B_2$-module $E_2$. Assuming $B_1$ is flat over $A$, we show that the Atiyah class morphism $F \to \LL_{B/B_2} \otimes^{\bL} F[1]$ in the derived category $D(B)$ factors naturally through (the shift of) a morphism $\beta : \Ker f \to \LL_{B/B_2} \otimes^{\bL} F$. We describe the obstruction to lifting $f$ over a (square zero) extension $B_1' \to B_1$ in terms of $\beta$ and the class of the extension. As an application, we use the reduced Atiyah class to construct a perfect obstruction theory on the Quot scheme of a vector bundle on a smooth curve (and more generally).


title: A simple model of cerebral blood flow dependence on arterial blood pressure

abstract: It is shown that the dependence of the cerebral blood flow (CBF) on mean arterial blood pressure (MABP) can be described with a simple model having the following assumptions. Below certain MABP (denoted as MABP1) there are no autoregulatory or feedback mechanisms influencing CBF. Between MABP1 and MABP2 (MABP at which breakthrough accurs) there is a linear (on MABP) dependent feedback with a sloap depending very much on the individual considered. The classical autoregulation model with a plateau in between MABP1 and MABP2 is a particular case of this model. The model describes well the experiments performed on dogs (Harper 1966), for which the individual feedback sloap parameter varied to great extent, indicating the importance of mesurments on individuals against averaged mesurments (or measurments on diffent individuals) which superficially support the classical autoregulation. New effect of decreased CBF, while increasing MABP, was observed.


title: HD 135344B: a young star has reached its rotational limit

abstract: We search for periodic variations in the radial velocity of the young Herbig star HD 135344B with the aim to determine a rotation period. We analyzed 44 high-resolution optical spectra taken over a time range of 151 days. The spectra were acquired with FEROS at the 2.2m MPG/ESO telescope in La Silla. The stellar parameters of HD 135344B are determined by fitting synthetic spectra to the stellar spectrum. In order to obtain radial velocity measurements, the stellar spectra have been cross-correlated with a theoretical template computed from determined stellar parameters. We report the first direct measurement of the rotation period of a Herbig star from radial-velocity measurements. The rotation period is found to be 0.16 d (3.9 hr), which makes HD 135344B a rapid rotator at or close to its break-up velocity. The rapid rotation could explain some of the properties of the circumstellar environment of HD 135344B such as the presence of an inner disk with properties (composition, inclination), that are significantly different from the outer disk.


title: Modeling the secular evolution of migrating planet pairs

abstract: The subject of this paper is the secular behaviour of a pair of planets evolving under dissipative forces. In particular, we investigate the case when dissipative forces affect the planetary semi-major axes and the planets move inward/outward the central star, in a process known as planet migration. To perform this investigation, we introduce fundamental concepts of conservative and dissipative dynamics of the three-body problem. Based on these concepts, we develop a qualitative model of the secular evolution of the migrating planetary pair. Our approach is based on analysis of the energy and the orbital angular momentum exchange between the two-planet system and an external medium; thus no specific kind of dissipative forces is invoked. We show that, under assumption that dissipation is weak and slow, the evolutionary routes of the migrating planets are traced by the Mode I and Mode II stationary solutions of the conservative secular problem. The ultimate convergence and the evolution of the system along one of these secular modes of motion is determined uniquely by the condition that the dissipation rate is sufficiently smaller than the proper secular frequency of the system. We show that it is possible to reassemble the starting configurations and migration history of the systems on the basis of their final states and consequently to constrain the parameters of the physical processes involved.


title: Classification of p-adic 6-dimensional filiform Leibniz algebras by solution of x^q=a

abstract: In this paper we study the $p$-adic equation $x^q=a$ over the field of $p$-adic numbers. We construct an algorithm of calculation of criteria of solvability in the case of $q=p^m$ and present a computer program to compute the criteria for fixed value of $m \leq p-1$. Moreover, using this solvability criteria for $q=2,3,4,5,6$, we classify $p$-adic 6-dimensional filiform Leibniz algebras.


title: The imprint of a symbiotic binary progenitor on the properties of Kepler's supernova remnant

abstract: We present a model for the Type Ia supernova remnant (SNR) of SN 1604, also known as Kepler's SNR. We find that its main features can be explained by a progenitor model of a symbiotic binary consisting of a white dwarf and an AGB donor star with an initial mass of 4-5 M_sun. The slow, nitrogen rich wind emanating from the donor star has partially been accreted by the white dwarf, but has also created a circumstellar bubble. Based on observational evidence, we assume that the system moves with a velocity of 250 km/s. Due to the systemic motion the interaction between the wind and the interstellar medium has resulted in the formation of a bow shock, which can explain the presence of a one-sided, nitrogen rich shell. We present two-dimensional hydrodynamical simulations of both the shell formation and the SNR evolution. The SNR simulations show good agreement with the observed kinematic and morphological properties of Kepler's SNR. Specifically, the model reproduces the observed expansion parameters (m=V/(R/t)) of m=0.35 in the north and m=0.6 in the south of Kepler's SNR. We discuss the variations among our hydrodynamical simulations in light of the observations, and show that part of the blast wave may have traversed through the one-sided shell completely. The simulations suggest a distance to Kepler's SNR of 6 kpc, or otherwise require that SN 1604 was a sub-energetic Type Ia explosion. Finally, we discuss the possible implications of our model for Type Ia supernovae and their remnants in general.


title: Population persistence under advection-diffusion in river networks

abstract: An integro-differential equation on a tree graph is used to model the evolution and spatial distribution of a population of organisms in a river network. Individual organisms become mobile at a constant rate, and disperse according to an advection-diffusion process with coefficients that are constant on the edges of the graph. Appropriate boundary conditions are imposed at the outlet and upstream nodes of the river network. The local rates of population growth/decay and that by which the organisms become mobile, are assumed constant in time and space. Imminent extinction of the population is understood as the situation whereby the zero solution to the integro-differential equation is stable. Lower and upper bounds for the eigenvalues of the dispersion operator, and related Sturm-Liouville problems are found, and therefore sufficient conditions for imminent extinction are given in terms of the physical variables of the problem.


title: Passive Scalar Structures in Supersonic Turbulence

abstract: We conduct a systematic numerical study of passive scalar structures in supersonic turbulent flows. We find that the degree of intermittency in the scalar structures increases only slightly as the flow changes from transonic to highly supersonic, while the velocity structures become significantly more intermittent. This difference is due to the absence of shock-like discontinuities in the scalar field. The structure functions of the scalar field are well described by the intermittency model of She and L\'{e}v\^{e}que [Phys. Rev. Lett. 72, 336 (1994)], and the most intense scalar structures are found to be sheet-like at all Mach numbers.